二维线段树（单点更新）

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <vector>
#include <queue>
using namespace std;
const int INF=1<<30;
const int maxn=2000+10;
struct INtervalTree2D{  //2维线段树x轴y轴的下标均从1开始
int Max[maxn][maxn],Min[maxn][maxn],Sum[maxn][maxn],n,m;
int xo,xleaf,x1,y1,x2,y2,x,y,v,vmax,vmin,vsum;
void query1D(int o,int L,int R){
if(y1<=L&&R<=y2){
vmax=max(Max[xo][o],vmax);vmin=min(Min[xo][o],vmin);
vsum+=Sum[xo][o];
}else{
int M=L+(R-L)/2;
if(y1<=M)query1D(o*2,L,M);
if(M<y2)query1D(o*2+1,M+1,R);
}
}
void query2D(int o,int L,int R){
if(x1<=L&&R<=x2){xo=o;query1D(1,1,m);}
else{
int M=L+(R-L)/2;
if(x1<=M)query2D(o*2,L,M);
if(M<x2)query2D(o*2+1,M+1,R);
}
}
void modify1D(int o,int L,int R){
if(L==R){
if(xleaf){Sum[xo][o]=Max[xo][o]=Min[xo][o]=v;return ;}
Max[xo][o]=max(Max[xo*2][o],Max[xo*2+1][o]);
Min[xo][o]=min(Min[xo*2][o],Min[xo*2+1][o]);
Sum[xo][o]=Sum[xo*2][o]+Sum[xo*2+1][o];
}else {
int M=L+(R-L)/2;
if(y<=M){modify1D(o*2,L,M);}
else {modify1D(o*2+1,M+1,R);}
Max[xo][o]=max(Max[xo][o*2],Max[xo][o*2+1]);
Min[xo][o]=min(Min[xo][o*2],Min[xo][o*2+1]);
Sum[xo][o]=Sum[xo][o*2]+Sum[xo][o*2+1];
}
}
void modify2D(int o,int L,int R){
if(L==R){xo=o;xleaf=1;modify1D(1,1,m);}
else{
int M=L+(R-L)/2;
if(x<=M)modify2D(o*2,L,M);
else modify2D(o*2+1,M+1,R);
xo=o;xleaf=0;modify1D(1,1,m);
}
}
void init_size(int h,int w)//设置矩阵的行数和列数
{
n=h;m=w;
}
void modify(int x,int y,int v){  //将(x,y)的值修改为v（初始化线段树也是用该函数，调用m*n次）
this->x=x;this->y=y;this->v=v;
modify2D(1,1,n);
}
void query(int x1,int y1,int x2,int y2){   //询问x1,y1,x2,y2间的最值，调用该函数后，vmax,vmin即为结果
this->x1=x1;this->y1=y1;this->x2=x2;this->y2=y2;
vmax=-INF;vmin=INF;vsum=0;
query2D(1,1,n);
}
}xdtree;
int main() {
//  freopen("in.txt","r",stdin);
int n,m;
int x,y,v;
cin>>n>>m;
xdtree.init_size(n,m);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
scanf("%d",&v);
xdtree.modify(i,j,v);
}
int q,x1,x2,y1,y2;
cin>>q;
char cas[3];
while(q--)
{
scanf("%s",cas);
if(cas[0]=='q')
{
cin>>x1>>y1>>x2>>y2;
xdtree.query(x1,y1,x2,y2);
cout<<xdtree.vmax<<' '<<xdtree.vmin<<' '<<xdtree.vsum<<endl;
}
else
{
cin>>x>>y>>v;
xdtree.modify(x,y,v);
}
}
return 0;
}