[SHOI2016]黑暗前的幻想乡

一、题目

点此看题

二、解法

好久没有做过这么版的题了

首先你要知道指数级容斥，也就是我们枚举一个公司集合，表示我们只在这些公司中选边，那么我们就用全部公司-少一个公司+少两个公司…来计算就行。

知道边算生成树，不就是矩阵树定理么？

#include <cstdio>
#include <vector>
#include <cstring>
#include <iostream>
using namespace std;
#define int long long
const int M = 20;
const int MOD = 1e9+7;
int read()
{
    int num=0,flag=1;
    char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
    return num*flag;
}
int n,m,a[M][M],rt,ans;
vector<pair<int,int> > b[M];
int qkpow(int a,int b)
{
    int r=1;
    while(b>0)
    {
        if(b&1) r=r*a%MOD;
        a=a*a%MOD;
        b>>=1;
    }
    return r;
}
void guass()
{
    for(int i=2;i<=n;i++)
    {
        for(int j=i+1;j<=n;j++)
            if(!a[i][i] && a[j][i])
            {
            	ans=-ans;
                swap(a[i],a[j]);
                break;
            }
        int inv=qkpow(a[i][i],MOD-2);
        for(int j=i+1;j<=n;j++)
        {
            int tmp=a[j][i]*inv%MOD;
            for(int k=i;k<=n;k++)
                a[j][k]=(a[j][k]-a[i][k]*tmp%MOD)%MOD;
        }
    }
    for(int i=2;i<=n;i++)
    	ans=ans*a[i][i]%MOD;
    ans=(ans+MOD)%MOD;
}
signed main()
{
	n=read();
	for(int i=0;i<n-1;i++)
	{
		int k=read();
		while(k--)
		{
			int u=read(),v=read();
			b[i].push_back(make_pair(u,v));
		}
	}
	m=1<<(n-1);
	for(int i=1;i<m;i++)
	{
		int cnt=n-1;
		memset(a,0,sizeof a);
		for(int j=0;j<n-1;j++)
			if(i&(1<<j))
			{
				cnt--;
				for(int k=0;k<b[j].size();k++)
				{
					int u=b[j][k].first,v=b[j][k].second;
					a[u][u]++;a[v][v]++;
					a[u][v]--;a[v][u]--;
				}
			}
		ans=1;
		guass();
		rt=(rt+(cnt&1?-1:1)*ans)%MOD;
	}
	printf("%lld\n",(rt+MOD)%MOD);
}