[BZOJ4653 区间]

Description

在数轴上有 n个闭区间 [l1,r1],[l2,r2],...,[ln,rn]。现在要从中选出 m 个区间，使得这 m个区间共同包含至少一个位置。换句话说，就是使得存在一个 x，使得对于每一个被选中的区间 [li,ri]，都有 li≤x≤ri。

对于一个合法的选取方案，它的花费为被选中的最长区间长度减去被选中的最短区间长度。区间 [li,ri] 的长度定义为 ri−li，即等于它的右端点的值减去左端点的值。

求所有合法方案中最小的花费。如果不存在合法的方案，输出 −1。

Input

第一行包含两个正整数 n,m用空格隔开，意义如上文所述。保证 1≤m≤n

接下来 n行，每行表示一个区间，包含用空格隔开的两个整数 li 和 ri 为该区间的左右端点。

N<=500000,M<=200000,0≤li≤ri≤10^9

Output

只有一行，包含一个正整数，即最小花费。

Sample Input

6 3
3 5
1 2
3 4
2 2
1 5
1 4

Sample Output

2

 

双指针扫一下即可，难得的noi水题

代码：

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #define M 1000010
 #define int long long
 #define ls node<<1
 #define rs node<<1|1
 using namespace std;
 int read()
 {
     char ch=getchar();int x=;
     while(ch>''||ch<'') ch=getchar();
     while(ch>=''&&ch<='') x=x*+ch-'',ch=getchar();
     return x;
 }

 int l,r,n,m,tmp,ans=1e18;
 int b[M<<],cnt[M<<],tag[M<<];
 struct point{int l,r;}a[M];
 bool cmp(point a1,point a2) {return a1.r-a1.l<a2.r-a2.l;}

 void getans(int l,int r)
 {
     ans=min(ans,b[a[r].r]-b[a[r].l]+b[a[l].l]-b[a[l].r]);
 }

 int get(int x)
 {
     int l=,r=tmp;
     while(l<=r)
     {
         int mid=(l+r)/;
         if(b[mid]==x) return mid;
         if(b[mid]<x) l=mid+;
         else r=mid-;
     }
 }

 void push(int node)
 {
     if(tag[node]!=)
     {
         tag[ls]+=tag[node];
         tag[rs]+=tag[node];
         cnt[ls]+=tag[node];
         cnt[rs]+=tag[node];
         tag[node]=;
     }
 }

 void change(int node,int l,int r,int l1,int r1,int v)
 {
     if(l1<=l&&r1>=r)
     {
         cnt[node]+=v;
         tag[node]+=v;
         return;
     }
     if(l1>r||r1<l) return;
     int mid=(l+r)/; push(node);
     change(ls,l,mid,l1,r1,v);
     change(rs,mid+,r,l1,r1,v);
     cnt[node]=max(cnt[ls],cnt[rs]);
 }

 main()
 {
     n=read(); m=read();
     for(int i=;i<=n;i++)
     {
         a[i].l=read(),a[i].r=read();
         b[++tmp]=a[i].l,b[++tmp]=a[i].r;
     }
     sort(a+,a++n,cmp);
     sort(b+,b++tmp);
     tmp=unique(b+,b++tmp)-b-;
     for(int i=;i<=n;i++)
     {
         a[i].l=get(a[i].l);
         a[i].r=get(a[i].r);
     }
     while()
     {
         while(cnt[]>=m)
         {
             getans(l,r);
             change(,,tmp,a[l].l,a[l].r,-);
             l++;
         }
         if(r==n) break;
         r++;
         change(,,tmp,a[r].l,a[r].r,);
     }
     if(ans==1e18) puts("-1");
     else printf("%lld",ans);
     return ;
 }