3294: [Cqoi2011]放棋子

Description

Input

输入第一行为两个整数n, m, c，即行数、列数和棋子的颜色数。第二行包含c个正整数，即每个颜色的棋子数。所有颜色的棋子总数保证不超过nm。
Output

输出仅一行，即方案总数除以 1,000,000,009的余数。
题解：
g[i][j]表示i行j列放k这种颜色的棋子的合法方案数。
f[k][i][j]表示前k种棋子占据i行j列有多少种合法的放法 。
那么，g[i][j]=总方案数-不合法的方案(即实际上有没被占的行或列的方案)

g[i][j]=$C_{i*j}^{a[k]}-\sum_{ii=1}^{ii<=i}\sum_{jj=1}^{jj<=j}g[ii][jj]*C_{i}^{ii}*C_{j}^{jj}$
对于每种颜色都求一次对应的g。

f[k][i][j]=$\sum_{ii=1}^{ii<=i}\sum_{jj=1}^{jj<=j}f[k-1][i-ii][j-jj]*g[ii][jj]*C_{n-(i-ii)}^{ii}*C_{m-(j-jj)}^{jj}$
代码：

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cstring>

using namespace std;
typedef long long ll;
int n,m,c1;
int a[20];
ll g[35][35],f[15][35][35],c[920][920];
ll mod=1000000009;
void prepare()
{
    c[0][0]=1ll;
    for(int i=1;i<=910;i++)
    {
        c[i][0]=1ll;
        for(int j=1;j<=910;j++)
        {
           c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
        }
    }

}
int main()
{
    scanf("%d%d%d",&n,&m,&c1);
    prepare();
    for(int i=1;i<=c1;i++) scanf("%d",&a[i]);

    memset(f,0,sizeof(f));
    f[0][0][0]=1ll;
    ll ans=0;
    for(int k=1;k<=c1;k++)
    {
        memset(g,0,sizeof(g));
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=m;j++)
            {
                if(i*j>=a[k])
                {
                    g[i][j]=c[i*j][a[k]];
                    for(int ii=1;ii<=i;ii++)
                    for(int jj=1;jj<=j;jj++)
                    {
                        if(ii==i&&jj==j) continue;
                        g[i][j]-=g[ii][jj]*c[i][ii]%mod*c[j][jj]%mod;
                        g[i][j]=(g[i][j]+mod)%mod;
                    }
                }
            }
        }

        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=m;j++)
            {
                for(int ii=1;ii<=i;ii++)
                {
                    for(int jj=1;jj<=j;jj++)
                    {
                        if(ii*jj>=a[k])
                        {
                            f[k][i][j]+=f[k-1][i-ii][j-jj]*g[ii][jj]%mod*c[n-(i-ii)][ii]%mod*c[m-(j-jj)][jj]%mod;
                            f[k][i][j]%=mod;

                        }
                    }
                }
                if(k==c1)
                {
                    ans+=f[k][i][j]%mod;
                    ans%=mod;
                }
            }
        }
    }
    printf("%lld",ans);
}
//蒟蒻300题留念