洛谷P4169 天使玩偶 CDQ分治

还是照着CDQ的思路来。

但是有一些改动：

• 要求4个方向的，但是可爱的CDQ分治只能求在自己一个角落方向上的。怎么办？旋转！做4次就好了。

• 统计的不是和，而是——max！理由如下：

　　　　设当前点是(x,y)，目标点是(x',y')，那么所求的|x-x'|+|y-y'|首先用旋转大法化为x-x'+y-y'，然后我们发现这个东西其实就是x+y-x'-y'=(x+y)-(x'+y')，而x+y我们是已知的。所以我们求一下max(x'+y')即可。具体实现是对树状数组魔改。

然后交上去发现狂T不止...

疯狂优化！首先把树状数组的max从作死的M改成了所有坐标中最大的mm，发现不行，跑去看题解。

发现不能每次sort，而且还要删点。

然后学习了一波操作之后，成功的到达了91分...还是有3个点死活过不了。把题解拿下来一交，发现居然也T了一个点，哈哈哈。

听从雨菲的建议，果断开O2，A了。

看代码理解删点优化和魔改树状数组。

 // luogu-judger-enable-o2
 #include <cstdio>
 #include <algorithm>
 using namespace std;
 const int M = ,N = ,INF=0x3f3f3f3f;
 inline int read()
 {
     int ans=;char ch=getchar();
     while(ch<''||ch>'') ch=getchar();
     while(ch>=''&&ch<='') ans=(ans<<)+(ans<<)+ch-'',ch=getchar();
     return ans;
 }
 inline int lowbit(int x){return x&(-x);}
 int ta[M+],n,mm,ans[N];int p,rx,ry;
 inline void add(int x,int a)
 {
     if(x<=) return;
     for(int i=x;i<=ry;i+=lowbit(i)) ta[i]=max(ta[i],a);
     return;
 }
 inline int getsum(int x)
 {
     if(x<=)return ;int ans=-INF;
     for(int i=x;i>;i-=lowbit(i)) ans=max(ans,ta[i]);
     return ans;
 }
 inline void del(int x)
 {
     if(x<=) return;
     for(int i=x;i<=ry;i+=lowbit(i)) ta[i]=-INF;
     return;
 }
 struct Node
 {
     int x,y,k=-INF,t,type;
 }node[N],temp[N],f[N];
 inline bool cmp_t(Node a,Node b){return a.t<b.t;}
 inline void rotate()
 {
     for(int i=;i<=n;i++)
     {
         int t=f[i].x;
         f[i].x=mm-f[i].y;
         f[i].y=t;
     }
     return;
 }
 void CDQ(int l,int r)
 {
     if(l==r) return;
     int mid=(l+r)>>;
     CDQ(l,mid);CDQ(mid+,r);
     int i=l,j=mid+,t=;
     while(i<=mid||j<=r)
     {
         if(j>r||(i<=mid&&node[i].x<=node[j].x))
         {
             if(node[i].type==) add(node[i].y,node[i].x+node[i].y);
             temp[++t]=node[i++];
         }
         else
         {
             if(node[j].type>) {node[j].k=max(node[j].k,getsum(node[j].y));
             ans[node[j].t]=min(ans[node[j].t],node[j].x+node[j].y-node[j].k);}
             ///if(node[j].type>1) node[j].k=min(node[j].k,(node[j].x+node[j].y-getsum(node[j].y)));
             temp[++t]=node[j++];
         }
     }
     t=;
     for(j=l;j<i;j++) if(node[j].type==) del(node[j].y);
     for(i=l;i<=r;i++) node[i]=temp[++t];
     return;
 }

 inline void Delete()
 {
     rx=ry=p=;
     for(int i=;i<=n;i++) if(f[i].type==) rx=max(rx,f[i].x),ry=max(ry,f[i].y);
     for(int i=;i<=n;i++) if(f[i].x<=rx&&f[i].y<=ry) node[++p]=f[i];
     return;
 }
 int main()
 {
     fill(ta,ta+M,-INF);
     int n2=read(),n3=read(),x,y,flag;
     for(int i=;i<=n2;i++)
     {
         x=read()+;y=read()+;
         f[++n].x=x;
         f[n].y=y;
         f[n].type=;
         f[n].t=n;
         mm=max(mm,x);
         mm=max(mm,y);
     }
     for(int i=;i<=n3;i++)
     {
         flag=read();x=read()+;y=read()+;
         f[++n].x=x;
         f[n].y=y;
         f[n].type=flag;
         f[n].t=n;
         mm=max(mm,x);
         mm=max(mm,y);
     }
     mm++;
     fill(ans+,ans+n+,INF);
     for(y=;y<=;y++)
     {
         Delete();
         ry++;
         //for(int i=1;i<=p;i++) printf("%d %d %d\n",node[i].type,node[i].x,node[i].y);
         CDQ(,p);
         //for(int i=1;i<=n;i++)if(f[i].type==2)printf("%d ",ans[i]);
         //printf("\n");
         rotate();
     }

     for(int i=;i<=n;i++)
     {
         if(f[i].type==) printf("%d\n",ans[i]);
     }

     return ;
 }
 /**
 1 3
 3 1
 2 5 5
 1 6 4
 2 5 5
 */

伪AC代码：