bzoj5028 小Z的加油店 差分+数论+线段树/树状数组

原题: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;
}