[go: up one dir, main page]

login
A349129
a(n) = Sum_{d|n} A003958(d) * A003959(n/d), where A003958 is fully multiplicative with a(p) = (p-1), and A003959 is fully multiplicative with a(p) = (p+1).
6
1, 4, 6, 13, 10, 24, 14, 40, 28, 40, 22, 78, 26, 56, 60, 121, 34, 112, 38, 130, 84, 88, 46, 240, 76, 104, 120, 182, 58, 240, 62, 364, 132, 136, 140, 364, 74, 152, 156, 400, 82, 336, 86, 286, 280, 184, 94, 726, 148, 304, 204, 338, 106, 480, 220, 560, 228, 232, 118, 780, 122, 248, 392, 1093, 260, 528, 134, 442, 276
OFFSET
1,2
COMMENTS
Dirichlet convolution of A003958 with A003959.
LINKS
FORMULA
Multiplicative with a(p^e) = ((p+1)^(e+1) - (p-1)^(e+1))/2. - Amiram Eldar, Nov 09 2021
For all n >= 1, A349130(n) <= a(n) <= A349170(n).
MATHEMATICA
f[p_, e_] := ((p + 1)^(e + 1) - (p - 1)^(e + 1))/2; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 09 2021 *)
PROG
(PARI)
A003958(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1]--); factorback(f); };
A003959(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1]++); factorback(f); };
A349129(n) = sumdiv(n, d, A003958(d)*A003959(n/d));
KEYWORD
nonn,mult
AUTHOR
Antti Karttunen, Nov 09 2021
STATUS
approved