[差分约束] HDU1384 Intervals

差分约束第一题，赶快学了下，这题算是入门题吧，对差分约束有了一点主观印象。

对于这题，没看差分约束的知识还是很难想到的，居然能和最长路联系起来，太玄妙了。

题意：先输入一个n，然后是n行，一行三个数表示一个范围[ai, bi]和ci，要求一个数集Z，要求Z和[ai, bi]至少有ci个公共数。

这题如何构造成差分约束系统呢？

差分约束系统有个特点，要包含减法。于是可以利用区间的减法，设d[i]表示[0,i]内包含在Z中的数。

然后就可以写出不等式了。

d[bi] - d[ai-1] >= ci

除了题目的条件外还有一些条件要加上，就是Z集合的关系。

d[i] - d[i-1] >= 0

d[i] - d[i-1] <= 1    =>   d[i-1] - d[i] >= -1

然后题目要求最小的Z，要像上面这样把不等号统一乘大于等于，跑一遍最长路。

至于为什么是大于等于，还不知道，试了一下改成小于等于，然后跑最短路，出来的结果是[min, max]内Z的最大值。

还有些细节，比如ai最小是0，如果减到-1的话数组就越界了，所以把所有条件加一，结果是不会变的。

再比如spfa+STLqueue+STLvector邻接表，我是跑超时的。

最近写网络流的教训，已经很少用vector了，有时候queue都手写，STL效率确实不够高，被卡时间卡到哭啊。

#include<cstring>
#include<cstdio>
#include<queue>
#include<iostream>
#include<algorithm>
#define pb push_back
#define pii pair<int,int>
#define LL __int64
#define INF 0x7fffffff
#define LLINF 0x7fffffffffffffff
using namespace std;
const int N =55005;
struct eg{
    int u, v, w;
    eg(){}
    eg(int a, int b, int c) { u = a, v = b, w = c; }
}edg[N<<4];
int fir[N<<4], nex[N<<4], ecnt;
void add(int a, int b, int c){
    edg[ecnt] = eg(a, b, c);
    nex[ecnt] = fir[a], fir[a] = ecnt++;
}
int dis[N];
bool vis[N];
int spfa(int s, int n){
    for(int i = 0 ; i <= n; ++i) dis[i] = -INF, vis[i] = 0;
    queue<int>q;
    dis[s] = 0; vis[s] = 1;
    q.push(s);
    while( !q.empty() ){
        int u = q.front(); q.pop(); vis[u] = 0;
        for(int k = fir[u]; k != -1; k = nex[k]){
            int v = edg[k].v, w = edg[k].w;
            if(dis[v] < dis[u] + w){
                dis[v] = dis[u] + w;
                if( !vis[v] ){
                    vis[v] = 1;
                    q.push(v);
                }
            }
        }
    }
    return dis[n];
}
int main(){
    int a, b, c, n, l, r;
    while( scanf("%d", &n) != EOF ){
        memset(fir, -1, sizeof(fir));
        l = 50005, r = 0, ecnt = 0;
        for(int i = 0; i < n; ++i){
            scanf("%d %d %d", &a, &b, &c);
            a += 1, b += 1;
            add(a-1, b, c);             // d[b] - d[a-1] >= c
            l = min (a, l);
            r = max (b, r);
        }
        for(int i = l; i <= r; ++i){
            add(i-1, i, 0);             //  d[i] - d[i-1] >= 0
            add(i, i-1, -1);            //  d[i-1] - d[i] >= -1
        }
        printf("%d\n", spfa(l-1, r) );
    }
}