OFFSET
0,3
COMMENTS
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..140, flattened
G. Castiglione, A. Frosini, E. Munarini, A. Restivo and S. Rinaldi, Combinatorial aspects of L-convex polyominoes, European J. Combin. 28 (2007), no. 6, 1724-1741.
FORMULA
G.f.: G(t,z) = (1-z)^2*(1+z)^2/(1-5z^2+2z^4-2tz).
The g.f. of column k is (2z)^k*(1-z^2)^2/(1-5z^2+2z^4)^{k+1} (we have a Riordan array).
The g.f. H(t,s,z), where z marks size and t (s) marks number of columns with an odd (even) sum, is H=(1-z^2)^2/(1-2z^2+z^4-2tz-3sz^2+sz^4).
EXAMPLE
T(2,2) = 4 because we have (1,0/0,1), (0,1/1,0), (1,1/0,0), and (0,0/1,1) (the 2-compositions are written as (top row/bottom row)).
Triangle starts:
1;
0, 2;
3, 0, 4;
0, 16, 0, 8;
14, 0, 52, 0, 16;
MAPLE
G := (1-z^2)^2/(1-5*z^2+2*z^4-2*t*z): Gser := simplify(series(G, z = 0, 15)): for n from 0 to 11 do P[n] := sort(coeff(Gser, z, n)) end do; for n from 0 to 11 do seq(coeff(P[n], t, k), k = 0 .. n) end do; # yields sequence in triangular form
# second Maple program:
b:= proc(n) option remember; `if`(n=0, 1,
expand(add(add(`if`(i=0 and j=0, 0, b(n-i-j)*
`if`(irem(i+j, 2)=1, x, 1)), i=0..n-j), j=0..n)))
end:
T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(n)):
seq(T(n), n=0..15); # Alois P. Heinz, Mar 16 2014
MATHEMATICA
b[n_] := b[n] = If[n == 0, 1, Expand[Sum[Sum[If[i == 0 && j == 0, 0, b[n-i-j]* If[Mod[i+j, 2] == 1, x, 1]], {i, 0, n-j}], {j, 0, n}]]]; T[n_] := Function[{p}, Table[Coefficient[p, x, i], {i, 0, Exponent[p, x]}]][b[n]]; Table[T[n], {n, 0, 15}] // Flatten (* Jean-François Alcover, Feb 19 2015, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Emeric Deutsch, Oct 13 2010
STATUS
approved