HDU 4819 Mosaic --二维线段树(树套树)
题意： 给一个矩阵，每次查询一个子矩阵内的最大最小值，然后更新子矩阵中心点为(Max+Min)/2.
解法： 由于是矩阵，且要求区间最大最小和更新单点，很容易想到二维的线段树，可是因为之前没写过二维的线段树，所以没跳出来。后来熟悉了一下，原来很多细节地方都没有考虑到。
这里build,update,query都分为两个函数，第一个为Y轴的(sub_update)，第二个为X轴的(update)，不仅每个sub_update或sub_build后面要加Y轴的pushup函数，而且每个update或build函数最后也要进行pushup，只不过是X轴的pushup，写成 sub_build(rt,-1,1,n,1) 或 sub_update(rt,-1,1,n,1);
二维线段树就是这里比较费解一点，我感觉， 其他地方还比较好理解。
代码：
#include <iostream> #include <cstdio> #include <cstring> #include <cstdlib> #include <cmath> #include <algorithm> #define Mod 1000000007 #define lll __int64 using namespace std; #define N 100007 struct sub_seg{ lll maxi,mini; }; struct seg{ sub_seg st[810<<2]; }tt[810<<2]; int n; lll Max,Min; void pushup(int xrt,int rt) { tt[xrt].st[rt].maxi = max(tt[xrt].st[rt<<1].maxi,tt[xrt].st[rt<<1|1].maxi); tt[xrt].st[rt].mini = min(tt[xrt].st[rt<<1].mini,tt[xrt].st[rt<<1|1].mini); } void sub_build(int xrt,int posY,int l,int r,int rt) { if(l == r) { if(posY != -1) { scanf("%I64d",&tt[xrt].st[rt].maxi); tt[xrt].st[rt].mini = tt[xrt].st[rt].maxi; } else { tt[xrt].st[rt].maxi = max(tt[xrt<<1].st[rt].maxi,tt[xrt<<1|1].st[rt].maxi); tt[xrt].st[rt].mini = min(tt[xrt<<1].st[rt].mini,tt[xrt<<1|1].st[rt].mini); } return; } int mid = (l + r) >> 1; sub_build(xrt,posY, l,mid, 2*rt); sub_build(xrt,posY, mid+1,r, 2*rt+1); pushup(xrt,rt); //pushup Y } void build(int l, int r, int rt) { if(l == r) { sub_build(rt,l,1,n,1); return; } int mid = (l + r) >> 1; build(l,mid,2*rt); build(mid+1,r,2*rt+1); sub_build(rt,-1,1,n,1); //pushup X } void sub_update(int xrt,int pos,int y,int l,int r,int val,int rt) { if(l == r) { if(pos!=-1) tt[xrt].st[rt].maxi = tt[xrt].st[rt].mini = val; else { tt[xrt].st[rt].maxi = max(tt[xrt<<1].st[rt].maxi,tt[xrt<<1|1].st[rt].maxi); tt[xrt].st[rt].mini = min(tt[xrt<<1].st[rt].mini,tt[xrt<<1|1].st[rt].mini); } return; } int mid = (l+r) >> 1 ; if(y <= mid) sub_update(xrt, pos, y, l,mid, val, rt<<1) ; else sub_update(xrt, pos, y, mid+1,r, val, rt<<1|1) ; pushup(xrt,rt); //pushup Y维 } void update(int x,int y,int l,int r,int val,int rt) //单点更新 { if(l == r) { sub_update(rt,l,y,1,n,val,1); return; } int mid = (l+r) >> 1; if(x<=mid) update(x,y, l,mid, val, rt<<1); else update(x,y, mid+1,r, val, rt<<1|1); sub_update(rt,-1,y,1,n,val,1); //pushup X维 } void sub_query(int xrt,int aa,int bb,int l,int r,int rt) { if(aa <= l && bb >= r) { Max = max(Max,tt[xrt].st[rt].maxi); Min = min(Min,tt[xrt].st[rt].mini); return; } int mid = (l+r)>>1; if(aa<=mid) sub_query(xrt,aa,bb,l,mid, rt<<1); if(bb>mid) sub_query(xrt,aa,bb,mid+1,r, rt<<1|1); } void query(int l,int r,int aa,int bb,int L1,int R1,int rt) { if(aa <= l && bb >= r) { sub_query(rt,L1,R1,1,n,1); return; } int mid = (l+r)>>1; if(aa <= mid) query(l,mid,aa,bb,L1,R1, rt<<1); if(bb > mid) query(mid+1,r,aa,bb,L1,R1, rt<<1|1); } int main() { int t,i,j,m,x,y,L,cs = 1; scanf("%d",&t); while(t--) { scanf("%d",&n); build(1,n,1); scanf("%d",&m); printf("Case #%d:\n",cs++); while(m--) { scanf("%d%d%d",&x,&y,&L); int Lx = max(1,x-L/2); int Ly = max(1,y-L/2); int Rx = min(n,x+L/2); int Ry = min(y+L/2,n); Max = -1, Min = Mod; query(1,n,Lx,Rx,Ly,Ry,1); lll mid = (Max+Min)/2LL; printf("%I64d\n",mid); update(x,y,1,n,mid,1); } } return 0; }