正题

题目链接:https://www.luogu.com.cn/problem/P4169

题目大意

nnn个点，然后每次操作

加一个新的点
询问一个点更近的点

解题思路

定义tit_iti表示第几个操作，开始就有的点tit_iti为000
假设最近的点在左上角，那么有要求ti<tj,xi<xj,yi<yjt_i<t_j,x_i<x_j,y_i<y_jti<tj,xi<xj,yi<yj构成了一个三维偏序。
对于一个点对(i,j)(i,j)(i,j)，权值是xi−xj+yi−yj⇒xi+yi−(xj+yj)x_i-x_j+y_i-y_j\Rightarrow x_i+y_i-(x_j+y_j)xi−xj+yi−yj⇒xi+yi−(xj+yj)。用树状数组维护xi+yix_i+y_ixi+yi的前缀maxmaxmax即可

然后旋转一下坐标就行了，时间复杂度O(nlog⁡2n)O(n\log^2n)O(nlog2n)

codecodecode

/*#pragma GCC optimize(2)
%:pragma GCC optimize(3)
%:pragma GCC optimize("Ofast")*/
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
#define lowbit(x) (x&-x)
using namespace std;
const int N=6e5+10,K=1e6+1;
struct node{int x,y,id,t;
}a[N],b[N];
int n,m,tot,cnt,t[K+10],ans[N];
int read() {int x=0,f=1; char c=getchar();while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}while(isdigit(c)) x=(x<<1)+(x<<3)+c-48,c=getchar();return x*f;
}
void print(int x){if (x>9) print(x/10); putchar(x%10+48); return;
}
void Change(int x,int val){x++;while(x<=K){t[x]=max(t[x],val);x+=lowbit(x);}return;
}
void Clear(int x){x++;while(x<=K){t[x]=-1e7;x+=lowbit(x);} return;
}
int Ask(int x){ int ans=-1e7;x++;while(x){ans=max(ans,t[x]);x-=lowbit(x);}return ans;
}
bool cmp(node x,node y)
{return x.t<y.t;}
void merge(int l,int mid,int r){int L=l,R=mid+1,z=l;while(z<=r){if((a[L].x<=a[R].x||R>r)&&L<=mid)b[z++]=a[L++];else if(a[L].x>a[R].x||L>mid)b[z++]=a[R++];}for(int i=l;i<=r;i++)a[i]=b[i];return;
}
void cdq(int l,int r){if(l==r)return;int mid=(l+r)>>1;cdq(l,mid);cdq(mid+1,r);int z=l;for(int i=mid+1;i<=r;i++){while(z<=mid&&a[z].x<=a[i].x){if(a[z].id==0)Change(a[z].y,a[z].x+a[z].y);z++;}if(a[i].id!=0)ans[a[i].id]=min(ans[a[i].id],a[i].x+a[i].y-Ask(a[i].y));}for(int i=l;i<z;i++)Clear(a[i].y);merge(l,mid,r);return;
}
int main()
{n=read();m=read(); for(int i=1;i<=n;i++){a[i].x=read()+1;a[i].y=read();a[i].t=i;}for(int i=1;i<=m;i++){int op=read();a[++n].t=n;a[n].x=read()+1;a[n].y=read(); if(op==2)a[n].id=++cnt;}sort(a+1,a+1+n,cmp);memset(t,0xcf,sizeof(t));memset(ans,0x3f,sizeof(ans));cdq(1,n);for(int i=1;i<=n;i++)a[i].x=K-a[i].x;sort(a+1,a+1+n,cmp);cdq(1,n);for(int i=1;i<=n;i++)a[i].y=K-a[i].y;sort(a+1,a+1+n,cmp);cdq(1,n);for(int i=1;i<=n;i++)a[i].x=K-a[i].x;sort(a+1,a+1+n,cmp);cdq(1,n);for(int i=1;i<=cnt;i++)print(ans[i]),putchar('\n');
}

P4169-[Violet]天使玩偶/SJY摆棋子【CDQ分治】相关推荐

P4169 [Violet]天使玩偶/SJY摆棋子 [CDQ分治]
天使玩偶天使玩偶天使玩偶 Ayu 在七年前曾经收到过一个天使玩偶,当时她把它当作时间囊埋在了地下.而七年后的今天,Ayu 却忘了她把天使玩偶埋在了哪里,所以她决定仅凭一点模糊的记忆来寻找它. 我们把 ...
[Luogu4169] [Violet]天使玩偶/SJY摆棋子 [cdq分治/k-d tree]
[Link\frak{Link}Link] 考虑用 cdq 分治解决? 显然依旧分成修改跟查询,考虑一维排序二维cdq?啊但是那显然不太对,cdq能解决三维偏序, 但是她不能解决最近点对啊.平面最近点 ...
P4169 [Violet]天使玩偶/SJY摆棋子（CDQ分治、暴力）
解析之前用KDtree做的一道题由于懒不想再码一遍了考虑CDQ分治关键就是如何拿掉绝对值如果只维护左下角的,显然就是一个经典的三维偏序问题了但是本题不一定在左下角,也可能在左上.右下.右上 ...
题解洛谷 P4169 [Violet]天使玩偶/SJY摆棋子【CDQ分治】
题目链接题意:平面上原有一些点,支持动态加点,动态查询与某个点曼哈顿距离最小的点的距离. 这题据说是有 KD-Tree 和 CDQ 分治两种做法,又据说 KDT 会被卡,于是我采用 CDQ (其实是 ...
[Violet]天使玩偶/SJY摆棋子
传送门题目大意就是先给定n个点,然后m个操作,每次添加一个点或者查询离一个点最近的一个点是哪个点. 因为要求的是曼哈顿距离,所以我们先把式子转化一下,就是\((x_i + y_i) - (x_j + ...
【洛谷P4169】天使玩偶/SJY摆棋子【CDQ分治】
传送门题意:动态加点,给定点询问曼哈顿距离最近的点 N,M≤3e5,x,y≤1e6N,M \leq 3e5,x,y \leq 1e6N,M≤3e5,x,y≤1e6 经(kan)过(le)分(ti)析 ...
[Luogu4169][Violet]天使玩偶/SJY摆棋子
luogu 题意一个平面上有\(n\)个点,\(m\)次操作,每次新增一个点,或者是询问离某个点最近的点的距离.这里的距离是曼哈顿距离. \(n,m\le3*10^5\) sol 写一发\(CDQ\ ...
天使玩偶/SJY摆棋子
P4169 [Violet]天使玩偶/SJY摆棋子 CDQ分治的题目. 我们发现题目要我们求的\(|A_x-B_x|+|A_y-B_y|\)的绝对值号比较恶心. 试想一下怎么去掉如果所有的点都在我们 ...
[BZOJ2716/2648][Violet 3]天使玩偶/SJY摆棋子[KDtree]
KDtree干这个复杂度是不对的,重构不一定有作用解释一下的话,因为复杂度是跟size相关的,所以重构作用不大,KDtree在查询最近点对中的作用仅仅是剪枝,可以构造数据使得他遍历O(n)个节点 h ...

P4169-[Violet]天使玩偶/SJY摆棋子【CDQ分治】

正题

题目大意

解题思路

codecodecode

P4169-[Violet]天使玩偶/SJY摆棋子【CDQ分治】相关推荐

最新文章

热门文章