login
Number of distinct primes dividing 2n (i.e., omega(2n)).
4

%I #21 Sep 21 2024 08:42:04

%S 1,1,2,1,2,2,2,1,2,2,2,2,2,2,3,1,2,2,2,2,3,2,2,2,2,2,2,2,2,3,2,1,3,2,

%T 3,2,2,2,3,2,2,3,2,2,3,2,2,2,2,2,3,2,2,2,3,2,3,2,2,3,2,2,3,1,3,3,2,2,

%U 3,3,2,2,2,2,3,2,3,3,2,2,2,2,2,3,3,2,3,2,2,3,3,2,3,2,3,2,2,2,3,2,2,3,2,2,4

%N Number of distinct primes dividing 2n (i.e., omega(2n)).

%C Bisection of A001221.

%H G. C. Greubel, <a href="/A099812/b099812.txt">Table of n, a(n) for n = 1..1000</a>

%F From _Amiram Eldar_, Sep 21 2024: (Start)

%F a(n) = A001221(2*n).

%F a(n) = omega(n) + 1 if n is odd, and a(n) = omega(n) if n is even.

%F Sum_{k=1..n} a(k) = n * (log(log(n)) + B + 1/2) + O(n/log(n)), where B is Mertens's constant (A077761). (End)

%e a(6) = 2 because 12 = 2*2*3 has 2 distinct prime divisors.

%e a(15) = 3 because 30 = 2*3*5 has 3 distinct prime divisors.

%p with(numtheory): omega:=proc(n) local div,A,j: div:=divisors(n): A:={}: for j from 1 to tau(n) do if isprime(div[j])=true then A:=A union {div[j]} else A:=A fi od: nops(A) end: seq(omega(2*n),n=1..130); # _Emeric Deutsch_, Mar 10 2005

%t Table[PrimeNu[2*n], {n,1,50}] (* _G. C. Greubel_, May 21 2017 *)

%o (PARI) for(n=1,50, print1(omega(2*n), ", ")) \\ _G. C. Greubel_, May 21 2017

%o (Magma) [#PrimeDivisors(2*n): n in [1..100]]; // _Vincenzo Librandi_, Jul 26 2017

%Y Cf. A001221, A077761, A092523.

%K nonn,easy

%O 1,3

%A _N. J. A. Sloane_, Nov 19 2004

%E More terms from _Emeric Deutsch_, Mar 10 2005