这一道题的话就我而言可能更像数学里的分类讨论？
就是讨论一下修改矩形的四个角在守护者的哪边就好了，只是我的代码可能有点麻烦，如果有神犇看到的话能给点建议吗，拜谢！
然后很多东西我会放在代码里面，大佬们就看代码吧
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <vector>
#include <cstring>
#include <queue>
#include <complex>
#include <stack>
#define LL long long int
#define dob double
#define FILE "chessa"
//#define FILE "魔幻棋盘"
using namespace std;
const int N = 500010;
int n,m,X,Y,q,rt,tot;
struct Array{//来自BOBOYANG的黑科技
LL t[N];
LL *operator [](int x){return &t[(x-1)*m];}
}A;
struct Tree{int l,r,ls,rs,tree;LL val;}T[N<<4];
inline LL gcd(LL a,LL b){return b?gcd(b,a%b):a;}
inline LL gi(){
LL x=0,res=1;char ch=getchar();
while(ch>'9'||ch<'0'){if(ch=='-')res*=-1;ch=getchar();}
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*res;
}
namespace Seg{
inline void build(int &x,int l,int r,int id){
x=++tot,T[x].l=l,T[x].r=r;int mid=(l+r)>>1;
if(l==r){T[x].val=A[id][l];return;}
build(T[x].ls,l,mid,id);build(T[x].rs,mid+1,r,id);
T[x].val=gcd(T[T[x].ls].val,T[T[x].rs].val);
}
inline LL query(int x,int xl,int xr){
int l=T[x].l,r=T[x].r,mid=(l+r)>>1;
if(l==xl && r==xr)return T[x].val;
if(xr<=mid)return query(T[x].ls,xl,xr);
if(xl>mid)return query(T[x].rs,xl,xr);
return gcd(query(T[x].ls,xl,mid),query(T[x].rs,mid+1,xr));
}
inline void update(int x,int xp,LL val){
int l=T[x].l,r=T[x].r,mid=(l+r)>>1;
if(l==r){T[x].val+=val;return;}
if(xp<=mid)update(T[x].ls,xp,val);
else update(T[x].rs,xp,val);
T[x].val=gcd(T[T[x].ls].val,T[T[x].rs].val);
}
}
namespace TwoD_Segment_Tree{
inline void merge(int &x,int r1,int r2,int s,int t){
if(!x)x=++tot,T[x].l=s,T[x].r=t;int mid=(s+t)>>1;
if(s==t){T[x].val=gcd(T[r1].val,T[r2].val);return;}
merge(T[x].ls,T[r1].ls,T[r2].ls,s,mid);
merge(T[x].rs,T[r1].rs,T[r2].rs,mid+1,t);
T[x].val=gcd(T[T[x].ls].val,T[T[x].rs].val);
}
inline void build(int &x,int s,int t){
x=++tot;T[x].l=s;T[x].r=t;int mid=(s+t)>>1;
if(s==t){Seg::build(T[x].tree,1,m,s);return;}
build(T[x].ls,s,mid);build(T[x].rs,mid+1,t);
merge(T[x].tree,T[T[x].ls].tree,T[T[x].rs].tree,1,m);
}
inline LL query(int x,int x1,int x2,int y1,int y2){
int l=T[x].l,r=T[x].r,mid=(l+r)>>1;
if(l==x1 && r==x2)return Seg::query(T[x].tree,y1,y2);
if(x2<=mid)return query(T[x].ls,x1,x2,y1,y2);
if(x1>mid)return query(T[x].rs,x1,x2,y1,y2);
return gcd(query(T[x].ls,x1,mid,y1,y2),query(T[x].rs,mid+1,x2,y1,y2));
}
inline void update(int x,int xp,int y,LL val){
int l=T[x].l,r=T[x].r,mid=(l+r)>>1;
if(l==r){Seg::update(T[x].tree,y,val);return;}
if(xp<=mid)update(T[x].ls,xp,y,val);
else update(T[x].rs,xp,y,val);
LL lval=Seg::query(T[T[x].ls].tree,y,y);
LL rval=Seg::query(T[T[x].rs].tree,y,y);
LL mval=Seg::query(T[x].tree,y,y);
Seg::update(T[x].tree,y,gcd(lval,rval)-mval);
}
}
using namespace TwoD_Segment_Tree;
int main()
{
n=gi();m=gi();X=gi();Y=gi();q=gi();
for(int i=1;i<=n;++i)
for(int j=1;j<=m;++j)
A[i][j]=gi();
for(int i=1;i<=n;++i){
for(int j=1;j<Y;++j)A[i][j]-=A[i][j+1];
for(int j=m;j>Y;--j)A[i][j]-=A[i][j-1];
}
for(int j=1;j<=m;++j){
for(int i=1;i<X;++i)A[i][j]-=A[i+1][j];
for(int i=n;i>X;--i)A[i][j]-=A[i-1][j];
}
build(rt,1,n);
for(int t=1,type,x1,y1,x2,y2;t<=q;++t){
scanf("%d%d%d%d%d",&type,&x1,&y1,&x2,&y2);
if(type==0){
x1=X-x1;y1=Y-y1;x2=X+x2;y2=Y+y2;
printf("%lld
",abs(query(rt,x1,x2,y1,y2)));
}
if(type==1){
LL val=gi();
//讨论修改矩形的四个角在守护者的哪边。
//左上角(x1,y1)
/*在守护者左上*/if(x1<=X && y1<=Y && x1>1 && y1>1)update(rt,x1-1,y1-1,val);
/*在守护者右上*/if(x1<=X && y1>Y && x1>1)update(rt,x1-1,y1,-val);
/*在守护者左下*/if(x1>X && y1<=Y && y1>1)update(rt,x1,y1-1,-val);
/*在守护者右下*/if(x1>X && y1>Y)update(rt,x1,y1,val);
//右上角(x1,y2)
/*在守护者左上*/if(x1<=X && y2<Y && x1>1)update(rt,x1-1,y2,-val);
/*在守护者右上*/if(x1<=X && y2>=Y && x1>1 && y2<m)update(rt,x1-1,y2+1,val);
/*在守护者左下*/if(x1>X && y2<Y)update(rt,x1,y2,val);
/*在守护者右下*/if(x1>X && y2>=Y && y2<m)update(rt,x1,y2+1,-val);
//左下角(x2,y1)
/*在守护者左上*/if(x2<X && y1<=Y && y1>1)update(rt,x2,y1-1,-val);
/*在守护者右上*/if(x2<X && y1>Y)update(rt,x2,y1,val);
/*在守护者左下*/if(x2>=X && y1<=Y && x2<n && y1>1)update(rt,x2+1,y1-1,val);
/*在守护者右下*/if(x2>=X && y1>Y && x2<n)update(rt,x2+1,y1,-val);
//右下角(x2,y2)
/*在守护者左上*/if(x2<X && y2<Y)update(rt,x2,y2,val);
/*在守护者右上*/if(x2<X && y2>=Y && y2<m)update(rt,x2,y2+1,-val);
/*在守护者左下*/if(x2>=X && y2<Y && x2<n)update(rt,x2+1,y2,-val);
/*在守护者右下*/if(x2>=X && y2>=Y && x2<n && y2<m)update(rt,x2+1,y2+1,val);
//讨论修改矩形的四条边跨守护者十字情况。
//矩形竖边对守护者横线的影响
if(x1<=X && x2>=X){
//左边(x1,y1) -> (x2,y1)
/*跨左横线*/if(y1<=Y && y1>1)update(rt,X,y1-1,-val);
/*跨右横线*/if(y1>Y)update(rt,X,y1,val);
//右边(x1,y2) -> (x2,y2)
/*跨左横线*/if(y2<Y)update(rt,X,y2,val);
/*跨右横线*/if(y2>=Y && y2<m)update(rt,X,y2+1,-val);
}
//矩形横边对守护者竖线的影响
if(y1<=Y && y2>=Y){
//上边(x1,y1) -> (x1,y2)
/*跨上竖线*/if(x1<=X && x1>1)update(rt,x1-1,Y,-val);
/*跨下竖线*/if(x1>X)update(rt,x1,Y,val);
//下边(x2,y1) -> (x2,y2)
/*跨上竖线*/if(x2<X)update(rt,x2,Y,val);
/*跨下竖线*/if(x2>=X && x2<n)update(rt,x2+1,Y,-val);
}
if(x1<=X && x2>=X && y1<=Y && y2>=Y)update(rt,X,Y,val);
}
}
return 0;
}