OFFSET
1,1
COMMENTS
Primes of the form (25*k^2 + 25*k + 2)/2.
Numbers k such that (25*k^2 + 25*k + 2)/2 is prime: 3, 4, 7, 12, 15, 19, 28, 31, 36, 39, 47, 48, 51, 55, 56, 60, 67, 71, 72, 80, 83, 87, 91, ...
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
OEIS Wiki, Centered polygonal numbers
Eric Weisstein's World of Mathematics, Centered Polygonal Number
MAPLE
select(isprime, [seq((25*k^2+25*k+2)/2, k=1..200)]); # Robert Israel, Sep 01 2016
MATHEMATICA
Intersection[Table[(25 k^2 + 25 k + 2)/2, {k, 0, 1000}], Prime[Range[28000]]]
PROG
(PARI) lista(nn) = for(n=1, nn, if(isprime(p=(25*n^2 + 25*n + 2)/2), print1(p, ", "))); \\ Altug Alkan, Aug 26 2016
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Ilya Gutkovskiy, Aug 26 2016
STATUS
approved