OFFSET
1,1
LINKS
PROG
(PARI) do(lim)=my(v=List()); for(n=9, lim+1, if(ispseudoprime(10^n\3-2), next); forprime(p=7, , if(Mod(10, p)^n==7, listput(v, p); next(2)))); Vec(v) \\ tests numbers with A089017(n) <= m; Charles R Greathouse IV, Oct 23 2013
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Lekraj Beedassy, Nov 04 2003
EXTENSIONS
More terms from Ray Chandler, Nov 04 2003
STATUS
approved
