|
%I #20 Sep 04 2023 02:01:39
%S 1,3,4,7,6,12,8,15,13,18,12,28,16,24,24,31,23,39,31,42,32,36,43,60,31,
%T 48,40,56,57,72,78,63,48,69,48,91,102,93,64,90,136,96,177,84,78,129,
%U 232,124,57,93,92,112,298,120,72,120,124,171,386,168,491,234,104,127
%N Inverse Möbius transform of A003964.
%H N. J. A. Sloane, <a href="/transforms.txt">Transforms</a>.
%F Multiplicative with a(p^e) = (q^(e+1)-1)/(q-1) where q = partition(pi(p)+1) = A000041(A000720(p)+1). - _David W. Wilson_, Sep 01 2001
%t f[p_, e_] := Module[{q = PartitionsP[PrimePi[p] + 1]}, (q^(e+1)-1)/(q-1)]; a[n_] := Times @@ f @@@ FactorInteger[n]; a[1] = 1; Array[a, 100] (* _Amiram Eldar_, Sep 04 2023 *)
%Y Cf. A000041, A000720, A003964, A003978.
%K nonn,easy,mult
%O 1,2
%A _Marc LeBrun_
%E More terms from _Vladeta Jovovic_, Dec 01 2001
|
