|
|
A084137
|
|
Binomial transform of A084136.
|
|
2
|
|
|
1, 2, 8, 32, 144, 672, 3200, 15360, 73984, 356864, 1722368, 8314880, 40144896, 193830912, 935886848, 4518838272, 21818834944, 105350561792, 508677324800, 2456111022080, 11859152338944, 57261051346944, 276480810549248
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (1-4*x)/((1-2*x)*(1-4*x-4*x^2)).
E.g.f.: exp(2*x)*cosh(sqrt(2)*x)^2 = (exp(x)*cosh(sqrt(2)*x))^2.
a(n) = ((2+sqrt(8))^n + (2-sqrt(8))^n + 2^(n+1))/4.
a(n) = 6*a(n-1) - 4*a(n-2) - 8*a(n-3). - G. C. Greubel, Oct 13 2022
|
|
MATHEMATICA
|
|
|
PROG
|
(PARI) a(n)=if(n<0, 0, polsym(4+4*x-x^2, n)[n+1]/4+2^(n-1))
(Magma)
A002203:= func< n | Round((1+Sqrt(2))^n + (1-Sqrt(2))^n) >;
(SageMath) [2^(n-2)*(2+lucas_number2(n, 2, -1)) for n in range(41)] # G. C. Greubel, Oct 13 2022
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|