bzoj2716-天使玩偶

开始给出平面上(n)个点，(m)个操作，每次加一个点或者查询一个点的曼哈顿距离最近点的曼哈顿距离。（坐标为整数）。

(n,mle 5 imes 10^5)。

分析

首先KD-Tree的做法是错误的，一个经典的例子是给出(n)个点的圆，询问圆心。

如果这个问题是只统计一个角方向（一个点左下方，右下方，右上方，左上方的曼哈顿距离最近点）的，那么可以排序加线段树解决。又因为这个问题满足加法，所以我们对每个方向做一次就好啦。

那动态插入怎么办呢？使用cdq分治，每次统计前一半的点对后一半的询问的贡献，由于满足加法，所以可以这样做。

复杂度为(O(nlog^2n))。

由于是角方向，所以可以发现查询总是查询([1,r])或者([l,max])之间的最小值，所以可以维护两个相反的树状数组来代替常数大的线段树。

代码

没有办法再优化了，还是过不了啊！！精神上过了。

#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#include<sys/mman.h>
#include<sys/stat.h>
#define lowbit(x) (x&-x)
struct BufferedInputStream {
	char *buf,*p;
	int size;
	inline void init() {
		register int fd=fileno(stdin);
		struct stat sb;
		fstat(fd, &sb);
		size=sb.st_size;
		buf=reinterpret_cast<char *>(mmap(0, size, PROT_READ, MAP_PRIVATE, fileno(stdin), 0));
		p=buf;													    
	}
	inline char gchar() {
		return (p==buf+size || *p==-1)?-1:*p++;			    
	}
} str;
using std::sort;
using std::fill;
inline int read() {
	register int x=0,f=1;
	register char c=str.gchar();
	for (;!isdigit(c);c=str.gchar()) if (c=='-') f=-1;
	for (;isdigit(c);c=str.gchar()) x=x*10+c-'0';
	return x*f;
}
inline void write(int &x) {
	if (!x) putchar('0'); else {
		static char s[15];
		register int tot=0;
		while (x) s[++tot]=x%10+'0',x/=10;
		while (tot) putchar(s[tot--]);
	}
	putchar('
');
}
const int maxn=2.1e6+10;
const int maxp=1e6+10;
const int inf=1e8+10;
struct node {
	int x,y;
	node (int x=0,int y=0):x(x),y(y) {}
	inline bool operator == (const node &a) const {return x==a.x && y==a.y;}
	inline bool operator < (const node &a) const {return y==a.y?x<a.x:y<a.y;}
};
struct Q {
	int tp,id;
	node p;
} a[maxn],b[maxn];
inline bool cmp1(const Q &a,const Q &b) {
	return a.p==b.p?a.tp<b.tp:a.p<b.p;
}
inline bool cmp2(const Q &a,const Q &b) {
	return a.p==b.p?a.tp>b.tp:a.p<b.p;
}
inline int min(const int &x,const int &y) {return x<y?x:y;}
inline void Min(int &x,const int y) {
	x=min(x,y);
}
int ans[maxn],n,m;
struct BIT {
	struct array {
		int a[maxp],t[maxp],clc;
		inline int & operator [] (const int &x) {
			if (t[x]!=clc) t[x]=clc,a[x]=inf;
			return a[x];
		}
	} a,b;
	BIT () {clr();}
	inline void clr() {
		++a.clc;
		++b.clc;
	}
	inline void insert(int p,const int d) {
		for (register int i=p;i<maxp;i+=lowbit(i)) Min(a[i],d);
		p=maxp-p;
		for (register int i=p;i<maxp;i+=lowbit(i)) Min(b[i],d); 
	}
	inline int query(int l,const int r) {
		register int ret=inf;
		if (l==1) {
			for (register int i=r;i;i-=lowbit(i)) Min(ret,a[i]);
		} else if (r==maxp-1) {
			l=maxp-l;
			for (register int i=l;i;i-=lowbit(i)) Min(ret,b[i]);
		}
		return ret;
	}
} sgt;
void solve(const int l,const int r) {
	if (l==r) return;
	register int mid=(l+r)>>1,all=0;
	solve(l,mid);
	solve(mid+1,r);
	for (register int i=l;i<=mid;++i) if (!a[i].tp) b[++all]=a[i];
	for (register int i=mid+1;i<=r;++i) if (a[i].tp) b[++all]=a[i];


	sort(b+1,b+all+1,cmp1);
	sgt.clr();
	for (register int i=1,j;i<=all;i=j+1) {
		for (j=i;j<all && b[j+1].p.y==b[i].p.y;++j);
		for (register int k=i;k<=j;++k) if (!b[k].tp) sgt.insert(b[k].p.x,-b[k].p.x-b[k].p.y); else 
		Min(ans[b[k].id],b[k].p.x+b[k].p.y+sgt.query(1,b[k].p.x));
	}

	sgt.clr();
	for (register int i=all,j;i;i=j-1) {
		for (j=i;j>1 && b[j-1].p.y==b[i].p.y;--j);
		for (register int k=j;k<=i;++k) if (!b[k].tp) sgt.insert(b[k].p.x,b[k].p.y-b[k].p.x); else
		Min(ans[b[k].id],b[k].p.x-b[k].p.y+sgt.query(1,b[k].p.x));
	}


	sort(b+1,b+all+1,cmp2);
	sgt.clr();
	for (register int i=1,j;i<=all;i=j+1) {
		for (j=i;j<all && b[j+1].p.y==b[i].p.y;++j);
		for (register int k=j;k>=i;--k) if (!b[k].tp) sgt.insert(b[k].p.x,b[k].p.x-b[k].p.y); else
		Min(ans[b[k].id],b[k].p.y-b[k].p.x+sgt.query(b[k].p.x,maxp-1));
	}

	sgt.clr();
	for (register int i=all,j;i;i=j-1) {
		for (j=i;j>1 && b[j-1].p.y==b[i].p.y;--j);
		for (register int k=i;k>=j;--k) if (!b[k].tp) sgt.insert(b[k].p.x,b[k].p.x+b[k].p.y); else 
		Min(ans[b[k].id],-b[k].p.x-b[k].p.y+sgt.query(b[k].p.x,maxp-1));
	}
}
int main() {
#ifndef ONLINE_JUDGE
	freopen("test.in","r",stdin);
	freopen("my.out","w",stdout);
#endif
	str.init();
	n=read(),m=read();
	fill(ans+1,ans+m+1,inf);
	for (register int i=1;i<=n;++i) {
		int x=read()+1,y=read()+1;
		a[i]=(Q){0,0,node(x,y)};
	}
	for (register int i=1;i<=m;++i) {
		int tp=read()-1,x=read()+1,y=read()+1;
		a[i+n]=(Q){tp,i,node(x,y)};
	}
	solve(1,n+m);
	for (register int *i=ans+1,*ed=ans+m+1;i!=ed;++i) if (*i<=(int)2e6) write(*i);
	return 0;
}