OFFSET
1,1
COMMENTS
Primes of the form 9*k^2 - 9*k + 1.
LINKS
OEIS Wiki, Figurate numbers
Eric Weisstein's World of Mathematics, Centered Polygonal Number
Eric Weisstein's World of Mathematics, Prime Number
PROG
(PARI) for(n=1, 1e3, if(isprime(k=9*n^2-9*n+1), print1(k, ", "))) \\ Altug Alkan, Nov 26 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 26 2015
STATUS
approved