【BZOJ5297】【CQOI2018】社交网络（有向图生成树计数）

最新推荐文章于 2024-01-11 16:08:06 发布

stargazer.

最新推荐文章于 2024-01-11 16:08:06 发布

阅读量231

点赞数

分类专栏：生成树

本文链接：https://blog.csdn.net/qq_42555009/article/details/87875944

版权

生成树专栏收录该内容

4 篇文章 0 订阅

订阅专栏

传送门

板子题+1
套一个有向图矩阵树就完事了

#include<bits/stdc++.h>
using namespace std;
#define ll long long
inline int read(){
	char ch=getchar();
	int res=0,f=1;
	while(!isdigit(ch)){if(ch=='-')f=-f;ch=getchar();}
	while(isdigit(ch))res=res*10+(ch^48),ch=getchar();
	return res*f;
}
const int mod=10007;
int n,m,tr[310][310];
int change(){
	int ans=1,tmp=1;
	for(int i=2;i<=n;i++){
		for(int j=i+1;j<=n;j++){
			int x=i,y=j;
			while(tr[y][i]!=0){
				int t=tr[x][i]/tr[y][i];
				for(int k=i;k<=n;k++){
					tr[x][k]=((tr[x][k]-tr[y][k]*t)%mod+mod)%mod;
				}
				swap(x,y);
			}
			if(x!=i){
				for(int k=2;k<=n;k++){
					swap(tr[i][k],tr[x][k]);
				}
				tmp=-tmp;
			}
		}
		ans=((ans*tr[i][i])%mod+mod%mod);
	}
	ans=(ans*tmp%mod+mod)%mod;
	return ans;
}
int main(){
	n=read(),m=read();
	for(int i=1;i<=m;i++){
		int x=read(),y=read();
		tr[x][x]++,tr[y][x]--;
	}
	for(int i=2;i<=n;i++){
		for(int j=2;j<=n;j++){
			tr[i][j]=(tr[i][j]%mod+mod)%mod;
		}
	}
	cout<<change();
}