[go: up one dir, main page]

login
A354792
Nonprime terms in A354790 in order of appearance.
5
1, 6, 35, 22, 39, 85, 14, 33, 65, 34, 133, 69, 319, 155, 481, 82, 119, 57, 253, 1247, 235, 403, 74, 287, 901, 177, 1159, 737, 1633, 2117, 215, 611, 2449, 166, 3293, 679, 697, 159, 1121, 671, 1541, 2059, 365, 4343, 1339, 5029, 3379, 158, 3071, 623, 1649, 4633
OFFSET
1,2
LINKS
Rémy Sigrist, C program
PROG
(Python)
from math import lcm, gcd
from itertools import count, islice
from collections import deque
from sympy import factorint
def A354792_gen(): # generator of terms
aset, aqueue, c, b, f = {1}, deque([1]), 2, 1, True
yield 1
while True:
for m in count(c):
if m not in aset and gcd(m, b) == 1 and all(map(lambda n:n<=1, fs:=factorint(m).values())):
if len(fs) > 1:
yield m
aset.add(m)
aqueue.append(m)
if f: aqueue.popleft()
b = lcm(*aqueue)
f = not f
while c in aset:
c += 1
break
A354792_list = list(islice(A354792_gen(), 30)) # Chai Wah Wu, Jul 17 2022
(C) See Links section.
CROSSREFS
Sequence in context: A230941 A359571 A145000 * A209182 A275407 A362306
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Rémy Sigrist, Jul 17 2022
STATUS
approved