bzoj 4653: [Noi2016]区间

额，是不是一到了晚上IQ就--；

这个题一开始完全没有思路。（貌似脑子就没动一下）

%了一下题解。

大概是决策是有单调性的，因为要去区间长度差最小，所以接排个序，然后扫描右端点，找出满足有点被覆盖m次的最右的左端点就好。

然后判断是不是有覆盖m个点的用线段树维护一下。

（23333，吐槽，为什么离散化二分的时候，把左端点搞成0，右端点搞成L就会RE？？excuse me？！！）

 1 #include<bits/stdc++.h>
 2 #define LL long long
 3 #define inf 0x7fffffff
 4 using namespace std;
 5 inline int ra()
 6 {
 7     int x=0,f=1; char ch=getchar();
 8     while (ch<‘0‘ || ch>‘9‘) {if (ch==‘-‘) f=-1; ch=getchar();}
 9     while (ch>=‘0‘ && ch<=‘9‘) {x=x*10+ch-‘0‘; ch=getchar();}
10     return x*f;
11 }
12
13 const int maxn=500005;
14
15 struct seg{int l,r,tag,mx;}t[maxn<<3];
16 struct data{int x,y,len;}a[maxn];
17
18 int L,n,m,b[maxn<<1],tot,ans=inf;
19 bool cmp_len(data a, data b){return a.len<b.len;}
20 int find(int x)
21 {
22     int l=1,r=L-1;
23     while (l<=r)
24     {
25         int mid=l+r>>1;
26         if (b[mid]==x) return mid;
27         else if (b[mid]<x) l=mid+1;
28         else r=mid-1;
29     }
30 }
31 void build(int k, int l, int r)
32 {
33     t[k].l=l; t[k].r=r;
34     if (l==r) return;
35     int mid=l+r>>1;
36     build(k<<1,l,mid); build(k<<1|1,mid+1,r);
37 }
38 void pushdown(int k)
39 {
40     int tmp=t[k].tag; t[k].tag=0;
41     t[k<<1].tag+=tmp; t[k<<1|1].tag+=tmp;
42     t[k<<1].mx+=tmp; t[k<<1|1].mx+=tmp;
43 }
44 void update(int k) {t[k].mx=max(t[k<<1].mx,t[k<<1|1].mx);}
45 void change(int k, int x, int y, int val)
46 {
47     int l=t[k].l,r=t[k].r;
48     if (x==l && y==r)
49     {
50         t[k].tag+=val; t[k].mx+=val;
51         return;
52     }
53     if (t[k].tag) pushdown(k);
54     int mid=l+r>>1;
55     if (y<=mid) change(k<<1,x,y,val);
56     else if (x>mid) change(k<<1|1,x,y,val);
57     else change(k<<1,x,mid,val),change(k<<1|1,mid+1,y,val);
58     update(k);
59 }
60 int main()
61 {
62 //    freopen("bzoj4653_data.in","r",stdin);
63 //    freopen("bzoj4653_data.out","w",stdout);
64     n=ra(); m=ra();
65     if (!m) {puts("0");return 0;}
66     for (int i=1; i<=n; i++)
67     {
68         a[i].x=ra(),a[i].y=ra(),a[i].len=a[i].y-a[i].x;
69         b[++tot]=a[i].x; b[++tot]=a[i].y;
70     }
71     sort(a+1,a+n+1,cmp_len); sort(b+1,b+tot+1);
72     L=unique(b+1,b+tot+1)-b;
73     for (int i=1; i<=n; i++) a[i].x=find(a[i].x),a[i].y=find(a[i].y);
74     int now=1; build(1,1,L);
75     for (int i=1; i<=n; i++)
76     {
77         change(1,a[i].x,a[i].y,1);
78         while (t[1].mx>=m)
79         {
80             ans=min(ans,a[i].len-a[now].len);
81             change(1,a[now].x,a[now].y,-1);
82             now++;
83         }
84     }
85     printf("%d\n",ans==inf?-1:ans);
86     return 0;
87 }