OFFSET
1,1
COMMENTS
Sum 1/a(n) appears to converge albeit slowly.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
FORMULA
MATHEMATICA
With[{pc14=Select[Prime[Range[100]], Mod[#, 4]==1&]}, Flatten[Table[ Prime[ Take[ pc14, {n}]], {n, 40}]]] (* Harvey P. Dale, Aug 12 2011 *)
PROG
(PARI) pip(n, m, r) = {sr=0; forprime(x=3, n, if(x%m == r, v=prime(x); sr+=1.0/v; print1(v" "); ) ); print(); print("m="m" r="r" sr="sr); }
(Magma) [NthPrime(n): n in PrimesUpTo(400) | n mod 4 eq 1]; // Bruno Berselli, Aug 14 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Cino Hilliard, Dec 14 2002
EXTENSIONS
Offset changed from 0 to 1 by Vincenzo Librandi, Aug 14 2012
STATUS
approved