[go: up one dir, main page]

login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

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