bzoj3294: [Cqoi2011]放棋子

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const int maxl=33;
const int maxC=13,sumC=1100;
const LL mod=1e9+9;
LL C[maxl*maxl][maxl*maxl];
void initC()
{
    C[0][0]=1;
    for(int i=1;i<maxl*maxl;i++)
    {
        C[i][0]=1;
        for(int j=1;j<maxl*maxl;j++)
            C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
    }
}
//~~~~~~~~~~~init~~~~~~~~~~~~~~~~

LL g[maxC][maxl][maxl];//用cn[i]个棋子占i行j列并且一定占满 ------> 的方案数 
LL f[maxC][maxl][maxl];//前c种棋子占i行j列 -----> 的方案数 
int main()
{
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
    initC();
    int n,m,CC;
    scanf("%d%d%d",&n,&m,&CC);
    for(int t=1;t<=CC;t++)
    {
        int c;
        scanf("%d",&c);
        for(int i=0;i<=n;i++)
            for(int j=0;j<=m;j++)
            {
                g[t][i][j]=C[i*j][c];
                for(int u=0;u<=i;u++)
                    for(int v=0;v<=j;v++)
                    {
                        if(u==i&&v==j)continue;
                        g[t][i][j]=(g[t][i][j]-g[t][u][v]*C[i][u]*C[j][v])%mod;
                    }
                g[t][i][j]=(g[t][i][j]+mod)%mod;
            }
    }
    
    f[0][0][0]=1;
    for(int t=1;t<=CC;t++)
        for(int i=0;i<=n;i++)
            for(int j=0;j<=m;j++)
                for(int u=0;u+i<=n;u++)    
                    for(int v=0;v+j<=m;v++)
                        f[t][i+u][j+v]+=f[t-1][i][j]*g[t][u][v]%mod*C[n-i][u]%mod*C[m-j][v]%mod;
    LL ans=0;
    for(int i=0;i<=n;i++) 
        for(int j=0;j<=m;j++)
            ans=(ans+f[CC][i][j])%mod;
    printf("%lld\n",ans);
    
    return 0;
}