The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 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

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified June 2 22:56 EDT 2023. Contains 363102 sequences. (Running on oeis4.)