Atcoder CF 2017 TR I

Atcoder CF 2017 TR I

给定一个有n个点，m条边的图，求为每条边定向，使得从1出发和2出发的两个人可以见面的方案数。

先把问题转换成求all-不能见面的方案数。那么可以把图划分成这样一个集合：

用\(f[0/1][i]\)表示包含0/1的集合i,从0/1出发刚好能到达所有点的方案数 。注意是“刚好“，所以统计的是闭合子图。利用容斥原理求出f，接下来的任务就是统计other集合中，与两个集合不相关的边的个数即可。

#include <cstdio> 
using namespace std;

const int maxn=20, maxm=120, maxs=4e4, p=1e9+7;
int n, m, a[maxm], b[maxm], c[maxs], d[maxs], mi2[maxm];
long long f[2][maxs], ans;  
//f[0/1][i]表示包含0/1的集合i,从0/1出发刚好能到达所有点的方案数 

int main(){
    scanf("%d%d", &n, &m);
    for (int i=mi2[0]=1; i<maxm; ++i)
        mi2[i]=mi2[i-1]*2%p;
    for (int i=0; i<m; ++i){
        scanf("%d%d", &a[i], &b[i]);  //a b 存边 
        --a[i]; --b[i]; }
    for (int i=0; i<(1<<n); ++i)   
        for (int j=0; j<m; ++j){
            if ((i>>a[j]&1)&&(i>>b[j]&1)) ++c[i];  //c:两边都在s中的点数 
            if ((i>>a[j]&1)||(i>>b[j]&1)) ++d[i];  //d:恰好一个点在s中的点数 
        }
    for (int i=0; i<2; ++i)
        for (int j=0; j<1<<n; ++j){
            if ((j>>i&1)==0) continue;
            f[i][j]=mi2[c[j]];
            for (int k=j; k>0; k=(k-1)&j)  //注意是“刚好能到所有点” 
                f[i][j]=(f[i][j]+p-f[i][j-k]*mi2[c[k]])%p;
        }
    for (int i=0; i<1<<n; ++i){
        for (int j=0; j<1<<n; ++j){
            if ((i&j)||c[i]+c[j]!=c[i|j]) continue;
            ans=(ans+f[0][i]*f[1][j]%p*mi2[m-d[i|j]])%p;  //不与集合相连的边随便选 
        }
    }
    printf("%lld\n", (mi2[m]-ans+p)%p);
    return 0;
}