cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A228183 Semiprimes generated by the Euler polynomial x^2 + x + 41.

This page as a plain text file.
%I A228183 #24 Sep 08 2022 08:46:05
%S A228183 1681,1763,2021,2491,3233,4331,5893,6683,6847,7181,7697,8051,8413,
%T A228183 9353,10547,10961,12031,13847,14803,15047,15293,16043,16297,17071,
%U A228183 18673,19223,19781,20633,21797,24221,25481,26123,26447,26773,27101,29111,29797,30143
%N A228183 Semiprimes generated by the Euler polynomial x^2 + x + 41.
%C A228183 This is a subsequence of A145292. The first numbers in A145292 but not in here are 176861, 186233, 241613, 242597, ...
%H A228183 Shyam Sunder Gupta and K. D. Bajpai, <a href="/A228183/b228183.txt">Table of n, a(n) for n = 1..10000</a> (first 4760 terms from Shyam Sunder Gupta)
%e A228183 The semiprime 1763 = 41^2 + 41 + 41 = 41*43 is in the sequence.
%t A228183 a = {}; Do[If[PrimeOmega[x^2 + x + 41] == 2, AppendTo[a, x^2 + x + 41]], {x, 1, 200}]; a
%t A228183 (* For the b-file: *) n = 0; Do[t = k^2 + k + 41; If[PrimeOmega[t] == 2, n++; Print[n, " ", t]], {k, 30000}] (* _K. D. Bajpai_, Apr 22 2014 *)
%o A228183 (Magma) IsSemiprime:=func<i | &+[d[2]: d in Factorization(i)] eq 2>; [s: x in [2..200] | IsSemiprime(s) where s is x^2+x+41]; // _Bruno Berselli_, Aug 15 2013
%Y A228183 Cf. A001358, A145292.
%K A228183 nonn,less
%O A228183 1,1
%A A228183 _Shyam Sunder Gupta_, Aug 15 2013