标签：
Bellman-Ford可以求有负权的图的最短路，也可以判断是否有负环存在。
单源最短路模板
//从s出发到所有点的最短距离 void BellmanFord(int s) { for (int i = 1; i <= n; ++i) d[i] = INF; d[s] = 0; while (true) { bool update = false; for (int i = 0; i < E; ++i) { edge e = es[i]; if (d[e.from] != INF && d[e.to] > d[e.from] + e.cost) { d[e.to] = d[e.from] + e.cost; update = true; } } if (!update) break; } }
如果没有环，最多更新n-1次，n是点的个数，如果更新了n次，证明图中有环。因为求的是最短路，所以一般判断的是负环，正环也可以判断的。
poj3259 Wormholes
求是否存在负环
#include <iostream> #include <vector> #include <cstdio> using namespace std; const int INF = 0x5f5f5f5f; const int MAX_V = 505; const int MAX_E = 2705; struct edge { int from, to, cost; }; vector<edge> G; int d[MAX_V]; int V, E, W; bool BellmanFord(int n) { for (int i = 1; i <= n; ++i) d[i] = INF; d[1] = 0; for (int i = 1; i <= n; ++i) { for (int j = 0; j < G.size(); ++j) { int u = G[j].from, v = G[j].to; if (d[u] != INF && d[v] > d[u]+G[j].cost) { d[v] = min(d[v], d[u]+G[j].cost); if (i == n) return true; } } } return false; } int main() { //freopen("in.txt", "r", stdin); int t; scanf("%d", &t); while (t--) { G.clear(); scanf("%d%d%d", &V, &E, &W); int u, v, w; for (int i = 0; i < E; ++ i) { scanf("%d%d%d", &u, &v, &w); G.push_back(edge{u, v, w}); G.push_back(edge{v, u, w}); } for (int i = 0; i < W; ++ i) { scanf("%d%d%d", &u, &v, &w); G.push_back(edge{u, v, -w}); } if (BellmanFord(V)) puts("YES"); else puts("NO"); } return 0; }
poj1860 Currency Exchange
求是否存在正环
//poj1860 #include <iostream> #include <cstdio> #include <vector> using namespace std; const int N = 105; const int INF = 0x5f5f5f5f; struct Edge { int from, to; double rate, com; }; vector<Edge> G; double d[N]; bool BellmanFord(int s, double money, int n) { for (int i = 1; i <= n; ++i) d[i] = 0; d[s] = money; for (int i = 1; i <= n; ++i) { for (int j = 0; j < G.size(); ++j) { int u = G[j].from, v = G[j].to; if (d[v] < (d[u]-G[j].com)*G[j].rate) { d[v] = (d[u]-G[j].com)*G[j].rate; if (i == n) return true; } } } return false; } int main() { //freopen("in.txt", "r", stdin); int n, m; int s; double money; while (~scanf("%d%d%d%lf", &n, &m, &s, &money)) { G.clear(); while (m--) { int u, v; double r1, c1, r2, c2; scanf("%d%d%lf%lf%lf%lf", &u, &v, &r1, &c1, &r2, &c2); G.push_back((Edge){u, v, r1, c1}); G.push_back((Edge){v, u, r2, c2}); } if (BellmanFord(s, money, n)) puts("YES"); else puts("NO"); } return 0; }
Bellman-Ford的复杂度很高，可以用队列来优化（也就变成了传说中的spfa）
模板：
const int INF = 0x5f5f5f5f; const int N = 105; struct Edge { int to, cost; }; vector<Edge> G[N]; int cnt[N], pre[N], d[N]; bool inq[N]; bool BellmanFord(int s, int n) { queue<int> q; memset(inq, 0, sizeof inq); memset(cnt, 0, sizeof cnt); for (int i = 1; i <= n; ++i) d[i] = INF; d[s] = 0; inq[s] = true; q.push(s); while (q.size()) { int u = q.front(); q.pop(); inq[u] = false; for (int i = 0; i < G[u].size(); ++i) { int v = G[u][i].to; int c = G[u][i].cost; if (d[u] < INF && d[v] > d[u] + c) { d[v] = d[u] + c; pre[v] = u; if (!inq[v]) { q.push(v); inq[v] = true; if (++cnt[v] > n) return false; // 存在负环 } } } } return true; }
poj1847 Tram
#include <iostream> #include <cstring> #include <cstdio> #include <queue> using namespace std; const int INF = 0x5f5f5f5f; const int N = 105; struct Edge { int to, cost; }; vector<Edge> G[N]; int d[N]; bool inq[N]; bool BellmanFord(int s, int n) { queue<int> q; memset(inq, 0, sizeof inq); for (int i = 1; i <= n; ++i) d[i] = INF; d[s] = 0; inq[s] = true; q.push(s); while (q.size()) { int u = q.front(); q.pop(); inq[u] = false; for (int i = 0; i < G[u].size(); ++i) { int v = G[u][i].to; int c = G[u][i].cost; if (d[u] < INF && d[v] > d[u] + c) { d[v] = d[u] + c; if (!inq[v]) { q.push(v); inq[v] = true; } } } } return true; } int main() { //freopen("in.txt", "r", stdin); int n, s, t; while (~scanf("%d%d%d", &n, &s, &t)) { for (int i = 0; i <= n; ++i) G[i].clear(); int cnt, v; for (int i = 1; i <= n; ++i) { scanf("%d", &cnt); for (int j = 0; j < cnt; ++j) { scanf("%d", &v); G[i].push_back(Edge{v, j!=0}); } } BellmanFord(s, n); printf("%d\n", d[t]==INF ? -1 : d[t]); } return 0; }
标签：
原文地址：http://www.cnblogs.com/wenruo/p/5831134.html
