原题:https://www.lydsy.com/JudgeOnline/problem.php?id=5028
题解:要求的其实是[l,r]的gcd。就是辗转相减法。但是修改很麻烦考虑差分,即
。很容易就能看出不管区间加法只会改变两个端点,所以差分的部分可以用线段树维护(区间查询,单点修改),a_l最简单可以用树状数组维护(区间修改,单点查询)
#include<bits/stdc++.h>
#define mid (l+r)/2
#define ls now<<1
#define rs now<<1|1
using namespace std;
const int N=1e5+10;
int n,a[N],m,b[N],tt[N<<2],t[N<<2],bit[N];
int gcd(int x,int y){
if(x==0) return y;
else return gcd(y%x,x);
}
inline int rd(){
int x=0;int f=1;char s=getchar();
while(!isdigit(s)) f=(s=='-'?-1:f),s=getchar();
while(isdigit(s)) x=(x<<1)+(x<<3)+s-'0',s=getchar();
return x*f;
}
inline int lowbit(int x){
return x&(-x);
}
void add(int x,int c){//维护原值
for(x;x<=n;x+=lowbit(x)) bit[x]+=c;
}
int sum(int x){
int ans=0;
for(x;x;x-=lowbit(x)) ans+=bit[x];
return ans;
}
inline void upd(int now){
tt[now]=gcd(tt[ls],tt[rs]);
}
void build(int now,int l,int r){
if(l==r){
tt[now]=b[l];return ;
}
build(ls,l,mid);
build(rs,mid+1,r);
upd(now);
}
int query(int now,int l,int r,int x,int y){//修改差分值
if(x>y || l>r) return 0;
if(x<=l && r<=y){
return tt[now];
}
int la=0,ra=0;
if(x<=mid) la=query(ls,l,mid,x,y);
if(mid<y) ra=query(rs,mid+1,r,x,y);
return gcd(la,ra);
}
void change(int now,int l,int r,int p,int c){
if(r>n || l>r || p>n || p<1) return ;
if(l==p && p==r) {
tt[now]+=c;
return;
}
if(p<=mid) change(ls,l,mid,p,c);
if(mid<p) change(rs,mid+1,r,p,c);
upd(now);
}
int main(){
// freopen("bzoj5028.in","r",stdin);
n=rd();m=rd();
for(int i=1;i<=n;i++) a[i]=rd();
b[1]=0;for(int i=2;i<=n;i++) b[i]=a[i]-a[i-1];
build(1,1,n);
for(int i=1,opt,x,y,c;i<=m;i++){
opt=rd();x=rd();y=rd();
if(x>y) continue;
if(opt==1){
int ans=abs(query(1,1,n,x+1,y));
int p=a[x]+sum(x);
printf("%d\n",gcd(ans,p));
}else {
c=rd();
change(1,1,n,x,c);
if(y+1<=n)change(1,1,n,y+1,-c);
add(x,c);add(y+1,-c);
}
}
return 0;
}