UVA 11324 强连通分量 最长链

UVA 11324

题目链接：

https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=25&page=show_problem&problem=2299

题意：

给一个有向图，问最大的子图使得任意两点有路径可达。

思路：

1A啊！

求强连通缩点后具有最多结点的链即可。

源码：

#include <cstdio>

#include <cstring>

#include <cmath>

#include <cstdlib>

#include <algorithm>

#include <iostream>

#include <queue>

#include <vector>

#include <stack>

using namespace std;

#define gmax(a,b) ((a) > (b) ? (a) : (b))

#define gmin(a,b) ((a) < (b) ? (a) : (b))

const int MAXN = 1000 + 5;

int pre[MAXN], low[MAXN], dfs_clock;

stack<int>sta;

queue<int>que;

vector<int>lin[MAXN], lv[MAXN];

int sccno[MAXN], scc_cnt;

int n, m;

int num[MAXN], in[MAXN], val[MAXN];

void dfs_scc(int u)

{

    pre[u] = low[u] = ++dfs_clock;

    sta.push(u);

    for(int i = 0 ; i < (int)lin[u].size() ; i++){

        int v = lin[u][i];

        if(!pre[v]){

            dfs_scc(v);

            low[u] = gmin(low[u], low[v]);

        }

        else if(!sccno[v])

            low[u] = gmin(low[u], pre[v]);

    }

    if(low[u] == pre[u]){

        scc_cnt++;

        while(!sta.empty()){

            int org = sta.top();    sta.pop();

            sccno[org] = scc_cnt;

            num[scc_cnt]++;

            if(org == u)

                break;

        }

    }

}

void Tarjan()

{

    memset(pre, 0, sizeof(pre));

    memset(sccno, 0, sizeof(sccno));

    memset(num, 0, sizeof(num));

    dfs_clock = 0;

    scc_cnt = 0;

    while(!sta.empty()) sta.pop();

    for(int i = 1 ; i <= n ; i++)

        if(!pre[i]) dfs_scc(i);

}

void topo()

{

    for(int i = 1 ; i <= n ; i++)

        lv[i].clear();

    memset(in, 0, sizeof(in));

    for(int i = 1 ; i <= n ; i++){

        for(int j = 0 ; j < (int)lin[i].size() ; j++){

            int u = i, v = lin[i][j];

            if(sccno[u] != sccno[v]){

                in[sccno[v]]++;

                lv[sccno[u]].push_back(sccno[v]);

            }

        }

    }

    memset(val, 0, sizeof(val));

    while(!que.empty()) que.pop();

    for(int i = 1 ; i <= scc_cnt ; i++)

        if(in[i] == 0)  que.push(i);

    while(!que.empty()){

        int org = que.front();  que.pop();

        for(int i = 0 ; i < (int)lv[org].size() ; i++){

            int v = lv[org][i];

            in[v]--;

            val[v] = gmax(val[v], val[org] + num[org]);

            if(!in[v])  que.push(v);

        }

    }

}

int main()

{

    int t;

    scanf("%d", &t);

    while(t--){

        scanf("%d%d", &n, &m);

        for(int i = 1 ; i <= n ; i++)

            lin[i].clear();

        int u, v;

        for(int i = 1 ; i <= m ; i++){

            scanf("%d%d", &u, &v);

            lin[u].push_back(v);

        }

        Tarjan();

        topo();

        int ans = 0;

//        for(int i = 1 ; i <= scc_cnt ; i++)

//            printf("num[%d] = %d\n", i, num[i]);

//        for(int i = 1 ; i <= scc_cnt ; i++)

//            printf("val[%d] = %d\n", i, val[i]);

//        for(int i = 1 ; i <= n ; i++)

//            printf("sccno[%d] = %d\n", i, sccno[i]);

        for(int i = 1 ; i <= scc_cnt ; i++)

            ans = gmax(ans, val[i] + num[i]);

        printf("%d\n", ans);

    }

    return 0;

}