A128428 Number of distinct prime factors of n^2+1.

1, 1, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 3, 1, 2, 1, 3, 2, 2, 1, 3, 2, 3, 1, 2, 1, 3, 2, 2, 2, 3, 2, 3, 2, 2, 1, 3, 2, 2, 1, 2, 2, 3, 2, 2, 2, 4, 2, 2, 2, 2, 2, 3, 1, 3, 1, 3, 2, 2, 2, 2, 2, 3, 2, 2, 1, 3, 2, 2, 2, 2, 3, 4, 1, 3, 2, 3, 2, 2, 2, 3, 2, 4, 1, 2, 2, 3, 2, 3, 1, 3, 2, 3, 1, 2, 2, 3, 3, 3, 2, 2, 2

Offset

1,3

Comments

a(n) is also the number of distinct prime factors, up to multiplication by units, of n + i in the ring of Gaussian integers. - Jason Kimberley, Dec 17 2011

Links

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000

Formula

a(n) = A001221(n^2+1).

Example

a(3) = 2 because 3^2+1 = 2*5.

Maple

a:= n-> nops(select(isprime, numtheory[divisors](n^2+1))):
seq(a(n), n=1..100);  # Alois P. Heinz, Dec 06 2020

Mathematica

a[n_]:=Length[FactorInteger[n^2 + 1]]

Prog

(PARI) a(n)=omega(n^2+1) \\ Charles R Greathouse IV, Jul 31 2011

Crossrefs

Cf. A193330 (counted with multiplicity).

Sequence in context: A329320 A316112 A317994 * A056171 A357328 A333749

Adjacent sequences: A128425 A128426 A128427 * A128429 A128430 A128431

Keyword

nonn

Author

Kent Horvath (kenthorvath(AT)gmail.com), May 10 2007

Status

approved