1. 程式人生 > >bzoj1087互不侵犯King(狀壓)

bzoj1087互不侵犯King(狀壓)

狀壓dp () int str color == namespace ble clas

題目:https://www.lydsy.com/JudgeOnline/problem.php?id=1087

簡單的狀壓dp。但是wa了好幾發。註意long long。

註意0和0的連邊。而且不能連成兩條,所以放在外面寫。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=15,M=1<<9;
int n,m,lm,list[M+5][M+5],cnt[M+5],zt[M+5],xnt,nm[M+5];//那一維是N*N! 
long long
ans,dp[N][N*N][M+5];// int num(int s) { int cnt=0; while(s)cnt++,s-=(s&-s); return cnt; } void init() { lm=(1<<n); for(int s=0;s<lm;s++) if(((s<<1)&s)==0) { zt[++xnt]=s; nm[xnt]=num(s); } list[1][++cnt[1]]=1;
for(int i=1;i<=xnt;i++) for(int j=i+1;j<=xnt;j++) if(((zt[i]<<1)&zt[j])==0&&((zt[i]>>1)&zt[j])==0&&(zt[i]&zt[j])==0) list[i][++cnt[i]]=j,list[j][++cnt[j]]=i; } int main() { scanf("%d%d",&n,&m); init();
for(int i=1;i<=xnt;i++)dp[1][nm[i]][i]=1; for(int i=2;i<=n;i++) for(int s=1;s<=xnt;s++) for(int j=nm[s];j<=m;j++) for(int k=1;k<=cnt[s];k++) dp[i][j][s]+=dp[i-1][j-nm[s]][list[s][k]]; // printf("dp[%d][%d][%d]=%d(dp[%d][%d][%d]=%d) %d %d\n",i,j,s,dp[i][j][s], // i-1,j-nm[s],list[s][k],dp[i-1][j-nm[s]][list[s][k]],zt[s],zt[list[s][k]]); for(int i=1;i<=xnt;i++)ans+=dp[n][m][i]; printf("%lld",ans); return 0; }

bzoj1087互不侵犯King(狀壓)