OFFSET
1,1
COMMENTS
Subsequence of A033210 (Primes of the form x^2+13*y^2).
LINKS
Jens Kruse Andersen, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1)=13 because the least prime of the form x^2+13*1^2 is 13 (at x=0).
a(100)=130121 because the least prime of the form x^2+13*100^2 is 130121 (at x=11).
a(200)=520361 because the least prime of the form x^2+13*200^2 is 520361 (at x=19).
PROG
(PARI) a(n) = {x = 0; while (!isprime(p=x^2+13*n^2), x++); p; } \\ Michel Marcus, Oct 06 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Oct 06 2014
STATUS
approved