Codeforces 1237F. Balanced Domino Placements

传送门

很妙的题

首先先考虑一个简化的问题，现在有一行格子让你填

你要么填一格 要么填两格 有的格子不让你填 问你填了 $a$ 个一格和填了 $b$ 个两格有多少种方案

那么显然先只考虑放两格的方案，这个可以简单 $dp$ 得到，设 $f[i][j]$ 表示前 $i$ 个格子放了 $j$ 个两格的方案数

那么如果 $i,i-1$ 都没障碍，那么 $f[i][j]=f[i-1][j]+f[i-2][j-1]$ ，否则 $f[i][j]=f[i-1][j]$

然后再来考虑填一格的，显然剩下的 $tot-2j$ 个位置都可以随便填，那么方案数为 $C[tot-2j][i]$ ，直接乘法原理乘起来即可

接下来可以考虑怎么把这道题简化到这个情况，假设放了 $a$ 个水平的多米诺，$b$ 个垂直的多米诺

对于行来说，相当于放 $a$ 个 $1$ ，$b$ 个 $2$，对于列就相当于放 $a$ 个 $2$ ，$b$ 个 $1$

注意到每个多米诺可以根据在第几行和第几列来唯一确定，所以我们对行列分别求一下之前那个东西然后乘起来再乘上 $a!b!$ 即可

乘上 $a!b!$ 就相当于把骨牌不同的放置顺序看成不同的放置方案，意思是强制行的第 $i$ 个放置和列的第 $i$ 个放置配对

就原本一维的方案我们让它放置有顺序，然后强制行和列两两匹配，这样才能确定二维平面上的具体位置

因为如果只是行列乘起来，那么没法确定某个骨牌的具体位置，考虑对于 $a$ ，行放了位置 $1,3$ ，列放了位置 $(1,2),(3,4)$

那么多米诺骨牌可以是 $((1,1),(1,2))$ 和 $((3,3),(3,4))$ ，但是也有可能是 $((1,3),(1,4))$ 和 $((3,1),(3,2))$

自己画画图就很容易理解了

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long ll;
inline int read()
{
    int x=,f=; char ch=getchar();
    while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
    while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
    return x*f;
}
const int N=,mo=;
inline int fk(int x) { return x>=mo ? x-mo : x; }
int n,m,K,Ans,fac[N],C[N][N];
int f[N][N],g[N][N];
bool px[N],py[N];
int main()
{
    n=read(),m=read(),K=read(); int mx=max(n,m);
    fac[]=; for(int i=;i<=mx;i++) fac[i]=1ll*fac[i-]*i%mo;
    for(int i=;i<=mx;i++)
    {
        C[i][]=;
        for(int j=;j<=i;j++)
            C[i][j]=fk(C[i-][j-]+C[i-][j]);
    }
    for(int i=;i<=K;i++)
    {
        int a=read(),b=read(),c=read(),d=read();
        px[a]=px[c]=; py[b]=py[d]=;
    }
    int cntx=,cnty=;
    for(int i=;i<=n;i++) cntx+=px[i];
    for(int i=;i<=m;i++) cnty+=py[i];
    for(int i=;i<=n;i++) f[i][]=;
    for(int i=;i<=n;i++)
        for(int j=;j<=mx;j++)
            if(px[i]||px[i-]) f[i][j]=f[i-][j];
            else f[i][j]=fk(f[i-][j]+f[i-][j-]);
    for(int i=;i<=m;i++) g[i][]=;
    for(int i=;i<=m;i++)
        for(int j=;j<=mx;j++)
            if(py[i]||py[i-]) g[i][j]=g[i-][j];
            else g[i][j]=fk(g[i-][j]+g[i-][j-]);
    for(int i=;i<=mx;i++)
        for(int j=;j<=mx;j++)
            if(i+j*<=n-cntx&&i*+j<=m-cnty)
                Ans=fk(Ans + 1ll*f[n][j]*C[n-cntx-j*][i]%mo *g[m][i]%mo *C[m-cnty-i*][j]%mo *fac[i]%mo *fac[j]%mo );
    printf("%d\n",Ans);
    return ;
}