HDU 5729 Rigid Frameworks（组合数学+dp）

Description
矩形是不稳定的，会变成平行四边形，但是可以在矩形对角线加边使这个矩形稳定下来。给一个n*m的矩形，可以在格里加两种对角线（主对角线或斜对角线），或者不加对角线，问使这个n*m的矩形稳定下来的方案数
Input
多组用例，每组用例占一行为两个整数n和m，以文件尾结束输入(1<=n,m<=10)
Output
对于每组用例，输出使这个n*m矩形稳定下来的方案数，结果模1e9+7
Sample Input
1 2
3 2
7 9
10 10
Sample Output
4
448
357533852
935300639
Solution
对这个n*m矩形，不管将其扭曲成什么样，每一行的所有竖边一定是平行的，同样的每一列的所有横边也一定是平行的，而在第i行第j列这个格子中加一条斜边会使得该格子的横边和竖边垂直，进而使得该格子的横边与该行所有竖边垂直，该格子的竖边与该列的所有横边垂直，将n行看作第一排点，m列看作第二排点，那么在第i行第j列放一个斜边相当于在第一排点的第i个点与第二排点的第j个点之间连一条边，而这个矩阵如果稳定，那么说明任意格子的横边与竖边都需要垂直，即这个二分图中第一排点中任一点都可以到第二排点中任一点，即这个二分图是连通的，所以问题转化为一个左边n个点右边m个点的连通二分图个数
令dp[i][j][k]表示左边i个点右边j个点连k条边的连通二分图个数
如果不考虑连通，这样的二分图个数为个
下面求不连通的二分图个数，考虑包含第一个点的连通块点数为(ii,jj)，即左边ii个点，右边jj个点，那么不连通二分图个数应该为
故有转移方程

最后的答案即为
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll mod=1000000007ll;
#define maxn 111
int n,m;
ll f[maxn],dp[11][11][maxn],C[maxn][maxn],ans;
void init()
{
    f[0]=1,C[0][0]=1;
    for(int i=1;i<=100;i++)
    {
        f[i]=f[i-1]*2%mod;
        C[i][i]=C[i][0]=1;
        for(int j=1;j<i;j++)C[i][j]=(C[i-1][j-1]+C[i-1][j])%mod;
    }
    memset(dp,-1,sizeof(dp));
}
int dfs(int n,int m,int k)
{
    if(n<m)swap(n,m);
    if(~dp[n][m][k])return dp[n][m][k];
    ll cnt=0;
    for(int i=1;i<=n;i++)
        for(int j=0;j<=m;j++)
        {
            if(i+j==n+m)break;
            for(int l=0;l<=min(k,i*j);l++)
                cnt=(cnt+dfs(i,j,l)*C[n-1][i-1]%mod*C[m][j]%mod*C[(n-i)*(m-j)][k-l]%mod)%mod;
        }
    return dp[n][m][k]=(C[n*m][k]-cnt+mod)%mod;
}
int main()
{
    init();
    while(~scanf("%d%d",&n,&m))
    {
        ans=0;
        if(n<m)swap(n,m);
        for(int i=0;i<=m*n;i++)ans=(ans+dfs(n,m,i)*f[i]%mod)%mod;
        printf("%I64d\n",ans);
    }
    return 0;
}