The OEIS is supported by the many generous donors to the OEIS Foundation. Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A144255 Semiprimes of the form k^2+1. 17
 10, 26, 65, 82, 122, 145, 226, 362, 485, 626, 785, 842, 901, 1157, 1226, 1522, 1765, 1937, 2026, 2117, 2305, 2402, 2501, 2602, 2705, 3365, 3482, 3601, 3722, 3845, 4097, 4226, 4762, 5042, 5777, 6085, 6242, 6401, 7226, 7397, 7745, 8465, 9026, 9217, 10001, 10202 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 COMMENTS Iwaniec proves that there are an infinite number of semiprimes or primes of the form n^2+1. Because n^2+1 is not a square for n>0, all such semiprimes have two distinct prime factors. Moreover, this implies that one prime factor p of n^2+1 is strictly smaller than n, and therefore also divisor of (the usually much smaller) m^2+1, where m = n % p (binary "mod" operation). - M. F. Hasler, Mar 11 2012 LINKS Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 Henryk Iwaniec, Almost-primes represented by quadratic polynomials, Inventiones Mathematicae 47 (2) (1978), pp. 171-188. FORMULA a(n) = A085722(n)^2 + 1. Equals { n^2+1 | A193432(n)=2 }. - M. F. Hasler, Mar 11 2012 MATHEMATICA Select[Table[n^2  + 1, {n, 100}], PrimeOmega[#] == 2&] (* Vincenzo Librandi, Sep 22 2012 *) PROG (PARI) select(n->bigomega(n)==2, vector(500, n, n^2+1)) \\ Zak Seidov Feb 24 2011 (Magma) IsSemiprime:= func; [s: n in [1..100] | IsSemiprime(s) where s is n^2 + 1]; // Vincenzo Librandi, Sep 22 2012 (Python) from sympy import primeomega from itertools import count, takewhile def aupto(limit):     form = takewhile(lambda x: x <= limit, (k**2+1 for k in count(1)))     return [number for number in form if primeomega(number)==2] print(aupto(10202)) # Michael S. Branicky, Oct 26 2021 CROSSREFS Cf. A001358, A085722, A069987, A193432. Subsequence of A134406. Sequence in context: A055710 A332596 A134420 * A259290 A072379 A005970 Adjacent sequences:  A144252 A144253 A144254 * A144256 A144257 A144258 KEYWORD nonn,easy AUTHOR T. D. Noe, Sep 16 2008 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 September 28 01:47 EDT 2022. Contains 357063 sequences. (Running on oeis4.)