bzoj2560串珠子——子集DP
阿新 • • 發佈:2018-06-05
-- lin () mes online 連通 劃分 bzoj2560 tdi
題目:https://www.lydsy.com/JudgeOnline/problem.php?id=2560
轉載:
很明顯的狀壓dp
一開始寫的dp可能會出現重復統計的情況 而且難以去重
假設 一個狀態s的隨意連邊集合是A;
那麽 A應該是 全部合法的方案(Ans)+sigma(某一部分合法(即某一部分是連通圖)的方案*其他任意連邊的方案);
那麽可以把最終答案設置為f[i], 隨意連邊(也可以完全連邊)設置成g[i];
先定一個基準點 x 和基準點相連的都是合法的, 其余集合 t=s^(1<<(x-1))可以隨便連;
f[i]=g[i]-sigma((t的所有子集i)f[i]*g[s^i]);
為什麽一個是f 一個是g 這樣其實是要保證不重不漏
而且必須註意劃分點一定是連通部分的,否則樣例輸出49...不知道為什麽。
代碼如下:
#include<iostream> #include<cstdio> #include<cstring> using namespace std; typedef long long ll; ll const mod=1000000007; ll n,a[20][20],f[1<<16],g[1<<16]; int main() { scanf("%d",&n); for(int i=1;i<=n;i++) for(intj=1;j<=n;j++) scanf("%d",&a[i][j]); int m=(1<<n)-1; for(int s=1;s<=m;s++) { g[s]=1; for(int i=1;i<=n;i++) if(s&(1<<(i-1))) for(int j=i+1;j<=n;j++) if(s&(1<<(j-1))) (g[s]*=(a[i][j]+1))%=mod; f[s]=g[s]; int nw; for(nw=n-1;;nw--) if(s&(1<<nw))break; nw=(s^(1<<nw));//去掉一個劃分點 for(int k=nw;k;k=((k-1)&nw))//枚舉nw的子集 ((f[s]-=g[k]*f[s^k])+=mod)%=mod;//f和g別反 } printf("%lld",f[m]); return 0; }
bzoj2560串珠子——子集DP