OFFSET
0,2
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (27, -242, 720).
FORMULA
If we define f(m,j,x)=sum(binomial(m,k)*stirling2(k,j)*x^(m-k),k=j..m) then a(n-2)=f(n,2,8), (n>=2). - Milan Janjic, Apr 26 2009
From Harvey P. Dale, Oct 25 2011: (Start)
a(0) = 1, a(1) = 27, a(2) = 487, a(n) = 27*a(n-1)-242*a(n-2)+720*a(n-3).
a(n) = 2^(2+3*(n+1))+2^(n+1)*5^(2+n)-9^(2+n). (End)
MAPLE
a:=seq(2^(2+3*(n+1))+2^(n+1)*5^(2+n)-9^(2+n), n=0..16); # Muniru A Asiru, Feb 10 2018
MATHEMATICA
CoefficientList[Series[1/((1-8x)(1-9x)(1-10x)), {x, 0, 20}], x] (* or *)
LinearRecurrence[{27, -242, 720}, {1, 27, 487}, 20] (* or *)
Table[2^(3n+2)+ 2^n 5^(n+1)- 9^(n+1), {n, 20}] (* Harvey P. Dale, Oct 25 2011 *)
PROG
(Magma) [2^(2+3*(n+1))+2^(n+1)*5^(n+2)-9^(n+2): n in [0..30]]; // Vincenzo Librandi, Oct 26 2011
(PARI) x='x+O('x^30); Vec(1/((1-8*x)*(1-9*x)*(1-10*x))) \\ G. C. Greubel, Feb 09 2018
(GAP) a:=[1, 27, 487];; for n in [4..17] do a[n]:=27*a[n-1]-242*a[n-2]+720*a[n-3]; od; # Muniru A Asiru, Feb 10 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
Typo in formula fixed by Vincenzo Librandi, Oct 26 2011
STATUS
approved