acdream1197 Points In Cuboid（hash树状数组）

题目链接：http://acdream.info/problem?pid=1197

题意：给出三维空间n个点，m个查询，每次查询某个立方体内的点的个数。

思路：按照一维排序，根据查询插入，其他两位用二位树状数组维护。由于这个坐标太大，二位数组开不出来。这时候就是hash，对于一个位置(x,y)，哈希成一个数，作为下标。查询的时候不存在的数字为0.

const int mod=4000007;

const int INF=1000000005;

const int N=100005;

const int M=20011;



struct node

{

	int x,y,z;

};



node a[N],b[N][2];

int n,m;



int ans[N];



pii Q[N];



int cmp(node a,node b)

{

	return a.z<b.z;

}



int A[mod];

int arr[mod];



inline int ha(int x)

{

	int k=x%mod;

	int i;

	for(i=k;;i++)

	{

		if(i==mod) i=0;

		if(A[i])

		{

			if(A[i]==x) return i;

		}

		else

		{

			A[i]=x;

			return i;

		}

	}

}



inline int find(int x)

{

	int k=x%mod;

	int i;

	for(i=k;;i++)

	{

		if(i==mod) i=0;

		if(A[i])

		{

			if(A[i]==x) return i;

		}

		else

		{

			return i;

		}

	}

	return 0;

}



int que(int x,int y)

{

	if(!x||!y) return 0;

	int ans=0;

	while(x)

	{

		int i;

		for(i=y;i;i-=i&-i) 

		{

			ans+=arr[find(x*M+i)];

		}

		x-=x&-x;

	}

	return ans;

}



void add(int x,int y)

{

	while(x<M)

	{

		int i=y;

		while(i<M) arr[ha(x*M+i)]++,i+=i&-i;

		x+=x&-x;

	}

}



int cal(int t)

{

	int ans=0;

	ans+=que(b[t][1].x,b[t][1].y);

	ans-=que(b[t][1].x,b[t][0].y-1);

	ans-=que(b[t][0].x-1,b[t][1].y);

	ans+=que(b[t][0].x-1,b[t][0].y-1);

	return ans;

}



void cal()

{

	clr(A,0);

	clr(arr,0);

	int i;

	for(i=1;i<=n;i++)

	{

		scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z);

		a[i].x+=10005;

		a[i].y+=10005;

	}

	scanf("%d",&m);

	int j;

	for(i=1;i<=m;i++)

	{

		for(j=0;j<2;j++)

		{

			scanf("%d%d%d",&b[i][j].x,&b[i][j].y,&b[i][j].z);

			b[i][j].x+=10005;

			b[i][j].y+=10005;

		}

		if(b[i][0].x>b[i][1].x) swap(b[i][0].x,b[i][1].x);

		if(b[i][0].y>b[i][1].y) swap(b[i][0].y,b[i][1].y);

		if(b[i][0].z>b[i][1].z) swap(b[i][0].z,b[i][1].z);



		Q[i*2-1]=MP(b[i][0].z-1,-i);

		Q[i*2]=MP(b[i][1].z,i);

	}

	sort(a+1,a+n+1,cmp);

	sort(Q+1,Q+2*m+1);

	j=1;

	for(i=1;i<=m;i++) ans[i]=0;

	for(i=1;i<=n;i++)

	{

		while(j<=2*m&&Q[j].first<a[i].z)

		{

			int x=Q[j].second;

			if(x>0) ans[x]+=cal(x);

			else ans[-x]-=cal(-x);

			j++;

		}

		add(a[i].x,a[i].y);

	}		

	while(j<=2*m)

	{

		int x=Q[j].second;

		if(x>0) ans[x]+=cal(x);

		else ans[-x]-=cal(-x);

		j++;

	}

	for(i=1;i<=m;i++) printf("%d\n",ans[i]);

}



int main()

{

	int num=0;

	while(scanf("%d",&n)!=-1)

	{

		printf("Case #%d:\n",++num);

		cal();

	}

}