1. 程式人生 > >bzoj2560串珠子——子集DP

bzoj2560串珠子——子集DP

-- 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(int
j=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