【NOI2016】区间

时间：2017-04-29 23:26:59 阅读：277 评论：0 收藏：0 [点我收藏+]

标签：tree struct ext input -o center 区间 nat blog

目链接：http://uoj.ac/problem/222

在数轴上有 $n$

对于一个合法的选取方案，它的花费为被选中的最长区间长度减去被选中的最短区间长度。区间 $[l_{i}, r_{i}]$

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

输入格式

第一行包含两个正整数 $n, m$

接下来 $n$

输出格式

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

样例一

input

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

output

explanation

技术分享

如图，当 $n = 6, m = 3$

限制与约定

所有测试数据的范围和特点如下表所示：

测试点编号	$n$	$m$	$l_{i}, r_{i}$
1	$20$	$9$	$0 \leq l_{i} \leq r_{i} \leq 100$
2	$20$	$10$	$0 \leq l_{i} \leq r_{i} \leq 100$
3	$199$	$3$	$0 \leq l_{i} \leq r_{i} \leq 100000$
4	$200$	$3$
5	$1000$	$2$
6	$2000$	$2$
7	$199$	$60$	$0 \leq l_{i} \leq r_{i} \leq 5000$
8	$200$	$50$	$0 \leq l_{i} \leq r_{i} \leq 5000$
9	$200$	$50$	$0 \leq l_{i} \leq r_{i} \leq 10^{9}$
10	$1999$	$500$	$0 \leq l_{i} \leq r_{i} \leq 5000$
11	$2000$	$400$	$0 \leq l_{i} \leq r_{i} \leq 5000$
12	$2000$	$500$	$0 \leq l_{i} \leq r_{i} \leq 10^{9}$
13	$30000$	$2000$	$0 \leq l_{i} \leq r_{i} \leq 100000$
14	$40000$	$1000$
15	$50000$	$15000$
16	$100000$	$20000$
17	$200000$	$20000$	$0 \leq l_{i} \leq r_{i} \leq 10^{9}$
18	$300000$	$50000$
19	$400000$	$90000$
20	$500000$	$200000$

时间限制： $3 s$

空间限制： $256 MB$

$256 MB$

首先对li,ri离散化。

如样例

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

经排序后——1 1 1 2 2 2 3 3 4 4 5 5。(3,5)就变成了(7,11)这样l,r就变成<=500000了。

将每条线段按长度排个序，变成！

2 2

1 2

3 4

3 5

1 4

1 5

接下来，维护一个双指针(two pointer)如

1 2 3 4 5 6(线段编号)

不断向右移动r直到有个点出现的次数>m为止如

1 2 3 4 5 6

r——>r

记录答案

再不断向右移动l，直到没有一个位置出现的次数>m如

l—->l

1 2 3 4 5 6

r——>r

记录答案，以此类推，直到r>n为止。

放代码跑

#include <cstdio>
#include <algorithm>
using namespace std;
struct data{int l,r,len;}a[500010];
struct tree{int l,r,max,lazy;}T[10001000];
int b[1000100],num,i,j,k,n,m,x,y,t;
bool cmp(const data&a,const data&b){return a.len<b.len;}
int find(int x){
    int l=1,r=num,ans;
    while (l<=r){
        int mid=(l+r)>>1;
        if (b[mid]>=x)ans=mid,r=mid-1;else l=mid+1;
    }
    return ans;
}
void build(int i,int l,int r){
    T[i].l=l;T[i].r=r;T[i].max=0;
    if (l==r)return;
    build(i*2,l,(l+r)>>1);build(i*2+1,((l+r)>>1)+1,r);
}
void pushdown(int i){
    if (T[i].lazy){
        T[i*2].max+=T[i].lazy;T[i*2+1].max+=T[i].lazy;
        T[i*2].lazy+=T[i].lazy;T[i*2+1].lazy+=T[i].lazy;
        T[i].lazy=0;
    }
}
void pushup(int i){T[i].max=max(T[i*2].max,T[i*2+1].max);}
void change(int i,int l,int r,int v){
    if (T[i].l==l&&T[i].r==r){T[i].max+=v;T[i].lazy+=v;return;}
    pushdown(i);
    int mid=(T[i].l+T[i].r)>>1;
    if (r<=mid)change(i*2,l,r,v);
    else if (l>mid)change(i*2+1,l,r,v);
    else {change(i*2,l,mid,v);change(i*2+1,mid+1,r,v);}
    pushup(i);
}
int main(){
    scanf("%d%d",&n,&m);
    for (i=1;i<=n;i++){
        scanf("%d%d",&x,&y);
        a[i].l=x;a[i].r=y;a[i].len=y-x;
        b[++num]=x;b[++num]=y;
    }
    sort(a+1,a+1+n,cmp);sort(b+1,b+1+num);
    for (i=1;i<=n;i++)a[i].l=find(a[i].l),a[i].r=find(a[i].r);
    build(1,1,num);
    int l=1,r=1,ans=1013687232;
    while(l<=n&&r<=n){
        while (T[1].max<m&&r<=n){change(1,a[r].l,a[r].r,1);r++;}
        if (T[1].max>=m)ans=min(ans,a[r-1].len-a[l].len);
        if (r>n)break;
        change(1,a[l].l,a[l].r,-1);l++;
    }
    printf("%d\n",ans==1013687232?-1:ans);
    return 0;
}