OFFSET
1,1
COMMENTS
Start with decimal expansion of Pi, 3,1,4,1,5,9,2,6,5,3,5,8,9,7,9,3,2,3,8,4,6,2,6,4,3... Make the sections with two successive digits: s={3,1},{1,4},{4,1},{1,5},{5,9},{9,2},{2,6},{6,5},{5,3},{3,5},{5,8},{8,9},{9,7},{7,9},{9,3},{3,2},{2,3},... Then A104822(n) = prime number from digits of [s(m)]: 31,41,59,53,89,97,79,23,...
Note that leading zeros are allowed, so the sequence of digits ...02...05...07...03... leads to the terms a(12)=2, a(18)=5, a(22)=7, a(24)=3, etc. - M. F. Hasler, Oct 25 2011
MATHEMATICA
Select[ FromDigits /@ Partition[ RealDigits[Pi, 10, 330][[1]], 2, 1], PrimeQ[ # ] &] (* Robert G. Wilson v, Mar 28 2005 *)
PROG
(PARI) for(i=1, 99, isprime(p=Pi\.1^i%100)&print1(p", ")) \\ M. F. Hasler, Oct 25 2011
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Zak Seidov, Mar 27 2005
EXTENSIONS
More terms from Robert G. Wilson v, Mar 28 2005
STATUS
approved