OFFSET
1,2
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..200 (computed by Anthony J. Guttmann)
Nicholas R. Beaton, Filippo Disanto, Anthony J. Guttmann and Simone Rinaldi, On the enumeration of column-convex permutominoes, in FPSAC 2011, Reykjavik, Iceland DMTCS proc. AO, 2011, 111-122.
MATHEMATICA
cecka = {}; Do[fxyz = Sum[If[h + d + w <= n, Subscript[c, h, d, w]*x^h*y^d*z^w, 0], {h, 0, n}, {d, 0, n}, {w, 0, n}]; fxyz = Simplify[fxyz /. cecka]; feq = Expand[Together[(x - y)*(z - y)*(y + y*z/(z - y)*(fxyz /. y -> z) - y^2/(z - y)*(fxyz /. z -> y) + x*y/(x - y)*(fxyz /. y -> x) - y^2/(x - y)*(fxyz /. x -> y))]]; fxyz1 = Expand[Simplify[(x - y)*(z - y)*fxyz]]; fxyz2 = 0; Do[p = fxyz1[[j]]; If[Total[Exponent[p, {x, y, z}]] <= n + 2, fxyz2 += p], {j, 1, Length[fxyz1]}]; poly2 = 0; Do[p = feq[[j]]; If[Total[Exponent[p, {x, y, z}]] <= n + 2, poly2 += p], {j, 1, Length[feq]}]; ceckanew = SolveAlways[fxyz2 == poly2, {x, y, z}][[1]]; cecka = Join[cecka, ceckanew]; gfxyz = fxyz /. cecka; , {n, 1, 15}]; Rest[CoefficientList[gfxyz /. {x -> y, z -> y}, y]] (* Vaclav Kotesovec, Jun 07 2021 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Sep 30 2011
EXTENSIONS
More terms from Vaclav Kotesovec, Jun 07 2021
STATUS
approved