[BZOJ1087][SCOI2005]互不侵犯King

最新推荐文章于 2019-08-14 11:30:17 发布

Zvezda_

最新推荐文章于 2019-08-14 11:30:17 发布

阅读量417

点赞数

分类专栏： BZOJ DP 状压DP

本文链接：https://blog.csdn.net/Zvezda_/article/details/46947771

版权

BZOJ 同时被 3 个专栏收录

159 篇文章 0 订阅

订阅专栏

31 篇文章 0 订阅

订阅专栏

状压DP

1 篇文章 0 订阅

订阅专栏

原题地址

经典的状压DP.

AC code:

#include <cstdio>
const int N=10;
const int M=1<<10;
typedef long long ll;
int n,m;
int tot[M];
ll ans;
ll f[N][M][N*N];
bool wrng[M];
bool ok[M][M];

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=0;i<(1<<n);i++){
        if(i&(i<<1)){
            wrng[i]=1;
            continue;
        }
        int t=i;
        while(t){
            tot[i]+=t&1;
            t>>=1;
        }
    }
    for(int i=0;i<(1<<n);i++){
        for(int j=0;j<(1<<n);j++){
            if(wrng[i]||wrng[j]) continue;
            if(!(i&j)&&!((i<<1)&j)&&!((i>>1)&j)) ok[i][j]=1;
        }
    }
    for(int i=0;i<(1<<n);i++){
        if(!wrng[i]) f[1][i][tot[i]]=1;
    }
    for(int i=2;i<=n;i++){
        for(int j=0;j<(1<<n);j++){
            if(wrng[j]) continue;
            for(int k=tot[j];k<=m;k++){
                for(int jj=0;jj<(1<<n);jj++){
                    if(ok[j][jj]) f[i][j][k]+=f[i-1][jj][k-tot[j]];
                }
            }
        }
    }
    for(int i=0;i<(1<<n);i++){
        if(wrng[i]) continue;
        ans+=f[n][i][m];
    }
    printf("%lld",ans);

    return 0;
}