原题链接:https://www.lydsy.com/JudgeOnline/problem.php?id=2716
Ayu 在七年前曾经收到过一个天使玩偶,当时她把它当作时间囊埋在了地下。而七年后 的今天,Ayu 却忘了她把天使玩偶埋在了哪里,所以她决定仅凭一点模糊的记忆来寻找它。
我们把 Ayu 生活的小镇看作一个二维平面坐标系,而 Ayu 会不定时地记起可能在某个点 (xmy) 埋下了天使玩偶;或者 Ayu 会询问你,假如她在 (x,y) ,那么她离近的天使玩偶可能埋下的地方有多远。
因为 Ayu 只会沿着平行坐标轴的方向来行动,所以在这个问题里我们定义两个点之间的距离为dist(A,B)=|Ax-Bx|+|Ay-By|。其中 Ax 表示点 A的横坐标,其余类似。
第一行包含两个整数n和m ,在刚开始时,Ayu 已经知道有n个点可能埋着天使玩偶, 接下来 Ayu 要进行m 次操作
接下来n行,每行两个非负整数 (xi,yi),表示初始n个点的坐标。
再接下来m 行,每行三个非负整数 t,xi,yi。
如果t=1 ,则表示 Ayu 又回忆起了一个可能埋着玩偶的点 (xi,yi) 。
如果t=2 ,则表示 Ayu 询问如果她在点 (xi,yi) ,那么在已经回忆出来的点里,离她近的那个点有多远
对于每个t=2 的询问,在单独的一行内输出该询问的结果。
2 3
1 1
2 3
2 1 2
1 3 3
2 4 2
1
2
n,m<=300 000
xi,yi<=1 000 000
我们把每个点看成三元组:,对于一个点当我们只考虑且的情况时,我们要找的点就是,且最大的点,那么就变成了一个三维偏序。
对于其他的情况,我们只需要略微变换一下坐标,就能变成上面的那种情况,做四遍即可。
这题真的卡常,博主上了都血。。。
代码放在这儿仅供参考,开才能,等博主学完再一雪前耻。。。
#include<bits/stdc++.h> #define M(v) mx[v]=max(mx[v<<1],mx[v<<1|1]) using namespace std; const int M=2e6,INF=1e9; struct sd{int op,x,y,t;}; bool operator <(sd a,sd b) { if(a.x!=b.x)return a.x<b.x; if(a.y!=b.y)return a.y<b.y; return a.op<b.op; } sd ope[M],tmp[M]; int n,m,ans[M],mx[M],MX,tot,base=1,i,r; char c; inline char nc() { static const int buflen=1e6; static char buf[buflen],*p1=buf,*p2=buf; return p1==p2&&(p2=(p1=buf)+fread(buf,1,buflen,stdin),p1==p2)?EOF:*p1++; } int read() { r=0; while(!isdigit(c))c=nc(); while(isdigit(c)){r=(r<<1)+(r<<3)+c-'0',c=nc();} return r; } void in() { int a,b,c; scanf("%d%d",&n,&m); for(int i=1;i<=n;++i)scanf("%d%d",&a,&b),++a,++b,ope[i]=(sd){0,a,b,i},MX=max(MX,max(a,b)); for(int i=n+1;i<=n+m;++i)scanf("%d%d%d",&a,&b,&c),++b,++c,ope[i]=(sd){a-1,b,c,i},MX=max(MX,max(b,c)); } void clear(int v){v+=base;for(;v;v>>=1)mx[v]=0;} void add(int v,int s){v+=base;mx[v]=max(mx[v],s);for(v>>=1;v;v>>=1)M(v);} int query(int ri) { int le=base,ans=0;ri+=base+1; for(;le^ri^1;le>>=1,ri>>=1) { if(le&1^1)ans=max(ans,mx[le+1]); if(ri&1)ans=max(ans,mx[ri-1]); } return ans?ans:-INF; } void cdq(int le,int ri) { if(le==ri)return; int mid=le+ri>>1,p1=le,p2=mid+1; for(i=le;i<=ri;++i) { if(ope[i].t<=mid&&!ope[i].op)add(ope[i].y,ope[i].x+ope[i].y); if(ope[i].t>mid&&ope[i].op)ans[ope[i].t]=min(ope[i].x+ope[i].y-query(ope[i].y),ans[ope[i].t]); } for(i=le;i<=ri;++i)if(ope[i].t<=mid&&!ope[i].op)clear(ope[i].y); for(i=le;i<=ri;++i) if(ope[i].t<=mid)tmp[p1++]=ope[i]; else tmp[p2++]=ope[i]; for(i=le;i<=ri;++i)ope[i]=tmp[i]; cdq(le,mid);cdq(mid+1,ri); } void ac() { MX++;while(base<MX)base<<=1;tot=n+m; fill(ans+n+1,ans+1+tot,INF); sort(ope+1,ope+1+tot); cdq(1,tot); for(i=1;i<=tot;++i)ope[i].x=MX-ope[i].x;sort(ope+1,ope+1+tot);cdq(1,tot); for(i=1;i<=tot;++i)ope[i].y=MX-ope[i].y;sort(ope+1,ope+1+tot);cdq(1,tot); for(i=1;i<=tot;++i)ope[i].x=MX-ope[i].x;sort(ope+1,ope+1+tot);cdq(1,tot); for(i=n+1;i<=tot;++i)if(ans[i]<INF)printf("%d\n",ans[i]); } int main() { in();ac(); return 0; }