BZOJ 3812主旋律

求一个图中强联通图的个数。

一看就是容斥啦，但这种二进制高端操作还是学习一下Candy?dalao

注释在代码里

好久没更了。。。

#include<bits/stdc++.h>
using namespace std;
const int N=(1<<16)+10,mod=1e9+7;
typedef long long ll;
int n,m,x,y,f[N],g[N],h[N],one[N],w[N],in[N],out[N];
ll bin[N];
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;++i)
    {
        scanf("%d%d",&x,&y);
        out[1<<x-1]|=(1<<y-1);in[1<<y-1]|=(1<<x-1);
    }
    bin[0]=1;
    for(int i=1;i<=m;++i)bin[i]=bin[i-1]*2%mod;
    for(int s=1;s<1<<n;s++)one[s]=1+one[s^(s&-s)];
    for(int s=1;s<1<<n;s++){
        int p=s&-s,ns=s^p;
        for(int t=ns;t;t=(t-1)&ns)
            g[s]=(g[s]-1ll*f[s^t]*g[t]%mod)%mod;
        h[s]=h[ns]+one[in[p]&ns]+one[out[p]&ns];//如何求一个诱导子图中边数 
        f[s]=bin[h[s]];
        for(int t=s;t;t=(t-1)&s){
            int p=(s^t)&-(s^t);
            if(t==s)w[t]=0;
            else w[t]=w[t^p]-one[out[p]&(s^t)]+one[in[p]&t];//求两个点集间的边数 
            f[s]=(f[s]-1ll*g[t]*bin[w[t]+h[s^t]]%mod)%mod;
        }
        if(f[s]<0)f[s]+=mod;
        g[s]=(g[s]+f[s])%mod;
    }
    printf("%d",f[(1<<n)-1]);
    return 0;
}