差分约束
每一个不等式约束xi - xj <=c,都可以变成xi <= c+xj;
和最短路径问题中的松弛操作对比:
if(dis[v]> dis[u]+w(u,v)){
dis[v] =dis[u] + w(u,v)
}
发现松弛操作之后就能满足这个不等式约束,而这也就是求最短路径的问题。
具体就是将xi, xj作为路径上的点,c作为权重。
在求解对端路径的时候,程序是按照等号来计算的,得到的是上界;如果要求下界,将<=变成>=,然后求最长路就可。
题目
给定一个数轴上的 n 个区间,要求在数轴上选取最少的点使得第 i 个区间 [ai, bi] 里至少有 ci 个点
使用差分约束系统的解法解决这道题
input
输入第一行一个整数 n 表示区间的个数,接下来的 n 行,每一行两个用空格隔开的整数 a,b 表示区间的左右端点。1 <= n <= 50000, 0 <= ai <= bi <= 50000 并且 1 <= ci <= bi - ai+1。
output
输出一个整数表示最少选取的点的个数
分析
如果记:sum[i]是0到i内选取点的个数。
区间至少有ci个点,则是:sum[bi]-sum[ai-1]>=ci
然后还需要满足的条件是:0<=sum[i]-sum[i-1]<=1 (加上一个整数单位,最多增加一个点–选的是整数点)。
不难分析需要的是下界,需要满足的不等式约束是:
1、sum[bi]-sum[ai-1]>=ci
2、sum[i]-sum[i-1]>=0
3、sum[i-1]-sum[i]>=-1
最开始我是从0开始进行最长路径算法的,但是WA了,然后我改成min(ai)-1就过了(难道sum(min(ai)-1)不是等于0吗?)
代码
#include<stdio.h>
#include<queue>
#include<math.h>
using namespace std;
int tot;
int a[100005];//点的权值
int head[100005],inq[100005],dis[100005];
queue<int> q;
//前向星
struct edge{
int to,next,w;
}e[1000005];
void add(int x,int y,int w){
e[++tot].to=y;
e[tot].next=head[x];
head[x]=tot;
e[tot].w=w;
}
void SPFA(int s){//求最长路
dis[s]=0;inq[s]=1;
q.push(s);
while(!q.empty()){
int u=q.front();q.pop();
inq[u]=0;
for(int i=head[u];i!=0;i=e[i].next){
int v=e[i].to;
if(dis[v]<dis[u]+e[i].w){
dis[v]=dis[u]+e[i].w;
if(!inq[v]){
q.push(v);
inq[v]=1;
}
}
}
}
}
int main(){
for(int i=0;i<50005;i++) dis[i]=-100000;
int T,A,B,ci;
scanf("%d",&T);//区间个数
int bmax=0,amax=50005;
for(int j=1;j<=T;j++){
scanf("%d %d %d",&A,&B,&ci);
add(A-1,B,ci);//sum【B】-sum【A-1】>=w
if(B>bmax) bmax=B;
if(A<amax) amax=A;
}
for(int i=1;i<=bmax;i++)
{
add(i-1,i,0);
add(i,i-1,-1);
}
SPFA(amax-1);
printf("%d\n",dis[bmax]);
}