A360997 Multiplicative with a(p^e) = e + 3.

1, 4, 4, 5, 4, 16, 4, 6, 5, 16, 4, 20, 4, 16, 16, 7, 4, 20, 4, 20, 16, 16, 4, 24, 5, 16, 6, 20, 4, 64, 4, 8, 16, 16, 16, 25, 4, 16, 16, 24, 4, 64, 4, 20, 20, 16, 4, 28, 5, 20, 16, 20, 4, 24, 16, 24, 16, 16, 4, 80, 4, 16, 20, 9, 16, 64, 4, 20, 16, 64, 4, 30, 4, 16

Offset: 1

Author

Vaclav Kotesovec, Feb 28 2023

Keywords

• nonn
• easy
• mult

Links

• Amiram Eldar, Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A005361 (multiplicative with a(p^e) = e), A000005 (e+1), A343443 (e+2), this sequence (e+3), A322327 (2*e), A048691 (2*e+1), A360908 (2*e-1), A226602 (3*e), A048785 (3*e+1), A360910 (3*e-1), A360909 (3*e+2), A360911 (3*e-2), A322328 (4*e), A360996 (5*e).

Cf. A000005, A007425, A007426, A074816.

Programs

Mathematica

g[p_, e_] := e+3; a[1] = 1; a[n_] := Times @@ g @@@ FactorInteger[n]; Array[a, 100]

PARI

for(n=1, 100, print1(direuler(p=2, n, (1+2*X-2*X^2)/(1-X)^2)[n], ", "))

Formula

Dirichlet g.f.: Product_{primes p} (1 + (4*p^s - 3)/(p^s - 1)^2).

Dirichlet g.f.: zeta(s)^4 * Product_{primes p} (1 - 5/p^(2*s) + 6/p^(3*s) - 2/p^(4*s)).

From Amiram Eldar, Sep 01 2023: (Start)

a(n) = A000005(A361264(n)).

a(n) = A074816(n)*A007426(n)/A007425(n). (End)