【题解】Bzoj2560串珠子

　　挺强的……容斥+状压DP。首先想到如果可以求出f[k]，f[k]代表联通状态为k的情况下的合法方案数，则f[k] = g[k] - 非法方案数。g[k]为总的方案数，这是容易求得的。那么非法方案数我们可以枚举 k 的子集 j，则 j 联通而剩下的则随意连（不与j联通）。可是做到这里以为自己做出来了，实际上并没有……

　　注意到枚举到子集 j 时，若 s' = k - j， 那如果 s' 中有一个联通的方案 s''，我们在这里减去一次，在之后枚举到s''时又会枚举到这个方案一次。实际上，这也就是说0111与1000这两个子集是对称的。所以我们为了避免这样的情况，就锁定一个点a，使得点 a 一定不出现在集合 j 中，可以使得 a 只能出现在集合 s' 中，也就避免了重复。

　　代码有参考，如有雷同，是我抄的 (o´ω`o)ﾉ

#include <bits/stdc++.h>
using namespace std;
#define maxn 20
#define maxm ((1 << 16) + 2)
#define mod 1000000007
#define int long long
int n, bin[maxn], a[maxn][maxn];
int f[maxm], g[maxm];

int read()
{
    int x = , k = ;
    char c;
    c = getchar();
    while(c < '' || c > '') { if(c == '-') k = -; c = getchar(); }
    while(c >= '' && c <= '') x = x *  + c - '', c = getchar();
    return x * k;
}

signed main()
{
    n = read(); bin[] = ;
    for(int i = ; i <= n; i ++)
        for(int j = ; j <= n; j ++)
            a[i][j] = read();
    for(int i = ; i <= n; i ++) bin[i] = bin[i - ] << ;
    for(int k = ; k < bin[n]; ++ k)
    {
        f[k] = ;
        for(int i = ; i < n; i ++)
            if(k & bin[i - ])
                for(int j = i + ; j <= n; j ++)
                    if(k & bin[j - ])
                        f[k] = f[k] * (a[i][j] + ) % mod;
        g[k] = f[k]; int K = (k ^ (k & -k));
        for(int j = K; j; j = (j - ) & K)
            f[k] = (f[k] - g[j] * f[k ^ j] % mod + mod) % mod;
    }
    printf("%lld\n", f[bin[n] - ]);
    return ;
}