|
%I #13 Mar 16 2017 04:19:41
%S 0,0,0,1,2,5,6,10,10,17,18,26,24,33,30,41,38,52,46,64,54,71,62,87,70,
%T 91,80,106,90,116,100,130,112,139,120,163,130,161,144,185,152,190,162,
%U 208,172,205,178,244,186,232,208,262,212,267,218,291,246,287,248,329,252
%N Sum(Omega(i)*Omega(j): i+j=n), with Omega=A001222.
%H Robert Israel, <a href="/A112967/b112967.txt">Table of n, a(n) for n = 1..10000</a>
%F G.f.: (Sum_{p prime, k>=1} x^(p^k)/(1 - x^(p^k)))^2. - _Ilya Gutkovskiy_, Mar 15 2017
%p X:= Vector(100, numtheory:-bigomega):
%p seq(add(X[i]*X[n-i],i=1..n-1),n=1..100); # _Robert Israel_, Mar 15 2017
%t Table[Sum[PrimeOmega[i] PrimeOmega[n - i],{i,1, n - 1} ], {n, 1, 61}] (* _Indranil Ghosh_, Mar 16 2017 *)
%o (PARI) for(n=1, 61, print1(sum(i=1, n - 1, bigomega(i) * bigomega(n - i)),", ")) \\ _Indranil Ghosh_, Mar 16 2017
%Y Cf. A112965, A112968, A022559.
%K nonn,look
%O 1,5
%A _Reinhard Zumkeller_, Oct 07 2005
|
