OFFSET
1,2
COMMENTS
A000244 is a subsequence.
PROG
(PARI) is(n) = my(p=11, q=1006003); Mod(n, p^2)^(p-1)==1 && Mod(n, q^2)^(q-1)==1
(Python)
def ok(b): return pow(b, 10, 121)==1 and pow(b, 1006002, 1006003**2)==1
print([k for k in range(10**6) if ok(k)]) # Michael S. Branicky, May 25 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Fröhlich, May 25 2022
STATUS
approved