login
A304117
If n = Product (p_j^k_j) then a(n) = Product (pi(p_j)*k_j), where pi() = A000720.
4
1, 1, 2, 2, 3, 2, 4, 3, 4, 3, 5, 4, 6, 4, 6, 4, 7, 4, 8, 6, 8, 5, 9, 6, 6, 6, 6, 8, 10, 6, 11, 5, 10, 7, 12, 8, 12, 8, 12, 9, 13, 8, 14, 10, 12, 9, 15, 8, 8, 6, 14, 12, 16, 6, 15, 12, 16, 10, 17, 12, 18, 11, 16, 6, 18, 10, 19, 14, 18, 12, 20, 12, 21, 12, 12, 16, 20, 12, 22, 12
OFFSET
1,3
FORMULA
a(n) = A005361(n)*A156061(n).
a(p^k) = A000720(p)*k where p is a prime.
a(A002110(m)^k) = k^m*m!.
As an example:
a(A000040(k)) = k.
a(A006450(k)) = A000040(k).
a(A001248(k)) = a(A031215(k)) = A005843(k).
a(A030078(k)) = a(A031336(k)) = A008585(k)
a(A061742(k)) = A000165(k).
a(A115964(k)) = A032031(k).
a(A002110(k)) = A000142(k).
a(A080696(k)) = A002110(k).
EXAMPLE
a(36) = 8 because 36 = 2^2*3^2 = prime(1)^2*prime(2)^2 and 1*2*2*2 = 8.
MATHEMATICA
a[n_] := Times @@ (PrimePi[#[[1]]] #[[2]] & /@ FactorInteger[n]); a[1] = 1; Table[a[n], {n, 1, 80}]
PROG
(PARI) a(n) = my(f=factor(n)); for (k=1, #f~, f[k, 1] = primepi(f[k, 1])*f[k, 2]; f[k, 2] = 1); factorback(f); \\ Michel Marcus, May 06 2018
KEYWORD
nonn,mult
AUTHOR
Ilya Gutkovskiy, May 06 2018
STATUS
approved