题解【bzoj4653 [NOI2016] 区间】
先按照长度排个序,然后依次添加区间。什么是添加?设这个区间是\([l,r]\),添加就是把\(a_l,a_{l+1},a_{l+2},{...},a_{r}\)都加上\(1\),其中\(a_i\)表示第\(i\)个位置被几个区间覆盖。拿走一个区间的含义就是把它们都减\(1\)。这个过程很显然可以用线段树维护。
如果在添加到一个区间 \(i\) 时,有一个点被区间覆盖了\(M\)次,那么先更新答案,再把前面的加入过的区间一直拿直到没有一个点被覆盖\(M\)次。如何判断有没有点被覆盖\(M\)次?因为是一个一个区间加的,所以只用维护一个\(a_i\)的最大值,看他是否\(=M\)就行了。
什么叫再把前面的加入过的区间一直拿直到没有一个点被覆盖\(M\)次?
比如你一直添加区间到第\(5\)个,此时有一个点被覆盖了\(M\)次。这时你就将第一个区间拿出,如果此时依然有有一个点被覆盖了\(M\)次,那么你就拿走第二个...
这个过程就好比一个队列,可以从后面添加区间达到一个点被覆盖了\(M\)次;从前面弹出区间直到没有一个点被覆盖了\(M\)次。
差不多就是这样,还有注意一下\(l_i,r_i \leq 10^9\),开线段树是要离散化的。上代码:
#include <bits/stdc++.h>
#define INF 1000000001
using namespace std;
const int N = 500500;
int n, m, cnt, tot, ans = INF;
struct Seg {
int l, r, len;
bool operator < (const Seg &x) const {
return len < x.len;
}
}a[N];
struct KEY {
int d, id, se;
}key[N * 2];
inline bool cmp1(KEY x, KEY y) {
return x.d < y.d;
}
inline bool cmp2(KEY x, KEY y) {
return x.id < y.id;
}
struct node {
int left, right, Max, lazy;
node *ch[2];
}pool[N * 4], *root;
inline void pushup(node *r) {
r->Max= max(r->ch[0]->Max, r->ch[1]->Max);
}
inline void pushdown(node *r) {
if(!r->lazy) return ;
r->Max += r->lazy;
if(r->ch[0]) r->ch[0]->lazy += r->lazy;
if(r->ch[1]) r->ch[1]->lazy += r->lazy;
r->lazy = 0; return ;
}
inline void build(node *r, int left, int right) {
r->left = left, r->right = right;
if(left == right) return ;
int mid = (left + right) >> 1;
node *lson = &pool[++cnt], *rson = &pool[++cnt];
r->ch[0] = lson, r->ch[1] = rson;
build(lson, left, mid), build(rson, mid + 1, right);
}
inline void change(node *r, int left, int right, int d) {
if(r->left == left && r->right == right) {
r->lazy += d; return ;
}
pushdown(r);
if(r->ch[0]->right >= right) change(r->ch[0], left, right, d);
else if(r->ch[1]->left <= left) change(r->ch[1], left, right, d);
else change(r->ch[0], left, r->ch[0]->right, d), change(r->ch[1], r->ch[1]->left, right, d);
pushdown(r->ch[0]), pushdown(r->ch[1]), pushup(r);
}
int main() {
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; i++) {
scanf("%d%d", &a[i].l, &a[i].r);
a[i].len = a[i].r - a[i].l;
key[++tot].d = a[i].l, key[tot].id = tot;
key[++tot].d = a[i].r, key[tot].id = tot;
}
sort(key + 1, key + tot + 1, cmp1);
key[0].d = -1; key[0].se = 0;
for(int i = 1; i <= tot; i++)
if(key[i].d == key[i - 1].d)
key[i].se = key[i - 1].se;
else key[i].se = key[i - 1].se + 1;
sort(key + 1, key + tot + 1, cmp2);
for(int i = 1; i <= n; i++)
a[i].l = key[i * 2 - 1].se, a[i].r = key[i * 2].se;
sort(a + 1, a + n + 1);
build(root = &pool[0], 1, 2 * n + 1);
int pos = 1;
change(root, a[1].l, a[1].r, 1);
if(m == 1) ans = 0;
for(int i = 2; i <= n; i++) {
change(root, a[i].l, a[i].r, 1);
while(root->Max >= m) {
change(root, a[pos].l, a[pos].r, -1);
ans = min(ans, a[i].len - a[pos].len);
pos++;
}
}
if(ans == INF) ans = -1;
printf("%d\n", ans);
return 0;
}