|
Search: id:A144848
|
|
|
| A144848 |
|
a(n) = number of distinct prime divisors (taken together) of numbers of the form x^2+1 for x<=10^n. |
|
+0 5
|
|
| 7, 70, 720, 7102, 70780, 704537, 7026559, 70122424, 700184485, 6993568566, 69870544960
(list; graph; listen)
|
|
|
OFFSET
|
1,1
|
|
|
LINKS
|
B. Helmes, Quadratisches Siebverfahren/Primzahlgenerator
|
|
MATHEMATICA
|
d = 10; l = 0; p = 1; 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[x == d, a = Union[a]; l = Length[a]; d = 10 d; Print[l]; AppendTo[c, l]], {x, 1, 10000}]; c (*Artur Jasinski*)
|
|
CROSSREFS
|
For primes of the form n^2+1 see A002496
Cf. A002496, A143835, A143868, A144848, A144849, A144850, A144851
Sequence in context: A063416 A043034 A015251 this_sequence A078246 A097184 A002296
Adjacent sequences: A144845 A144846 A144847 this_sequence A144849 A144850 A144851
|
|
KEYWORD
|
base,nonn
|
|
AUTHOR
|
Artur Jasinski & Bernhard Helmes (bhelmes(AT)gmx.de) (grafix(AT)csl.pl), Sep 22 2008, Sep 24 2008
|
|
|
Search completed in 0.002 seconds
|