CDQ分治,归并排序写错调了一天可海星
考虑哈夫曼距离特点:
$$dis(A,B)=|A_x-B_x|+|A_y-B_y|$$
若能够把绝对值去掉岂不是很妙!
$$dis(A,B)=(A_x+A_y)-(B_x+B_y)$$
由于A坐标已经定下来了,我们只要找B坐标使得$B_x+B_y$ 最大,就可以是距离最小了!
然后套上CDQ三维偏序求解就好了
因为不一定每个点都在当前询问点的左下方(也就是说绝对值不能直接去掉,我们把所有点转个方向再做3遍就吼了)
Tip:这道题卡常,请使用归并排序,并且求解区间最值时不要使用线段树,改用树状数组!!
(虽然某人说正常人用的都是树状数组,可我开始就是用线段树了..)
1 #include<bits/stdc++.h> 2 #define lowbit(x) ((x)&(-x)) 3 #define writeln(x) write(x),puts("") 4 #define writep(x) write(x),putchar(' ') 5 using namespace std; 6 inline int read(){ 7 int ans=0,f=1;char chr=getchar(); 8 while(!isdigit(chr)){if(chr=='-') f=-1;chr=getchar();} 9 while(isdigit(chr)){ans=(ans<<3)+(ans<<1)+chr-48;chr=getchar();} 10 return ans*f; 11 }void write(int x){ 12 if(x<0) putchar('-'),x=-x; 13 if(x>9) write(x/10); 14 putchar(x%10+'0'); 15 }const int M = 1e6+5; 16 int n,m,tot,s[M],maxx,maxy,maxn,ans[M]; 17 struct P{int x,y,id,type;}q[M],p[M],h[M]; 18 inline void cmax(int &x,int y){if(x<y)x=y;} 19 inline void cmin(int &x,int y){if(x>y)x=y;} 20 inline void Add(int x,int y) {for(;x<=maxn;x+=lowbit(x))cmax(s[x],y);} 21 inline void Clear(int x) {for(;x<=maxn;x+=lowbit(x))s[x]=0;} 22 inline int Max(int x){int ans=0;for(;x;x-=lowbit(x))cmax(ans,s[x]);return ans;} 23 void Solve(int l,int r){ 24 if(l==r) return; 25 int mid=l+r>>1; 26 Solve(l,mid),Solve(mid+1,r); 27 int j=l,ll=l,rr=mid+1,nn=l-1; 28 for(int i=mid+1;i<=r;i++){ 29 while(i<=r&&p[i].type!=2) ++i;if(i>r) break; 30 for(;j<=mid&&p[j].x<=p[i].x;++j)if(p[j].type==1)Add(p[j].y,p[j].y+p[j].x); 31 int t=Max(p[i].y); 32 if(t) cmin(ans[p[i].id],p[i].x+p[i].y-t); 33 } 34 for(int i=l;i<j;i++)if(p[i].type==1) Clear(p[i].y); 35 while(ll<=mid&&rr<=r){ 36 if(p[ll].x<=p[rr].x) h[++nn]=p[ll++]; 37 else h[++nn]=p[rr++]; 38 } 39 while(ll<=mid) h[++nn]=p[ll++]; 40 while(rr<=r) h[++nn]=p[rr++]; 41 for(int i=l;i<=r;i++) p[i]=h[i]; 42 } 43 int main(){ 44 n=read(),m=read(); 45 for(int x,y,i=1;i<=n;i++)x=read(),y=read(),q[++tot]=(P){x+1,y+1,i,1},cmax(maxx,x+1),cmax(maxy,y+1); 46 for(int x,y,opt,i=1;i<=m;i++){ 47 opt=read(),x=read(),y=read(); 48 cmax(maxx,x+1),cmax(maxy,y+1); 49 if(opt==1) q[++tot]=(P){x+1,y+1,i+n,1}; 50 else q[++tot]=(P){x+1,y+1,i+n,2}; 51 }n+=m;memset(ans,0x3f,sizeof(ans)); 52 maxn=max(maxx,maxy)+1; 53 for(int i=1;i<=n;i++) p[i]=q[i]; 54 Solve(1,n); 55 for(int i=1;i<=n;i++) p[i]=q[i],p[i].x=maxn-p[i].x; 56 Solve(1,n); 57 for(int i=1;i<=n;i++) p[i]=q[i],p[i].y=maxn-p[i].y; 58 Solve(1,n); 59 for(int i=1;i<=n;i++) p[i]=q[i],p[i].x=maxn-p[i].x,p[i].y=maxn-p[i].y; 60 Solve(1,n); 61 for(int i=1;i<=n;i++)if(q[i].type==2) printf("%d\n",ans[i]); 62 return 0; 63 }