poj3169 差分约束系统

题意：

• 从1到n，n个数，从左向右一次排列。给定两种形式的约束条件：
  1.$x_i$与$y_i$的最大距离为$d_k$
  2.$x_i$与$y_i$的最小距离为$d_k$
  问满足这些限定条件的情况下，数1和n的最大距离是多少？（若约束条件相互矛盾则输出-1，若最大距离可以为无穷大则输出-2）

  知识补充：

  • 差分约束系统的概念：由n个变量和m个约束条件（实数）组成，且都是形如： $$x_i-y_j \le b_k（x,y为变量，b为实数)$$ 的形式。
  • 用Bellman-Ford算法求解差分约束系统：由于最短路三角不等式:$d[v] - d[u] \le e[u, v]$与差分约束的不等式形式一样，故构建j到i长度为bk的边来建成一个图，再加一个起点s，到每个顶点的距离都为0，由于可能存在负边所以用Bellman-Ford算法来求解最短路，最终得到的d[i]数组是满足该差分约束系统的一个可行解。
  • 注：若{d[1], d[2], ….. ,d[n]}是差分约束系统的一个可行解，那么{d[1] + x, d[2] + x, ….., d[n] + x}也是可行解。

  思路：

• 题目为有三个限制条件的差分约束系统，即：

  $$d[i] - d[i + 1] \le 0$$ $$d[BL] - d[AL] \le d[AL]$$ $$d[AD] - DD \le d[BD]$$ 根据这三个不等式建立图，这里核心难点是：为什么d[n] - d[1]（即起点1到终点n的最短距离）就是对于约束条件下，1到n的最长距离？（这里尚未思考明白…..)

代码：

#include <cstdio>
#include <iostream>
#include <vector>
using namespace std;
const int INF = 0x3fffffff;
struct edge{int from, to, cost;}E[100009];
int v, n, m, size;


void input(void) {
    int x, y, z;
    for (int i = 0; i < n; i++) {
        scanf("%d%d%d", &x, &y, &z);
        E[size++] = (edge){x - 1, y - 1, z};
    }
    for (int i = 0; i < m; i++) {
        scanf("%d%d%d", &x, &y, &z);
        E[size++] = (edge){y - 1, x - 1, -z};
    }
}

int bellman_ford(void) {
    int d[v];
    fill(d, d + v, INF);
    d[0] = 0;
    for (int k = 0; k < v; k++) {    //由于可能存在负圈会无限更新的情况要注意设置更新次数上限为顶点个数。
        bool update = false;
        for (int i = 0; i < size; i++) {
            if(d[E[i].to] > d[E[i].from] + E[i].cost && d[E[i].from] < INF) {
                update = true;
                d[E[i].to] = d[E[i].from] + E[i].cost;
            }
        }
        if(!update) break;
    }
    if(d[0] < 0) return -1;
    if(d[v - 1] == INF) return -2;
    return d[v - 1];
}

int main(void)
{
    while (~scanf("%d%d%d", &v, &n, &m)) {
        size = 0;
        for (int i = 0; i < v - 1; i++) {
            E[size++] = (edge){i + 1, i, 0};    //这个形式很有趣
        }
        input();
        printf("%d\n", bellman_ford());
    }
    return 0;
}