OFFSET
0,1
COMMENTS
32 is the first prime power > 1 not in this sequence.
PROG
(Python)
from itertools import count, islice
from sympy import factorint
def A062559_gen(): # generator of terms
aset2, alist = set(), []
for k in count(0):
if len(f:=factorint(k).values()) == 1 and max(f) > 1:
bset2 = set()
for a in alist:
if (b:=a+k) in aset2:
break
bset2.add(b)
else:
yield k
alist.append(k)
aset2.update(bset2)
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Jul 03 2001
EXTENSIONS
More terms from Sean A. Irvine, Apr 03 2023
STATUS
approved