A144852 - OEIS

The OEIS is supported by the many generous donors to the OEIS Foundation.

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Other ways to Give

A144852

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

9, 87, 836, 8000, 78124, 766585, 7556731, 74771106, 741554656, 7366252759, 73261462211, 729280694469

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,1

COMMENTS

Primes of the form 4x^2+1 see A121326(n) = A002496(n+1).

LINKS

Table of n, a(n) for n=1..12.

Bernhard Helmes, Prime sieving on the polynomial f(n)=4n^2+1.

MATHEMATICA

d = 10; l = 0; p = 4; 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

Cf. A002383, A121326, A143835, A143868, A144848, A144850, A144851.

Sequence in context: A100814 A055725 A219123 * A153191 A223277 A267265

Adjacent sequences: A144849 A144850 A144851 * A144853 A144854 A144855

KEYWORD

nonn

AUTHOR

Artur Jasinski & Bernhard Helmes (bhelmes(AT)gmx.de), Sep 22 2008

EXTENSIONS

Fixed broken link, corrected and extended to agree with website. - Ray Chandler, Jun 30 2015

STATUS

approved