[go: up one dir, main page]

login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A105050
a(0)=2, a(n) is the smallest prime of the form a(n-1)*k^2 + 1.
0
2, 3, 13, 53, 30529, 122117, 17584849, 70339397, 22789964629, 11030342880437, 9927308592393301, 91490075987496662017, 23421459452799145476353, 918214896387537699254943013, 528891780319221714770847175489, 8665362928750128574805560123211777, 311953065435004628693000164435623973
OFFSET
0,1
EXAMPLE
a(5) = 30529 = 53*24^2 + 1 is the smallest prime of the form 53*k^2 + 1.
PROG
(PARI) lista(nn) = my(k, p=2); print1(p); for(n=1, nn, k=1; while(!isprime(p*k^2+1), k++); print1(", ", p=p*k^2+1)); \\ Jinyuan Wang, Aug 02 2021
CROSSREFS
Sequence in context: A203985 A164511 A184256 * A184179 A100102 A208202
KEYWORD
nonn
AUTHOR
John L. Drost, Apr 04 2005
EXTENSIONS
a(15)-a(16) from Jinyuan Wang, Aug 02 2021
STATUS
approved