OFFSET
0,1
LINKS
Max Alekseyev, Table of n, a(n) for n = 0..31 (shortened by N. J. A. Sloane, Jan 13 2019)
EXAMPLE
The corresponding prime powers are 2 + 3 = 5^1, 2 + 3 + 11 = 2^4, 2 + 3 + 11 + 13 = 29^1, etc.
MAPLE
a := [2, 3] ; while true do as := add(i, i=a) ; p := nextprime(op(-1, a)) ; while nops(numtheory[factorset](p+as)) > 1 do p := nextprime(p) ; od; a := [op(a), p] ; print(a) ; od: # R. J. Mathar, Apr 28 2008
PROG
(PARI) { printA139021() = my(a=2, s=2); print1(2, ", "); for(n=2, 100, if( s%2==0, until(isprimepower(s+a), a=nextprime(a+1)), t=log(s+a)\log(2) + 1; while( !ispseudoprime(2^t-s), t++); a=2^t-s; ); s+=a; print1(a, ", "); ); } /* Max Alekseyev, Oct 17 2015 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Apr 06 2008
EXTENSIONS
9 more terms from R. J. Mathar, Apr 28 2008
a(14)-a(19) from Donovan Johnson, Nov 26 2008
a(20)-a(21) from Max Alekseyev, Oct 14 2012
a(22)-a(31) in b-file from Max Alekseyev, Oct 17 2015
STATUS
approved