发现自己的离散化姿势一直有问题

今天终于掌握了正确的姿势

虽然这并不能阻挡我noip退役爆零的历史进程

还是先来看看离散化怎么写吧,我以前都是这么写的

for(std::set<int>::iterator it=s.begin();it!=s.end();it++)
ma[*it]=++tot;

这是使用\(set+map\)的离散化,但是显然巨大的常数是极大的劣势

正确的操作应该是这个样子

std::sort(a+1,a+n+1)
int tot=unique(a+1,a+n+1)-a-1;
for(re int i=1;i<=tot;i++)
ma[a[i]]=i;

\(unique\)能将一个有序数组去重,返回值是去重之后的尾地址,减去首地址就可以得到去重之后的数量了

之后再来看这道题,也是一道非常神的题

这道题告诉我们暴力通向错误的解,错误的解往往跟正解很接近了

首先\(O(n^2logn)\)还是比较好想的做法,我们先将所有的区间按照长度排序,之后我们定住一个区间作为长度最小的区间,强行套用线段树暴力覆盖之后的比它大的区间,直到这个区间有一个点被覆盖了\(m\)次,那么最后覆盖上去的区间的长度减去定住的最小值就是答案了

对所有的答案取一个\(min\)就好了

这样是显然不行的呀,我们得想个办法

于是我就想出来一种显然错误的做法

我大胆的猜测答案是单调的

少年谁给你的勇气

于是这个显然错误的做法是这样的,先将最短的区间一直覆盖,直到覆盖到符合条件为止,之后由于所谓的"答案单调",那么使下一个次短的区间符合条件的区间一定在前面的答案的后面

这样的复杂度均摊下来是对的

答案单调很有道理,但是凉凉了

很容易就能找到反例了

我们再覆盖第一个最短的区间的时候,由于我们只是在判断最短的区间是否满足条件,所以可能这个时候之后的某个区间就突然满足条件了,所以这个答案显然不是单调的

但是这个错误的思路和正解已经非常接近了,差别只有一点,我们判断的时候判断的不是最短的区间是否符合条件,而是判断整个区间内是否有点被覆盖了\(m\)次

这有什么道理呢,其实联想一下尺取法,和尺取法差不多

我们从最短的区间开始覆盖,可能覆盖的过程中满足条件的点并没有来自当前的区间,但是没有关系,我们直接用当前的区间作为最小的区间就行了,这样显然只会导致答案偏大,而真正的最小区间我们在后面也会取到

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstring>
#include<set>
#include<algorithm>
#include<map>
#define maxn 500005
#define re register
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
std::set<int> s;
std::map<int,int> ma;
struct node
{
int ll,rr,len;
int L,R;
}a[maxn];
int b[maxn<<1];
int n,m,N;
int ans=19999999999;
int l[4000005],r[4000005],d[4000005],tag[4000005];
inline int read()
{
char c=getchar();
int x=0;
while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9')
x=(x<<3)+(x<<1)+c-48,c=getchar();
return x;
}
inline int cmp(node K,node M)
{
return K.len<M.len;
}
void build(int x,int y,int i)
{
l[i]=x;
r[i]=y;
if(x==y) return;
int mid=x+y>>1;
build(x,mid,i<<1),build(mid+1,y,i<<1|1);
}
inline void pushdown(int i)
{
if(!tag[i]) return;
tag[i<<1]+=tag[i];
tag[i<<1|1]+=tag[i];
d[i<<1|1]+=tag[i];
d[i<<1]+=tag[i];
tag[i]=0;
}
void change(int x,int y,int v,int i)
{
if(x<=l[i]&&y>=r[i])
{
d[i]+=v;
tag[i]+=v;
return;
}
pushdown(i);
int mid=l[i]+r[i]>>1;
if(y<=mid) change(x,y,v,i<<1);
else if(x>mid) change(x,y,v,i<<1|1);
else change(x,y,v,i<<1),change(x,y,v,i<<1|1);
d[i]=max(d[i<<1],d[i<<1|1]);
}
int main()
{
n=read();
m=read();
for(re int i=1;i<=n;i++)
a[i].ll=read(),a[i].rr=read(),a[i].len=a[i].rr-a[i].ll,b[++N]=a[i].ll,b[++N]=a[i].rr;
std::sort(a+1,a+n+1,cmp);
std::sort(b+1,b+N+1);
int tot=std::unique(b+1,b+N+1)-b-1;
for(re int i=1;i<=tot;i++)
ma[b[i]]=i;
build(1,tot,1);
for(re int i=1;i<=n;i++) a[i].L=ma[a[i].ll],a[i].R=ma[a[i].rr];
int now=-1;
for(re int i=1;i<=n;i++)
{
change(a[i].L,a[i].R,1,1);
if(d[1]==m)
{
ans=a[i].len-a[1].len;
now=i;
break;
}
}
if(now==-1)
{
puts("-1");
return 0;
}
for(re int i=2;i<=n;i++)
{
change(a[i-1].L,a[i-1].R,-1,1);
while(d[1]<m)
{
if(now==n)
{
printf("%d\n",ans);
return 0;
}
now++;
change(a[now].L,a[now].R,1,1);
}
ans=min(ans,a[now].len-a[i].len);
}
printf("%d\n",ans);
return 0;
}

