A354792 Nonprime terms in A354790 in order of appearance.

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

Author

Michael De Vlieger and N. J. A. Sloane, Jul 17 2022

Keywords

• nonn

Links

• Rémy Sigrist, Table of n, a(n) for n = 1..5715
• Rémy Sigrist, C program

Crossrefs

Cf. A354790, A354791.

Programs

C

See Links section.

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

Extensions

More terms from Rémy Sigrist, Jul 17 2022