HDU 5729 Rigid Frameworks(连通性DP)

forever97 posted @ 2016年7月22日 13:26 in 算法-连通性DP with tags 连通性DP , 354 阅读

 

【题目链接】 http://acm.hdu.edu.cn/showproblem.php?pid=5729

 

【题目大意】

    给出一个n*m的方格框,可以在单位矩形中添加两种对角线的线,使得其变得稳定,问使得其变成稳定图形的方案数。

 

【题解】

    稳定状态指的是在n*m范围内每行每列都有一个固定的格子,并且联动,计算合法的情况非常的复杂,难以枚举,因此我们可以枚举非法情况,从组合数中减去即可。非法情况即将图形划分为合法部分和非法部分,注意枚举全面。

 

【代码】

#include <cstdio>
using namespace std;
#define rep(i,n) for(int i=1;i<=n;i++)
const long long mod=1e9+7;
long long dp[15][15][105],c[105][105];
int n,m;
long long pow(long long a,long long b){
    long long tmp=1; a%=mod;
    while(b){
        if(b&1)tmp=tmp*a%mod;
        a=a*a%mod;
        b>>=1;
    }return tmp;
}
int main(){
    rep(i,100){
        c[i][0]=c[i][i]=1;
        rep(j,i-1)c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
    }c[0][0]=dp[1][0][0]=dp[0][1][0]=1;
    rep(i,10)rep(j,10){
        if(i==1||j==1){dp[i][j][i+j-1]=1;continue;}
        rep(k,i*j){
            dp[i][j][k]=c[i*j][k];
            rep(ii,i)for(int jj=0;jj<=j;jj++){
                if(i+j==ii+jj)continue;
                for(int kk=0;kk<=k;kk++){
                    long long 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(~scanf("%d%d",&n,&m)){
        long long ans=0;
        rep(i,n*m)ans=(ans+pow(2,i)*dp[n][m][i]%mod)%mod;
        printf("%lld\n",ans);
    }return 0;
}

 

  • 无匹配
  • 无匹配

登录 *


loading captcha image...
(输入验证码)
or Ctrl+Enter