OFFSET
1,1
COMMENTS
Primes in this sequence include 271639843 and 7245669763.
EXAMPLE
a(1) = 781 = g(5,4) = g(7,3) + g(19,2)
a(2) = 1464 = g(11,3) = g(37,2) + g(7,2)
a(3) = 5403 = g(73,2) = g(17,3) + g(13,2)
a(4) = 36673 = g(191,2) = g(2,14) + g(5,5)
MAPLE
N:= 10^10: # for terms <= N
S:= {}:
for p from 1 do
p:= nextprime(p);
if 1+p+p^2 > N then break fi;
t:= 1+p;
for k from 2 do
t:= t+p^k;
if t > N then break fi;
S:= S union {t};
od
od:
sort(convert(`union`(seq(map(`+`, S, S[i])) intersect S, i=1..nops(S)), list));
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Dec 06 2020
STATUS
approved