LA 4287 Proving Equivalences(tarjan+强连通缩点)

最新推荐文章于 2019-09-26 18:17:22 发布

克莉丝汀娜

最新推荐文章于 2019-09-26 18:17:22 发布

阅读量454

点赞数

CC 4.0 BY-SA版权

分类专栏：图论

本文链接：https://blog.youkuaiyun.com/u014679804/article/details/46648899

图论专栏收录该内容

98 篇文章

订阅专栏

本文介绍了一个关于图论的问题，通过构建有向图并利用深度优先搜索求解强连通分量，来确定为了证明一组命题全部等价所需的最少额外推导次数。文章提供了完整的算法实现代码。

题意：n个命题，已经做了m次推导，问为了证明这n个命题全部等价，至少还需要几次推导。

分析：n个命题相当于n个顶点，m次推导相当于m条有向边。全部等价相当于图强连通。

求出图中的强连通分量，若只有一个，即图本身就是强连通，则无需推导，答案为0。

否则进行缩点，得到一个DAG图，分别统计该图中点的入度和出度总数a和b，答案就是max(a,b)。

#include<bits/stdc++.h>
using namespace std;
#define maxn 20005
vector<int> G[maxn];
int pre[maxn],low[maxn],sccno[maxn],dfs_clock,scc_cnt;
stack<int> S;

void dfs(int u){
    pre[u]=low[u]=++dfs_clock;
    S.push(u);
    for(int i=0;i<G[u].size();++i){
        int v=G[u][i];
        if(!pre[v]){
            dfs(v);
            low[u]=min(low[u],low[v]);//用后代的low函数更新自身
        }
        else if(!sccno[v]){
            low[u]=min(low[u],pre[v]);//用反向边更新
        }
    }
    if(low[u]==pre[u]){
        ++scc_cnt;
        for(;;){
            int x=S.top();S.pop();
            sccno[x]=scc_cnt;
            if(x==u) break;
        }
    }
}

void find_scc(int n){
    dfs_clock=scc_cnt=0;
    memset(sccno,0,sizeof(sccno));
    memset(pre,0,sizeof(pre));
    for(int i=0;i<n;++i)
        if(!pre[i]) dfs(i);
}
int in[maxn],out[maxn];
int main()
{
    int i,t,a,b,j,n,m;
    cin>>t;
    while(t--)
    {
        scanf("%d%d",&n,&m);
        for(i=0;i<n;++i) G[i].clear();
        for(i=1;i<=m;++i){
            scanf("%d%d",&a,&b);
            --a,--b;
            G[a].push_back(b);
        }
        find_scc(n);
        if(scc_cnt==1) {puts("0");continue;}
        for(i=1;i<=scc_cnt;++i) in[i]=out[i]=0;
        for(i=0;i<n;++i){
            for(j=0;j<G[i].size();++j)
            {
                int v=G[i][j];
                if(sccno[v]!=sccno[i]) out[sccno[i]]=in[sccno[v]]=1;
            }
        }
        a=b=0;
        for(i=1;i<=scc_cnt;++i)
        {
            if(!in[i])++a;
            if(!out[i]) ++b;
        }
        printf("%d\n",max(a,b));
    }
    return 0;
}