A090685 Primes of the form 8*k^2 + 1.

73, 1153, 1801, 2593, 3529, 8713, 10369, 18433, 20809, 34849, 48673, 52489, 69193, 83233, 103969, 115201, 139393, 180001, 187273, 209953, 259201, 267913, 362953, 405001, 415873, 426889, 438049, 449353, 472393, 484129, 508033, 622729, 636193

Offset

1,1

Links

Daniel Starodubtsev, Table of n, a(n) for n = 1..10000

Mathematica

Select[8 Range[300]^2+1, PrimeQ]  (* Harvey P. Dale, Apr 02 2011 *)

Prog

(PARI) A090685(n) = { for(x=1, n, y = 8*x^2+1; if(isprime(y), print1(y", ")) ) }
(Magma) [ a: n in [0..350] | IsPrime(a) where a is 8*n^2+1 ]; // Vincenzo Librandi, Apr 04 2011

Crossrefs

Cf. A188549.

Sequence in context: A320214 A346642 A008400 * A232297 A008392 A221479

Adjacent sequences: A090682 A090683 A090684 * A090686 A090687 A090688

Keyword

easy,nonn

Author

Cino Hilliard, Dec 18 2003

Status

approved