OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..1510
P. Z. Chinn, R. Grimaldi and S. Heubach, Tiling with Ls and Squares, J. Int. Sequences 10 (2007) #07.2.8.
S. Heubach, Tiling with Ls and Squares, 2005.
Index entries for linear recurrences with constant coefficients, signature (6,5,-44,-39,2,-29,4,-4).
FORMULA
G.f.: 3x(1-x)^2(1+6x+3x^2)/(1-3x-7x^2+x^3-2x^4)^2.
EXAMPLE
a(2) = 30 because the 3 X 2 board can be tiled in one way with only square tiles, in 8 ways using one L-tile and 3 square tiles and in 2 ways with 2 L-tiles, so there are altogether 6 + 8*3 = 30 square tiles in all of the 3x2 tilings.
MAPLE
f:= gfun:-rectoproc({a(n) - 6*a(n-1)-5*a(n-2)+44*a(n-3)+39*a(n-4)-2*a(n-5)+29*a(n-6)-4*a(n-7)+4*a(n-8), a(0) = 0, a(1) = 3, a(2) = 30, a(3) = 171, a(4) = 1044, a(5) = 5691, a(6) = 30678, a(7) = 159891}, a(n), remember):
seq(f(n), n=1..40); # Robert Israel, Dec 22 2015
MATHEMATICA
Table[Coefficient[Normal[Series[3x(1-x)^2(1+6x+3x^2)/(1-3x-7x^2+x^3-2x^4)^2, {x, 0, 30}]], x, n], {n, 0, 30}]
LinearRecurrence[{6, 5, -44, -39, 2, -29, 4, -4}, {3, 30, 171, 1044, 5691, 30678, 159891, 821100}, 25] (* Vincenzo Librandi, Dec 23 2015 *)
PROG
(PARI) my(x='x+O('x^100)); Vec(3*x*(1-x)^2*(1+6*x+3*x^2)/(1-3*x-7*x^2+x^3-2*x^4)^2) \\ Altug Alkan, Dec 22 2015
(Magma) I:=[3, 30, 171, 1044, 5691, 30678, 159891, 821100]; [n le 8 select I[n] else 6*Self(n-1)+5*Self(n-2)-44*Self(n-3)-39*Self(n-4)+2*Self(n-5)-29*Self(n-6)+4*Self(n-7)-4*Self(n-8): n in [1..30]]; // Vincenzo Librandi, Dec 23 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Silvia Heubach (sheubac(AT)calstatela.edu), Feb 03 2007
STATUS
approved