BZOJ4653: [Noi2016]区间(线段树 双指针)

题意

题目链接

Sol

按照dls的说法，一般这一类的题有两种思路，一种是枚举一个点\(M\)，然后check它能否成为答案。但是对于此题来说好像不好搞

另一种思路是枚举最小的区间长度是多少，这样我们把所有区间按长度排序后可以二分出满足条件的最短的区间长度

观察后不难发现，较长区间的长度一定是随着短区间长度的增加而单调递增的。

直接用双指针维护即可。

判断是否可行也就是是否有一个点被覆盖了\(m\)次，离散化后线段树维护。。

经验：

1. 该类问题的两种思路
2. 最大值的单调性

#include<bits/stdc++.h>
#define ls k << 1
#define rs k << 1 | 1
using namespace std;
const int MAXN = 4e6 + 10, INF = 1e9 + 10;
inline int read() {
	char c = getchar(); int x = 0, f = 1;
	while(c < '0' || c > '9') {if(c == '-')f =- 1; c = getchar();}
	while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
	return x * f;
}
int N, M;
struct Node {
	int l, r, mx, f, siz; 
}T[MAXN];
struct Line {
	int l, r;
	bool operator < (const Line &rhs) const {
		return abs(r - l) <= abs(rhs.r - rhs.l);
	}
}a[MAXN];
int date[MAXN], num;
void add(int k, int val) {
	T[k].mx += val;
	T[k].f += val;
} 
void pushdown(int k) {
	if(!T[k].f) return ;
	add(ls, T[k].f); add(rs, T[k].f);
	T[k].f = 0;
}
void update(int k) {
	T[k].mx = max(T[ls].mx, T[rs].mx);
}
void Build(int k, int ll, int rr) {
	T[k].l = ll; T[k].r = rr;
	if(ll == rr) return ;
	int mid = ll + rr >> 1;
	Build(ls, ll ,mid); Build(rs, mid + 1, rr);
}
void Add(int k, int ll, int rr, int val) {
	if(ll <= T[k].l && T[k].r <= rr) {add(k, val); return ;}
	pushdown(k);
	int mid = T[k].l + T[k].r >> 1;
	if(ll <= mid) Add(ls, ll, rr, val);
	if(rr >  mid) Add(rs, ll, rr, val);
	update(k);
}
main() {
//	freopen("testdata.in", "r", stdin);
	N = read(); M = read();
	for(int i = 1; i <= N; i++) 
		a[i].l = read(), a[i].r = read(), date[++num] = a[i].l, date[++num] = a[i].r;
//	puts("-1");
	//sort(a + 1, a + N + 1);
	stable_sort(a + 1, a + N + 1);
//	puts("-1");
	sort(date + 1, date + num + 1);
	num = unique(date + 1, date + num + 1) - date - 1;
	for(int i = 1; i <= N; i++) a[i].l = lower_bound(date + 1, date + num + 1, a[i].l) - date,
								a[i].r = lower_bound(date + 1, date + num + 1, a[i].r) - date;
	Build(1, 1, num);
	
	int ans = INF, now = 0;
	for(int i = 1; i <= N; i++) {
		while(T[1].mx < M && (now < N)) 
			now++, Add(1, a[now].l, a[now].r, 1); 
		if(T[1].mx >= M) ans = min(ans, date[a[now].r] - date[a[now].l] - (	date[a[i].r] - date[a[i].l]));
		Add(1, a[i].l, a[i].r, -1);
	}
	printf("%d\n", ans == INF ? -1 : ans);
}
/*
*/