luogu P1712 [NOI2016]区间 贪心 尺取法 线段树 二分

LINK：区间

没想到尺取法.

先说暴力 可以发现答案一定可以转换到端点处 所以在每个端点从小到大扫描线段就能得到答案 复杂度(ncdot m)

再说我的做法 想到了二分 可以进行二分答案 从左到右加入线段 加到线段的每个端点的时候 将所有加入的线段 插到主席树里面 考虑判定当前ans合法.

每加入一条线段权值为v 那么意味着 权值在 v-ans这个权值区间中加1 且看一下之前是否存在相同权值的线段加入 更新其本身.

这样查一下全局最大值即可完成判定 实现起来非常繁琐 复杂度(nlog^2n) 真的是人傻 什么垃圾算法都能想出来就是想不出来正解.

考虑正解 考虑答案ans 其权值线段的最大最小值为l,r 那么意味着 l-r这个权值区间中的端点被覆盖次数>=m.

将所以线段由小到大排序 此时找到这个l,r 可以利用尺取法轻松求出 利用线段树判定即可.

复杂度(nlogn)

code

//#include<bits/stdc++.h>
#include<iostream>
#include<cstdio>
#include<ctime>
#include<cctype>
#include<queue>
#include<deque>
#include<stack>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<cctype>
#include<cstdlib>
#include<queue>
#include<deque>
#include<stack>
#include<vector>
#include<algorithm>
#include<utility>
#include<bitset>
#include<set>
#include<map>
#define ll long long
#define db double
#define INF 1000000001
#define ldb long double
#define pb push_back
#define put_(x) printf("%d ",x);
#define get(x) x=read()
#define gt(x) scanf("%d",&x)
#define gi(x) scanf("%lf",&x)
#define put(x) printf("%d
",x)
#define putl(x) printf("%lld
",x)
#define rep(p,n,i) for(RE int i=p;i<=n;++i)
#define go(x) for(int i=lin[x],tn=ver[i];i;tn=ver[i=nex[i]])
#define fep(n,p,i) for(RE int i=n;i>=p;--i)
#define vep(p,n,i) for(RE int i=p;i<n;++i)
#define pii pair<int,int>
#define mk make_pair
#define RE register
#define P 1000000007ll
#define gf(x) scanf("%lf",&x)
#define pf(x) ((x)*(x))
#define uint unsigned long long
#define ui unsigned
#define EPS 1e-4
#define sq sqrt
#define S second
#define F first
#define mod 998244353
#define l(p) t[p].l
#define r(p) t[p].r
#define sum(p) t[p].sum
#define zz p<<1
#define yy p<<1|1
using namespace std;
char *fs,*ft,buf[1<<15];
inline char gc()
{
	return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++;
}
inline int read()
{
	RE int x=0,f=1;RE char ch=gc();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=gc();}
	return x*f;

}
const int MAXN=500010;
int n,m,cnt,num;
struct wy
{
	int l,r;
	int sum;
}t[MAXN<<3],s[MAXN];
inline int cmp(wy a,wy b){return a.sum<b.sum;}
int a[MAXN],b[MAXN],tag[MAXN<<3],c[MAXN<<1];
inline void build(int p,int l,int r)
{
	l(p)=l;r(p)=r;
	if(l==r)return;
	int mid=(l+r)>>1;
	build(zz,l,mid);
	build(yy,mid+1,r);
}
inline void pushdown(int p)
{
	tag[zz]+=tag[p];
	tag[yy]+=tag[p];
	sum(zz)+=tag[p];
	sum(yy)+=tag[p];
	tag[p]=0;
}
inline void change(int p,int l,int r,int x)
{
	if(l<=l(p)&&r>=r(p))return sum(p)+=x,tag[p]+=x,void();
	int mid=(l(p)+r(p))>>1;
	if(tag[p])pushdown(p);
	if(l<=mid)change(zz,l,r,x);
	if(r>mid)change(yy,l,r,x);
	sum(p)=max(sum(zz),sum(yy));
}
int main()
{
	//freopen("1.in","r",stdin);
	get(n);get(m);
	rep(1,n,i)
	{
		get(s[i].l);get(s[i].r);
		c[++num]=s[i].l;
		c[++num]=s[i].r;
	}
	sort(c+1,c+1+num);
	rep(1,num,i)if(i==1||c[i]!=c[i-1])c[++cnt]=c[i];
	rep(1,n,i)
	{
		s[i].sum=s[i].r-s[i].l;
		s[i].r=lower_bound(c+1,c+1+cnt,s[i].r)-c;
		s[i].l=lower_bound(c+1,c+1+cnt,s[i].l)-c;
	}
	sort(s+1,s+1+n,cmp);
	build(1,1,cnt);
	int L=1,R=1,ans=INF;
	while(R<=n)
	{
		while(sum(1)<m&&R<=n)
		{
			change(1,s[R].l,s[R].r,1);
			//put(sum(1));
			++R;
		}
		while(sum(1)>=m)
		{
			ans=min(ans,s[R-1].sum-s[L].sum);
			change(1,s[L].l,s[L].r,-1);
			++L;
		}
	}
	put(ans==INF?-1:ans);
	return 0;
}