OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
MATHEMATICA
s = Select[Prime@ Range[10^6], PrimeQ[# + 2] &]; Select[Range@ 1300, s[[#]] + s[[# + 3]] == s[[# + 1]] + s[[# + 2]] &] (* after Robert G. Wilson v at A001359 *)
PROG
(PARI) t(n, p=3) = { while( p+2 < (p=nextprime( p+1 )) || n-->0, ); p-2}
b(n) = t(n) + t(n+3) - t(n+1) - t(n+2);
for(n=1, 2000, if(b(n) == 0, print1(n, ", ")));
(PARI) list(lim) = {my(k = 0, p1 = 2, t = [0, 0, 0, 0]); forprime(p2 = 3, lim, if(p2 - p1 == 2, k++; t = concat(t[2..4], p1); if(t[1] + t[4] == t[2] + t[3], print1(k-3, ", "))); p1 = p2); } \\ Amiram Eldar, Feb 22 2025
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Altug Alkan, Mar 18 2016
STATUS
approved