HDU5729 Rigid Frameworks

题目连接：http://acm.hdu.edu.cn/showproblem.php?pid=5729

【题意】对于一个n*m的矩形，矩形由1*1的正方形组成，求添加斜线使得整个矩形固定的方案数。

【分析】对于i行j列，假如在格子(i,j)上添加斜线，那么i行j列将联动。于是可以得到一个左边n个点，右边m个点的二分图，题目转换为求使该二分图连通的方案数。注意每个格子都有两种添加方法，最后解为方案数*2^(连线数)。

【代码】

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#define LL long long
const int mod=1e9+7;
LL dp[15][15][105];
LL c[105][105];
int pow[111];
int main(){
        int n,m;
        c[0][0]=pow[0]=1;
        for(int i=1;i<=100;++i){
                pow[i]=(pow[i-1]*2)%mod;
                c[i][0]=c[i][i]=1;
                for(int j=1;j<i;++j)
                        c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
        }
        dp[0][1][0]=dp[1][0][0]=1;
        for(int i=1;i<=10;++i)
                for(int j=1;j<=10;++j){
                        if(i==1 || j==1){
                                dp[i][j][i+j-1]=1;
                                continue;
                        }
                        for(int k=1;k<=i*j;++k){
                                dp[i][j][k]=c[i*j][k];
                                for(int ii=1;ii<=i;++ii)
                                        for(int jj=0;jj<=j;++jj){
                                                if(i+j==ii+jj)
                                                        continue;
                                                for(int kk=0;kk<=k;++kk){
                                                        LL tmp=c[i-1][ii-1]*c[j][jj]%mod;
                                                        tmp=tmp*dp[ii][jj][kk]%mod;
                                                        tmp=tmp*c[(i-ii)*(j-jj)][k-kk]%mod;
                                                        dp[i][j][k]=(dp[i][j][k]+mod-tmp)%mod;
                                                }
                                        }
                        }
                }

        while(cin>>n>>m){
                LL ans=0;
                for(int i=0;i<=n*m;++i)
                        ans=(ans+pow[i]*dp[n][m][i])%mod;
                cout<<ans<<endl;;
        }
}