login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


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 | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 28 13:24 EDT 2020. Contains 337393 sequences. (Running on oeis4.)