BZOJ3659: Which Dreamed It

题目大意：给定一个有向图，求从1号节点出发回到1号节点并经过所有边的方案数*1的度数

.....

为什么会有这种题出现啊

定理题...

1.有向图以i为根的树形图的数目=基尔霍夫矩阵去掉第i行和第i列的主子式的行列式的值(即Matrix-Tree定理不仅适用于求无向图生成树数目,也适用于求有向图树形图数目)
2.以某个点为起点的欧拉回路数=该点为根的树形图数*所有点出度-1的乘积(本名BEST theorem)

http://blog.csdn.net/creationaugust/article/details/51074919

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 110
using namespace std;
long long n,nn;
long long a[N][N],du[N][N],lj[N][N];
long long rd[N],cd[N];
bool vis[N];
long long ni[1000010],pre[200010],mod=1000003;
long long ksm(long long d,long long c)
{
	long long ret=1;
	while(c)
	{
		if(c&1) ret=ret*d%mod;
		d=d*d%mod;c/=2;
	}
	return ret;
}
void dfs(long long x)
{
	long long i,j;
	vis[x]=true;
	for(i=1;i<=n;i++)
	if(lj[x][i]&&!vis[i])
	dfs(i);
}
long long sj[N];
void gaosi()
{
	long long i,j,k,t;
	for(i=1;i<=n;i++)
	{
		j=i;
		while(!a[j][i]) j++;
		for(k=1;k<=n;k++)
		swap(a[i][k],a[j][k]);
		for(j=i+1;j<=n;j++)
		{
			t=ni[a[i][i]]*a[j][i]%mod;
			for(k=i;k<=n;k++)
			{
				a[j][k]-=a[i][k]*t;
				a[j][k]=(a[j][k]%mod+mod)%mod;
			}
		}
	}
	long long ans=1;
	for(i=1;i<=n;i++)
	ans=ans*a[i][i]%mod;
	for(i=1;i<=nn;i++)
	ans=ans*pre[cd[i]-1]%mod;
	ans=ans*cd[1]%mod;
	printf("%lld\n",ans);
}
void doit()
{
	long long i,j,x,y;
	memset(lj,0,sizeof(lj));
	memset(vis,0,sizeof(vis));
	memset(rd,0,sizeof(rd));
	for(i=1;i<=n;i++)
	{
		scanf("%lld",&x);
		cd[i]=x;
		for(j=1;j<=x;j++)
		{
			scanf("%lld",&y);
			rd[y]++;
			lj[i][y]++;
		}
	}
	dfs(1);
	for(i=1;i<=n;i++)
	if(rd[i]!=cd[i]||(cd[i]!=0&&!vis[i])) {puts("0");return;}
	for(i=1;i<=n;i++)
	du[i][i]=cd[i];
	nn=0;
	for(i=2;i<=n;i++)
	if(cd[i])
	{
		nn++;
		sj[nn]=i;
	}
	swap(n,nn);
	for(i=1;i<=n;i++)
	for(j=1;j<=n;j++)
	a[i][j]=((du[sj[i]][sj[j]]-lj[sj[i]][sj[j]])%mod+mod)%mod;
	gaosi();
}
int main()
{
	long long i;
	for(i=1;i<mod;i++)
	ni[i]=ksm(i,mod-2);
	pre[0]=1;
	for(i=1;i<=200000;i++)
	pre[i]=pre[i-1]*i%mod;
	while(scanf("%lld",&n)==1&&n!=0)
	doit();
}