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.

A144850 a(n) = number of distinct prime divisors (taken together) of numbers of the form x^2+x+1 for x<=10^n.

This page as a plain text file.
%I A144850 #16 Jun 30 2015 11:06:26
%S A144850 8,74,734,7233,71653,712026,7090655,70686855,705173825,7038475146,
%T A144850 70278276834,701910715473
%N A144850 a(n) = number of distinct prime divisors (taken together) of numbers of the form x^2+x+1 for x<=10^n.
%H A144850 Bernhard Helmes, <a href="http://www.devalco.de/quadr_Sieb_x%5E2+x+1.php#4a">Prime sieving on the polynomial f(n)=n^2+n+1</a>.
%t A144850 d = 10; l = 0; p = 1; c = {}; a = {}; Do[k = p x^2 + x + 1; b = Divisors[k]; Do[If[PrimeQ[b[[n]]], AppendTo[a, b[[n]]]], {n, 1, Length[b]}]; If[x == d, a = Union[a]; l = Length[a]; d = 10 d; Print[l]; AppendTo[c, l]], {x, 1, 10000}]; c (*_Artur Jasinski_*)
%Y A144850 Cf. A002383, A143835, A143868, A144848, A144851.
%K A144850 nonn
%O A144850 1,1
%A A144850 _Artur Jasinski_ & Bernhard Helmes (bhelmes(AT)gmx.de), Sep 22 2008
%E A144850 Fixed broken link, corrected and extended to agree with website. - _Ray Chandler_, Jun 30 2015