OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
GP:= [0, seq(op([m*(3*m-1)/2, m*(3*m+1)/2]), m=1..50)]:
N:= GP[-1]:
V:= Array(0..N, datatype=integer[4]):
for i from 1 to nops(GP) do
for j from 1 to i do
r:= GP[i]+GP[j];
if r > N then break fi;
V[r]:= V[r]+1
od od:
select(t -> V[t] = 0, [$0..N]); # Robert Israel, Feb 26 2025
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Jon Perry, Mar 29 2004
EXTENSIONS
Definition clarified by Robert Israel, Feb 26 2025
STATUS
approved