bzoj 4596: [Shoi2016]黑暗前的幻想乡

这个题貌似奶了一口trump，，，

本来看到一个公司一条边，而且N是那么科学，肯定要容斥了，然后蛋疼的想应该怎么求生成树。。

本来一开始乱说了个矩阵树定理，然而忘了，就像奇奇怪怪的东西了。。然并卵，正解就是矩阵树定理，IQ--

需要注意的是，求解行列式和一般的线性方程是不一样的，行列式是消去本行，然后交换，变号。。。

#include <bits/stdc++.h>
#define LL long long
using namespace std;
inline int ra()
{
    int x=0,f=1; char ch=getchar();
    while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}

const int mod=1e9+7;

LL g[20][20];
int n;
vector<pair<int , int > > orz[20];

int Gauss()
{
    LL ans=1;
    for (int i=1; i<n; i++)
    {
        for (int j=i+1; j<n; j++)
            while (g[j][i])
            {
                LL t=g[i][i]/g[j][i];
                for (int k=i; k<n; k++)
                    g[i][k]=(g[i][k]-g[j][k]*t)%mod;
                swap(g[i],g[j]);
                ans=-ans;
            }
        ans=ans*g[i][i]%mod;
        if (!ans) return 0;
    }
    return (ans+mod)%mod;
}

int main(int argc, char const *argv[])
{
    n=ra();
    for (int i=1; i<n; i++)
    {
        int m=ra();
        for (int j=1; j<=m; j++)
        {
            int x=ra(),y=ra();
            orz[i].push_back(make_pair(x,y));
        }
    }
    int ans=0;
    for (int i=0; i<1<<(n-1); i++)
    {
        int cnt=0; memset(g,0,sizeof(g));
        for (int j=1; j<n; j++)
            if ((i>>(j-1))&1)
            {
                for (int k=0; k<orz[j].size(); k++)
                {
                    int x=orz[j][k].first,y=orz[j][k].second;
                    g[x][x]++; g[y][y]++;
                    g[x][y]--; g[y][x]--;
                }
                cnt++;
            }
        if ((n-cnt)&1) ans=(ans+Gauss())%mod; else ans=(ans-Gauss()+mod)%mod;
    }
    printf("%d
",ans);
    return 0;
}