Tiling n X m rectangles with 1 X 1 and s X s squares
  We consider tilings of a rectangle which is n units wide and m units long by non-overlapping 1 X 1 squares and s X s squares. Bivariate generating functions are computed with the Transfer Matrix Method for moderately large but fixed widths n as a function of the parameter m and of the number of s X s squares in the rectangle. 
Richard J. Mathar
math.CO