OFFSET
2,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 2..10000
MAPLE
P:= select(isprime, {seq(i, i=3..10000, 2)}):
TP:= sort(convert(P intersect map(`-`, P, 2), list)):
f:= proc(p) local k;
for k from 6 by 6 do if isprime(p + k^2) and isprime(p + k^2 + 2) then return k fi od
end proc:
map(f, TP[2..-1]);
MATHEMATICA
With[{s = Select[Prime@ Range[3, 332], PrimeQ[# + 2] &]}, Array[Block[{k = 1}, While[! AllTrue[s[[#]] + k^2 + {0, 2}, PrimeQ], k++]; k] &, Length@ s]] (* Michael De Vlieger, Sep 03 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Sep 03 2019
STATUS
approved