http://www.lydsy.com/JudgeOnline/problem.php?id=2716

样例输入

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

样例输出

1

2

————————————————————————————————

整整一个上午的时间,终于DEBUG出来了,原来是因为我数组开小了……

首先我们很容易想到一个三元组(t,x,y),其中t为操作时间。

默认最开始给定的一些点的操作为插入,且时间优先级均高于其他的点。

那么显然就是三维偏序,CDQ可以上了。

给t排个序,边归并x边树状数组记录y,然后查即可……

但是我们怎么求dis啊?

我们有种玄学的方法叫做分情况讨论。

分成四种情况:

对于(x,y)(x1,y1)的dis:

x+y-(x1+y1)      (x>x1,y>y1)

x-y-(x1-y1)        (x>x1,y<y1)

-x+y-(y1-x1)      (x<x1,y>y1)

-x-y-(-x1-y1)      (x<x1,y<y1)

其中括号内的东西就是我们要用树状数组来维护的东西,显然当括号内的东西最大的时候我们的距离有最小值。

那么就变成了维护区间最大值的问题了,显然树状数组是能够胜任这项工作的。

具体的实现方法看我代码吧……不好描述。

#include<cstdio>
#include<queue>
#include<cctype>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
typedef long long ll;
const int M=;
const int N=;
const int INF=*N;
inline int read(){
int X=,w=;char ch=;
while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();
return w?-X:X;
}
struct rem{
int t;
int x;
int y;
int pos;
}q[*M],tmp[*M];
int m,n,cnt=,ans[M],tree[N],maxy=-;
bool cmp(rem a,rem b){
return (a.t<b.t||(a.t==b.t&&a.x<b.x)||(a.t==b.t&&a.x==b.x&&a.y<b.y)||(a.t==b.t&&a.x==b.x&&a.y==b.y&&a.pos<b.pos));
}
inline int lowbit(int t){return t&(-t);}
void add(int x,int y){
for(int i=x;i<=maxy;i+=lowbit(i))tree[i]=max(tree[i],y);
return;
}
int query(int x){
int res=-INF;
for(int i=x;i>;i-=lowbit(i))res=max(res,tree[i]);
return res;
}
void cdq(int l,int r){
if(l>=r)return;
int mid=(l+r)>>;
cdq(l,mid);cdq(mid+,r);
for(int i=l,j=l,p=mid+;i<=r;i++){
if(j<=mid&&(p>r||q[j].x<=q[p].x))tmp[i]=q[j++];
else tmp[i]=q[p++];
}
///////////
for(int i=l;i<=r;i++){
q[i]=tmp[i];
if(q[i].t<=mid&&!q[i].pos)add(q[i].y,q[i].x+q[i].y);
if(q[i].t>mid&&q[i].pos)ans[q[i].pos]=min(ans[q[i].pos],q[i].x+q[i].y-query(q[i].y));
}
for(int i=l;i<=r;i++){
if(q[i].t<=mid&&!q[i].pos){
for(int j=q[i].y;j<=maxy;j+=lowbit(j))tree[j]=-INF;
}
}
///////////
for(int i=l;i<=r;i++){
if(q[i].t<=mid&&!q[i].pos)add(maxy-q[i].y,q[i].x-q[i].y);
if(q[i].t>mid&&q[i].pos)ans[q[i].pos]=min(ans[q[i].pos],q[i].x-q[i].y-query(maxy-q[i].y));
}
for(int i=l;i<=r;i++){
if(q[i].t<=mid&&!q[i].pos){
for(int j=maxy-q[i].y;j<=maxy;j+=lowbit(j))tree[j]=-INF;
}
}
///////////
for(int i=r;i>=l;i--){
if(q[i].t<=mid&&!q[i].pos)add(q[i].y,q[i].y-q[i].x);
if(q[i].t>mid&&q[i].pos)ans[q[i].pos]=min(ans[q[i].pos],-q[i].x+q[i].y-query(q[i].y));
}
for(int i=l;i<=r;i++){
if(q[i].t<=mid&&!q[i].pos){
for(int j=q[i].y;j<=maxy;j+=lowbit(j))tree[j]=-INF;
}
}
///////////
for(int i=r;i>=l;i--){
if(q[i].t<=mid&&!q[i].pos)add(maxy-q[i].y,-q[i].x-q[i].y);
if(q[i].t>mid&&q[i].pos)ans[q[i].pos]=min(ans[q[i].pos],-q[i].x-q[i].y-query(maxy-q[i].y));
}
for(int i=l;i<=r;i++){
if(q[i].t<=mid&&!q[i].pos){
for(int j=maxy-q[i].y;j<=maxy;j+=lowbit(j))tree[j]=-INF;
}
}
///////////
return;
}
void clear(){
maxy++;
for(int i=;i<=cnt;i++)ans[i]=INF;
for(int i=;i<=maxy;i++)tree[i]=-INF;
return;
}
int main(){
n=read();
m=read();
for(int i=;i<=n;i++){
q[i].x=read()+;
q[i].y=read()+;
q[i].t=i;
maxy=max(maxy,q[i].y);
}
for(int i=n+;i<=m+n;i++){
int w=read();
if(w==){
q[i].x=read()+;
q[i].y=read()+;
q[i].t=i;
maxy=max(maxy,q[i].y);
}else{
q[i].x=read()+;
q[i].y=read()+;
q[i].t=i;
q[i].pos=++cnt;
maxy=max(maxy,q[i].y);
}
}
sort(q+,q+m+n+,cmp);
clear();
cdq(,m+n);
for(int i=;i<=cnt;i++)printf("%d\n",ans[i]);
return ;
}

