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

最新推荐文章于 2019-06-27 21:51:20 发布

puck_just_me

最新推荐文章于 2019-06-27 21:51:20 发布

阅读量498

点赞数

分类专栏：连通分量

本文链接：https://blog.csdn.net/Doom_bringer/article/details/50625896

版权

连通分量专栏收录该内容

2 篇文章 0 订阅

订阅专栏

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 1000001

queue<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;
}