BZOJ2716: [Violet 3]天使玩偶

题解：

kdtree裸题？

感觉判断是否进入左右子树真的是个估价函数。。。

inline int get(int k)
 {
     if(!k)return inf;
    int ret=0;
    for0(i,1)ret+=max(0,t[k].mi[i]-now[i]);
    for0(i,1)ret+=max(0,now[i]-t[k].mx[i]);
    return ret;
 }

我们来看这个函数

按我的理解，这就是在估计以k为根的子树中与now的距离最小是多少

然后就根据这个看看需不需要进入左右子树。。。

好没理由啊。。。怎么保证复杂度的。。。

代码：

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<string>
#define inf 1000000000
#define maxn 500000+5
#define maxm 500+100
#define eps 1e-10
#define ll long long
#define pa pair<int,int>
#define for0(i,n) for(int i=0;i<=(n);i++)
#define for1(i,n) for(int i=1;i<=(n);i++)
#define for2(i,x,y) for(int i=(x);i<=(y);i++)
#define for3(i,x,y) for(int i=(x);i>=(y);i--)
#define for4(i,x) for(int i=head[x],y;i;i=e[i].next)
#define mod 1000000007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
    return x*f;
}
int n,m,cur,tot,rt,ans;
struct rec
{
    int d[2],mi[2],mx[2],l,r;
    int& operator[](int i){return d[i];}
}p[maxn],now,t[maxn<<1];
bool operator<(rec a,rec b){return a[cur]<b[cur];}
inline int dis(rec a,rec b){return abs(a[0]-b[0])+abs(a[1]-b[1]);}
inline void pushup(int k)
 {
    int l=t[k].l,r=t[k].r;
    for0(i,1)
    {
        t[k].mi[i]=min(t[k].mi[i],min(t[l].mi[i],t[r].mi[i]));
        t[k].mx[i]=max(t[k].mx[i],max(t[l].mx[i],t[r].mx[i]));
    }
 }
int build(int l,int r,int dir)
 {
    cur=dir;int mid=(l+r)>>1;
    nth_element(p+l,p+mid,p+r+1);
    t[mid]=p[mid];
    for0(i,1)t[mid].mi[i]=t[mid].mx[i]=t[mid][i];
    if(l<mid)t[mid].l=build(l,mid-1,dir^1);
    if(mid<r)t[mid].r=build(mid+1,r,dir^1);
    pushup(mid);
    return mid;
 }
void insert(int &k,int dir)
 {
    if(!k)
    {
      k=++tot;t[tot]=now;
      for0(i,1)t[tot].mi[i]=t[tot].mx[i]=t[tot][i];
      return;
    }
    if(now[dir]<t[k][dir])insert(t[k].l,dir^1);
    else insert(t[k].r,dir^1);
    pushup(k);
 }
inline int get(int k)
 {
     if(!k)return inf;
    int ret=0;
    for0(i,1)ret+=max(0,t[k].mi[i]-now[i]);
    for0(i,1)ret+=max(0,now[i]-t[k].mx[i]);
    return ret;
 }
void query(int k,int dir)
 {
    int dl=get(t[k].l),dr=get(t[k].r),d0=dis(t[k],now);
    if(d0<ans)ans=d0;
    if(dl<dr)
    {
        if(dl<ans)query(t[k].l,dir^1);
        if(dr<ans)query(t[k].r,dir^1);
    }else
    {
        if(dr<ans)query(t[k].r,dir^1);
        if(dl<ans)query(t[k].l,dir^1);
    }
 }
int a[20];
inline void print(int x)
{
    if(!x)putchar('0');
    for(a[0]=0;x;x/=10)a[++a[0]]=x%10;
    for3(i,a[0],1)putchar('0'+a[i]);
    putchar('
');
}   
int main()
{
    freopen("input.txt","r",stdin);
    freopen("output.txt","w",stdout);
    n=read();m=read();tot=n;
    for0(i,1)t[0].mi[i]=inf,t[0].mx[i]=-inf;
    for1(i,n)p[i][0]=read(),p[i][1]=read();
    rt=build(1,n,0);
    while(m--)
    {
        int op=read();now[0]=read();now[1]=read();ans=inf;
        if(op==1)insert(rt,0);else query(rt,0),print(ans);
    }
    return 0;
}

加了输出优化还是没进第一版，好桑心QAQ

2716: [Violet 3]天使玩偶

Time Limit: 80 Sec  Memory Limit: 128 MB
Submit: 569  Solved: 242
[Submit][Status]

Description

Input

Output