A143868 a(n)= Number of distinct prime divisors (taken together) of numbers of the form 2x^2-1 for x<=10^n.
8, 84, 815, 7922, 77250, 759077, 7492588, 74198995, 736401956, 7319543972, 72834161468, 725344237597
Offset: 1
Keywords
Links
- Bernhard Helmes, Prime sieving on the polynomial f(n)=2n^2-1.
Crossrefs
Cf. A143835.
Programs
-
Mathematica
l = 0; p = 2; c = {}; a = {}; Do[k = p x^2 - 1; b = Divisors[k]; Do[If[PrimeQ[b[[n]]], AppendTo[a, b[[n]]]], {n, 1, Length[b]}]; If[N[Log[x]/Log[10]] == Round[N[Log[x]/Log[10]]], a = Union[a]; l = Length[a]; Print[l]; AppendTo[c, l]], {x, 1, 10000}]; c (*Artur Jasinski*)
Extensions
Fixed broken link and corrected terms to agree with website - Ray Chandler, Jun 30 2015