1093: [ZJOI2007]最大半连通子图 tarjan+拓扑最长链

首先一个强连通分量一定是一个半联通子图，所以我们可以先tarjan缩一下环。然后在得到的图中跑最长链即可。注意重边，用的map还是不慢的。。

#include<iostream>
#include<cstdio>
#include<map>
#define N 100005
using namespace std;
int head[N],head0[N],dfn[N],low[N],belong[N],num[N],ind[N],f[N],g[N],q[N],stack[N];
int n,m,M,scc,tot,cnt,top,ans1,ans2;
map <int,int> mp[100005];
bool inset[N];
int next[2000005],list[2000005],next0[2000005],list0[2000005];
inline int read()
{
    int a=0,f=1; char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
    return a*f;
}
inline void insert(int x,int y)
{
    next[++cnt]=head[x];
    head[x]=cnt;
    list[cnt]=y;
}
inline void insert0(int x,int y)
{
    next0[++cnt]=head0[x];
    head0[x]=cnt;
    list0[cnt]=y;
    ind[y]++;
}
void tarjan(int x)
{
    dfn[x]=low[x]=++tot;
    stack[++top]=x;
    inset[x]=1;
    for (int i=head[x];i;i=next[i])
        if (!dfn[list[i]])
        {
            tarjan(list[i]);
            low[x]=min(low[x],low[list[i]]);
        }
        else if (inset[list[i]]) low[x]=min(low[x],dfn[list[i]]);
    if (low[x]==dfn[x])
    {
        int i=-1; scc++;
        while (i!=x)
        {
            i=stack[top--];
            inset[i]=0;
            belong[i]=scc;
            num[scc]++;
        }
    }
}
inline void rebuild()
{
    cnt=0;
    for (int x=1;x<=n;x++)
        for (int i=head[x];i;i=next[i])
            if (belong[x]!=belong[list[i]]&&!mp[belong[x]][belong[list[i]]])
            {
                mp[belong[x]][belong[list[i]]]=1;
                insert0(belong[x],belong[list[i]]);
            }
}
inline void dp()
{
    int t=0,w=0;
    for (int i=1;i<=scc;i++)
    {
        if (!ind[i]) q[++w]=i;
        f[i]=num[i]; g[i]=1;
    }
    while (t!=w)
    {
        int x=q[++t];
        for (int i=head0[x];i;i=next0[i])
        {
            ind[list0[i]]--;
            if (!ind[list0[i]]) q[++w]=list0[i];
            if (f[x]+num[list0[i]]>f[list0[i]])
            {
                f[list0[i]]=f[x]+num[list0[i]];
                g[list0[i]]=g[x];
            }
            else if (f[x]+num[list0[i]]==f[list0[i]]) 
                g[list0[i]]=(g[list0[i]]+g[x])%M;
        }
    }
}
int main()
{
    n=read(); m=read(); M=read();
    for (int i=1,u,v;i<=m;i++)
        u=read(),v=read(),insert(u,v);
    for (int i=1;i<=n;i++) if (!dfn[i]) tarjan(i);
    rebuild();
    dp();
    for (int i=1;i<=scc;i++)
    {
        if (f[i]>ans1) ans1=f[i],ans2=g[i];
        else if (f[i]==ans1) ans2=(ans2+g[i])%M;
    }
    cout << ans1 << endl << ans2 << endl;
    return 0;
}