//博主自己的板子,有问题可以私聊
#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e3 + 5;
const int inf = 0x3f3f3f3f;
int n, m;
int mp[maxn][maxn];
int ve[maxn], vl[maxn], ee[maxn], el[maxn];
int into[maxn], start;
int Stack[maxn];
struct node {
int u, v, w;
} e[maxn];
/*
ve 是指从始点开始到顶点Vj的最大路径长度
vl 是指在不推迟整个工期的前提下,事件vj允许的最晚发生时间
ee 若活动ai由弧<vk,vj>表示,则活动ai的最早开始时间应该等于事件vk的最早发生时间
el 若活动ai由弧<vk,vj>表示,则ai的最晚开始时间要保证事件vj的最迟发生时间不拖后
*/
void topo() {
int cnt = 0;
for(int i = 1; i <= n; i++) {
int k = -1;
for(int j = 1; j <= n; j++) {
if(into[j] == 0) {
Stack[++cnt] = j;
k = j;
into[j] = -1;
break;
}
}
for(int j = 1; j <= n; j++) {
if(mp[k][j] != inf) {
ve[j] = max(ve[j], ve[k] + mp[k][j]);
into[j]--;
}
}
}
// for(int i=1;i<=n;i++) printf("ve[%d]=%d\n",i,ve[i]);
}
void solve() {
topo();
memset(vl, inf, sizeof(vl));
vl[Stack[n]] = ve[Stack[n]];
for(int i = n; i >= 1; i--) {
for(int j = 1; j <= n; j++) {
if(mp[Stack[i]][j] != inf) {
vl[Stack[i]] = min(vl[j] - mp[Stack[i]][j], vl[Stack[i]]);
}
}
}
// for(int i = 1; i <= n; i++) printf("vl[%d]=%d\n", i, vl[i]);
//计算ee
for(int i = 1; i <= m; i++) {
ee[i] = ve[e[i].u];
}
//计算el
for(int i = 1; i <= m; i++) {
el[i] = vl[e[i].v] - e[i].w;
}
printf("\n");
for(int i = 1; i <= m; i++) {
if(ee[i] == el[i]) printf("%d %d %d\n", e[i].u, e[i].v, e[i].w);
}
}
int main() {
scanf("%d%d", &n, &m);
memset(mp, inf, sizeof(mp));
for(int i = 1; i <= m; i++) {
int u, v, w;
scanf("%d%d%d", &u, &v, &w);
e[i] = {u, v, w};
// e[i].u = u;
// e[i].v = v;
// e[i].w = w;
mp[u][v] = w;
into[v]++;
}
solve();
return 0;
}
/*
7 10
1 2 3
1 4 6
1 3 2
3 4 1
2 5 4
2 4 2
4 5 1
5 7 3
6 7 4
3 6 3
*/
求解关键路径代码
猜你喜欢
转载自blog.csdn.net/yiqzq/article/details/80823014
今日推荐
周排行