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,m,X,cnt
,low
,dfn
,sta
,tot[2],id,ans,size
,bj
,mx
,tt
;
bool s
;
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;
}