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

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一下即可

code：

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
#define maxn 100005
#define maxm 1000005
using namespace std;
char ch;
bool ok;
void read(int &x){
    for (ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if (ch=='-') ok=1;
    for (x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar());
    if (ok) x=-x;
}
int n,m,mod,a,b;
int idx,dfn[maxn],low[maxn],top,stack[maxn],cnt,belong[maxn],sum[maxn];
bool in[maxn];
int num;
struct Edge{
    int u,v;
    inline bool operator==(Edge b){return u==b.u&&v==b.v;}
}edge[maxm];
inline bool cmp(Edge a,Edge b){
    if (a.u==b.u) return a.v<b.v;
    return a.u<b.u;
}
int head,tail,list[maxn],f[maxn],g[maxn],deg[maxn];
struct Graph{
    int tot,now[maxn],pre[maxm],son[maxm];
    void put(int a,int b){pre[++tot]=now[a],now[a]=tot,son[tot]=b;}
    void dfs(int u){
        dfn[u]=low[u]=++idx,stack[++top]=u,in[u]=1;
        for (int p=now[u],v=son[p];p;p=pre[p],v=son[p])
            if (!dfn[v]) dfs(v),low[u]=min(low[u],low[v]);
            else if (in[v]) low[u]=min(low[u],dfn[v]);
        if (dfn[u]==low[u]){
            int v; ++cnt;
            do{v=stack[top--],in[v]=0,belong[v]=cnt,sum[cnt]++;}while(v!=u);
        }
    }
    void get(){
        for (int u=1;u<=n;u++)
            for (int p=now[u],v=son[p];p;p=pre[p],v=son[p])
                if (belong[u]!=belong[v]) edge[++num]=(Edge){belong[u],belong[v]};
        sort(edge+1,edge+num+1,cmp),num=unique(edge+1,edge+num+1)-edge-1;
    }
    void calc(){
        head=0,tail=0;
        for (int i=1;i<=cnt;i++) if (!deg[i]) list[++tail]=i,g[i]=1;
        while (head<tail){
            int u=list[++head];
            f[u]+=sum[u];
            for (int p=now[u],v=son[p];p;p=pre[p],v=son[p]){
                if (f[v]<f[u]) f[v]=f[u],g[v]=g[u];
                else if (f[v]==f[u]) g[v]+=g[u],g[v]%=mod;
                deg[v]--;
                if (!deg[v]) list[++tail]=v;
            }
        }
        int ans1=0,ans2=0;
        for (int i=1;i<=cnt;i++) ans1=max(ans1,f[i]);
        printf("%d\n",ans1);
        for (int i=1;i<=cnt;i++) if (ans1==f[i]) ans2+=g[i],ans2%=mod;
        printf("%d\n",ans2);
    }
}G1,G2;
int main(){
    read(n),read(m),read(mod);
    for (int i=1;i<=m;i++) read(a),read(b),G1.put(a,b);
    for (int i=1;i<=n;i++) if (!dfn[i]) G1.dfs(i);
    G1.get();
    for (int i=1;i<=num;i++) G2.put(edge[i].u,edge[i].v),deg[edge[i].v]++;
    G2.calc();
    return 0;
}