区间选点
给定一个数轴上的 n 个区间,要求在数轴上选取最少的点使得第 i 个区间 [ai, bi] 里至少有 ci 个点
Input
输入第一行一个整数 n 表示区间的个数,接下来的 n 行,每一行两个用空格隔开的整数 a,b 表示区间的左右端点。1 <= n <= 50000, 0 <= ai <= bi <= 50000 并且 1 <= ci <= bi - ai+1
Output
输出一个整数表示最少选取的点的个数
Sample Input
5
3 7 3
8 10 3
6 8 1
1 3 1
10 11 1
Sample Output
6
问题分析
差分约束
差分约束系统 是一种特殊的 n 元一次不等式组,它包含 n个变量以及m个约束条件。
求解差分约束系统,都可以转化为图论中单元最短路问题 ,对于差分约束中的每一个不等式约束𝑥i−𝑥j≤𝑐k 都可以移项变形为 𝑥i≤𝑐k+𝑥j。我们可以把变量𝑥i看作图中的一个结点,对于每个不等式约束𝑥i− 𝑥j≤𝑐k,从节点𝑗到节点𝑖连一条长度为ck的有向边,于是求解差分约束问题变成了求解最短路问题。
解题思路
在本题中,要求在数轴上选取最少的点使得第 i 个区间 [ai, bi] 里至少有 ci 个点。若记𝑠𝑢𝑚[𝑖]表示数轴上[0,𝑖]之间选点的个数 ,对于第i个区间 [𝑎i,𝑏i] 需要满足𝑠𝑢𝑚[bi+1]−𝑠𝑢𝑚[𝑎i]≥𝑐i。
若sum[i] 为变量,则我们可以得到m个不等式,此题可以用差分约束解决。
将不等式变形可得 sum[bi+1] ≥ ci +sum[ai]。即从节点sum[ai] 到sum[bi+1]有一条有向边ci。
因此,题目输入的是一些特定的边,我们只需要加入这些特定的边即可。
除了这些特定的边,还需关注到,区间[0,n] 中 每一个子区间 [i,i+1] 选的点数大于等于0,而 [i+1,i] 区间不符合规定,不能选点,即为-1。向图中加入这些边。
添加完边之后,使用SPFA算法跑最长路,最终结果即为dis[maxb]。
代码实现
#include<stdio.h>
#include<queue>
#include<math.h>
#include<string.h>
#define N 50500
using namespace std;
struct edge
{
int v,net;
int w;
};
int head[N],dis[N],vis[N];
edge e[160000];
queue<int> q;
int index=0;
void add(int u, int v, int w)
{
e[index].v = v;
e[index].w = w;
e[index].net = head[u];
head[u] = index;
index++;
}
void sqfa(int s)
{
memset(vis,0,sizeof(vis));
memset(dis,-1,sizeof(dis));
dis[s] = 0;
vis[s] = 1;
q.push(s);
while(!q.empty())
{
// printf("!!!");
int a=q.front();
q.pop();
vis[a] = 0;
for(int i=head[a]; i!=-1; i=e[i].net)
{
int b = e[i].v;
int w = e[i].w;
if(dis[b] < dis[a]+w)
{
dis[b] = dis[a]+w;
if(vis[b] == 0)
{
vis[b]=1;
q.push(b);
}
}
}
}
}
int main()
{
int t=0;
scanf("%d",&t);
memset(head,-1,sizeof(head));
int maxb=0, mina=52000;
for(int k=1; k<=t; k++)
{
int a,b,c;
scanf("%d%d%d",&a,&b,&c);
if(b+1>maxb) maxb = b+1;
if(a<mina) mina = a;
add(a,b+1,c);
}
for(int k=mina; k<maxb; k++)
{
add(k, k+1, 0);
add(k+1, k, -1);
}
sqfa(mina);
printf("%d",dis[maxb]);
}