OFFSET
1,1
PROG
(PARI) is(n) = forprime(p=1, 19, if(Mod(n, p^2)^(p-1)!=1, return(0))); 1
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Felix Fröhlich, Dec 08 2020
STATUS
approved