bzoj3294 [Cqoi2011]放棋子

最新推荐文章于 2019-07-10 15:10:00 发布

Euryale_

最新推荐文章于 2019-07-10 15:10:00 发布

阅读量822

点赞数 1

分类专栏：其他dp 文章标签：组合数动态规划

本文链接：https://blog.csdn.net/u010336344/article/details/70061843

版权

其他dp 专栏收录该内容

8 篇文章 0 订阅

订阅专栏

题目链接：bzoj3294
题目大意：
n行m列的棋盘放棋子，使每个格子最多放一个棋子，且不同颜色的棋子不能在同一行或者同一列。问有多少种方法。方案总数除以 1,000,000,009的余数。
N,M<=30 C<=10 总棋子数<=250

题解：
dp+组合数
%%%HYC
设g[i][j]表示选i行j列出来填k个棋子的方案数，要保证每行每列都要有棋子。
那么g[i][j]=总方案数-不合法的方案(即实际上有没被占的行或列的方案)
即

g [i] [j] = C k i \times j - \sum i i = 1 i \sum j j = 1 j g [i i] [j j] \times C i i i \times C j j j

$g[i][j]=C_{i\times j}^k-\sum_{ii=1}^i\sum_{jj=1}^jg[ii][jj]\times C_{i}^{ii}\times C_{j}^{jj}$
对于每种颜色都求一次对应的g。
然后设f[k][i][j]表示前k种颜色填了i行j列的方案数。
枚举第k种颜色选了ii行jj行来占领，然后转移就好了。
即

f [k] [i] [j] = \sum i i = 1 i \sum j j = 1 j f [k - 1] [i - i i] [j - j j] \times g [i i] [j j] \times C i i n - (i - i i) \times C j j m - (j - j j)

$f[k][i][j]=\sum_{ii=1}^{i}\sum_{jj=1}^{j}f[k-1][i-ii][j-jj]\times g[ii][jj]\times C_{n-(i-ii)}^{ii}\times C_{m-(j-jj)}^{jj}$
于是

ans=∑ni=1∑mj=1f[c][i][j] $ans=\sum_{i=1}^n\sum_{j=1}^mf[c][i][j]$

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;

const LL mod=1000000009;
int a[20];
LL C[910][910],f[15][35][35],g[35][35];
void init()
{
    for (int i=0;i<=910;i++) C[i][0]=1LL;
    for (int i=1;i<=910;i++)
     for (int j=1;j<=910;j++)
      C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
}
int main()
{
    //freopen("a.in","r",stdin);
    //freopen("a.out","w",stdout);
    int n,m,c,i,j,k,ii,jj;LL ans=0;
    scanf("%d%d%d",&n,&m,&c);
    init();
    for (i=1;i<=c;i++) scanf("%d",&a[i]);
    memset(f,0,sizeof(f));
    f[0][0][0]=1LL;
    for (k=1;k<=c;k++)
    {
        memset(g,0,sizeof(g));
        for (i=1;i<=n;i++)
         for (j=1;j<=m;j++)
          if (i*j>=a[k])
          {
              g[i][j]=C[i*j][a[k]];
              for (ii=1;ii<=i;ii++)
               for (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 (i=1;i<=n;i++)
         for (j=1;j<=m;j++)
         {
             for (ii=1;ii<=i;ii++)
              for (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==c) 
            {
                ans+=f[k][i][j]%mod;
                ans%=mod;
            }
         }
    }
    printf("%lld\n",ans);
    return 0;
}