OFFSET
1,3
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[n_] := n (n + 1)^2/2; Table[k = 0; While[! PrimeQ[f[n] + f[k]], k++]; k, {n, 100}] (* T. D. Noe, Apr 03 2012 *)
PROG
(PARI) a(n)=my(K=n*(n+1)^2/2, m); while(!isprime(K+m*(m+1)^2/2), m++); m \\ Charles R Greathouse IV, Aug 03 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Gerasimov Sergey, Mar 27 2012
EXTENSIONS
Corrected by R. J. Mathar, Mar 31 2012
STATUS
approved