login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A304662
Total number of domino tilings of Ferrers-Young diagrams summed over all partitions of 2n.
11
1, 2, 6, 16, 42, 106, 268, 650, 1580, 3750, 8862, 20598, 47776, 109248, 248966, 562630, 1264780, 2823958, 6282198, 13884820, 30590124, 67051982, 146463790, 318588916, 690882926, 1492592450, 3215372064, 6904561416, 14786529836, 31574656096, 67261524262
OFFSET
0,2
FORMULA
a(n) = Sum_{k=0..A304790(n)} k * A304789(n,k).
a(n) = Sum_{k=0..n} A304718(n,k).
a(n) = A296625(n) for n < 7.
EXAMPLE
a(2) = 6:
._. .___. ._._. .___. ._.___. .___.___.
| | |___| | | | |___| | |___| |___|___|
|_| | | |_|_| |___| |_|
| | |_|
|_|
MAPLE
h:= proc(l, f) option remember; local k; if min(l[])>0 then
`if`(nops(f)=0, 1, h(map(x-> x-1, l[1..f[1]]), subsop(1=[][], f)))
else for k from nops(l) while l[k]>0 by -1 do od;
`if`(nops(f)>0 and f[1]>=k, h(subsop(k=2, l), f), 0)+
`if`(k>1 and l[k-1]=0, h(subsop(k=1, k-1=1, l), f), 0)
fi
end:
g:= l-> `if`(add(`if`(l[i]::odd, (-1)^i, 0), i=1..nops(l))=0,
`if`(l=[], 1, h([0$l[1]], subsop(1=[][], l))), 0):
b:= (n, i, l)-> `if`(n=0 or i=1, g([l[], 1$n]), b(n, i-1, l)
+b(n-i, min(n-i, i), [l[], i])):
a:= n-> b(2*n$2, []):
seq(a(n), n=0..12);
MATHEMATICA
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]]];
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];
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]]];
a[n_] := b[2n, 2n, {}];
Table[a[n], {n, 0, 12}] (* Jean-François Alcover, Aug 29 2021, after Alois P. Heinz *)
CROSSREFS
Row sums of A304718.
Bisection (even part) of A304680.
Sequence in context: A102699 A266124 A217194 * A296625 A156664 A025169
KEYWORD
nonn
AUTHOR
Alois P. Heinz, May 16 2018
STATUS
approved