BZOJ 1093 tajan

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

Time Limit: 30 Sec  Memory Limit: 162 MB
Submit: 3363  Solved: 1324
[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

 

//http://hzwer.com/4633.html

 

#include<iostream>
#include<cstdio>
#include<cstdlib>
using namespace std;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int mx,ans;
int ind,cnt,scc,top;
int n,m,X;
int last[100005],last2[100005];
int dfn[100005],low[100005],hav[100005],belong[100005],q[100005];
int r[100005],f[100005],g[100005],vis[100005];
bool inq[100005];
struct edge{int to,next;}e[2000005],ed[2000005];
void insert(int u,int v)
{
    e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;
}
void ins(int u,int v)
{
    ed[++cnt].to=v;ed[cnt].next=last2[u];last2[u]=cnt;
    r[v]++;
}
void tarjan(int x)
{
    dfn[x]=low[x]=++ind;
    q[++top]=x;inq[x]=1;
    for(int i=last[x];i;i=e[i].next)
        if(!dfn[e[i].to])
            tarjan(e[i].to),low[x]=min(low[x],low[e[i].to]);
        else if(inq[e[i].to])low[x]=min(low[x],dfn[e[i].to]);
    int now=0;
    if(low[x]==dfn[x])
    {
        scc++;
        while(now!=x)
        {
            now=q[top];top--;
            inq[now]=0;
            hav[scc]++;
            belong[now]=scc;
        }
    }
}
void rebuild()
{
    cnt=0;
    for(int x=1;x<=n;x++)
    {
        for(int i=last[x];i;i=e[i].next)
            if(belong[x]!=belong[e[i].to])
                ins(belong[x],belong[e[i].to]);
    }
}
void dp()
{
    int head=0,tail=0;
    for(int i=1;i<=scc;i++)
 {
  if(!r[i])q[tail++]=i;
  f[i]=hav[i];g[i]=1;
 }
    while(head!=tail)
    {
        int now=q[head];head++;
        for(int i=last2[now];i;i=ed[i].next)
  {
   r[ed[i].to]--;
   if(!r[ed[i].to])q[tail++]=ed[i].to;
   if(vis[ed[i].to]==now)continue;//重构边的时候可能会重复加边 例如加了两次1->2
   if(f[now]+hav[ed[i].to]>f[ed[i].to])
   {
    f[ed[i].to]=f[now]+hav[ed[i].to];
    g[ed[i].to]=g[now];
   }
   else if(f[now]+hav[ed[i].to]==f[ed[i].to])
    g[ed[i].to]=(g[ed[i].to]+g[now])%X;
   vis[ed[i].to]=now;
  }
    }
}
int main()
{
    n=read();m=read();X=read();
    for(int i=1;i<=m;i++)
    {
        int 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]>mx)mx=f[i],ans=g[i];
  else if(f[i]==mx)ans=(ans+g[i])%X;
 }
    printf("%d\n%d\n",mx,ans);
    return 0;
}