本文主要是介绍HDU 1384(差分约束系统)，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！
题目要求的是求的最短路，
则对于 不等式 f(b)-f(a)>=c,建立 一条 a 到 b 的边 权值为 c(因为当前点b由源点a与值c来判断),则求的最长路 即为 最小值(集合)
并且有隐含条件：0<=f(a)-f(a-1)<=1 则有边权关系(a,a-1,0)以及(a-1,a,-1);
将源点到各点的距离初始化为INF(无穷大)，其中之1为0，最终求出的最短路满足 它们与该点之间相互差值最大
差分约束
在实际的应用中，一般使用SPFA(Shortest Path Fast Algorithm)算法来实现。
差分约束系统中源点到每个点的距离确定
关于Dist[]的初始化化
1.如果将源点到各点的距离初始化为0，最终求出的最短路满足 它们之间相互最接近了
2.如果将源点到各点的距离初始化为INF(无穷大)，其中之1为0，最终求出的最短路满足 它们与该点之间相互差值最大。
3.差分约束系统的确立要根据自己确定的约束条件，从约束点走向被约束点
连边一般有两种方法，第一种是连边后求最长路的方法，第二种是连边后求最短路的方法。
例：d[x]-d[y]>=Z
如果想连边后求最长路 那么将不等式变形为这种形式 d[x]>=d[y]+z y---x连一条权值为z的边
求最短路则变形成d[y]<=d[x]-z x---y连一条权值为-z的边。
如果是别的不等式，也可以根据情况变形。但是要保证的是 两个变量（x,y)的系数一定要是正的。而常量则不一定。
第一：
感觉难点在于建图
第二：
①：对于差分不等式，a - b <= c ，建一条 b 到 a 的权值为 c 的边，求的是最短路，得到的是最大值
②：对于不等式 a - b >= c ，建一条 b 到 a 的权值为 c 的边，求的是最长路，得到的是最小值
③：存在负环的话是无解
④：求不出最短路（dist[ ]没有得到更新）的话是任意解
第三：
一种建图方法：
设x[i]是第i位置（或时刻）的值（跟所求值的属性一样），那么把x[i]看成数列，前n项和为s[n]，则x[i] = s[i] - s[i-1]；
那么这样就可以最起码建立起类似这样的一个关系：0 <= s[i] - s[i-1] <= 1;
其他关系就要去题目探索了
题目意思就是: 给出一系列三元组 R ( x, y, v ) , 表示在闭区间 [x,y] 中有 v 个不同的整数出现,
现在要你求一个最小的集合, 使得集合中的元素满足所有的元组, 且集合元素最少.
定义函数 S( x ) , 表示[0,x] 中出现元素的个数, 那么对于元组 R ( x, y, v ) 有 S( y ) – S( x-1 ) >= v.
因为这样的元组又很多, 所以就构成了一组约束序列, 也就是差分约束了 ( 大概是因为它们的差值S( y ) – S( x-1 ) >= v.满足一些约束 ) .
我们看看 SPFA 的一个松弛操作 : if ( R[y] > R[x] + v[y] ) R[y] = R[x] + v[y]; 可以发现, 它是在新增条件
比原定条件更优时, 对原定条件进行优化( 也就是松弛 ) . 这样进行了一系列的松弛以后, 结果就是
最优的了( 详细可以看SPFA的证明 ). 对比一个元组的关系 S( y ) – S( x-1 ) >= v, 以及最终关系 :
S( MAX ) – S ( MIN ) >= RES, (其中 MAX 为序列中的出现的最大值, MIN 为最小值). 那么必然存在一个
一个min ( RES ), 使得最终关系成立, 怎样求这个最小值呢? 看到前面的SPFA的松弛操作就应该明白
了, 这也是一步步对结果进行松弛的过程, 直到不能松弛.
另外, 仅仅使用题给的元组进行构图, 那么这个图不是强连通的, S 存在隐藏关系 :
S(i+1) – S (i) >= 0
S ( i ) – S ( i+1 ) >= -1
加上这2个关系整个图就是强连通的了.
add（a, b+1, c）由a指向b+1，权值为c。。a走到b+1要加上值c。
按最小值最长路：
#include <iostream>
#include <cstring>
#include <queue>
#include <cstdio>
#include <climits>
using namespace std;
#define maxn 50010
#define maxedges maxn*3
struct edge
{
int v, w, next;
} edges[maxedges];
int head[maxn], cnt;
void add(int u, int v, int w)
{
edges[cnt].v = v;
edges[cnt].w = w;
edges[cnt].next = head[u];
head[u] = cnt++;
}
int SPFA(int start, int target)
{
queue<int> q;
bool inque[maxn];
int d[maxn];
for (int i = 0; i < maxn; ++i) d[i] = -maxn;
memset(inque, false, sizeof inque);
d[start] = 0;
q.push(start);
while (!q.empty())
{
int u = q.front(); q.pop();
inque[u] = false;
for (int e = head[u]; e != -1; e = edges[e].next)
if (d[edges[e].v] < d[u] + edges[e].w)
{
d[edges[e].v] = d[u] + edges[e].w;
if (!inque[edges[e].v])
{
inque[edges[e].v] = true;
q.push(edges[e].v);
}
}
}
return d[target];
}
int main()
{
int n, a, b, c, i, _min, _max;
while (cin >> n)
{
_min = maxn, _max = 0;
memset(head, -1, sizeof head);
cnt = 0;
for (i = 0; i < n; ++i)
{
cin >> a >> b >> c;
add(a, b + 1, c);
_min = _min > a ? a : _min;
_max = _max < b + 1 ? b + 1 : _max;
}
for (i = _min + 1; i <= _max; ++i)
{
add(i - 1, i, 0);
add(i, i - 1, -1);
}
printf("%d\n", SPFA(_min, _max));
}
}
按最短路最大值：
#include <iostream>
#include <cstring>
#include <queue>
#include <cstdio>
#include <climits>
using namespace std;
#define maxn 50010
#define maxedges maxn*3
struct edge
{
int v, w, next;
} edges[maxedges];
int head[maxn], cnt;
void add(int u, int v, int w)
{
edges[cnt].v = v;
edges[cnt].w = w;
edges[cnt].next = head[u];
head[u] = cnt++;
}
int SPFA(int start, int target)
{
queue<int> q;
bool inque[maxn];
int d[maxn];
for (int i = 0; i < maxn; ++i) d[i] = maxn;
memset(inque, false, sizeof inque);
d[start] = 0;
q.push(start);
while (!q.empty())
{
int u = q.front(); q.pop();
inque[u] = false;
for (int e = head[u]; e != -1; e = edges[e].next)
if (d[edges[e].v] > d[u] + edges[e].w)
{
d[edges[e].v] = d[u] + edges[e].w;
if (!inque[edges[e].v])
{
inque[edges[e].v] = true;
q.push(edges[e].v);
}
}
}
return d[target];
}
int main()
{
int n, a, b, c, i, _min, _max;
while (scanf("%d", &n) != EOF)
{
_min = maxn, _max = 0;
memset(head, -1, sizeof head);
cnt = 0;
for (i = 0; i < n; ++i)
{
cin >> a >> b >> c;
add(b + 1, a, -c);
_min = _min > a ? a : _min;
_max = _max < b + 1 ? b + 1 : _max;
}
for (i = _min + 1; i <= _max; ++i)
{
add(i - 1, i, 1);
add(i, i - 1, 0);
}
printf("%d\n", -SPFA(_max, _min)); //SPFA(_min, _max)
}
}
最小值最长路，add（a，b+1，c）换成add（b+1，a，c），由大走向小。：
对比上面与最小值最长路、最大值最短路的差别：
#include <iostream>
#include <cstring>
#include <queue>
#include <cstdio>
#include <climits>
using namespace std;
#define maxn 50010
#define maxedges maxn*3
struct edge
{
int v, w, next;
} edges[maxedges];
int head[maxn], cnt;
void add(int u, int v, int w)
{
edges[cnt].v = v;
edges[cnt].w = w;
edges[cnt].next = head[u];
head[u] = cnt++;
}
int SPFA(int start, int target)
{
queue<int> q;
bool inque[maxn];
int d[maxn];
for (int i = 0; i < maxn; ++i) d[i] = -maxn;
memset(inque, false, sizeof inque);
d[start] = 0;
q.push(start);
while (!q.empty())
{
int u = q.front(); q.pop();
inque[u] = false;
for (int e = head[u]; e != -1; e = edges[e].next)
if (d[edges[e].v] < d[u] + edges[e].w)
{
d[edges[e].v] = d[u] + edges[e].w;
if (!inque[edges[e].v])
{
inque[edges[e].v] = true;
q.push(edges[e].v);
}
}
}
return d[target];
}
int main()
{
int n, a, b, c, i, _min, _max;
while (scanf("%d", &n) != EOF)
{
_min = maxn, _max = 0;
memset(head, -1, sizeof head);
cnt = 0;
for (i = 0; i < n; ++i)
{
cin >> a >> b >> c;
add(b + 1, a, c);
_min = _min > a ? a : _min;
_max = _max < b + 1 ? b + 1 : _max;
}
for (i = _min + 1; i <= _max; ++i)
{
add(i - 1, i, -1);
add(i, i - 1, 0);
}
printf("%d\n", SPFA(_max, _min)); //SPFA(_min, _max)
}
}
最各节点进行排序，去重，建边，对比与上面代码的区别：
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <iostream>
using namespace std;
const int N = 55555;
const int INF = 0x55555555;
struct Edge {
int id, nx, val;
Edge() {}
Edge(int id, int nx, int val) : id(id), nx(nx), val(val) {}
} edge[N << 2];
int eh[N], ec;
void init() {
ec = 0;
memset(eh, -1, sizeof(eh));
}
void addedge(int u, int v, int w) {
edge[ec] = Edge(v, eh[u], w);
eh[u] = ec++;
}
int rx[N << 1], dis[N];
bool vis[N];
queue<int> Q;
void spfa(int s) {
while (!Q.empty()) Q.pop();
memset(vis, 0, sizeof(vis));
for (int i = 0; i < N; i++) dis[i] = -INF;
Q.push(s);
vis[s] = true;
dis[s] = 0;
while (!Q.empty()) {
int cur = Q.front();
Q.pop();
vis[cur] = false;
for (int t = eh[cur]; ~t; t = edge[t].nx) {
if (dis[edge[t].id] < dis[cur] + edge[t].val) {
dis[edge[t].id] = dis[cur] + edge[t].val;
if (vis[edge[t].id]) continue;
Q.push(edge[t].id);
vis[edge[t].id] = true;
}
}
}
}
int main() {
// freopen("in", "r", stdin);
int n, x, y, z;
while (~scanf("%d", &n)) {
init();
int cnt = 0;
for (int i = 0; i < n; i++) {
scanf("%d%d%d", &x, &y, &z);
addedge(x, y + 1, z);
rx[cnt++] = x;
rx[cnt++] = y + 1;
}
sort(rx, rx + cnt);
cnt = unique(rx, rx + cnt) - rx;
for (int i = 1; i < cnt; i++) addedge(rx[i], rx[i - 1], rx[i - 1] - rx[i]), addedge(rx[i - 1], rx[i], 0);
spfa(rx[0]);
printf("%d\n", dis[rx[cnt - 1]]);
}
return 0;
}
这篇关于HDU 1384(差分约束系统)的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！
