cp's OEIS Frontend

A349712 a(n) = Sum_{d|n} sopf(d) * sopf(n/d).

%I A349712 #17 Nov 28 2021 12:35:40
%S A349712 0,0,0,4,0,12,0,8,9,20,0,32,0,28,30,12,0,42,0,48,42,44,0,52,25,52,18,
%T A349712 64,0,124,0,16,66,68,70,87,0,76,78,76,0,164,0,96,78,92,0,72,49,90,102,
%U A349712 112,0,72,110,100,114,116,0,234,0,124,102,20,130,244,0,144,138,236,0,132,0,148,110
%N A349712 a(n) = Sum_{d|n} sopf(d) * sopf(n/d).
%C A349712 Dirichlet convolution of A008472 with itself.
%H A349712 Antti Karttunen, <a href="/A349712/b349712.txt">Table of n, a(n) for n = 1..20000</a>
%F A349712 Dirichlet g.f.: ( zeta(s) * primezeta(s-1) )^2.
%F A349712 a(n) = Sum_{d|n} A061397(d) * A319131(n/d).
%F A349712 a(p) = 0 for p prime. - _Michael S. Branicky_, Nov 26 2021
%F A349712 a(p^k) = (k-1)*p^2 for p prime and k > 0. - _Chai Wah Wu_, Nov 28 2021
%t A349712 sopf[n_] := DivisorSum[n, # &, PrimeQ[#] &]; a[n_] := Sum[sopf[d] sopf[n/d], {d, Divisors[n]}]; Table[a[n], {n, 1, 75}]
%o A349712 (PARI) sopf(n) = vecsum(factor(n)[, 1]); \\ A008472
%o A349712 a(n) = sumdiv(n, d, sopf(d)*sopf(n/d)); \\ _Michel Marcus_, Nov 26 2021
%o A349712 (Python)
%o A349712 from sympy import divisors, factorint
%o A349712 def sopf(n): return sum(factorint(n))
%o A349712 def a(n): return sum(sopf(d)*sopf(n//d) for d in divisors(n))
%o A349712 print([a(n) for n in range(1, 76)]) # _Michael S. Branicky_, Nov 26 2021
%Y A349712 Cf. A008472, A034761, A061397, A070288, A319131, A349711.
%K A349712 nonn
%O A349712 1,4
%A A349712 _Ilya Gutkovskiy_, Nov 26 2021