cp's OEIS Frontend

A372931 a(n) = Sum_{1 <= x_1, x_2, x_3, x_4 <= n} gcd(x_1, x_2, x_3, x_4, n)^4.

%I A372931 #19 May 22 2024 01:57:47
%S A372931 1,31,161,736,1249,4991,4801,15616,19521,38719,29281,118496,57121,
%T A372931 148831,201089,311296,167041,605151,260641,919264,772961,907711,
%U A372931 559681,2514176,1170625,1770751,2106081,3533536,1414561,6233759,1847041,5963776,4714241,5178271,5996449
%N A372931 a(n) = Sum_{1 <= x_1, x_2, x_3, x_4 <= n} gcd(x_1, x_2, x_3, x_4, n)^4.
%H A372931 Amiram Eldar, <a href="/A372931/b372931.txt">Table of n, a(n) for n = 1..10000</a>
%F A372931 a(n) = Sum_{d|n} mu(n/d) * d^4 * tau(d), where mu is the Moebius function A008683.
%F A372931 From _Amiram Eldar_, May 21 2024: (Start)
%F A372931 Multiplicative with a(p^e) = (e - e/p^4 + 1) * p^(4*e).
%F A372931 Dirichlet g.f.: zeta(s-4)^2/zeta(s).
%F A372931 Sum_{k=1..n} a(k) ~ (n^5/(5*zeta(5))) * (log(n) + 2*gamma - 1/5 - zeta'(5)/zeta(5)), where gamma is Euler's constant (A001620). (End)
%t A372931 f[p_, e_] := (e - e/p^4 + 1) * p^(4*e); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, May 21 2024 *)
%o A372931 (PARI) a(n) = sumdiv(n, d, moebius(n/d)*d^4*numdiv(d));
%Y A372931 Column k=4 of A372938.
%Y A372931 Cf. A343498, A372926, A372929, A372937.
%Y A372931 Cf. A000005, A008683.
%Y A372931 Cf. A001620, A013663.
%K A372931 nonn,mult
%O A372931 1,2
%A A372931 _Seiichi Manyama_, May 17 2024