|
|
A347136
|
|
a(n) = Sum_{d|n} d * A003961(n/d), where A003961 shifts the prime factorization of its argument one step towards larger primes.
|
|
5
|
|
|
1, 5, 8, 19, 12, 40, 18, 65, 49, 60, 24, 152, 30, 90, 96, 211, 36, 245, 42, 228, 144, 120, 52, 520, 109, 150, 272, 342, 60, 480, 68, 665, 192, 180, 216, 931, 78, 210, 240, 780, 84, 720, 90, 456, 588, 260, 100, 1688, 247, 545, 288, 570, 112, 1360, 288, 1170, 336, 300, 120, 1824, 128, 340, 882, 2059, 360, 960, 138
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Dirichlet convolution of the identity function (A000027) with the prime shifted identity (A003961). Multiplicative because both A000027 and A003961 are.
Dirichlet convolution of Euler phi (A000010) with the prime shifted sigma (A003973).
Dirichlet convolution of sigma (A000203) with the prime shifted phi (A003972).
Inverse Möbius transform of A347137.
|
|
LINKS
|
|
|
FORMULA
|
For all primes p, a(p) = p + A003961(p).
|
|
MATHEMATICA
|
f[p_, e_] := ((np = NextPrime[p])^(e + 1) - p^(e + 1))/(np - p); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Aug 24 2021 *)
|
|
PROG
|
(PARI)
A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,mult
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|