OFFSET
0,1
COMMENTS
LINKS
Stanislav Sykora and Amiram Eldar, Table of n, a(n) for n = 0..64 (terms 0..44 from Stanislav Sykora)
Clifford A. Pickover, Belphegor's Prime: 1000000000000066600000000000001
Wikipedia, Belphegor's prime
MATHEMATICA
Table[FactorInteger[(10^(n + 3) + 666)*10^(n + 1) + 1][[-1, 1]], {n, 20}] (* T. D. Noe, Nov 25 2013 *)
PROG
(PARI) default(factor_proven, 1);
Belphegor(k)=(10^(k+3)+666)*10^(k+1)+1;
LargestPrimeFactor(k)={local(f); f=factor(k); return(f[#f[, 1], 1])};
nmax=40; v=vector(nmax);
for (n=0, #v-1, v[n+1]=LargestPrimeFactor(Belphegor(n)); print(v[n+1]))
CROSSREFS
KEYWORD
nonn
AUTHOR
Stanislav Sykora, Nov 24 2013
STATUS
approved