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.

A342420 a(n) = Sum_{k=1..n} (n/gcd(k,n))^(n/gcd(k,n)).

This page as a plain text file.
%I A342420 #9 Mar 11 2021 17:32:21
%S A342420 1,5,55,517,12501,93371,4941259,67109381,2324522989,40000012505,
%T A342420 2853116706111,35664401886907,3634501279107037,66672040958289359,
%U A342420 3503151123046887555,147573952589743522309,13235844190181388226833,236078448451781550068849,35611553801885644604231623
%N A342420 a(n) = Sum_{k=1..n} (n/gcd(k,n))^(n/gcd(k,n)).
%F A342420 a(n) = Sum_{d|n} phi(d^(d+1)) = Sum_{d|n} phi(d) * d^d.
%F A342420 G.f.: Sum_{k>=1} phi(k^(k+1))*x^k/(1 - x^k).
%t A342420 a[n_] := Sum[(n/GCD[k, n])^(n/GCD[k, n]), {k, 1, n}]; Array[a, 20] (* _Amiram Eldar_, Mar 11 2021 *)
%o A342420 (PARI) a(n) = sum(k=1, n, (n/gcd(k, n))^(n/gcd(k, n)));
%o A342420 (PARI) a(n) = sumdiv(n, d, eulerphi(d^(d+1)));
%o A342420 (PARI) a(n) = sumdiv(n, d, eulerphi(d)*d^d);
%o A342420 (PARI) my(N=20, x='x+O('x^N)); Vec(sum(k=1, N, eulerphi(k^(k+1))*x^k/(1-x^k)))
%Y A342420 Cf. A000010, A226459, A342411.
%K A342420 nonn
%O A342420 1,2
%A A342420 _Seiichi Manyama_, Mar 11 2021