当前位置:首页 > 随笔 > 正文内容

Tarjan算法及其应用

访客 随笔 2026年6月3日 1
为了更好地理解Tarjan算法,避免常见错误,以下是对一个典型问题的分析与解决方案。

初始代码(部分正确):

//缩点
#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int MAXN = 1e6 + 5, MAXM = 1e4 + 5;
int n, m, a[MAXN], result, distance[MAXN], componentID[MAXN], predecessor[MAXN];
int edge[MAXN], nextEdge[MAXN], head[MAXN], weight[MAXN], indexCounter;
int discoveryTime[MAXN], lowLink[MAXN], sizeComponent[MAXN], stackNode[MAXN];
int timeCounter, topStack, sccCount;
bool visited[MAXN];

void addEdge(int from, int to) {
    edge[indexCounter] = to, nextEdge[indexCounter] = head[from], head[from] = indexCounter++;
}

void tarjanAlgorithm(int node) {
    lowLink[node] = discoveryTime[node] = ++timeCounter;
    stackNode[++topStack] = node, visited[node] = true;
    for (int i = head[node]; ~i; i = nextEdge[i]) {
        int adjNode = edge[i];
        if (!discoveryTime[adjNode]) tarjanAlgorithm(adjNode), lowLink[node] = min(lowLink[node], lowLink[adjNode]);
        else if (visited[adjNode]) lowLink[node] = min(lowLink[node], discoveryTime[adjNode]);
    }
    if (lowLink[node] == discoveryTime[node]) {
        int y;
        ++sccCount;
        do {
            y = stackNode[topStack--];
            visited[y] = false, componentID[y] = sccCount;
            if (node == y) break;
            a[node] += a[y];
        } while (y != node);
    }
}

void topologicalSort() {
    queue<int> q;
    for (int i = 1; i <= n; i++) if (componentID[i] == i && !predecessor[i]) q.push(i), distance[i] = a[i];
    while (!q.empty()) {
        int current = q.front(); q.pop();
        for (int i = head[current]; ~i; i = nextEdge[i]) {
            int nextNode = edge[i];
            predecessor[nextNode]--;
            distance[nextNode] = max(distance[nextNode], distance[current] + a[nextNode]);
            if (!predecessor[nextNode]) q.push(nextNode);
        }
    }
}

修正后的代码:

//使用Tarjan算法将图转换为DAG并计算最长路径
#include<bits/stdc++.h>
using namespace std;
const int MAXN = 1e5 + 5;
int n, m, result, predecessor[MAXN], distance[MAXN], value[MAXN];
int edge[MAXN], nextEdge[MAXN], head[MAXN], index;
vector<int> graph[MAXN];
bool visited[MAXN];
int discovery[MAXN], lowLink[MAXN], stackNode[MAXN], componentID[MAXN];
int sccCount, timeCounter, topStack;

void addEdge(int u, int v) {
    edge[index] = v, nextEdge[index] = head[u], head[u] = index++;
}

void tarjanAlgorithm(int u) {
    lowLink[u] = discovery[u] = ++timeCounter;
    stackNode[++topStack] = u, visited[u] = true;
    for (int i = head[u]; ~i; i = nextEdge[i]) {
        int v = edge[i];
        if (!discovery[v]) tarjanAlgorithm(v), lowLink[u] = min(lowLink[u], lowLink[v]);
        else if (visited[v]) lowLink[u] = min(lowLink[u], discovery[v]);
    }
    if (lowLink[u] == discovery[u]) {
        int v;
        ++sccCount;
        do {
            v = stackNode[topStack--];
            visited[v] = false, value[sccCount] += value[v], componentID[v] = sccCount;
        } while (v != u);
    }
}

void topologicalSort() {
    queue<int> q;
    for (int i = 1; i <= sccCount; i++) if (!predecessor[i]) q.push(i), distance[i] = value[i];
    while (!q.empty()) {
        int current = q.front(); q.pop();
        for (int i = 0; i < graph[current].size(); i++) {
            int nextNode = graph[current][i];
            predecessor[nextNode]--;
            distance[nextNode] = max(distance[nextNode], distance[current] + value[nextNode]);
            if (!predecessor[nextNode]) q.push(nextNode);
        }
    }
}

