[子集DP] BZOJ 2560 串珠子

子集DP的裸题

采用补集转化的思想 减去不连通的 

我们可以 任选一个点 枚举这个点的联通块 也就是当前状态的子集

子集枚举 http://www.cnblogs.com/jffifa/archive/2012/01/16/2323999.html

for (int x = S; x; x = (x-1)&S)

称S中的1所在位为有效位，0所在位为无效位，则x中的无效位必为0，有效位为0或1，比如S=1011，x=1001（有效位加下划线）。-1就是加上-1补码1111…，可以想成把无效位的1先加上去，比如x=1001变成1101，再加有效位的1。由于无效位加完肯定是1，会把有效位的进位“传递”下去，然后再位与S使得无效位变成0，实际就相当于有效位加上1111…，也就是有效位-1。

#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef long long ll;

const int N=20;
const int P=1000000007;

int n,a[N][N];
ll f[1<<N],g[1<<N];
int bin[1<<N];

int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  scanf("%d",&n);
  bin[1]=1;
  for (int i=2;i<=(1<<n);i<<=1) bin[i]=bin[i>>1]+1;
  for (int i=1;i<=n;i++) for (int j=1;j<=n;j++) scanf("%d",&a[i][j]);
  for (int t=1;t<(1<<n);t++){
    g[t]=1;
    for (int i=t;i;i-=i&-i)
      for (int j=t;j;j-=j&-j)
	if (bin[i&-i]<bin[j&-j])
	  (g[t]*=(a[bin[i&-i]][bin[j&-j]]+1))%=P;
    f[t]=g[t];
    int k=t-(t&-t);
    for (int i=k;i;i=(i-1)&k)
      (f[t]+=P-g[i]*f[t^i]%P)%=P;
  }
  printf("%lld\n",f[(1<<n)-1]);
  return 0;
}