bzoj2716 2716: [Violet 3]天使玩偶

2716: [Violet 3]天使玩偶

Time Limit: 80 Sec  Memory Limit: 128 MB
Submit: 2943  Solved: 1292
[Submit][Status][Discuss]

Description

Input

Output

同： 2648 ，需要插入的kd树。  

#include<bits/stdc++.h>
using namespace std;  
int const maxn=1000000+100; 
int const inf=1e9;    
int n,m,cur,root;   
struct P
{
    int d[2],mx[2],mn[2],lch,rch;
    int& operator[](int x) {return d[x];}
    friend bool operator<(P x,P y) {return x[cur]<y[cur];}
    friend int dis(P x,P y) {return abs(x[0] -y[0])+abs(x[1]-y[1]);}
}p[maxn]; 
struct kdtree{
  P t[maxn],T;
  int ans;
  void update(int k){
        int l=t[k].lch;  
        int r=t[k].rch;   
        if(l){
            t[k].mn[0]=min(t[k].mn[0],t[l].mn[0]);  
            t[k].mx[0]=max(t[k].mx[0],t[l].mx[0]);  
            t[k].mn[1]=min(t[k].mn[1],t[l].mn[1]);  
            t[k].mx[1]=max(t[k].mx[1],t[l].mx[1]); 
        } 
        if(r){
            t[k].mn[0]=min(t[k].mn[0],t[r].mn[0]); 
            t[k].mx[0]=max(t[k].mx[0],t[r].mx[0]);  
            t[k].mn[1]=min(t[k].mn[1],t[r].mn[1]);  
            t[k].mx[1]=max(t[k].mx[1],t[r].mx[1]);  
        }  
    } 
  int build(int now,int l,int r)
  {
      cur=now;
      int mid=(l+r)/2;
      nth_element(p+l,p+mid,p+r+1);
      t[mid]=p[mid];
      for (int i=0;i<2;i++) t[mid].mn[i]=t[mid].mx[i]=t[mid][i];
      if (l<mid) t[mid].lch=build(now^1,l,mid-1);
      if (r>mid) t[mid].rch=build(now^1,mid+1,r);
      update(mid);
      return mid;
  }
  int getmn(P c){
        int ret=0;  
        ret+=max(0,c.mn[0]-T[0]);  
        ret+=max(0,T[0]-c.mx[0]);
        ret+=max(0,c.mn[1]-T[1]);  
        ret+=max(0,T[1]-c.mx[1]);  
        return ret;  
    }   
    void querymn(int k){
        ans=min(ans,dis(t[k],T)); 
        int l=t[k].lch;  
        int r=t[k].rch; 
        int vl=inf,vr=inf;  
        if(l) vl=getmn(t[l]);  
        if(r) vr=getmn(t[r]); 
        if(vl<vr){
            if(vl<ans) querymn(l);  
            if(vr<ans) querymn(r);  
        }else{
            if(vr<ans) querymn(r); 
            if(vl<ans) querymn(l); 
        } 
    }     

    void insert(int k,int now){
        if(T[now]>=t[k][now]){
            if(t[k].rch) insert(t[k].rch,now^1);  
            else {
                t[k].rch=++n;  t[n]=T;  
                //t[n].ch[0]=t[n].ch[1]=0;  
                t[n].mx[0]=t[n].mn[0]=T[0];  
                t[n].mx[1]=t[n].mn[1]=T[1]; 
            } 
        }else {
            if(t[k].lch) insert(t[k].lch,now^1);  
            else {
                t[k].lch=++n; t[n]=T;  
            //    t[n].ch[0]=t[n].ch[1]=0; 
                t[n].mx[0]=t[n].mn[0]=T[0]; 
                t[n].mx[1]=t[n].mn[1]=T[1];   
            } 
        }
        update(k); 
    }     
    int query(int x,int y)
    {
        ans=inf;
        T[0]=x;T[1]=y;T.lch=0;T.rch=0;
        querymn(root);
        return ans;
    }
    void insert1(int x,int y)
    {
        T[0]=x;T[1]=y;T.lch=0;T.rch=0;insert(root,0);
    }
}kd;    
int main(){   
    scanf("%d%d",&n,&m);  
    for(int i=1;i<=n;i++){
        scanf("%d%d",&p[i][0],&p[i][1]);    
    }  
    root=kd.build(0,1,n); 
    while (m--){
        int t,x,y;  
        scanf("%d%d%d",&t,&x,&y);  
        if(t==1) kd.insert1(x,y);  
        else printf("%d\n",kd.query(x,y)); 
    } 
    return 0; 
} 

 

本题还可以进行cdq分治。  

 

 

#include<bits/stdc++.h>
using namespace std;  
int const N=500000+10;  
int const M=1000000+10;  
int const inf=1e9;  
struct node{
    int x,y,opt,id;  
}q[N<<1],t[N<<1],tmp[N<<1];   
int n,m,ans[N],len,c[M];   
void update(int x,int v){
    for(int i=x;i<=len;i+=(i&-i))  
        c[i]=max(c[i],v);  
}
void clear(int x){
    for(int i=x;i<=len;i+=(i&-i))  
        c[i]=0;  
}
int get(int x){
    int s=0;  
    for(int i=x;i;i-=(i&-i)) 
        s=max(s,c[i]); 
    return s? s:-inf;  
}    
void cdq(int l,int r){
    if(l==r) return;   
    int mid=(l+r)/2;  
    cdq(l,mid);    
    cdq(mid+1,r);   
    int i=l,j=mid+1,k=l;  
    while (i<=mid && j<=r){
        if(t[i].x<=t[j].x){
            if(t[i].opt==1)  
                update(t[i].y,t[i].x+t[i].y); 
            tmp[k++]=t[i++]; 
        }else {
            if(t[j].opt==2)   
            {
                ans[t[j].id]=min(ans[t[j].id],t[j].x+t[j].y-get(t[j].y));  
            }
            tmp[k++]=t[j++];  
        }
    }
    while (i<=mid) {
        if(t[i].opt==1)  update(t[i].y,t[i].x+t[i].y); 
        tmp[k++]=t[i++];  
    }
    while (j<=r) {
        if(t[j].opt==2)  ans[t[j].id]=min(ans[t[j].id],t[j].x+t[j].y-get(t[j].y)); 
        tmp[k++]=t[j++];  
    }
    for(int i=l;i<=mid;i++)  
        if(t[i].opt==1) clear(t[i].y);   
    for(int i=l;i<=r;i++) t[i]=tmp[i];  
}
void solve(int a,int b){
    for(int i=1;i<=n+m;i++){ 
        t[i]=q[i];  
        if(a) t[i].x=len-t[i].x;  
        if(b) t[i].y=len-t[i].y;  
    }
    cdq(1,n+m); 
}
int main(){ 
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        q[i].opt=1;  
        scanf("%d%d",&q[i].x,&q[i].y);   
        q[i].x++; q[i].y++;   
        len=max(q[i].x,len); 
        len=max(q[i].y,len);  
    } 
    for(int i=n+1;i<=n+m;i++){
        scanf("%d%d%d",&q[i].opt,&q[i].x,&q[i].y);  
        q[i].id=i-n;  
        q[i].x++;  
        q[i].y++;  
        len=max(q[i].x,len);  
        len=max(q[i].y,len); 
    }
    for(int i=1;i<=m;i++) ans[i]=inf;  
    len++;  
    solve(0,0);  
    solve(0,1);      
    solve(1,0);  
    solve(1,1);  
    for(int i=n+1;i<=n+m;i++)  
        if(q[i].opt==2)  printf("%d\n",ans[q[i].id]); 
    return 0; 
}