Uva11324——最大团

最新推荐文章于 2018-11-15 16:34:15 发布

cool_Fires

最新推荐文章于 2018-11-15 16:34:15 发布

阅读量686

点赞数

CC 4.0 BY-SA版权

文章标签： Graph theory Dynamic programming

本文链接：https://blog.youkuaiyun.com/cool_Fires/article/details/9283175

本文详细介绍了一种基于强连通分量(SCC)的算法实现，包括预处理阶段、深度优先搜索进行SCC划分及逆向图求解最长路径等步骤。通过具体代码展示了如何解决涉及图论的经典问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

经典问题，没什么好说的。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <stack>
#include <vector>
using namespace std;

const int maxn = 1000 + 10;
vector<int> g[maxn], g1[maxn];
int n, m;
int pre[maxn], lowlink[maxn], sccno[maxn], dfs_clock, scc_cnt;
int scc_num[maxn], top[maxn], vis[maxn], dp[maxn], cnt;
stack<int> s;

void prework()
{
    memset(pre, 0, sizeof(pre));
    memset(lowlink, 0, sizeof(lowlink));
    memset(sccno, 0, sizeof(sccno));
    memset(vis, 0, sizeof(vis));
    dfs_clock = scc_cnt = 0;
    cnt = 1;
    scanf("%d %d", &n, &m);
    for(int i = 1; i <= n; ++i) { g[i].clear(); g1[i].clear(); }
    while(!s.empty()) s.pop();
    for(int i = 0; i < m; ++i)
    {
        int a, b;
        scanf("%d %d", &a, &b);
        g[a].push_back(b);
    }
}

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++;
        scc_num[scc_cnt] = 0;
        while(1)
        {
            int x = s.top(); s.pop();
            sccno[x] = scc_cnt;
            scc_num[scc_cnt]++;
            if(x == u) break;
        }
    }
}

void dfs1(int u)
{
    vis[u] = 1;
    for(int i = 0; i < g1[u].size(); ++i)
    {
        int v = g1[u][i];
        if(!vis[v]) dfs1(v);
    }
    top[cnt++] = u;
}

void slove()
{
    for(int i = 1; i <= n; ++i)
        if(!pre[i]) dfs(i);
    for(int i = 1; i <= n; ++i)
    {
        for(int j = 0; j < g[i].size(); ++j)
        {
            int v = g[i][j];
            if(sccno[i] != sccno[v])
                g1[sccno[i]].push_back(sccno[v]);
        }
    }
    for(int i = 1; i <= scc_cnt; ++i)
        if(!vis[i]) dfs1(i);
    for(int i = 1; i <= scc_cnt; ++i)
    {
        int u = top[i];
        dp[u] = scc_num[u];
        for(int j = 0; j < g1[u].size(); ++j)
        {
            int v = g1[u][j];
            dp[u] = max(dp[u], scc_num[u] + dp[v]);
        }
    }
    int ans = 0;
    for(int i = 1; i <= scc_cnt; ++i)
        ans = max(ans, dp[i]);
    printf("%d\n", ans);
}

int main()
{
    freopen("in.txt", "r", stdin);
    int t; cin >> t;
    while(t--)
    {
        prework();
        slove();
    }
    return 0;
}