A304662 - OEIS

%I #37 Aug 29 2021 12:01:47

%S 1,2,6,16,42,106,268,650,1580,3750,8862,20598,47776,109248,248966,

%T 562630,1264780,2823958,6282198,13884820,30590124,67051982,146463790,

%U 318588916,690882926,1492592450,3215372064,6904561416,14786529836,31574656096,67261524262

%N Total number of domino tilings of Ferrers-Young diagrams summed over all partitions of 2n.

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/FerrersDiagram.html">Ferrers Diagram</a>

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Domino_(mathematics)">Domino</a>

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Domino_tiling">Domino tiling</a>

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Ferrers_diagram">Ferrers diagram</a>

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Mutilated_chessboard_problem">Mutilated chessboard problem</a>

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Partition_(number_theory)">Partition (number theory)</a>

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Young_tableau#Diagrams">Young tableau, Diagrams</a>

%H Gus Wiseman, <a href="/A304662/a304662.png">All 42 domino tilings of integer partitions of 8.</a>

%H Gus Wiseman, <a href="/A304662/a304662_1.png">All 106 domino tilings of integer partitions of 10.</a>

%H <a href="/index/Do#domino">Index entries for sequences related to dominoes</a>

%F a(n) = Sum_{k=0..A304790(n)} k * A304789(n,k).

%F a(n) = Sum_{k=0..n} A304718(n,k).

%F a(n) = A296625(n) for n < 7.

%e a(2) = 6:

%e ._.  .___.  ._._.  .___.  ._.___.  .___.___.

%e | |  |___|  | | |  |___|  | |___|  |___|___|

%e |_|  | |    |_|_|  |___|  |_|

%e | |  |_|

%e |_|

%p h:= proc(l, f) option remember; local k; if min(l[])>0 then

%p      `if`(nops(f)=0, 1, h(map(x-> x-1, l[1..f[1]]), subsop(1=[][], f)))

%p     else for k from nops(l) while l[k]>0 by -1 do od;

%p         `if`(nops(f)>0 and f[1]>=k, h(subsop(k=2, l), f), 0)+

%p         `if`(k>1 and l[k-1]=0, h(subsop(k=1, k-1=1, l), f), 0)

%p       fi

%p     end:

%p g:= l-> `if`(add(`if`(l[i]::odd, (-1)^i, 0), i=1..nops(l))=0,

%p         `if`(l=[], 1, h([0$l[1]], subsop(1=[][], l))), 0):

%p b:= (n, i, l)-> `if`(n=0 or i=1, g([l[], 1$n]), b(n, i-1, l)

%p                   +b(n-i, min(n-i, i), [l[], i])):

%p a:= n-> b(2*n$2, []):

%p seq(a(n), n=0..12);

%t h[l_, f_] := h[l, f] = Module[{k}, If[Min[l]>0, If[Length[f] == 0, 1, h[l[[1 ;; f[[1]]]]-1, ReplacePart[f, 1 -> Nothing]]], For[k = Length[l], l[[k]]>0, k--]; If[Length[f]>0 && f[[1]] >= k, h[ReplacePart[l, k -> 2], f], 0] + If[k>1 && l[[k-1]] == 0, h[ReplacePart[l, {k -> 1, k-1 -> 1}], f], 0]]];

%t g[l_] := If[Sum[If[OddQ[l[[i]]], (-1)^i, 0], {i, 1, Length[l]}] == 0, If[l == {}, 1, h[Table[0, {l[[1]]}], ReplacePart[l, 1 -> Nothing]]], 0];

%t b[n_, i_, l_] := If[n == 0 || i == 1, g[Join[l, Table[1, {n}]]], b[n, i-1, l] + b[n-i, Min[n-i, i], Append[l, i]]];

%t a[n_] := b[2n, 2n, {}];

%t Table[a[n], {n, 0, 12}] (* _Jean-François Alcover_, Aug 29 2021, after _Alois P. Heinz_ *)

%Y Row sums of A304718.

%Y Bisection (even part) of A304680.

%Y Cf. A004003, A296625, A300060, A300789, A304677, A304710, A304790.

%K nonn

%O 0,2

%A _Alois P. Heinz_, May 16 2018