login
A346845
E.g.f.: log(1 + x) / (1 - x)^3.
3
1, 5, 29, 186, 1374, 11352, 105048, 1070640, 11978640, 145558080, 1914027840, 27035890560, 408891369600, 6585851059200, 112656894336000, 2038285492992000, 38915729475840000, 781515776369664000, 16475855040820224000, 363685261902133248000, 8391522945839007744000
OFFSET
1,2
LINKS
FORMULA
a(n) = n! * Sum_{k=1..n} (-1)^(k+1) * binomial(n-k+2,2) / k.
a(n) ~ log(2) * n^2 * n! / 2. - Vaclav Kotesovec, Aug 06 2021
D-finite with recurrence: -(n + 1)*(3 + n)^2*a(n) + (n^2 + 3*n - 1)*a(n + 1) + (6 + n)*a(n + 2) - a(3 + n) = 0. - Robert Israel, Nov 18 2025
MAPLE
f:= proc(m) option remember; -(m - 2)*m^2*procname(m - 3) + (m^2 - 3*m - 1)*procname(m - 2) + (m + 3)*procname(m - 1) end proc:
f(0):= 0: f(1):= 1: f(2):= 5:
map(f, [$0..50]); # Robert Israel, Nov 18 2025
MATHEMATICA
nmax = 21; CoefficientList[Series[Log[1 + x]/(1 - x)^3, {x, 0, nmax}], x] Range[0, nmax]! // Rest
Table[n! Sum[(-1)^(k + 1) Binomial[n - k + 2, 2]/k , {k, 1, n}], {n, 1, 21}]
Table[n!*(((-1)^n*(2*n + 5) - 4*n - 5)/8 + (n+1)*(n+2)*(Log[2] - (-1)^n * LerchPhi[-1, 1, 1 + n])/2), {n, 1, 21}] // Simplify (* Vaclav Kotesovec, Aug 06 2021 *)
PROG
(PARI) my(x='x+O('x^25)); Vec(serlaplace(log(1+x)/(1-x)^3)) \\ Michel Marcus, Aug 06 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Aug 05 2021
STATUS
approved