|
|
A064549
|
|
a(n) = n * Product_{primes p|n} p.
|
|
46
|
|
|
1, 4, 9, 8, 25, 36, 49, 16, 27, 100, 121, 72, 169, 196, 225, 32, 289, 108, 361, 200, 441, 484, 529, 144, 125, 676, 81, 392, 841, 900, 961, 64, 1089, 1156, 1225, 216, 1369, 1444, 1521, 400, 1681, 1764, 1849, 968, 675, 2116, 2209, 288, 343, 500, 2601, 1352
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
Multiplicative with a(p^k)=p^(k+1) when k>0.
G.f.: Sum_{k>=1} mu(k)^2*phi(k)*k*x^k/(1 - x^k)^2. - Ilya Gutkovskiy, Nov 02 2018
Dirichlet g.f.: zeta(s-2) * zeta(s-1) * Product_{primes p} (1 + p^(3-2*s) - p^(4-2*s) - p^(1-s)).
Sum_{k=1..n} a(k) ~ c * n^3, where c = A065463/3 = A065464*Pi^2/18 = 0.234814...
(End)
Sum_{k>=1} (-1)^(k+1)/a(k) = zeta(2)*zeta(3)/(3*zeta(6)) = (1/3) * A082695. - Amiram Eldar, Nov 18 2020
|
|
EXAMPLE
|
a(12) = 72 since 12 = 2^2*3 and 12*2*3 = 72.
|
|
MAPLE
|
a:= n -> n * convert(numtheory:-factorset(n), `*`):
|
|
MATHEMATICA
|
Table[n*Product[If[PrimeQ[d], d, 1], {d, Divisors[n]}], {n, 1, 100}] (* Vaclav Kotesovec, Jun 15 2019 *)
|
|
PROG
|
(PARI) popf(n)= { local(f, p=1); f=factor(n); for(i=1, matsize(f)[1], p*=f[i, 1]); return(p) } { for (n=1, 1000, write("b064549.txt", n, " ", n*popf(n)) ) } \\ Harry J. Smith, Sep 18 2009
(PARI) A064549(n) = { my(f=factor(n)); for (i=1, #f~, f[i, 2]++); factorback(f); }; \\ Antti Karttunen, Aug 30 2018
(PARI) for(n=1, 100, print1(direuler(p=2, n, (1 - p*X + p^2*X)/(1 - p*X))[n], ", ")) \\ Vaclav Kotesovec, Jun 24 2020
(Haskell)
(Magma) [n^2/( (&+[Floor(k^n/n)-Floor((k^n - 1)/n) : k in [1..n]]) ): n in [1..50]]; // G. C. Greubel, Nov 02 2018
|
|
CROSSREFS
|
A permutation of the powerful numbers A001694.
|
|
KEYWORD
|
mult,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|