题意:给出n个区间,让你选出m个区间,使这些区间共同包含至少一个位置,一个区间的花费就是右端点减去左端点,,让你求复合题意的最小花费是多少。
思路:花费就是区间的长度,求总长度最小,先按长度排序,然后用一个线段树维护区间长度,使线段树里的区间数为m,若有线段覆盖,则加一,若该线段被抛弃,则减一,用尺取法选一遍就好了。统计最大的值若为m的话则这m个区间共同包含一个位置,最大值就是直接看tree[1]就好。然后看数据范围是1e9,所以要先对左右端点进行离散化。
ps:我觉得我的思路和题解一样,不过不知道为啥,一直有三个点过不去,可能是我离散化比较慢。
#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
#define ll long long
#define mem(ar,num) memset(ar,num,sizeof(ar))
#define me(ar) memset(ar,0,sizeof(ar))
using namespace std;
const int N = 5e5 + 10;
int n, m, cnt = 1, L[N << 2], R[N << 2], zuo[N], you[N], he[N << 1];
int tree[N << 2], tag[N << 2];
struct no {
int len, id;
} a[N];
int cmp(no a, no b) {
return a.len < b.len;
}
void pushdown(int rt) {
if(tag[rt]) {
tree[rt << 1] += tag[rt];
tree[rt << 1 | 1] += tag[rt];
tag[rt << 1] += tag[rt];
tag[rt << 1 | 1] += tag[rt];
tag[rt] = 0;
}
}
void update(int x, int y, int rt, int l, int r, int v) {
if(x > r || y < l)
return;
if(x <= l && y >= r) {
tree[rt] += v;
tag[rt] += v;
return;
}
int m = (l + r) >> 1;
pushdown(rt);
update(x, y, rt << 1, l, m, v);
update(x, y, rt << 1 | 1, m + 1, r, v);
tree[rt] = max(tree[rt << 1], tree[rt << 1 | 1]);
}
int main() {
ios::sync_with_stdio(0);
cin.tie(0);
cout.tie(0);
cin >> n >> m;
for(int l, r, i = 1; i <= n; i++) {
cin >> l >> r;
a[i] = no({
r - l, i
});
he[cnt] = l, zuo[cnt++] = l;
he[cnt] = r, you[cnt++] = r;
}
sort(he + 1, he + cnt + 1);
int ty = unique(he + 1, he + cnt + 1) - he - 1;
int li = 1, ri = 1, rig = 0;
for(int i = 1; i < cnt; i++) {
if(i & 1)
L[li++] = lower_bound(he + 1, he + 1 + ty, zuo[i]) - he - 1;
else
R[ri++] = lower_bound(he + 1, he + 1 + ty, you[i]) - he - 1;
rig = max(rig, L[li - 1]);
rig = max(rig, R[ri - 1]);
}
sort(a + 1, a + n + 1, cmp);
int ans = INF, le = 0, re = 0;
while(1) {
while(tree[1] < m && re <= n) {
re++;
int s = a[re].id, u = L[s], v = R[s];
update(u, v, 1, 1, rig, 1);
}
if(tree[1] < m)
break;
while(tree[1] >= m && le <= n) {
le++;
int s = a[le].id, u = L[s], v = R[s];
update(u, v, 1, 1, rig, -1);
}
ans = min(ans, a[re].len - a[le].len);
}
cout << ((ans == INF) ? -1 : ans) << endl;
return 0;
}