int main() {
    cin >> n >> m;
    memset(head, -1, sizeof head);
    for (int i = 1; i <= n; i++) cin >> value[i];
    for (int i = 1; i <= m; i++) {
        int x, y;
        cin >> x >> y;
        addEdge(x, y);
    }
    for (int i = 1; i <= n; i++) if (!discovery[i]) tarjanAlgorithm(i);
    for (int i = 1; i <= n; i++)
        for (int j = head[i]; ~j; j = nextEdge[j]) {
            int k = edge[j];
            int x = componentID[i], y = componentID[k];
            if (x != y) graph[x].push_back(y), predecessor[y]++;
        }
    topologicalSort();
    for (int i = 1; i <= n; i++) result = max(result, distance[i]);
    cout << result << endl;
    return 0;
}

相关文章

可以按小时收费的VPS

很多 VPS 提供商都支持 按小时计费(hourly billing),想短期试用 / 临时搭建节点、测试网络、短期项目等场景非常合适。下面是当前最主流且靠谱的按小时 VPS 选项,分别按不同需求场景整理: 1. Vultr(全球节点,包括日本) 按小时计费 可选机房:东京 / 大阪 / 洛杉矶 / 法兰克福 / 伦敦 … 支持 PayPal(部分情况),但更常用信用卡/PayPal+卡价格参考$...

在 iPhone 上下载国外App

地区/国家限制App Store 会根据 Apple ID 的国家或地区限制应用下载。如果你的 Apple ID 绑定的是中国大陆,就可能无法下载 OpenAI 官方的 ChatGPT 应用,因为它在大陆 App Store 不上架。解决办法:换成美国、加拿大、香港等地区的 Apple ID。或者在现有 Apple ID 上更改地区。注册一个国外 Apple ID(推荐)比如注册 美国区 Appl...

Node.js 中的异步编程:回调与 Promise

Node.js 是一个基于 JavaScript 构建的单线程、非阻塞运行环境,它通过异步编程机制来高效处理多个操作。在执行如文件读取、API 请求或数据库查询等任务时,Node.js 不会等待这些操作完成,而是使用回调函数和 Promise 来避免阻塞主线程。 回调方式实现异步 那么当异步操作完成后,Node.js 如何知道接下来要做什么呢?这就要用到 回调函数(callback)。 回调本质上...

Selenium自动化测试入门指南

Selenium自动化测试入门指南

什么是自动化测试? 自动化测试是指利用软件工具自动执行测试用例,模拟用户操作,如打开网页、点击链接、输入文本等,并验证结果是否符合预期。 其主要优点包括: 大幅减少人工成本 测试速度快 可以在非工作时间运行 支持持续集成和交付 然而,它也存在一些局限性,例如开发成本较高、不适合快速变化的项目、依赖稳定的UI界面等。 自动化测试的应用条件 适合引入自动化测试的情况包括: 手动测试耗时且需要大量...

MariaDB Galera集群故障快速恢复指南

OpenStack控制节点采用三节点MariaDB Galera集群架构。当数据库集群因故障重启时,有时会出现Galera集群无法正常启动的问题。虽然有多种方法可以恢复数据库服务,但如何实现快速启动同时确保数据完整性呢? 通过分析日志发现,MariaDB Galera集群节点宕机时会在日志中输出以下信息: [Note] WSREP: 新集群视图:全局状态: 874d8e7e-5980-11e8-8...

Android 中 EventBus 的通信机制与实现原理深度解析

EventBus 核心设计思想 EventBus 是一个基于观察者模式的事件总线框架,广泛应用于 Android 平台以实现组件解耦。它通过中心化的消息分发机制,使不同层级、不同线程的对象能够以"发布-订阅"方式通信,避免了传统接口回调或广播带来的强依赖问题。 核心角色说明 事件(Event):任意 Java 对象,作为数据载体,如网络状态变更通知、用户登录信息等。 发布者(Publi...

发表评论

访客

◎欢迎参与讨论,请在这里发表您的看法和观点。