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

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

Time Limit: 30 Sec  Memory Limit: 162 MB
Submit: 1986  Solved: 802
[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。

Source

题解：缩点dp。注意处理好大点间的重边。

你问我写哪种dp？首推topo啊！记忆化什么的常数大的。。。（雾

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<stack>
#include<queue>
#include<cstring>
#define PAU putchar(' ')
#define ENT putchar('
')
#define MSE(a,b) memset(a,b,sizeof(a))
#define REN(x) for(ted*e=fch[x];e;e=e->nxt)
#define REN2(x) for(ted*e=fch2[x];e;e=e->nxt)
#define TIL(x) for(int i=1;i<=x;i++)
#define ALL(x) for(int j=1;j<=x;j++)
using namespace std;
const int maxn=100000+10,maxm=1000000+10,inf=1e9;
struct ted{int x,y;ted*nxt;}adj[maxm],adj2[maxm],*fch[maxn],*fch2[maxn],*ms=adj,*ms2=adj2;
void add(int x,int y){*ms=(ted){x,y,fch[x]};fch[x]=ms++;return;}int in[maxn];
void add2(int x,int y){*ms2=(ted){x,y,fch2[x]};fch2[x]=ms2++;in[y]++;return;}
int dfn[maxn],low[maxn],beg[maxn],siz[maxn],scc,cz,f[maxn],g[maxn],vis[maxn],mod;bool ins[maxn];stack<int>S;
void tarjan(int x){
    dfn[x]=low[x]=++cz;ins[x]=true;S.push(x);REN(x){int v=e->y;
        if(!dfn[v])tarjan(v),low[x]=min(low[x],low[v]);else if(ins[v])low[x]=min(low[x],dfn[v]);
    }if(low[x]==dfn[x]){scc++;int t;
        do beg[t=S.top()]=scc,ins[t]=false,S.pop(),siz[scc]++;while(t!=x);
    }return;
}
int n,m;
void rebuild(){
    TIL(n)REN(i)if(beg[i]!=beg[e->y])add2(beg[i],beg[e->y]);return;
}
void topodp(){
    queue<int>Q;TIL(scc)if(!in[i])Q.push(i);TIL(scc)f[i]=siz[i],g[i]=1;
    while(!Q.empty()){
        int x=Q.front();Q.pop();REN2(x){int v=e->y;
            if(--in[v]==0)Q.push(v);
            if(vis[v]!=x){vis[v]=x;
                if(f[x]+siz[v]>f[v])f[v]=f[x]+siz[v],g[v]=g[x];
                else if(f[x]+siz[v]==f[v])(g[v]+=g[x])%=mod;
            }
        }
    }return;
}
inline int read(){
    int x=0;bool sig=true;char ch=getchar();
    for(;!isdigit(ch);ch=getchar())if(ch=='-')sig=false;
    for(;isdigit(ch);ch=getchar())x=10*x+ch-'0';return sig?x:-x;
}
inline void write(int x){
    if(x==0){putchar('0');return;}if(x<0)putchar('-'),x=-x;
    int len=0;static int buf[20];while(x)buf[len++]=x%10,x/=10;
    for(int i=len-1;i>=0;i--)putchar(buf[i]+'0');return;
}
int main(){
    n=read();m=read();mod=read();int x,y;
    TIL(m)x=read(),y=read(),add(x,y);
    TIL(n)if(!dfn[i])tarjan(i);rebuild();topodp();
    int mx=-inf,ans;
    TIL(scc){
        if(mx<f[i])mx=f[i],ans=g[i];
        else if(f[i]==mx)(ans+=g[i])%=mod;
    }write(mx);ENT;write(ans);
    return 0;
}