原题链接: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
我们把每个点看成三元组: (时间,x,y) ( 时 间 , x , y ) ,对于一个点 (ti,xi,yi) ( t i , x i , y i ) 当我们只考虑 xj≤xi x j ≤ x i 且 yj≤yi y j ≤ y i 的情况时,我们要找的点就是, xj≤xi,yj≤yi,tj≤ti x j ≤ x i , y j ≤ y i , t j ≤ t i 且 xj+yj x j + y j 最大的点,那么就变成了一个三维偏序。
对于其他的情况,我们只需要略微变换一下坐标,就能变成上面的那种情况,做四遍 cdq c d q 即可。
这题真 TM T M 的卡常,博主上了 buff b u f f 都血 TLE T L E 。。。
代码放在这儿仅供参考,开 O2 O 2 才能 A A ,等博主学完 KD−Tree K D − T r e e 再一雪前耻。。。
#include
#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.xif(a.y!=b.y)return a.yreturn a.opint 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(base1;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]printf("%d\n",ans[i]);
}
int main()
{
in();ac();
return 0;
}