此文出处:https://www.hrwhisper.me/algorithm-graph-dijkstra-spfa-bellmanford-prim-kruskal/ 第四大点
差分约束系统就是给出一些形如x-y<=b不等式的约束,问你是否有满足问题的解,或者求最小,最大解。
这个问题的神奇之处是可以转化为图论的最短路问题。
(一)原理
对于图论的最短路径,有:对于d(v) <= d(u) + w(u, v) ,而差分约束系统的解法利用到了单源最短路径问题中的三角形不等式。
移项得:d(v) – d(u) <= w(u, v),是不是和上面的x-y<=b的一样?
是的,这就是转化为最短路径算法的原理。
(二)建图
差分约束最难的地方在于建图。
1.求最小的解
有一个序列,题目用n个整数组合 [ai,bi,ci]来描述它,[ai,bi,ci]表示在该序列中处于[ai,bi]这个区间的整数至少有ci个。如果存在这样的序列,请求出满足题目要求的最短的序列长度是多少。
我们可以设s[i]为从[1,i]闭区间内的整数个数。
那么显然对于组合[a,b,c] 我们有 s[b] – s[a-1] >= c (a-1是因为我们需要算上a)
还有:0<= s[b]-s[b-1] <=1
变形得
S[b]-S[b-1] >=0
S[b-1]-S[b] >=-1
接着,对于求最小值的来说,我们用最长路来求解。为什么是最长路?一开始把距离d设为无穷小,之后最长路的更新公式为:if(d[v] < d[u]+w(u,v)) d[v]=d[u]+w(u,v); 可以看到d[v]不断的增大,即当满足条件的时候为最小的解。
所以我们在建图的时候,需要满足最长路的三角不等式: d(v) >= d(u) + w(u, v) =>d(v) – d(u) >=w(u,v) 的形式=>从u到v连接一条权重为w(u,v)的边。
故本题中三个不等式都是如下的形式:
s[b] – s[a-1] >= c
S[b]-S[b-1] >=0
S[b-1]-S[b] >=-1
题解见笔者CSDN:http://blog.csdn.net/murmured/article/details/18793199
2.求最大的解
一些母牛按升序排成一条直线。有两种要求,A和B距离(A < B)不得超过X,还有一种是C和D(C < D)距离不得少于Y,问可能的最大距离。如果没有输出-1,如果可以随便排输出-2,否则输出最大的距离。
题目按升序排,所以有:
s[i] – s[i-1] >=0
由题目的约束条件得到:
B – A <= X
D – C >= Y
因为题目求的是最大距离,故用最短路来求解
最短路三角不等式为 d(v) <= d(u) + w(u,v) => d(v) – d(u) <= w(u,v) =>从u到v连接一条权重为w(u,v)的边
所以我们都转化为小于等于的符号,所以有:
s[i-1] – s[i] <=0
B – A <= X
C – D <= -Y
建完图后SPFA即可。(有负环说明无解输出-1 , 1与n不连通说明可以随意摆放,没有约束嘛。输出-2,否则输出dis [n])
题解见笔者CSDN:http://blog.csdn.net/murmured/article/details/18819955
(三)注意点
1.不等式取>=还是<=
当你在纠结用小于号大于号的时候,看看题目求的是最大还是最小,如果只是判断有木有解,那么大于号小于都可以,只不过要注意全部等式要统一。
当求解最大解,使用最短路径求解, 由图最短路径的三角不等式:
d(v) <= d(u) + w(u,v) => d(v) – d(u) <= w(u,v) =>从u到v连接一条权重为w(u,v)的边
求最大的解用最短路,为d(v) – d(u) <= w(u,v) 的形式
求最小的解用最长路,为d(v) – d(u) >= w(u,v) 的形式
2.图不连通
原图可能不是连通图,故需要加一个超级源点S,从S到任意的顶点边权为0,然后从该点出发。为什么?添加从虚点S到每个顶点的权为0的边.这是为了保证构造出来的图是连通的.由于虚点本身并不引入负圈,所以设置虚点以后最短路仍然存在,并且每个约束仍然满足.
或者差分约束不用什么附加顶点, 附加顶点的唯一用处就是保证图的连通性, 不让你有负环判不到的情况, 解决这种问题的最佳途径就是初始把所有顶点都加入队列, 并且将所有dis
置0, 这就相当于加了一个不存在的附加顶点, 它与所有的顶点的直连长度都是0.
当然推荐第二种,效率也高。