bzoj1093(强联通分量+缩点+dp记忆化搜索（算是乱搞找方法吧！））

先强连通缩点，然后重新构图，然后找出包含点数最多的链，统计个数即可，可以用拓扑排序搞~

 

pS：重新构图时有重边，然后导致统计方案数的重复

注意缩点构图时，重边的情况处理

#include<cstdio>
#include<algorithm>
#include<cstring>

using namespace std;
const int N=100009;
struct aa
{
	int to,w,pre;
}edge[2][N*20];
int top,head[2][N],n,m,X,cnt[N],low[N],dfn[N],sta[N],tot[2],id,ans,size[N],bj[N],mx[N],tt[N];
bool s[N];
void add(int h,int x,int y)
{
	tot[h]++;
	edge[h][tot[h]].to=y;
	edge[h][tot[h]].pre=head[h][x];
	head[h][x]=tot[h];
}
void dfs(int u)
{
	dfn[u]=low[u]=++id;
	s[u]=true;
	sta[++top]=u;
	for (int i=head[0][u];i;i=edge[0][i].pre)
	{
		int v=edge[0][i].to;
		if (!dfn[v]) dfs(v),low[u]=min(low[u],low[v]);
		else if (s[v]) low[u]=min(low[u],low[v]);
	}
	if (dfn[u]==low[u])
	{
		ans++;
		while (sta[top+1]!=u)
		{
			int v=sta[top];
			s[v]=false;
			bj[v]=ans;
			size[ans]++;
			top--;
		}//求出所有联通分量
	}
}
void work()
{
	for (int i=1;i<=n;i++) if (!dfn[i]) dfs(i);
}
bool pan(int u,int v)//处理重边！！
{
	for (int i=head[1][u];i;i=edge[1][i].pre)
	if (edge[1][i].to==v) return false;
	return true;
}
void suodian()//缩点
{
	for (int i=1;i<=n;i++)
	for (int j=head[0][i];j;j=edge[0][j].pre)
	{
		int v=edge[0][j].to;
		if (bj[i]!=bj[v]&&pan(bj[i],bj[v])) add(1,bj[i],bj[v]),cnt[bj[v]]++;
	}
}
void dfs2(int u)
{
	if (mx[u]) return ; //记忆化，如果之前已经算过，直接返回就行，第一遍时没加这个就超时了
	mx[u]=size[u]%X;
	tt[u]=1;
	for (int j=head[1][u];j;j=edge[1][j].pre)
	{
		int v=edge[1][j].to;
		dfs2(v);
		if (mx[v]+size[u]==mx[u]) tt[u]=(tt[u]+tt[v])%X;
		if (mx[v]+size[u]>mx[u]) mx[u]=mx[v]+size[u],tt[u]=tt[v];//类似于树形dp
	}
}
int main()
{
	scanf("%d%d%d",&n,&m,&X);
	int x,y;
	for (int i=1;i<=m;i++) scanf("%d%d",&x,&y),add(0,x,y);
	work();
	suodian();
	int k=0,z=0;
	for (int i=1;i<=ans;i++)
	if (cnt[i]==0)//这个双联通分量起点一定是没有入度的否则可以更大
	{
		dfs2(i);
		if (mx[i]==k) z=(z+tt[i])%X;
		if (mx[i]>k) k=mx[i],z=tt[i];		
	}
	printf("%d\n%d",k,z);
	return 0;
}