%I #17 Sep 04 2023 01:59:25
%S 1,1,2,2,3,2,5,4,6,3,6,4,8,5,6,8,9,6,11,6,10,6,14,8,12,8,18,10,15,6,
%T 18,16,12,9,15,12,20,11,16,12,21,10,23,12,18,14,26,16,30,12,18,16,29,
%U 18,18,20,22,15,30,12,33,18,30,32,24,12,35,18,28,15,36,24,39,20,24,22,30
%N Möbius transform of A003962.
%F Multiplicative with a(p^e) = ((q-1)/2)((q+1)/2)^(e-1) where q = nextPrime(p). - _David W. Wilson_, Sep 01 2001
%t f[p_, e_] := Module[{q = NextPrime[p]}, ((q-1)/2)((q+1)/2)^(e-1)]; a[n_] := Times @@ f @@@ FactorInteger[n]; a[1] = 1; Array[a, 100] (* _Amiram Eldar_, Sep 04 2023 *)
%Y Cf. A003962, A003962, A003975, A151800.
%K nonn,easy,mult
%O 1,3
%A _Marc LeBrun_
%E Corrected and extended by _David W. Wilson_, Aug 29 2001
|