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

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

Time Limit: 30 Sec   Memory Limit: 162 MB
Submit: 4070   Solved: 1639
[ Submit][ Status][ Discuss]

Description

　　一个有向图G=(V,E)称为半连通的(Semi-Connected)，如果满足：?u,v∈V，满足u→v或v→u，即对于图中任意
两点u，v,存在一条u到v的有向路径或者从v到u的有向路径。若G'=(V',E')满足V'?V，E'是E中所有跟V'有关的边，
则称G'是G的一个导出子图。若G'是G的导出子图，且G'半连通，则称G'为G的半连通子图。若G'是G所有半连通子图
中包含节点数最多的，则称G'是G的最大半连通子图。给定一个有向图G，请求出G的最大半连通子图拥有的节点数K
，以及不同的最大半连通子图的数目C。由于C可能比较大，仅要求输出C对X的余数。

Input

　　第一行包含两个整数N，M，X。N，M分别表示图G的点数与边数，X的意义如上文所述接下来M行，每行两个正整
数a, b，表示一条有向边(a, b)。图中的每个点将编号为1,2,3…N，保证输入中同一个(a,b)不会出现两次。N ≤1
00000, M ≤1000000；对于100%的数据， X ≤10^8

Output

　　应包含两行，第一行包含一个整数K。第二行包含整数C Mod X.

Sample Input

6 6 20070603
1 2
2 1
1 3
2 4
5 6
6 4

Sample Output

3

3

思路：
先缩点，然后每个叶子节点到根的路径所经过的点就是一个
半联通图。拓扑排序从叶子节点更新整张图。
代码：

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+10;
int dfn[maxn],low[maxn],fa[maxn],instack[maxn];
vector<int>G[maxn],mp[maxn];
stack<int>P;
int n,m,mod,tol,num;
int hav[maxn],f[maxn],g[maxn],in[maxn],vis[maxn];
void tarjan(int v)
{
    dfn[v]=low[v]=++tol;
    P.push(v); instack[v]=1;
    for(int i=0;i<G[v].size();i++)
    {
        int to=G[v][i];
        if(!dfn[to])
        {
            tarjan(to);
            low[v]=min(low[v],low[to]);
        }
        else if(instack[to]) low[v]=min(low[v],dfn[to]);
    }
    if(dfn[v]==low[v])
    {
        num++;int t;
        do
        {
            t=P.top();P.pop();
            fa[t]=num;
            instack[t]=0;
            hav[num]++;
        }
        while(t!=v);
    }
}
void pp()
{
    queue<int>Q;
    for(int i=1;i<=num;i++)
    {
        if(!in[i])Q.push(i);
        f[i]=hav[i];g[i]=1;
    }
    while(!Q.empty())
    {
        int v=Q.front();Q.pop();
        for(int i=0;i<mp[v].size();i++)
        {
            int to=mp[v][i];
            in[to]--;
            if(in[to]==0)Q.push(to);
            if(vis[to]==v)continue;
            if(f[to]<f[v]+hav[to])
            {
                f[to]=f[v]+hav[to];
                g[to]=g[v];
            }
            else if(f[to]==f[v]+hav[to])
                g[to]=(g[to]+g[v])%mod;
            vis[to]=v;
        }
    }
}
int main()
{
    scanf("%d%d%d",&n,&m,&mod);
    while(m--)
    {
        int x,y;scanf("%d%d",&x,&y);
        G[x].push_back(y);
    }
    for(int i=1;i<=n;i++)
        if(!dfn[i]) tarjan(i);
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<G[i].size();j++)
        {
            int to=G[i][j];
            if(fa[i]!=fa[to])
                mp[fa[i]].push_back(fa[to]),in[fa[to]]++;
        }
    }
    pp();
    int ma=0,ans=0;
    for(int i=1;i<=num;i++)
    {
        if(f[i]>ma) ma=f[i],ans=g[i];
        else if(f[i]==ma) ans=(ans+g[i])%mod;
    }
    printf("%d %d\n",ma,ans);
    return 0 ;
}