OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
15+21+25 = 61, 9+15+21+25+27 = 97.
MAPLE
N:= 1000: # for terms <= N
tmax:= 2*ceil(N/6)+7:
T:= remove(isprime, [seq(i, i=3..tmax, 2)]):
S:= ListTools:-PartialSums([0, op(T)]):
R:= {}:
for m from 3 by 2 while S[m] <= N do
for j from 1 do
v:= S[m+j]-S[j];
if v > N then break fi;
if isprime(v) then R:= R union {v} fi
od od:
sort(convert(R, list)); # Robert Israel, Apr 17 2026
MATHEMATICA
lst={}; Do[If[PrimeQ[m], Continue[]]; s=m; Do[If[PrimeQ[n], Continue[]]; s+=n; If[PrimeQ[s], If[s<=2917, AppendTo[lst, s]]], {n, m+2, 2*6!, 2}], {m, 1, 2*6!, 2}]; lst=Take[Union@lst, 200]
PROG
(PARI) N=1000; v=vector(N); L=listcreate(); n=9; while(n<N, if(!isprime(n), listput(L, n)); n=n+2); vv=Vec(L); forstep(k=3, #vv, 2, for(offs=1, #vv-k+1, s=sum(i=offs, offs+k-1, vv[i]); if(isprime(s)&&s<=N, v[s]=1))); for(n=60, #v, if(v[n], print1(n, ", "))) \\ Ralf Stephan, Nov 26 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Sep 13 2009
EXTENSIONS
Edited by Ralf Stephan, Nov 26 2013
Corrected by Robert Israel, Apr 17 2026
STATUS
approved
