BZOJ 1093 [ZJOI2007]最大半连通子图

以前做过poj的一个判断图是否为弱连通的题，然后，这个题和poj那个差不多。

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

pS：重新构图时有重边，然后导致统计方案数的重复。。wa了好久。。还是wzc神犇告诉我这个蒟蒻的。。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
 
#define N 200000
#define M 5000000
#define BUG system("pause")
 
using namespace std;
 
int head[N],to[M],next[M];
int dfn[N],low[N];
int st[M],ed[M];
int n,m,cnt,ans,ansnum,mod;
int divg,belong[N],t,p,stk[N],val[N];
bool fg[N];
int num[N],in[N],dp[N],q[M],vis[N];
 
inline void add(int u,int v)
{
    to[cnt]=v; next[cnt]=head[u]; head[u]=cnt++;
}
 
inline void read()
{
    memset(head,-1,sizeof head); cnt=0;
    scanf("%d%d%d",&n,&m,&mod);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&st[i],&ed[i]);
        add(st[i],ed[i]);
    }
}
 
inline void dfs(int u)
{
    low[u]=dfn[u]=++t;
    stk[++p]=u; fg[u]=true;
    for(int i=head[u];~i;i=next[i])
    {
        if(!dfn[to[i]])
        {
            dfs(to[i]);
            low[u]=min(low[u],low[to[i]]);
        }
        else if(fg[to[i]]) low[u]=min(low[u],dfn[to[i]]);
    }
    if(dfn[u]==low[u])
    {
        divg++;
        int tmp=-1;
        while(tmp!=u)
        {
            tmp=stk[p--];
            belong[tmp]=divg;
            val[divg]++;
            fg[tmp]=false;
        }
    }
}
 
inline void topsort()
{
    int h=1,t=1,u;
    for(int i=1;i<=divg;i++)
        if(in[i]==0)
        {
            q[t++]=i;
            dp[i]=val[i];
            num[i]=1;
        }
    while(h<t)
    {
        u=q[h++];
        for(int i=head[u];~i;i=next[i])
        {
            in[to[i]]--;
            if(in[to[i]]==0) q[t++]=to[i];
            if(vis[to[i]]==u) continue;//有重边!! 
            if(dp[to[i]]<dp[u]+val[to[i]])
            {
                dp[to[i]]=dp[u]+val[to[i]];
                num[to[i]]=num[u];
            }
            else if(dp[to[i]]==dp[u]+val[to[i]])
            {
                num[to[i]]=(num[to[i]]+num[u])%mod;
            }
            vis[to[i]]=u;
        }
    }
}
 
inline void go()
{
    for(int i=1;i<=n;i++)
        if(!dfn[i]) dfs(i);
    memset(head,-1,sizeof head); cnt=0;
    for(int i=1;i<=m;i++)
        if(belong[st[i]]!=belong[ed[i]])
        {
            add(belong[st[i]],belong[ed[i]]);
            in[belong[ed[i]]]++;
        }
    topsort();
    for(int i=1;i<=divg;i++)
    {
        if(dp[i]>ans)
        {
            ans=dp[i];
            ansnum=num[i];
        }
        else if(dp[i]==ans) ansnum=(ansnum+num[i])%mod;
    }
    printf("%d\n%d\n",ans,ansnum);
}
 
int main()
{
    read();
    go();
    return 0;
}