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;
}