OFFSET
1,1
COMMENTS
See A355061 for further details.
PROG
(Python)
from sympy import primefactors
from itertools import count, islice
def agen(): # generator of terms
an1, an, f1, f, pset = 2, 6, {2}, {2, 3}, {2, 12}
for n in count(4):
if an == 6: yield n-1
an2, an1, an, f2, f1 = an1, an, 6, f1, f
f = set(primefactors(an))
while an*an1 in pset or f1&f == set() or f2&f != set() or f <= f1:
an += 1; f = set(primefactors(an))
pset.add(an*an1)
print(list(islice(agen(), 54))) # Michael S. Branicky, Jun 20 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Scott R. Shannon, Jun 20 2022
STATUS
approved