差分约束第一题,赶快学了下,这题算是入门题吧,对差分约束有了一点主观印象。
对于这题,没看差分约束的知识还是很难想到的,居然能和最长路联系起来,太玄妙了。
题意:先输入一个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) );
}
}