【BZOJ3659】Which Dreamed It

Description

有n个房间，每个房间有若干把钥匙能够打开特定房间的门。
你会做这么件事情：
最初你在房间1。
每当你到达一个房间，你可以选择该房间的一把钥匙，前往该钥匙对应的房间，并将该钥匙丢到垃圾桶中。
你希望：最终回到房间1，且垃圾桶中有所有的钥匙。
求方案数。两组方案不同，当且仅当使用钥匙的顺序不同。注意，每
把钥匙都是不同的。
Input

有多组数据。
对于每组数据第一行输入一个数n，表示房间数。
接下来n行依次描述每个房间：
首先一个数s，表示这个房间的钥匙数目，接下来s个数，分别描述每把
钥匙能够打开的房间的门。
输入以n-0结尾。
Output

对于每组数据，输出方案数，为了方便你的输出，请将答案对1000003取模。
Sample Input

1

0

2

1 1

1 2

0

Sample Output

1

0

HINT

在第一组样例中，没有钥匙，则方案数为1。

在第二组样例中，你不可能使用第二个房间的钥匙，所以方案数为0。

房间数小于等于100，钥匙数小于等于200000。

数据组数也不是特别多。

Source

这题是去年wangxz学长给我们出的胡策题真怀念以前的时光T_T虽然傻逼但是欢乐
题意就是求以1为起点的欧拉回路的个数乘1的度数
然后这里有两个结论
1.有向图以i为根的树形图的数目=基尔霍夫矩阵去掉第i行和第i列的主子式的行列式的值(即Matrix-Tree定理不仅适用于求无向图生成树数目,也适用于求有向图树形图数目)
2.以某个点为起点的欧拉回路数=该点为根的树形图数*所有点出度-1的乘积(本名BEST theorem)
关于BEST theorem可以提供一个wiki上的讲解:about BEST theorem
所以高斯消元求一遍行列式的值然后乘一下就没了
感觉这个题出的没什么意思..就是考你知不知道定理

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 110
#define LL long long
#define GET (ch>='0'&&ch<='9')
#define P 1000003
using namespace std;
int n,x,u,v;
int deg[MAXN];
LL a[MAXN][MAXN],fac[P+1],ans;
inline bool in(int &x)
{
    char ch=getchar();x=0;
    while (!GET)    {   if (ch==EOF)    return 0;ch=getchar();  }
    while (GET) x=x*10+ch-'0',ch=getchar();return x!=0;
}
int gauss(int n)
{
    if (n==0)   return 1;
    LL ret=1;int flag=1;
    for (int i=1;i<=n;i++)  for (int j=1;j<=n;j++)  ((a[i][j]%=P)+=P)%=P;
    for (int i=1;i<=n;i++)
    {
        for (int j=i+1;j<=n;j++)
            while (a[j][i])
            {
                LL tmp=a[i][i]/a[j][i];
                for (int k=i;k<=n;k++)  ((a[i][k]-=a[j][k]*tmp%P)+=P)%=P;
                swap(a[i],a[j]);flag*=-1;
            }
        if (!a[i][i])   return 0;
        (ret*=a[i][i])%=P;
    }
    return flag==-1?P-ret:ret;
}
int main()
{
    fac[0]=1;for (int i=1;i<=P;i++) fac[i]=fac[i-1]*i%P;
    while (in(n))
    {
        ans=1;
        for (int i=0;i<=n;i++)
        {
            deg[i]=0;
            for (int j=0;j<=n;j++)  a[i][j]=0;
        }
        for (int i=1;i<=n;i++)
        {
            in(deg[i]);
            for (int j=1;j<=deg[i];j++)
            {
                in(x);
                if (i!=x)   a[i][x]--,a[i][i]++;
            }
        }
        if (n==1&&!deg[1])  {   puts("1");continue; }
        for (int i=1;i<=n;i++)  (ans*=fac[deg[i]-1])%=P;
        (ans*=gauss(n-1))%=P;(ans*=deg[1])%=P;
        printf("%lld\n",ans);
    }
}