标签：open getchar zoj through esc 网上 ++ ack 大神
标准输入。输入第一行包含一个正整数n，表示珠子的个数。接下来n行，每行包含n个非负整数，用空格隔开。这n行中，第i行第j个数为ci,j。
标准输出。输出一行一个整数，为连接方案数对1000000007取模的结果。
对于100%的数据，n为正整数，所有的ci,j为非负整数且不超过1000000007。保证ci,j=cj,i。每组数据的n值如下表所示。
编号 1 2 3 4 5 6 7 8 9 10
n 8 9 9 10 11 12 13 14 15 16
正解：状压dp。
这道题很玄学，感觉还是懵懵懂懂。。今天的出题人太丧病了。。
考虑状压dp，我们把n个点压成二进制数。
我们设两个数组，g[s]表示s状态下的所有情况，即s状态下的点两两之间任意连边（包括不连边的情况），f[s]表示s状态下的合法情况，即使得s状态下所有点连通的合法情况。那么答案就是f[2^n-1]。
然后我们可以很容易地求出g[s]，我们求出g[s]后，考虑如何求f[s]，f[s]就是g[s]减去所有的不合法情况。那么我们可以枚举s的所有子集，设子集为i，那么不合法的情况就是g[i]*f[s^i]，我们减去这些情况，就能求出f[s]了。
枚举子集很玄学。。我反正不知道这是怎么回事，看了网上大神的博客。。
http://www.cnblogs.com/jffifa/archive/2012/01/16/2323999.html
for (int x = S; x; x = (x-1)&S)
大概是这个鬼东西。。大神的证明：
x = (x-1)&S实际上是把S中的0全部忽略，并不断减1的结果，比如S=1011，则x分别为：1011, 1010, 1001, 1000, 0011, 0010, 0001。忽略S中第二位的0其实就是111, 110, 101, 100, 011, 010, 001。
称S中的1所在位为有效位，0所在位为无效位，则x中的无效位必为0，有效位为0或1，比如S=1011，x=1001（有效位加下划线）。-1就是加上-1补码1111…，可以想成把无效位的1先加上去，比如x=1001变成1101，再加有效位的1。由于无效位加完肯定是1，会把有效位的进位“传递”下去，然后再位与S使得无效位变成0，实际就相当于有效位加上1111…，也就是有效位-1。
于是我们就能解决这题了。其实我还没搞懂为什么要这么搞。
1 //It is made by wfj_2048~ 2 #include <algorithm> 3 #include <iostream> 4 #include <cstring> 5 #include <cstdlib> 6 #include <cstdio> 7 #include <vector> 8 #include <cmath> 9 #include <queue> 10 #include <stack> 11 #include <map> 12 #include <set> 13 #define inf (1<<30) 14 #define il inline 15 #define RG register 16 #define ll long long 17 #define rhl 1000000007 18 #define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout) 19 20 using namespace std; 21 22 ll f[1<<16],g[1<<16],bin[20],a[20][20],n; 23 24 il int gi(){ 25 RG int x=0,q=1; RG char ch=getchar(); while ((ch<‘0‘ || ch>‘9‘) && ch!=‘-‘) ch=getchar(); 26 if (ch==‘-‘) q=-1,ch=getchar(); while (ch>=‘0‘ && ch<=‘9‘) x=x*10+ch-48,ch=getchar(); return q*x; 27 } 28 29 il void work(){ 30 n=gi(),bin[0]=1; 31 for (RG int i=1;i<=n;++i) 32 for (RG int j=1;j<=n;++j) a[i][j]=gi(); 33 for (RG int i=1;i<=n;++i) bin[i]=bin[i-1]<<1; 34 for (RG int k=0;k<bin[n];++k){ 35 f[k]=1; 36 for (RG int i=1;i<n;++i) if (k&bin[i-1]) 37 for (RG int j=i+1;j<=n;++j) if (k&bin[j-1]) f[k]=f[k]*(a[i][j]+1)%rhl; 38 g[k]=f[k]; RG int i=k^(k&-k); //枚举真子集？！ 39 for (RG int j=i;j;j=(j-1)&i) f[k]=(f[k]-g[j]*f[k^j]%rhl+rhl)%rhl; 40 } 41 printf("%lld\n",f[bin[n]-1]); return; 42 } 43 44 int main(){ 45 File("bunch"); 46 work(); 47 return 0; 48 }
标签：open getchar zoj through esc 网上 ++ ack 大神
原文地址：http://www.cnblogs.com/wfj2048/p/6501856.html
