OFFSET
1,12
LINKS
Robert Israel, Table of n, a(n) for n = 1..1000
MAPLE
N:= 1000: # to get terms for primes <= N
P, NP:= selectremove(isprime, [seq(i, i=1..N, 2)]):
V:= Vector(N):
for ri from 1 while NP[ri] <= N/3 do
for si from ri+1 while NP[ri]+2*NP[si] <= N do
for ti from si+1 do
x:= NP[ri]+NP[si]+NP[ti];
if x > N then break fi;
V[x]:= V[x]+1;
od od od:
0, seq(V[P[i]], i=1..nops(P)); # Robert Israel, Nov 03 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved