BZOJ4653 [NOI2016]区间 [线段树，离散化]

　　题目传送门

区间

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

---

　　分析：

　　一开始直接贪心+离散线段树，然后WA得天花乱坠，T飞到了九霄云外。。。还是太naive了。。。

　　首先，我们能想到这样一个思路：首先对区间按照长度进行排序，这个贪心应该是显然的；然后依次将加入区间，这里加入区间是指将该区间$[l,r]$内的所有权值+1，这样就可以得到，只要有一个点的权值大于或等于$m$，那么就可以更新答案。

　　维护权值不难想到用权值线段树，但是数据范围太大需要离散化（一开始还在离散卡了好久。。。太菜了。。。）。

　　更新答案的时候依次将前面添加的区间减掉，直到所有点的权值都小于$m$，然后就可以找到该情况下更新的答案。这是用到尺取法的思想。

　　Code：

//It is made by HolseLee on 23rd July 2018
//BZOJ 4653
#include<bits/stdc++.h>
using namespace std;

const int N=5e5+;
int n,m,ans,l[N],r[N],inf=-;

struct Seg{
    int id,len;
    bool operator < (const Seg x) const {
        return len<x.len;
    }
}a[N],p[N<<];

inline int Max(int x,int y)
{
    return x>y?x:y;
}

struct segment{
    int s[N<<],sign[N<<];

    void ready()
    {
        memset(s,,sizeof(s));
        memset(sign,,sizeof(sign));
    } 

    void pushup(int rt)
    {
        s[rt]=Max(s[rt<<],s[rt<<|]);
    }

    void pushdown(int rt)
    {
        if(!sign[rt])return;
        s[rt<<]+=sign[rt];
        s[rt<<|]+=sign[rt];
        sign[rt<<]+=sign[rt];
        sign[rt<<|]+=sign[rt];
        sign[rt]=;
    }

    void update(int l,int r,int rt,int L,int R,int C)
    {
        if(l>R||r<L)return;
        if(L<=l&&r<=R){
            s[rt]+=C;sign[rt]+=C;return;}
        int mid=(l+r)>>;
        pushdown(rt);
        if(L<=mid)update(l,mid,rt<<,L,R,C);
        if(R>mid)update(mid+,r,rt<<|,L,R,C);
        pushup(rt);
    }
}T;

inline int read()
{
    char ch=getchar();int num=;bool flag=false;
    while(ch<''||ch>''){if(ch=='-')flag=true;ch=getchar();}
    while(ch>=''&&ch<=''){num=num*+ch-'';ch=getchar();}
    return flag?-num:num;
}

int main()
{
    n=read();m=read();int x,y,z,cnt=,tot=;
    for(int i=;i<=n;i++){
        x=read();y=read();
        a[i].len=y-x;a[i].id=i;
        p[++cnt].len=x;p[cnt].id=i;
        p[++cnt].len=y;p[cnt].id=i;
    }
    sort(p+,p+cnt+);
    for(int i=;i<=cnt;i++){
        x=p[i].id;tot++;
        if(!l[x]) l[x]=tot;
        else r[x]=tot;
    }
    sort(a+,a+n+);
    inf=tot;T.ready();
    int le=,ri=;
    ans=;
    while(){
        while(T.s[]<m&&ri<=n){
            z=a[++ri].id;x=l[z];y=r[z];
            T.update(,inf,,x,y,);
        }
        if(T.s[]<m)break;
        while(T.s[]>=m&&le<=n){
            z=a[++le].id;x=l[z];y=r[z];
            T.update(,inf,,x,y,-);
        }
        ans=min(ans,a[ri].len-a[le].len);
    }
    if(ans==)ans=-;
    printf("%d",ans);
    return ;
}