OFFSET
1,1
COMMENTS
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, ", ")));
CROSSREFS
KEYWORD
nonn
AUTHOR
Altug Alkan, Mar 18 2016
STATUS
approved