OFFSET
1,1
COMMENTS
The first term is the Hardy-Ramanujan number.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..109 from Chai Wah Wu)
MATHEMATICA
Select[10*Range[0, 3*10^7] + 9, CompositeQ[#] && Divisible[# - 1, CarmichaelLambda[#]] &] (* Amiram Eldar, May 28 2022 *)
PROG
(Python)
from itertools import islice
from sympy import factorint, nextprime
def A352970_gen(): # generator of terms
p, q = 3, 5
while True:
for n in range(p+11-((p+2) % 10), q, 10):
f = factorint(n)
if max(f.values()) == 1 and not any((n-1) % (p-1) for p in f):
yield n
p, q = q, nextprime(q)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Omar E. Pol, Apr 12 2022
STATUS
approved