|
|
|
|
0, 1, 2, 16, 48, 256, 928, 4288, 16896, 73728, 301056, 1283072, 5318656, 22446080, 93634560, 393543680, 1646034944, 6906380288, 28918677504, 121248940032, 507934408704, 2129004593152, 8920531730432, 37385660268544, 156658714017792
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,3
|
|
LINKS
|
|
|
FORMULA
|
T(n, k) = [x^k] ((x+sqrt(x+5))^n - (x-sqrt(x+5))^n)/(2*sqrt(x+5)).
a(n) = (1/2)*Sum_{k=0..n-1} T(n, k)*Fibonacci(k-1).
a(n) = 2*a(n-1) + 12*a(n-2) - 8*a(n-3) - 16*a(n-4).
G.f.: x^2/(1-2*x-12*x^2+8*x^3+16*x^4).
|
|
EXAMPLE
|
The first five polynomials p(n,x) and their reductions are as follows:
p(0,x) = 1 -> 1
p(1,x) = 2*x -> 2*x
p(2,x) = 3 + x + 3*x^2 -> 8 + 4*x
p(3,x) = 12*x + 4*x^2 + 4*x^3 -> 8 + 32*x
p(4,x) = 9 + 6*x + 31*x^2 + 10*x^3 + 5*x^4 -> 96 + 96*x.
From these, read a(n) = (0, 2, 4, 32, 96, ...)/2 = (0, 1, 2, 16, 48, ...).
|
|
MATHEMATICA
|
LinearRecurrence[{2, 12, -8, -16}, {0, 1, 2, 16}, 40] (* G. C. Greubel, Jul 10 2023 *)
|
|
PROG
|
(Magma)
R<x>:=PowerSeriesRing(Integers(), 41);
[0] cat Coefficients(R!( x^2/(1-2*x-12*x^2+8*x^3+16*x^4) )); // G. C. Greubel, Jul 10 2023
(SageMath)
@CachedFunction
if (n<5): return (0, 0, 1, 2, 16)[n]
else: return 2*a(n-1) +12*a(n-2) -8*a(n-3) -16*a(n-4)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|