|
| |
|
|
A156226
|
|
Primes of the form 9*n^2 + 1.
|
|
3
|
|
|
|
37, 577, 1297, 2917, 4357, 7057, 8101, 14401, 15877, 22501, 24337, 32401, 41617, 44101, 57601, 69697, 72901, 90001, 93637, 147457, 156817, 176401, 197137, 224677, 324901, 331777, 352837, 404497, 427717, 476101, 484417, 509797, 562501
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
1,1
|
|
|
COMMENTS
|
9*n^2 + 1 can be a prime only for n's of the form n=10m or n=10m+-2.
Primes in this sequence must end with 1 or 7 and will have to be 1 modulo 30 or 7 modulo 30.
|
|
|
LINKS
|
Vincenzo Librandi, Table of n, a(n) for n = 1..1400
|
|
|
EXAMPLE
|
a(5) = 4357 = 9*22^2 + 1.
|
|
|
MAPLE
|
A156226:=n->`if`(isprime(9*n^2+1), 9*n^2+1, NULL): seq(A156226(n), n=1..500); # Wesley Ivan Hurt, Sep 19 2014
|
|
|
MATHEMATICA
|
Union[Select[9#^2+1&/@Flatten[Table[{10m, 10m+2, 10m-2}, {m, 0, 50}]], PrimeQ]] (* Harvey P. Dale, Dec 16 2010 *)
Select[Table[9 n^2 + 1, {n, 0, 2000}], PrimeQ] (* Vincenzo Librandi, Sep 20 2014 *)
|
|
|
PROG
|
(Magma) [a: n in [0..250] | IsPrime(a) where a is 9*n^2+1]; // Vincenzo Librandi, Dec 13 2010
(PARI) for(n=1, 10^3, if(isprime(9*n^2+1), print1(9*n^2+1, ", "))) \\ Derek Orr, Sep 19 2014
|
|
|
CROSSREFS
|
Intersection of A002496 and A016777.
Sequence in context: A220938 A221292 A231514 * A217501 A133998 A231381
Adjacent sequences: A156223 A156224 A156225 * A156227 A156228 A156229
|
|
|
KEYWORD
|
nonn
|
|
|
AUTHOR
|
Avik Roy (avik_3.1416(AT)yahoo.co.in), Feb 06 2009
|
|
|
EXTENSIONS
|
Corrected, extended, comments added by Rick L. Shepherd and Zak Seidov, Feb 08 2009
|
|
|
STATUS
|
approved
|
| |
|
|
