|
%I #17 Dec 30 2015 05:22:09
%S 1,2,5,7,16,20,64,140,304,352,1104,1248,2640,5568,22848,25152,77760,
%T 84672,260928,542592,1126656,1209600,4921344,10174464,21012480,
%U 64364544,197074944,209018880,639000576,674832384,3409993728
%N a(n) = Product_{i=2..n} A001222(i) * Sum_{i=2..n} 1/A001222(i).
%H Robert Israel, <a href="/A067580/b067580.txt">Table of n, a(n) for n = 2..2374</a>
%e A001222(2) = 1, A001222(3) = 1, A001222(4) = 2, so a(4) = (1*1*2)*(1+1+1/2)=5. - _Robert Israel_, Dec 29 2015
%p P:= 1: S:= 0:
%p for n from 2 to 100 do
%p b:= numtheory:-bigomega(n);
%p P:= P*b; S:= S + 1/b;
%p A[n]:= P*S;
%p od:
%p seq(A[n], n=2..100); # _Robert Israel_, Dec 29 2015
%t Table[Product[PrimeOmega[k], {k,2,n}] * Sum[1/PrimeOmega[k], {k,2,n}], {n,2,50}] (* _Vaclav Kotesovec_, Dec 30 2015 *)
%Y Cf. A001222.
%K nonn
%O 2,2
%A _Benoit Cloitre_, Jan 30 2002
%E Name and offset corrected by _Robert Israel_, Dec 29 2015
|
