OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000
FORMULA
G.f.: (x^9 - x^8 + x^7 - 4*x^6 + x^5 - 3*x^4 - x^3 - 2*x^2 - 1) / (x^19 - x^18 - x^16 + 2*x^12 + x^10 - x^9 - 2*x^7 + x^3 + x - 1).
EXAMPLE
a(4) = 9 because there are 9 ways to divide a 4 X 5 rectangle into subsquares, considering only the list of parts: [20(1 X 1)], [16(1 X 1), 1(2 X 2)], [12(1 X 1), 2(2 X 2)], [11(1 X 1), 1(3 X 3)], [8(1 X 1), 3(2 X 2)], [7(1 X 1), 1(2 X 2), 1(3 X 3)], [4(1 X 1), 4(2 X 2)], [4(1 X 1), 1(4 X 4)], [3(1 X 1), 2(2 X 2), 1(3 X 3)]. There is no way to divide this rectangle into [2(1 X 1), 2(3 X 3)].
MAPLE
gf:= (x^9-x^8+x^7-4*x^6+x^5-3*x^4-x^3-2*x^2-1)/
(x^19-x^18-x^16+2*x^12+x^10-x^9-2*x^7+x^3+x-1):
a:= n-> coeff(series(gf, x, n+1), x, n):
seq(a(n), n=0..60);
PROG
(Magma) m:=50; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!((1+2*x^2+x^3+3*x^4-x^5+4*x^6-x^7+x^8-x^9)/((1-x)^5*(1+x)^2*(1+x^2)*(1-x +x^2)*(1+x+x^2)^2*(1+x+x^2+x^3+x^4)))); // Bruno Berselli, Apr 17 2013
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Alois P. Heinz, Apr 17 2013
STATUS
approved