cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A344195 a(n) = Sum_{k=1..n} tau(gcd(k,n))^(n/gcd(k,n)), where tau(n) is the number of divisors of n.

This page as a plain text file.
%I A344195 #18 May 11 2021 10:35:16
%S A344195 1,3,4,9,6,26,8,49,25,140,12,240,14,782,156,1215,18,3349,20,5130,800,
%T A344195 20498,24,19558,151,98324,3148,111492,30,270624,32,551091,20520,
%U A344195 2097176,924,1716189,38,9437210,98348,8630496,42,25362724,44,43714584,266346,184549406,48,137141048,813,671096867
%N A344195 a(n) = Sum_{k=1..n} tau(gcd(k,n))^(n/gcd(k,n)), where tau(n) is the number of divisors of n.
%F A344195 a(n) = Sum_{d|n} phi(n/d) * tau(d)^(n/d).
%F A344195 If p is prime, a(p) = 1 + p.
%t A344195 a[n_] := DivisorSum[n, EulerPhi[n/#] * DivisorSigma[0, #]^(n/#) &]; Array[a, 50] (* _Amiram Eldar_, May 11 2021 *)
%o A344195 (PARI) a(n) = sum(k=1, n, numdiv(gcd(k, n))^(n/gcd(k, n)));
%o A344195 (PARI) a(n) = sumdiv(n, d, eulerphi(n/d)*numdiv(d)^(n/d));
%Y A344195 Cf. A000005, A279789, A342424, A344140, A344194.
%K A344195 nonn
%O A344195 1,2
%A A344195 _Seiichi Manyama_, May 11 2021