一道关于行列式的经典题

【问题描述】

给定 n 点 m 边有向无环图，其中没有入度的点被视为源点，没有出度的点被视为

汇点。保证源点和汇点数目相同。

考虑所有把源汇点两两配对，用两两点不相交的路径把它们两两连接的所有方案

如果这个方案中，把源点按标号排序后，得到的对应汇点序列的逆序数对的个数

是奇数，那么 A 给 B 一块钱，否则 B 给 A 一块钱。 

问最后 A 的收益，对一个 p 取模。n，m，p 均为整数，p 为质数。

【输入格式】

第一行读入三个数，n，m，p，接下来 m 行，读入有向边。

【输出格式】

一行，输出答案。

【样例输入】

4 4 1000003

2 1

2 4

3 1

3 4

【样例输出】

0

【数据规模与约定】

对于 30%的数据，n<=100

对于 100%的数据，1<=n<=600,1<=m<=50000

看一下行列式讲解就会了

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>

using namespace std;

typedef long long ll;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=610,M=50100;

int last[N],ecnt;
struct EDGE{int to,nt;}e[M];
inline void add(int u,int v)
{e[++ecnt]=(EDGE){v,last[u]};last[u]=ecnt;}

int n,m,mod,num;

int pos[N];
int ind[N],outd[N];

int A[N][N];
int vis[N],f[N];

void topo(int u)
{
	if(vis[u]==num) return ;
	for(int i=last[u],v;i;i=e[i].nt)
	{
		topo(v=e[i].to);
		(f[u]+=f[v])%=mod;
	}
	vis[u]=num;
}

int det()
{
	register int i,j,k,f=1,x,y,t,res=1;
	for(i=1;i<=num;++i)
	{
		for(j=i+1;j<=num;++j)
		{
			x=A[i][i],y=A[j][i];
			while(y)
			{
				t=x/y;x%=y;swap(x,y);
				for(k=i;k<=num;++k)
					A[i][k]=((A[i][k]-1ll*t*A[j][k]%mod)%mod+mod)%mod;
				for(k=i;k<=num;++k) swap(A[i][k],A[j][k]);
				f=-f;
			}
		}
		if(!A[i][i]) return 0;
		res=1ll*res*A[i][i]%mod;
	}
	res=res*f;
	return (res%mod+mod)%mod;
}

int main()
{
	n=read();m=read();mod=read();
	register int i,j,u,v;
	for(i=1;i<=m;++i)
		u=read(),v=read(),
		ind[v]++,outd[u]++,
		add(v,u);
	for(i=1;i<=n;++i)
		if(!ind[i])
			pos[i]=++num;
	num=0;
	for(i=1;i<=n;++i)
		if(!outd[i])
			pos[i]=++num;
	num=0;
	for(i=1;i<=n;++i)
		if(!ind[i])
		{
			vis[i]=++num;
			memset(f,0,sizeof(f));f[i]=1;
			for(j=1;j<=n;++j)
				topo(j);
			for(j=1;j<=n;++j)
				if(!outd[j])
					A[pos[i]][pos[j]]=f[j];
		}
	cout<<det()<<endl;
	return 0;
}