OFFSET
1,1
COMMENTS
a(n) is divisible by 4, for all n.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
first and second 4k+1 primes are 5 and 13, so a(1)=13-5=8;
MATHEMATICA
k=0; m=4; r=1; Do[s=Mod[Prime[n], m]; If[Equal[s, r], rp=ep; k=k+1; ep=Prime[n]; Print[(ep-rp)]; ], {n, 1, 1000}]
Differences[Select[Prime[Range[200]], Mod[#, 4]==1&]] (* Harvey P. Dale, Feb 05 2020 *)
PROG
(PARI) p=5; forprime(q=7, 1e3, if(q%4==1, print1(q-p", "); p=q)) \\ Charles R Greathouse IV, May 13 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Labos Elemer, Apr 07 2003
STATUS
approved