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.
%I A350125 #23 Oct 04 2023 03:54:49 %S A350125 1,8,40,345,3303,50225,833569,17045934,388654659,10039636255, %T A350125 285508661853,8924967326015,302927979357701,11114722212099135, %U A350125 437913155876193839,18447871416712820782,827249276230172525622,39347009369000530723017 %N A350125 a(n) = Sum_{k=1..n} k^2 * floor(n/k)^n. %H A350125 Seiichi Manyama, <a href="/A350125/b350125.txt">Table of n, a(n) for n = 1..386</a> %F A350125 a(n) = Sum_{k=1..n} k^2 * Sum_{d|k} (d^n - (d - 1)^n)/d^2. %F A350125 a(n) = [x^n] (1/(1 - x)) * Sum_{k>=1} (k^n - (k - 1)^n) * x^k * (1 + x^k)/(1 - x^k)^3. %F A350125 a(n) ~ n^n. - _Vaclav Kotesovec_, Dec 16 2021 %t A350125 a[n_] := Sum[k^2 * Floor[n/k]^n, {k, 1, n}]; Array[a, 18] (* _Amiram Eldar_, Oct 04 2023 *) %o A350125 (PARI) a(n) = sum(k=1, n, k^2*(n\k)^n); %o A350125 (PARI) a(n) = sum(k=1, n, k^2*sumdiv(k, d, (d^n-(d-1)^n)/d^2)); %Y A350125 Cf. A332469, A350109, A350123, A350124, A350128. %K A350125 nonn %O A350125 1,2 %A A350125 _Seiichi Manyama_, Dec 15 2021