BZOJ4653: [Noi2016]区间

BZOJ4653

（蒟蒻只想到了$n^2$暴力。。）
如果确定了最大值与最小值，那么就是所有长度在这之间的区间中选取答案。
于是就可以先按区间长度排序，固定右端点（最大值）向右扫，然后类似于单调队列，更新答案。
要离散化。

【代码】

#include <cstdio>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <string>
#include <algorithm>
#include <cmath>
#define N 500005
#define INF 0x7fffffff
using namespace std;

int n,m,tot,ans=INF;
int hash[N<<1];

int read()
{
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

class Section{
    public:
        int L,R,Len;
}E[N]; 

bool operator <(Section A,Section B){
    return A.Len<B.Len;
}

class Seg_Tree{
    public:
        int l,r,mx,tag;
}e[N<<3];

int Find(int x)
{
    int l=1,r=tot,rtn;
    while(l<=r)
    {
        int mid=l+r>>1;
        if(hash[mid]>=x) rtn=mid,r=mid-1;
        else l=mid+1;
    }
    return rtn;
}

void Build(int p,int l,int r)
{
    e[p].l=l,e[p].r=r;
    if(l==r) return;
    int mid=l+r>>1;
    Build(p<<1,l,mid);Build(p<<1|1,mid+1,r);
}

void pushdown(int p)
{
    e[p<<1].mx+=e[p].tag;e[p<<1|1].mx+=e[p].tag;
    e[p<<1].tag+=e[p].tag;e[p<<1|1].tag+=e[p].tag;
    e[p].tag=0; 
}

void Update(int p,int x,int y,int z)
{
    int l=e[p].l,r=e[p].r,mid=l+r>>1;
    if(x==l&&y==r)
    {
        e[p].mx+=z;e[p].tag+=z;
        return;
    }
    if(e[p].tag) pushdown(p);
    if(y<=mid) Update(p<<1,x,y,z);
    else if(x>mid) Update(p<<1|1,x,y,z);
    else Update(p<<1,x,mid,z),Update(p<<1|1,mid+1,y,z);
    e[p].mx=max(e[p<<1].mx,e[p<<1|1].mx);
}

int main() 
{
    n=read(),m=read();
    for(int i=1;i<=n;i++)
    {
        E[i].L=read(),E[i].R=read();
        E[i].Len=E[i].R-E[i].L;
        hash[++tot]=E[i].L,hash[++tot]=E[i].R;
    }
    sort(E+1,E+1+n);
    sort(hash+1,hash+1+tot);
    tot=unique(hash+1,hash+1+tot)-hash-1;
    for(int i=1;i<=n;i++)
    {
        E[i].L=Find(E[i].L);
        E[i].R=Find(E[i].R);
    } 
    Build(1,1,tot);
    for(int i=1,j=1;i<=n;i++)
    {
        Update(1,E[i].L,E[i].R,1);
        while(e[1].mx>=m)
        {
            ans=min(ans,E[i].Len-E[j].Len);
            Update(1,E[j].L,E[j].R,-1);j++;
        }
    }
    if(ans==INF) ans=-1;
    printf("%d\n",ans);
    return 0;
}