BZOJ2716:[Violet 3]天使玩偶——题解的更多相关文章

  1. bzoj2648SJY摆棋子&&bzoj2716[Violet 3]天使玩偶*

    bzoj2648SJY摆棋子 bzoj2716[Violet 3]天使玩偶 题意: 棋盘上有n个棋子,现在有m个操作,一种是加棋子,一种是查询离某个点最近的棋子.n,m≤500000. 题解: 先将已 ...

  2. [BZOJ2716] [Violet 3]天使玩偶(CDQ分治)

    [BZOJ2716] [Violet 3]天使玩偶(CDQ分治) 题面 Ayu 在七年前曾经收到过一个天使玩偶,当时她把它当作时间囊埋在了地下.而七年后 的今天,Ayu 却忘了她把天使玩偶埋在了哪里, ...

  3. 【kd-tree】bzoj2716 [Violet 3]天使玩偶

    #include<cstdio> #include<cmath> #include<algorithm> using namespace std; #define ...

  4. BZOJ2716 [Violet 3]天使玩偶 【CDQ分治】

    题目 输入格式 输出格式 输入样例 //样例太长就不贴了.... 输出样例 //见原题 提示 题解 我们将曼哈顿距离式子中的绝对值去掉,每次只考虑x,y比当前点小的更新答案. 为了使所有点都对答案进行 ...

  5. bzoj2716: [Violet 3]天使玩偶

    #include <iostream> #include <cstdio> #include <cstring> #include <cmath> #i ...

  6. BZOJ2716: [Violet 3]天使玩偶(KD-Tree)

    Description Input Output Sample Input 100 100 81 23 27 16 52 58 44 24 25 95 34 2 96 25 8 14 97 50 97 ...

  7. BZOJ2648: SJY摆棋子&&2716: [Violet 3]天使玩偶

    BZOJ2648: SJY摆棋子 BZOJ2716: [Violet 3]天使玩偶 BZOJ氪金无极限... 其实这两道是同一题. 附上2648的题面: Description 这天,SJY显得无聊. ...

  8. BZOJ 2716: [Violet 3]天使玩偶

    2716: [Violet 3]天使玩偶 Time Limit: 80 Sec  Memory Limit: 128 MBSubmit: 1473  Solved: 621[Submit][Statu ...

  9. BZOJ 2716: [Violet 3]天使玩偶( CDQ分治 + 树状数组 )

    先cdq分治, 然后要处理点对答案的贡献, 可以以询问点为中心分成4个区域, 然后去掉绝对值(4种情况讨论), 用BIT维护就行了. --------------------------------- ...

随机推荐

  1. JS基础,课堂作业,计算器

    网页内的简单计算器 <script> var a = parseInt(prompt("请输入第一个数字:")); var b = parseInt(prompt(&q ...

  2. abp 关闭审计日志

    在WebProject.Web.Mvc项目中的Startup文件夹下中的WebProjectWebMvcModule.cs public override void PreInitialize() { ...

  3. Maven学习(六)-----Maven仓库的详细介绍

    Maven仓库的详细介绍 在Maven中,任何一个依赖.插件或者项目构建的输出,都可以称之为构件.Maven在某个统一的位置存储所有项目的共享的构件,这个统一的位置,我们就称之为仓库.(仓库就是存放依 ...

  4. 二、利用EnterpriseFrameWork快速开发Web系统(B/S)

    EnterpriseFrameWork框架实例源代码下载: 实例下载 本章通过一个开发实例来讲解Web系统的开发经过,以及会碰到的一些问题.实例功能就是业务系统中最常见的增.删.改.查来演示,用一个界 ...

  5. 技本功丨利用 Atomic 构建 React 项目工作流,so easy!

    近日刷微博,#2018年结婚率创新低#荣登热门话题NO.1,沪浙最不积极. 生活压力越大,缺爱的人也越来越多...据本萌的不完全观察,程序猿虽然是压力加成的职业,在袋鼠云还是有不少早早脱了单.至于,脱 ...

  6. Python基础灬函数补充(作用域,迭代器,生成器)

    变量作用域 函数里面操作外部变量时,作用域仅限于函数里面. var1 = 123 def func(): var1 = 456 print("函数里:", var1) func() ...

  7. redis与mysql性能对比、redis缓存穿透、缓存雪崩

    写在开始 redis是一个基于内存hash结构的缓存型db.其优势在于速读写能力碾压mysql.由于其为基于内存的db所以存储数据量是受限的. redis性能 redis读写性能测试redis官网测试 ...

  8. 通过Nrgok映射外网调试微信

    一.注册账号 注册地址:http://www.ngrok.cc/login 登录系统,新增域名 二.下载客户端,修改配置文件 修改ngrok.cfg auth_token值登录平台管理系统可查看 su ...

  9. 【最小生成树+LCA】Imperial roads

    http://codeforces.com/gym/101889 I 先跑一遍最小生成树,把经过的边和答案记录下来 对于每个询问的边,显然如果处于MST中,答案不变 如果不在MST中,假设这条边连上了 ...

  10. C语言之goto浅析

    1.  读代码时遇了的疑惑点: static int do_bind(const char *host, int port, int protocol, int *family) { int fd; ...