Proving Equivalences UVALive - 4287（强连通分量水题）

强连通分量与入出度统计

最新推荐文章于 2021-10-14 21:04:13 发布

转载最新推荐文章于 2021-10-14 21:04:13 发布 · 49 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/WTSRUVF/p/9361303.html

本文介绍了一种算法，用于统计图中各强连通分量的入度和出度，最终输出入度或出度为0的分量数量。如果图中仅有一个强连通分量，则输出0。该算法使用了深度优先搜索和栈来识别强连通分量。

就是统计入度为0 的点和出度为0 的点输出大的那一个，，

若图中只有一个强连通分量则输出0即可

和https://www.cnblogs.com/WTSRUVF/p/9301096.html 这题差不多 poj1236

#include <iostream>
#include <cstdio>
#include <sstream>
#include <cstring>
#include <map>
#include <set>
#include <vector>
#include <stack>
#include <queue>
#include <algorithm>
#include <cmath>
#define rap(a, n) for(int i=a; i<=n; i++)
#define MOD 2018
#define LL long long
#define ULL unsigned long long
#define Pair pair<int, int>
#define mem(a, b) memset(a, b, sizeof(a))
#define _  ios_base::sync_with_stdio(0),cin.tie(0)
//freopen("1.txt", "r", stdin);
using namespace std;
const int maxn = 20005, INF = 0x7fffffff;
vector<int> G[maxn];
int pre[maxn], lowlink[maxn], sccno[maxn], dfs_clock, scc_cnt, in[maxn], out[maxn];
stack<int> s;

void dfs(int u)
{
    pre[u] = lowlink[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);
            lowlink[u] = min(lowlink[u], lowlink[v]);
        }
        else if(!sccno[v])
            lowlink[u] = min(lowlink[u], pre[v]);
    }
    if(lowlink[u] == pre[u])
    {
        scc_cnt++;
        for(;;)
        {
            int x = s.top(); s.pop();
            sccno[x] = scc_cnt;
            if(x == u) break;
        }
    }
}

void init()
{
    dfs_clock = scc_cnt = 0;
    mem(sccno, 0);
    mem(pre, 0);
    mem(out, 0);
    mem(in, 0);
    for(int i=0; i<maxn; i++) G[i].clear();
}


int main()
{
    int T, n, m;
    scanf("%d", &T);
    while(T--)
    {
        init();
        scanf("%d%d", &n, &m);
        rap(1, m)
        {
            int u, v;
            scanf("%d%d", &u, &v);
            G[u].push_back(v);
   //         G[v].push_back(u);
        }
        rap(1, n)
            if(!pre[i])
                dfs(i);
        rap(1, n)
        {
            for(int j=0; j<G[i].size(); j++)
            {
                if(sccno[i] != sccno[G[i][j]])
                    out[sccno[i]]++, in[sccno[G[i][j]]]++;
            }
        }
        if(scc_cnt == 1)
        {
            cout<< '0' <<endl;
            continue;
        }
        int a = 0, b = 0;
        rap(1, scc_cnt)
        {
            if(!in[i])
                a++;
            if(!out[i])
                b++;
        }
        cout<< max(a, b) <<endl;

    }

    return 0;
}