【[NOI2016]区间】的更多相关文章

  1. BZOJ_4653_[Noi2016]区间_线段树+离散化+双指针

    BZOJ_4653_[Noi2016]区间_线段树+离散化+双指针 Description 在数轴上有 n个闭区间 [l1,r1],[l2,r2],...,[ln,rn].现在要从中选出 m 个区间, ...

  2. [Noi2016]区间[离散化+线段树维护+决策单调性]

    4653: [Noi2016]区间 Time Limit: 60 Sec  Memory Limit: 256 MBSubmit: 621  Solved: 329[Submit][Status][D ...

  3. [NOI2016]区间 线段树

    [NOI2016]区间 LG传送门 考虑到这题的代价是最长边减最短边,可以先把边按长度排个序,双指针维护一个尺取的过程,如果存在包含某个点的区间数\(\ge m\),就更新答案并把左指针右移,这样做的 ...

  4. [BZOJ4653][NOI2016]区间 贪心+线段树

    4653: [Noi2016]区间 Time Limit: 60 Sec  Memory Limit: 256 MB Description 在数轴上有 n个闭区间 [l1,r1],[l2,r2],. ...

  5. 【BZOJ4653】[Noi2016]区间 双指针法+线段树

    [BZOJ4653][Noi2016]区间 Description 在数轴上有 n个闭区间 [l1,r1],[l2,r2],...,[ln,rn].现在要从中选出 m 个区间,使得这 m个区间共同包含 ...

  6. [NOI2016]区间 题解(决策单调性+线段树优化)

    4653: [Noi2016]区间 Time Limit: 60 Sec  Memory Limit: 256 MBSubmit: 1593  Solved: 869[Submit][Status][ ...

  7. Luogu P1712 [NOI2016]区间(线段树)

    P1712 [NOI2016]区间 题意 题目描述 在数轴上有 \(N\) 个闭区间 \([l_1,r_1],[l_2,r_2],...,[l_n,r_n]\) .现在要从中选出 \(M\) 个区间, ...

  8. 【题解】P1712 [NOI2016]区间(贪心+线段树)

    [题解]P1712 [NOI2016]区间(贪心+线段树) 一个observe是,对于一个合法的方案,将其线段长度按照从大到小排序后,他极差的来源是第一个和最后一个.或者说,读入的线段按照长度分类后, ...

  9. 洛谷P1712 [NOI2016]区间 尺取法+线段树+离散化

    洛谷P1712 [NOI2016]区间 noi2016第一题(大概是签到题吧,可我还是不会) 链接在这里 题面可以看链接: 先看题意 这么大的l,r,先来个离散化 很容易,我们可以想到一个结论 假设一 ...

  10. BZOJ4653: [Noi2016]区间

    传送门 UOJ上卡掉一个点,COGS上卡掉两个点..弃疗,不改了,反正BZOJ上过啦hhh 先把区间按长度递增排序.然后每次用线段树维护区间最大覆盖次数,用一个指针随便扫扫就行了. //NOI 201 ...

随机推荐

  1. grunt-contrib-cssmin使用指南

    原文:http://riny.net/2014/grunt-cssmin/ grunt-contrib-cssmin v0.7.0 使用cssmin压缩css文件 Getting Started 这个 ...

  2. Memcached 查询stats及各项状态解释

    一.两个最常用状态查询(掌握第一个就完全OK了) 1)查看状态:printf “stats\r\n” |nc 127.0.0.1 11211      2)模拟top命令查看状态:watch “ech ...

  3. jQuery比较时间处理过程(日期转换为时间戳比较)

    不管在JS中还是PHP中,终会遇到时间比较,但是很尴尬的是,举个例子,2017-04-12和2017-5-21是不一样或者说不能比较的,这在PHP和JS都是这种情况,这种情况怎么办呢,换成统一的格式比 ...

  4. Firebird Connection pool is full

    今天在做Firebird V3.0.3  x64 版本内存测试,本地PC上,准备开启800个事务(保持不关闭),每个事务做些事,尽量不释放内存. 每次测试当事务数达到时,就提示Connection p ...

  5. openlayers 各种图层,持续更新

    /*高德地图*/ var vectorLayerLine = new ol.layer.Tile({ source: new ol.source.XYZ({ urls: [ "http:// ...

  6. [android] 练习使用ListView(一)

    练习使用ListView,BaseAdapter,先展示文字的,再练习图片的 MainActivity.java package com.android.test; import android.ap ...

  7. java——程序的导出与导入

    导出: 选择项目,右击选择 最下面的properties——Resource——Location,就是你的项目所在地, 找到文件所在,拷贝到你的U盘中(或者直接点击项目直接拖到桌面)完成复制 导入: ...

  8. 来看看javaweb的自定义标签

    1.为什么需要自定义标签? jsp的简单标签其实就是jsp的自定义标签,主要作用就是移除jsp页面中的java代码,使得jsp页面只有标签和EL表达式,而没有java代码.利用自定义标签,可以使软件开 ...

  9. 二 Channel

    Java NIO的通道类似流,但又有些不同 既可以从通道中读取数据,也可以写数据到通道.但是流的读写通常是单向的 通道可以异步读写 通道中的数据通常总是要先读到一个Buffer,或者总是从Buffer ...

  10. spring 代理

    java动态代理实现 1. Java自带的动态代理,反射生成字节码 2. Cglib调用asm生成子类 spring 中代理实现 1. 如果类实现了接口,使用java动态代理 2. 没有实现接口,使用 ...