[bzoj1093][ZJOI2007]最大半连通子图

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

Time Limit: 30 Sec Memory Limit: 162 MB
Submit: 2028 Solved: 819
[Submit][Status][Discuss]
Description

Input

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

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

对于100%的数据， N ≤100000, M ≤1000000；对于100%的数据， X ≤10^8。

首先先求一下强联通分量缩一下点，把图变成一个DAG。
然后这个题就变成了在DAG上找一条最长路和最长路条数。
这样我们做一遍拓扑，在拓扑的时候dp一下就好了（其实我感觉就是一个类似于SPFA的东西），路径的条数就在dp更新的时候顺便处理一下就行了。
tips:还要注意一下重边的问题。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=100010;
int n,m,D,point_pre[N],next_pre[N*10],pre[N],lowlink[N],sccno[N],sccn,tot,dfsn,stack[N];
int point_now[N],next_now[N*10],v[N],ru[N],l[N*10],h,t,sum[N],f[N],use[N],maxn,ans,top;
struct S{int st,en;}aa[N*10],bb[N*10];
void add_pre(int x,int y)
{
    tot+=1;next_pre[tot]=point_pre[x];point_pre[x]=tot;
    aa[tot].st=x;aa[tot].en=y;
}
void add_now(int x,int y)
{
    tot+=1;next_now[tot]=point_now[x];point_now[x]=tot;
    bb[tot].st=x;bb[tot].en=y;
}
void dfs(int x)
{
    int i,u;
    pre[x]=lowlink[x]=++dfsn;
    stack[++top]=x;
    for(i=point_pre[x];i;i=next_pre[i]){
        u=aa[i].en;
        if(!pre[u]){
            dfs(u);
            lowlink[x]=min(lowlink[u],lowlink[x]);
        }
        else if(!sccno[u]) lowlink[x]=min(lowlink[x],pre[u]);
    }
    if(lowlink[x]==pre[x]){
        sccn+=1;
        while(1){
            u=stack[top];
            top-=1;
            sccno[u]=sccn;
            if(x==u) break;
        }
    }
}
int main()
{
    int i,j,x,y,u;
    scanf("%d%d%d",&n,&m,&D);
    for(i=1;i<=m;++i){
        scanf("%d%d",&x,&y);
        add_pre(x,y);
    }
    for(i=1;i<=n;++i)
      if(!pre[i])
        dfs(i);
    tot=0;
    for(i=1;i<=n;++i){
        v[sccno[i]]+=1;
        for(j=point_pre[i];j;j=next_pre[j])
          if(sccno[i]!=sccno[aa[j].en]){
            add_now(sccno[i],sccno[aa[j].en]);
            ru[sccno[aa[j].en]]+=1;
          }
    }
    h=0;t=1;
    for(i=1;i<=sccn;++i)
      if(!ru[i]){
        t+=1;
        l[t]=i;
        f[i]=v[i];
        sum[i]=1;
      }
    while(h<=t){
        u=l[h];
        for(i=point_now[u];i;i=next_now[i]){
            ru[bb[i].en]-=1;
            if(!ru[bb[i].en]) t+=1,l[t]=bb[i].en;
            if(use[bb[i].en]==u) continue;
            if(f[bb[i].en]==f[u]+v[bb[i].en]) sum[bb[i].en]=(sum[bb[i].en]+sum[u])%D;
            else if(f[bb[i].en]<f[u]+v[bb[i].en]){
                f[bb[i].en]=f[u]+v[bb[i].en];
                sum[bb[i].en]=sum[u];
            }
            use[bb[i].en]=u;
        }
        h+=1;
    }
    for(i=1;i<=sccn;++i){
        if(f[i]==maxn) ans=(ans+sum[i])%D;
        if(f[i]>maxn){
            maxn=f[i];
            ans=sum[i]; 
        }
    }
    printf("%d\n%d\n",maxn,ans);
}