最大半连通子图(有向图缩点)

先tarjan缩点 再在拓扑图上找一条最长链

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<set>
using namespace std;
#define ll long long
const int N=1e5+10;
const int M=1e6+10;
int h[N],e[M],ne[M],idx;
int hc[N],ec[M],nec[M],idxc;
int low[N],dfn[N],c,st[N],top,is[N],Size[N],num[N],cnt;
int f[N],g[N];
void add(int a,int b)
{
	e[idx]=b;
	ne[idx]=h[a];
	h[a]=idx++;
}
void addc(int a,int b)
{
	ec[idxc]=b;
	nec[idxc]=hc[a];
	hc[a]=idxc++;
}
void dfs(int x)
{
	low[x]=dfn[x]=++c;
	st[++top]=x;is[x]=1;
	for(int i=h[x];i!=-1;i=ne[i])
	{
		int j=e[i];
		if(dfn[j]==0)
		{
			dfs(j);
			low[x]=min(low[x],low[j]);
		}
		else if(is[j])
		{
			low[x]=min(low[x],dfn[j]);
		}
	}
	if(dfn[x]==low[x])
	{
		int y;
		cnt++;
		do{
			y=st[top--];
			is[y]=0;
			num[y]=cnt;
			Size[cnt]++;
		}while(x!=y);
	}
}
int main()
{
	int n,m,x;
	scanf("%d%d%d",&n,&m,&x);
	memset(h,-1,sizeof h);
	memset(hc,-1,sizeof hc);
	for(int i=1;i<=m;i++)	
	{
		int a,b;
		scanf("%d%d",&a,&b);
		add(a,b);
	}
	for(int i=1;i<=n;i++)
	{
		if(!dfn[i])
		{
			dfs(i);
		}
	}
	set<ll>mp;
	for(int i=1;i<=n;i++)
	{
		for(int j=h[i];j!=-1;j=ne[j])
		{
			if(num[i]!=num[e[j]])
			{
				int a=num[i],b=num[e[j]];
				ll ans=a*1000000ll+b;//要计算数量所以不能有重边
				if(!mp.count(ans))
				addc(num[i],num[e[j]]);
				mp.insert(ans);
			}
		}
	}
	for(int i=cnt;i>=1;i--)
	{
		if(!f[i])
		{
			f[i]=Size[i];
			g[i]=1;
		}
		for(int j=hc[i];j!=-1;j=nec[j])
		{
			int k=ec[j];
			if(f[k]<f[i]+Size[k])
			{
				f[k]=f[i]+Size[k];
				g[k]=g[i];
			}
			else if(f[k]==f[i]+Size[k]) g[k]=(g[k]+g[i])%x;
		}
		
	}
	int maxx=0,sum=0;
	for(int i=1;i<=cnt;i++)
	{
		if(f[i]>maxx)
		{
			maxx=f[i];
			sum=g[i];
		}
		else if(f[i]==maxx)
		{
			sum=(sum+g[i])%x;
		}
	}
	printf("%d\n%d\n",maxx,sum);
}