SJY摆棋子

一、题目

点此看题

二、解法

也是板子题，不过要多写一个插入和替罪羊的拍扁重建，实现看代码好了。

#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
const int M = 1000005;
const int inf = 1e9;
const double al = 0.75;
int read()
{
    int x=0,flag=1;char c;
    while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
    while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
    return x*flag;
}
int n,m,k,opt,top,ans,rt,cnt,siz[M],b[M];
int mi[M][2],mx[M][2],ls[M],rs[M];
struct node
{
    int x[2];
}a[M],v[M];
int newnode()
{
    if(top) return b[top--];
    return ++cnt;
}
void up(int x)
{
    for(int i=0;i<2;i++)
        mi[x][i]=mx[x][i]=v[x].x[i];
    if(ls[x])
        for(int i=0;i<2;i++)
        {
            mx[x][i]=max(mx[x][i],mx[ls[x]][i]);
            mi[x][i]=min(mi[x][i],mi[ls[x]][i]);
        }
    if(rs[x])
        for(int i=0;i<2;i++)
        {
            mx[x][i]=max(mx[x][i],mx[rs[x]][i]);
            mi[x][i]=min(mi[x][i],mi[rs[x]][i]);
        }
    siz[x]=siz[ls[x]]+siz[rs[x]]+1;
}
int cmp(node a,node b)
{
    return a.x[opt]<b.x[opt];
}
void build(int &x,int l,int r,int wd)
{
    if(l>r) return ;
    x=newnode();
    int mid=(l+r)>>1;
    opt=wd;
    nth_element(a+l,a+mid,a+r+1,cmp);
    v[x]=a[mid];
    build(ls[x],l,mid-1,wd^1);
    build(rs[x],mid+1,r,wd^1);
    up(x);
}
int Abs(int x)
{
    return x>0?x:-x;
}
int mht(node a,node b)
{
    return Abs(a.x[0]-b.x[0])+Abs(a.x[1]-b.x[1]);
}
int dis(node a,int b)
{
    int r=0;
    for(int i=0;i<2;i++)
    {
        r+=max(0,a.x[i]-mx[b][i]);
        r+=max(0,mi[b][i]-a.x[i]);
    }
    return r;
}
void dfs(int &x)
{
    if(!x) return ;
    dfs(ls[x]);
    b[++top]=x;
    a[++n]=v[x];
    dfs(rs[x]);
    x=0;
}
void check(int &x,int wd)
{
    if(siz[ls[x]]>al*siz[x] || siz[rs[x]]>al*siz[x])
    {
        n=0;dfs(x);
        build(x,1,n,wd);
    }
}
void ins(int &x,int wd,node y)
{
    if(!x)
    {
        x=newnode();v[x]=y;up(x);
        return ;
    }
    if(v[x].x[wd]>=y.x[wd]) ins(ls[x],wd^1,y);
    else ins(rs[x],wd^1,y);
    up(x);check(x,wd);
}
void get(int x,node y)
{
    int dl=inf,dr=inf,t=mht(v[x],y);
    ans=min(ans,t);
    if(ls[x]) dl=dis(y,ls[x]);
    if(rs[x]) dr=dis(y,rs[x]);
    if(dl<dr)
    {
        if(dl<ans) get(ls[x],y);
        if(dr<ans) get(rs[x],y);
    }
    else
    {
        if(dr<ans) get(rs[x],y);
        if(dl<ans) get(ls[x],y);
    }
}
signed main()
{
    n=read();m=read();
    for(int i=1;i<=n;i++)
        a[i].x[0]=read(),a[i].x[1]=read();
    build(rt,1,n,0);
    for(int i=1;i<=m;i++)
    {
        int t=read(),x=read(),y=read();
        if(t==1) ins(rt,0,node{x,y});
        else
        {
            ans=inf;
            get(rt,node{x,y});
            printf("%d\n",ans);
        }
    }
}