bzoj4653 [Noi2016]区间

http://www.elijahqi.win/archives/2873
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
考虑将所有的线段按照长度排序 然后依次将我的线段覆盖进去 查找我所有的区间是否有点被覆盖了m次 那么这个时候用最大值减去最小值 至少是一个答案但是是不是最小的就不好说了但是如果满足条件的话 在我讲左端点的指针挪动的时候我整个区间仍然会满足条件但是答案会减小 所以直接设置两个指针扫描一遍即可 答案一定不重不漏

#include<cstdio>
#include<cctype>
#include<algorithm>
#define N 1000010
#define lc x<<1
#define rc x<<1|1
#define inf 0x7f7f7f7f
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
struct node1{
    int l,r,len;    
}line[N];
struct node{
    int max,add;
}tree[N<<2];
int a[N],nn,ans=inf,n,m;
inline bool cmp(const node1 &a,const node1 &b){return a.len<b.len;}
inline void update(int x){tree[x].max=max(tree[lc].max,tree[rc].max);}
inline void pushdown(int x){
    if(!tree[x].add) return;int ad=tree[x].add;tree[lc].add+=ad;
    tree[rc].add+=ad;tree[lc].max+=ad;tree[rc].max+=ad;tree[x].add=0;
}
inline void modify(int x,int l,int r,int l1,int r1,int v){
    if (l1<=l&&r1>=r){tree[x].max+=v;tree[x].add+=v;return;}
    int mid=l+r>>1;pushdown(x);
    if (l1<=mid) modify(lc,l,mid,l1,r1,v);
    if (r1>mid) modify(rc,mid+1,r,l1,r1,v);update(x);
}
inline int query(int x,int l,int r,int l1,int r1){
    if (l1<=l&&r1>=r) return tree[x].max;int mid=l+r>>1,tmp=0;pushdown(x); 
    if (l1<=mid) tmp=max(tmp,query(lc,l,mid,l1,r1));
    if (r1>mid) tmp=max(tmp,query(rc,mid+1,r,l1,r1));return tmp;
}
inline void print(int x,int l,int r){
    if (x>10) return;int mid=l+r>>1;if (l!=r) pushdown(x);
    print(lc,l,mid);
    printf("%d %d %d %d\n",l,r,tree[x].max,tree[x].add);
    print(rc,mid+1,r);
}
int main(){
    freopen("bzoj4653.in","r",stdin);
//  freopen("bzoj4653.out","w",stdout);
    n=read();m=read();
    for (int i=1;i<=n;++i){
        int l=read(),r=read();line[i].len=r-l;
        line[i].l=l,line[i].r=r;a[++nn]=l;a[++nn]=r;
    }sort(a+1,a+nn+1);nn=unique(a+1,a+nn+1)-a-1;
    for (int i=1;i<=n;++i){
        line[i].l=lower_bound(a+1,a+nn+1,line[i].l)-a;
        line[i].r=lower_bound(a+1,a+nn+1,line[i].r)-a;
    }int now=1;sort(line+1,line+n+1,cmp);
    for (int i=1;i<=n;++i){
        while(now<=n&&query(1,1,nn,1,a[nn])<m)
            modify(1,1,nn,line[now].l,line[now].r,1),++now;
        if(ans==inf&&now>n) break;
        if (query(1,1,nn,1,a[nn])>=m) ans=min(ans,line[now-1].len-line[i].len);else break;
        modify(1,1,nn,line[i].l,line[i].r,-1);
    }printf("%d\n",ans==inf?-1:ans);
    return 0;
}