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

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

HINT

Source

容易推出结论 将强连通分量缩点，要求的就是最长链（根据size计算长度）

于是就是一个强连通分量的问题了

基本就是个tarjan的板，写熟了就好了

#include<bits/stdc++.h>
 
using namespace std;
 
const int maxn=100010;
 
int n,m,mod;
 
struct Edge
{
    int to,nxt;
}e[maxn*10],edge[maxn*10];
 
int low[maxn],dfn[maxn],head[maxn];
 
int st[maxn],top,scc[maxn];
 
int f[maxn],g[maxn];
 
int hhead[maxn],cnt;
 
inline void addedge(int x,int y)
{
    edge[++cnt].to=y;
    edge[cnt].nxt=hhead[x];
    hhead[x]=cnt;
}
 
int tot,num,size[maxn];
 
void dfs(int x)
{
    int tmp=0;
    low[x]=dfn[x]=++tot;
    st[++top]=x;
    for(int i=hhead[x];i;i=edge[i].nxt)
    {
        int y=edge[i].to;
        if(!dfn[y])
        {
            dfs(y);
            low[x]=min(low[x],low[y]);
        }
        else if(!scc[y]) low[x]=min(low[x],dfn[y]);
    }
    if(low[x]==dfn[x])
    {
        num++;
        while(tmp!=x)
        {
            tmp=st[top--];
            scc[tmp]=num;
            size[num]++;
        }
    }
}
 
map <pair<int ,int >,bool > mp;
 
int d[maxn];
 
inline void rebuild()
{
    cnt=0;
    for(int i=1;i<=n;i++)
    {
        for(int j=hhead[i];j;j=edge[j].nxt)
        {
            int x=scc[i],y=scc[edge[j].to];
            if(x!=y&&!mp[make_pair(x,y)])
            {
                d[y]++;
                e[++cnt].to=y;
                e[cnt].nxt=head[x];
                head[x]=cnt;
                mp[make_pair(x,y)]=true;
            }
        }
    }
}
 
int s[maxn];
 
int main()  
{
    scanf("%d%d%d",&n,&m,&mod);
    for(int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        addedge(x,y);
    }
    for(int i=1;i<=n;i++) if(!dfn[i]) dfs(i);
    rebuild();
    queue <int> q;
    for(int i=1;i<=num;i++)
    {
        if(!d[i])
        {
            q.push(i);
            g[i]=size[i];
            s[i]=1;
        }
    }
    tot=0;
    while(!q.empty())
    {
        int tmp=q.front();q.pop();
        f[++tot]=tmp;
        for(int i=head[tmp];i;i=e[i].nxt)
        {
            int y=e[i].to;
            d[y]--;
            if(!d[y]) q.push(y);
        }
    }
    for(int i=1;i<=num;i++)
    {
        int x=f[i];
        for(int j=head[x];j;j=e[j].nxt)
        {
            int y=e[j].to;
            g[y]=max(g[y],g[x]+size[y]);
        }
    }
    for(int i=1;i<=num;i++)
    {
        int x=f[i];
        for(int j=head[x];j;j=e[j].nxt)
        {
            int y=e[j].to;
            if(g[y]==g[x]+size[y]) (s[y]+=s[x])%=mod;
        }
    }
    int mx=0;
    for(int i=1;i<=num;i++) mx=max(mx,g[i]);
    int ans=0;
    for(int i=1;i<=num;i++) if(g[i]==mx) (ans+=s[i])%=mod;
    printf("%d\n%d\n",mx,ans);
}