A356131 - OEIS

A356131

a(n) = Sum_{k=1..n} (k - 1)^n * binomial(floor(n/k)+1,2).

0, 1, 9, 100, 1302, 20648, 377022, 7921039, 186926431, 4916562309, 142373072781, 4506381442625, 154721361953489, 5729251983077521, 227585590018322461, 9654855432715969784, 435659531345223039702, 20836069677785611552293

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,3

LINKS

Table of n, a(n) for n=1..18.

FORMULA

a(n) = Sum_{k=1..n} k * (sigma_{n-1}(k) - floor(n/k)^n) = A356129(n) - A350109(n).

a(n) = Sum_{k=1..n} k * Sum_{d|k} (d - 1)^n / d.

a(n) = [x^n] (1/(1-x)) * Sum_{k>=1} (k - 1)^n * x^k/(1 - x^k)^2.

MATHEMATICA

a[n_] := Sum[(k - 1)^n * Binomial[Floor[n/k]+1, 2], {k, 1, n}]; Array[a, 18] (* Amiram Eldar, Jul 28 2022 *)

PROG

(PARI) a(n) = sum(k=1, n, (k-1)^n*binomial((n\k)+1, 2));

(PARI) a(n) = sum(k=1, n, k*(sigma(k, n-1)-(n\k)^n));

(PARI) a(n) = sum(k=1, n, k*sumdiv(k, d, (d-1)^n/d));

CROSSREFS

Cf. A350109, A356100, A356129.

Sequence in context: A056002 A060150 A202833 * A287039 A360348 A174509

Adjacent sequences: A356128 A356129 A356130 * A356132 A356133 A356134

KEYWORD

nonn

AUTHOR

Seiichi Manyama, Jul 27 2022

STATUS

approved