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);
}
}