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.

A321294 a(n) = Sum_{d|n} mu(n/d)*d*sigma_n(d).

This page as a plain text file.
%I A321294 #20 Mar 12 2021 09:30:50
%S A321294 1,9,83,1058,15629,282381,5764807,134480900,3486902505,100048836321,
%T A321294 3138428376731,107006403495850,3937376385699301,155572843119518781,
%U A321294 6568408661060858767,295150157013526773768,14063084452067724991025,708236697425777157039381
%N A321294 a(n) = Sum_{d|n} mu(n/d)*d*sigma_n(d).
%H A321294 Seiichi Manyama, <a href="/A321294/b321294.txt">Table of n, a(n) for n = 1..385</a>
%F A321294 a(n) = [x^n] Sum_{i>=1} Sum_{j>=1} mu(i)*j^(n+1)*x^(i*j)/(1 - x^(i*j))^2.
%F A321294 a(n) = Sum_{d|n} phi(n/d)*d^(n+1).
%F A321294 a(n) = Sum_{k=1..n} gcd(n,k)^(n+1).
%F A321294 a(n) ~ n^(n+1). - _Vaclav Kotesovec_, Nov 02 2018
%t A321294 Table[Sum[MoebiusMu[n/d] d DivisorSigma[n, d], {d, Divisors[n]}], {n, 18}]
%t A321294 Table[Sum[EulerPhi[n/d] d^(n + 1), {d, Divisors[n]}], {n, 18}]
%t A321294 Table[Sum[GCD[n, k]^(n + 1), {k, n}], {n, 18}]
%o A321294 (PARI) a(n) = sumdiv(n, d, moebius(n/d)*d*sigma(d, n)); \\ _Michel Marcus_, Nov 03 2018
%o A321294 (Python)
%o A321294 from sympy import totient, divisors
%o A321294 def A321294(n):
%o A321294     return sum(totient(d)*(n//d)**(n+1) for d in divisors(n,generator=True)) # _Chai Wah Wu_, Feb 15 2020
%Y A321294 Cf. A018804, A069097, A320940, A332517, A342432, A342433.
%K A321294 nonn
%O A321294 1,2
%A A321294 _Ilya Gutkovskiy_, Nov 02 2018