貌似想复杂了……
#include <iostream>
#include <cstring>
#include <cstdio>
#define mid ((l+r)>>1)
using namespace std;
const int maxn=;
const int maxm=;
int N,M,Q;
struct Array{
long long num[maxm];
long long *operator [](int x){
return &num[(x-)*M];
}
}a,b; long long ABS(long long x){return x>?x:-x;}
long long GCD(long long x,long long y){
return y?GCD(y,x%y):ABS(x);
} int cnt,ch[maxn][];
long long w[maxn]; void Merge(int &x,int p1,int p2,int l,int r){
if(!x)x=++cnt;
w[x]=GCD(w[p1],w[p2]);
if(l!=r){
Merge(ch[x][],ch[p1][],ch[p2][],l,mid);
Merge(ch[x][],ch[p1][],ch[p2][],mid+,r);
}
} void Build(int &x,int l,int r,long long t[]){
x=++cnt;
if(l==r){
w[x]=t[l];
return;
}
Build(ch[x][],l,mid,t);
Build(ch[x][],mid+,r,t);
w[x]=GCD(w[ch[x][]],w[ch[x][]]);
} int rt[maxm<<];
void Build(int x,int l,int r){
if(l==r){
Build(rt[x],,M,a[l]);
return;
}
Build(x<<,l,mid);
Build(x<<|,mid+,r);
Merge(rt[x],rt[x<<],rt[x<<|],,M);
} int x,y,tx,ty;
int tp,x1,y1,x2,y2;
long long d; void Change(int x,int l,int r,long long d){
if(l==r){w[x]+=d;return;}
if(mid>=ty)Change(ch[x][],l,mid,d);
else Change(ch[x][],mid+,r,d);
w[x]=GCD(w[ch[x][]],w[ch[x][]]);
} void Update(int x,int p1,int p2,int l,int r){
w[x]=GCD(w[p1],w[p2]);
if(l==r)return;
if(mid>=ty)Update(ch[x][],ch[p1][],ch[p2][],l,mid);
else Update(ch[x][],ch[p1][],ch[p2][],mid+,r);
} void Modify(int x,int l,int r,long long d){
if(l==r){
Change(rt[x],,M,d);
return;
}
if(mid>=tx)Modify(x<<,l,mid,d);
else Modify(x<<|,mid+,r,d);
Update(rt[x],rt[x<<],rt[x<<|],,M);
} long long Que(int x,int l,int r){
if(l>=y1&&r<=y2)
return w[x];
long long ret=;
if(mid>=y1)ret=Que(ch[x][],l,mid);
if(mid<y2)ret=GCD(ret,Que(ch[x][],mid+,r));
return ret;
} long long Query(int x,int l,int r){
if(l>=x1&&r<=x2)
return Que(rt[x],,M);
long long ret=;
if(mid>=x1)ret=Query(x<<,l,mid);
if(mid<x2)ret=GCD(ret,Query(x<<|,mid+,r));
return ret;
} int main(){
freopen("chessa.in","r",stdin);
freopen("chessa.out","w",stdout);
scanf("%d%d",&N,&M);
scanf("%d%d",&x,&y);
scanf("%d",&Q);
for(int i=;i<=N;i++)
for(int j=;j<=M;j++)
scanf("%lld",&b[i][j]);
for(int i=;i<=N;i++)
for(int j=;j<=M;j++){
a[i][j]=b[i][j];
if(i<x)a[i][j]-=b[i+][j];
if(i>x)a[i][j]-=b[i-][j];
if(j<y)a[i][j]-=b[i][j+];
if(j>y)a[i][j]-=b[i][j-];
if(i!=x&&j!=y)a[i][j]+=b[i+(i<x?:-)][j+(j<y?:-)];
} Build(,,N);
while(Q--){
scanf("%d",&tp);
scanf("%d%d",&x1,&y1);
scanf("%d%d",&x2,&y2);
if(tp==){
x1=x-x1;x2=x+x2;
y1=y-y1;y2=y+y2;
printf("%lld\n",Query(,,N));
}
if(tp==){
scanf("%lld",&d); {
if(x1<=x&&y1<=y){
tx=x1-;ty=y1-;
if(tx&&ty)Modify(,,N,d);
} if(x1<=x&&y1>y){
tx=x1-;ty=y1;
if(tx)Modify(,,N,-d);
} if(x1>x&&y1<=y){
tx=x1;ty=y1-;
if(ty)Modify(,,N,-d);
} if(x1>x&&y1>y){
tx=x1;ty=y1;
Modify(,,N,d);
}
} {
if(x1<=x&&y2<y){
tx=x1-;ty=y2;
if(tx)Modify(,,N,-d);
} if(x1<=x&&y2>=y){
tx=x1-;ty=y2+;
if(tx&&ty<=M)Modify(,,N,d);
} if(x1>x&&y2<y){
tx=x1;ty=y2;
Modify(,,N,d);
} if(x1>x&&y2>=y){
tx=x1;ty=y2+;
if(ty<=M)Modify(,,N,-d);
}
} {
if(x2<x&&y1<=y){
tx=x2;ty=y1-;
if(ty)Modify(,,N,-d);
} if(x2<x&&y1>y){
tx=x2;ty=y1;
Modify(,,N,d);
} if(x2>=x&&y1<=y){
tx=x2+;ty=y1-;
if(tx<=N&&ty)Modify(,,N,d);
} if(x2>=x&&y1>y){
tx=x2+;ty=y1;
if(tx<=N)Modify(,,N,-d);
}
} {
if(x2<x&&y2<y){
tx=x2;ty=y2;
Modify(,,N,d);
} if(x2<x&&y2>=y){
tx=x2;ty=y2+;
if(ty<=M)Modify(,,N,-d);
} if(x2>=x&&y2<y){
tx=x2+;ty=y2;
if(tx<=N)Modify(,,N,-d);
} if(x2>=x&&y2>=y){
tx=x2+;ty=y2+;
if(tx<=N&&ty<=M)Modify(,,N,d);
}
} if(x1<=x&&x2>=x){
if(y1<=y){
tx=x;ty=y1-;
if(ty)Modify(,,N,-d);
}
if(y1>y){
tx=x;ty=y1;
Modify(,,N,d);
} if(y2<y){
tx=x;ty=y2;
Modify(,,N,d);
}
if(y2>=y){
tx=x;ty=y2+;
if(ty<=M)Modify(,,N,-d);
}
} if(y1<=y&&y2>=y){
if(x1<=x){
tx=x1-;ty=y;
if(tx)Modify(,,N,-d);
}
if(x1>x){
tx=x1;ty=y;
Modify(,,N,d);
} if(x2<x){
tx=x2;ty=y;
Modify(,,N,d);
}
if(x2>=x){
tx=x2+;ty=y;
if(tx<=N)Modify(,,N,-d);
}
} if(x1<=x&&x2>=x&&y1<=y&&y2>=y){
tx=x;ty=y;
Modify(,,N,d);
}
}
}
return ;
}