BZOJ2716[Violet] 天使玩偶/SJY摆棋子

原题链接：https://www.lydsy.com/JudgeOnline/problem.php?id=2716

天使玩偶/SJY摆棋子

题目描述

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 的询问，在单独的一行内输出该询问的结果。

输入输出样例

输入样例#1：

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

输出样例#1：

1
2

说明

n,m<=300 000

xi,yi<=1 000 000

题解

我们把每个点看成三元组：$(时间,x,y)$，对于一个点$(t_i,x_i,y_i)$当我们只考虑$x_j\le x_i$且$y_j\le y_i$的情况时，我们要找的点就是，$x_j\le x_i,y_j\le y_i,t_j\le t_i$且$x_j+y_j$最大的点，那么就变成了一个三维偏序。

对于其他的情况，我们只需要略微变换一下坐标，就能变成上面的那种情况，做四遍$cdq$即可。

代码

这题真$^{TM}$的卡常，博主上了$\mathcal{buff}$都血$\mathcal{TLE}$。。。

代码放在这儿仅供参考，开$O2$才能$A$，等博主学完$KD-Tree$再一雪前耻。。。

#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;
}