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

题解

第一次写离散化，挂了好几次
一开始看到这道题，想法是按照右端点排序，从前往后依次加入每个区间，如果出现最后加入的区间的左端点大于最先加入的区间的右端点就将最先加入的区间删除，用线段树维护。结果写完后连样例都过不了。。究其原因：有时候虽然一个小区间中的点虽然已经被超过m个区间覆盖，但是无法从这些区间中选出恰好m个区间使得答案最小。
正确做法：
将所有端点离散化，按照每个区间的长度进行排序，每次还有未加入的区间并且当前没有点被覆盖过m次，就加入后面的区间。更新答案，将最开始的区间删除。更新答案只需要用最右边的区间长度减去最左边的区间长度就可以了。这些都用线段树维护。
在这些区间中，也许有一些暂时没有贡献的区间，但是这并不影响答案，因为加入的区间长度是单调递增的。

CODE：

#include<bits/stdc++.h>
using namespace std;
struct node
{
    int l,r;
}a[500001];
struct tree
{
    int plus,sum;
}t[4000001];
int s[1000001];
int n,m,tot,L=1,R,ans=0x7fffffff;
inline int max(int a,int b){return a>b?a:b;}
inline int min(int a,int b){return a<b?a:b;}
inline bool cmp(node a,node b)
{
    return (a.r-a.l)<(b.r-b.l);
}
void add(int L,int R,int l,int r,int now,int num)
{
    if(L<=l&&r<=R)
    {
        t[now].plus+=num,t[now].sum+=num;
        return;
    }
    int mid=(l+r)>>1;
    if(L<=mid) add(L,R,l,mid,now<<1,num);
    if(R>mid) add(L,R,mid+1,r,now<<1|1,num);
    t[now].sum=max(t[now<<1].sum,t[now<<1|1].sum)+t[now].plus;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%d%d",&a[i].l,&a[i].r);
        s[++tot]=a[i].l;
        s[++tot]=a[i].r;
    }
    sort(s+1,s+tot+1);
    tot=unique(s+1,s+tot+1)-s-1;
    sort(a+1,a+n+1,cmp);
    for(int i=1;i<=n;i++)
    {
        a[i].l=lower_bound(s+1,s+tot+1,a[i].l)-s;
        a[i].r=lower_bound(s+1,s+tot+1,a[i].r)-s;
    }
    while(L<=n)
    {
        while(R<n&&t[1].sum<m) R++,add(a[R].l,a[R].r,1,tot,1,1);
        if(t[1].sum>=m) ans=min(ans,s[a[R].r]-s[a[R].l]-s[a[L].r]+s[a[L].l]);
        add(a[L].l,a[L].r,1,tot,1,-1),L++;
    }
    if(ans==0x7fffffff) printf("-1");
    else printf("%d",ans);
    return 0;
}