【连通分量】BZOJ 1093: [ZJOI2007]最大半连通子图

来源:互联网 发布:北京美工培训班学费 编辑:程序博客网 时间:2024/03/29 13:29

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


Description

每一头牛的愿望就是变成一头最受欢迎的牛。现在有N头牛,给你M对整数(A,B),表示牛A认为牛B受欢迎。 这种关系是具有传递性的,如果A认为B受欢迎,B认为C受欢迎,那么牛A也认为牛C受欢迎。你的任务是求出有多少头牛被所有的牛认为是受欢迎的。


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。


Solution

tarjan缩点
然后topo dp
维护一个dis 一个count

Code

#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#include<stack>#include<queue>using namespace std;#define maxn 1000001queue<int>q;stack<int>s;struct edge{    int to,lst;}e[maxn];int last[maxn],tot,n,m,low[maxn],dfn[maxn],Tot,num[maxn],size,f[maxn],mod,in[maxn];bool ins[maxn];void add(int u,int v){    e[++tot]=(edge){v,last[u]};last[u]=tot;}void tarjan(int poi){    dfn[poi]=low[poi]=++Tot;    s.push(poi);ins[poi]=1;    for(int i=last[poi];i;i=e[i].lst)    {        int u=e[i].to;        if(!dfn[u])        {            tarjan(u);            low[poi]=min(low[poi],low[u]);        }        else if(ins[u])low[poi]=min(low[poi],low[u]);    }    if(dfn[poi]==low[poi])    {        size++;        int u;        do{            u=s.top();            s.pop();            ins[u]=0;            f[u]=size;            num[size]++;        }while(poi!=u);    }}struct ED{    int u,v;}E[maxn];bool cmp(const ED A,const ED B){    return f[A.u]==f[B.u]?f[A.v]<f[B.v]:f[A.u]<f[B.u];}int top=0,F[maxn],cot[maxn];void topo_dp(){    for(int i=1;i<=size;i++)if(!in[i])q.push(i),F[i]=num[i],cot[i]=1;    while(q.size())    {        int u=q.front();q.pop();        for(int i=last[u];i;i=e[i].lst)        {            int v=e[i].to;            in[v]--;            if(!in[v])q.push(v);            if(F[v]<F[u]+num[v])F[v]=F[u]+num[v],cot[v]=cot[u];            else if(F[v]==F[u]+num[v])cot[v]=(cot[v]+cot[u])%mod;        }    }    int mx=0,sum=0;    for(int i=1;i<=size;i++)    {        if(F[i]>F[mx])sum=cot[i]%mod,mx=i;        else if(F[i]==F[mx])sum=(cot[i]+sum)%mod;    }    printf("%d\n%d",F[mx],sum);}int main(){    freopen("semi.in","r",stdin);    freopen("semi.out","w",stdout);    scanf("%d%d%d",&n,&m,&mod);    for(int i=1;i<=m;i++)    {        int u,v;        scanf("%d%d",&u,&v);        add(u,v);        E[i]=(ED){u,v};    }    for(int i=1;i<=n;i++)if(!dfn[i])tarjan(i);    sort(1+E,1+E+m,cmp);    memset(last,0,sizeof(last));tot=0;    for(int i=1;i<=m;i++)    {        int A=f[E[i].u],B=f[E[i].v];        if(A==f[E[i-1].u]&&B==f[E[i-1].v])continue;        if(A==B)continue;        add(A,B);        in[B]++;    }    topo_dp();    return 0;}
0 0
原创粉丝点击