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,2

Links

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

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

Cf. A354790, A354791.

Sequence in context: A230941 A359571 A145000 * A209182 A275407 A362306

Adjacent sequences: A354789 A354790 A354791 * A354793 A354794 A354795

Keyword

nonn

Author

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

Extensions

More terms from Rémy Sigrist, Jul 17 2022

Status

approved