OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..650
Index entries for linear recurrences with constant coefficients, signature (0, 13, 2, -57, -12, 190, -10, -453, 396, -2, -88, 308, -160, -80).
FORMULA
G.f.: (20*x^12+40*x^11 +18*x^10+52*x^9 +35*x^8-26*x^7 +34*x^6-4*x^5 -21*x^4 +2*x^3 +9*x^2-1) / (-80*x^14-160*x^13 +308*x^12-88*x^11 -2*x^10+396*x^9 -453*x^8-10*x^7 +190*x^6-12*x^5 -57*x^4+2*x^3 +13*x^2-1).
EXAMPLE
a(2) = 4, because there are 4 tilings of a 5 X 2 rectangle using right trominoes and 2 X 2 tiles:
.___. .___. .___. .___.
| . | | . | | ._| |_. |
|___| |___| |_| | | |_|
| ._| |_. | |___| |___|
|_| | | |_| | . | | . |
|___| |___| |___| |___|
MAPLE
a:= n-> (Matrix(14, (i, j)-> `if`(i=j-1, 1, `if`(i=14, [-80, -160, 308, -88, -2, 396, -453, -10, 190, -12, -57, 2, 13, 0][j], 0)))^n. <<0, 1/4, 0, 1, 0, 4, 0, 16, 0, 136, 0, 1128, 384, 8120>>)[4, 1]: seq(a(n), n=0..30);
MATHEMATICA
a[n_] := (MatrixPower[ Table[ If[i == j-1, 1, If[i == 14, {-80, -160, 308, -88, -2, 396, -453, -10, 190, -12, -57, 2, 13, 0}[[j]], 0]], {i, 1, 14}, {j, 1, 14}], n] . {0, 1/4, 0, 1, 0, 4, 0, 16, 0, 136, 0, 1128, 384, 8120})[[4]]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Dec 05 2013, translated from Alois P. Heinz's Maple program *)
CROSSREFS
KEYWORD
easy,nice,nonn
AUTHOR
Alois P. Heinz, May 18 2011
STATUS
approved