|
|
A289213
|
|
a(n) = n! * Laguerre(n,-7).
|
|
3
|
|
|
1, 8, 79, 916, 12113, 179152, 2921911, 51988748, 1000578817, 20686611736, 456805020959, 10721879413252, 266382974861521, 6980304560060384, 192311632290456007, 5555079068684580988, 167822887344661475969, 5290815252203206305832, 173713426149927498289903
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
FORMULA
|
E.g.f.: exp(7*x/(1-x))/(1-x).
a(n) = n! * Sum_{i=0..n} 7^i/i! * binomial(n,i).
a(n) ~ exp(-7/2 + 2*sqrt(7*n) - n) * n^(n + 1/4) / (sqrt(2)*7^(1/4)) * (1 + 367/(48*sqrt(7*n))). - Vaclav Kotesovec, Nov 13 2017
Sum_{n>=0} a(n) * x^n / (n!)^2 = exp(x) * Sum_{n>=0} 7^n * x^n / (n!)^2. - Ilya Gutkovskiy, Jul 17 2020
|
|
MAPLE
|
a:= n-> n! * add(binomial(n, i)*7^i/i!, i=0..n):
seq(a(n), n=0..20);
|
|
MATHEMATICA
|
Table[n!*LaguerreL[n, -7], {n, 0, 20}] (* Indranil Ghosh, Jul 04 2017 *)
|
|
PROG
|
(Python)
from mpmath import *
mp.dps=100
def a(n): return int(fac(n)*laguerre(n, 0, -7))
(PARI) x = 'x + O('x^30); Vec(serlaplace(exp(7*x/(1-x))/(1-x))) \\ Michel Marcus, Jul 04 2017
(Magma) m:=25; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(Exp(7*x/(1-x))/(1-x))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, May 13 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|