[ZOJ-2859]（二维线段树）

头有点晕，就不写解释了，随缘看。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
#define mem(a,x) memset(a,x,sizeof(a));
#define s1(x) scanf("%d",&x);
#define s2(x,y) scanf("%d%d",&x,&y);
#define s3(x,y,z) scanf("%d%d%d",&x,&y,&z);
#define s4(x,y,z,k) scanf("%d%d%d%d",&x,&y,&z,&k);
#define ls 2*rt
#define rs 2*rt+1
#define lson ls,L,mid
#define rson rs,mid+1,R
#define ll long long
using namespace std;
typedef pair<int,int> pii;
const ll inf = 0x3f3f3f3f;
/*void dis(int a[], int n){
	printf("总数为%d个\n",n); 
	for(int i = 0; i < n; i++) 	cout<<a[i]<<", ";
	cout<<endl<<"------------------"<<endl;		
}*/

const int mx = 0x3f3f3f3f;
int no[1300][1300]; 
int n,m;
void buildy(int x,int rt, int L, int R){

	no[x][rt] = mx;    //放错位置 
	if(L == R){
		return;
	}
	
	int mid = (L+R)/2;
	buildy(x,lson);
	buildy(x,rson);
	
}
void buildx(int rt, int L, int R){
	buildy(rt,1,1,n);
	if(L == R)
		return;
	int mid = (L+R)/2;
	buildx(lson);
	buildx(rson);
}
void inserty(int x,int rt, int L, int R, int a, int d){
	if(L == R){
		no[x][rt] = min(no[x][rt], d);
		return; 
	} 
	int mid = (L + R) / 2;
	if(a <= mid)
		inserty(x,lson,a,d);
	else
		inserty(x,rson,a,d);
		
	no[x][rt] = min(no[x][ls],no[x][rs]);   //这句不能少，与其叫做插入不如叫做更新  	
}
void insertx(int rt, int L, int R, int h, int a, int d){
	inserty(rt,1,1,n,a,d);
	if(L == R)
		return ;
	int mid = (L + R) / 2;
	if(h <= mid)
		insertx(lson,h,a,d);
	else
		insertx(rson,h,a,d);
}

int quey(int x, int rt, int L, int R, int a, int a1){
	if(a <= L&& R<=a1)                 //写反了 
		return no[x][rt];
	int t1=mx,t2 = mx;	
	int mid = (L + R) / 2;
	if(a<=mid) t1 = quey(x,lson,a,a1);
	if(a1 >mid) t2 = quey(x,rson,a,a1);

	return min(t1,t2); 
	
}
int quex(int rt, int L, int R, int h, int h1, int a, int a1){
	
	if(h <= L&& R<=h1)
		return quey(rt,1,1,n,a,a1);
	int t1=mx,t2 = mx;	
	int mid = (L + R) / 2;
	if(h<=mid) t1 = quex(lson,h,h1,a,a1);
	if(h1 >mid) t2 = quex(rson,h,h1,a,a1);
	
	return min(t1,t2); 	
} 
int main(){
	int T=10,x1,x2,y1,y2,te;
	scanf("%d",&T);
	
	while(T--){
		s1(n);
		buildx(1,1,n);
		for(int i = 1; i <= n; i++)
			for(int j = 1; j<= n; j++){
				s1(te);
				insertx(1,1,n,i,j,te);
			}
		s1(m);
		for(int i = 0; i <m ;i++){
			s4(x1,y1,x2,y2);
			printf("%d\n",quex(1,1,n,x1,x2,y1,y2));
		}	

	}
	
	return 0;
}