bzoj 5028 线段树&裴蜀定理&差分

根据倒水问题易知，区间[l,r]的gcd就是区间[l,r]的解。且有 gcd(a,b,c)=gcd(a,a-b,b-c) ，故易知，可以将数列转换成差分数组，这样区间加法[l,r]可以变成对l、r+1两点的单点修改add(l,num),add(r+1,-num)，不然就要对区间的所有值全部更新，会退化成最差的暴力更新了。区间查询的时候，又由差分数组可知，gcd[l,r]=gcd(1 ~ l的前缀和，l+1 ~ r的gcd）。因为由上述等式可知，gcd[l,r]的l是没有进行过差分的，也就是原来的数，而l+1 ~ r是进行过差分的，那么为了保证答案的正确性，要求1~l的前缀和得到原来l上的数，然后求的l+1 ~ r的gcd，两者再次gcd即可。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<algorithm>
#include<sstream>
#include<bitset>
#define scand(a) scanf("%d",&a)
#define scandd(a,b) scanf("%d%d",&a,&b)
#define scanddd(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define mst(a,b) memset(a,b,sizeof(a))
#define lowbit(x) x&-x
#define lson rt<<1
#define rson rt<<1|1
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
const int inf=0x3f3f3f3f;
const int maxn=100005;
const int mod=1e9+7;
struct nod
{
	int l,r,x,d;
}tre[maxn<<2];
int a[maxn],cf[maxn],b,n,m;
void pushup(int rt)
{
	tre[rt].x=tre[lson].x+tre[rson].x;
	tre[rt].d=__gcd(tre[lson].d,tre[rson].d);
}
void build(int rt,int l,int r)
{
	tre[rt].l=l,tre[rt].r=r;
	if(l==r)
	{
		tre[rt].x=tre[rt].d=cf[l];
		return;
	}
	int mid=l+r>>1;
	build(lson,l,mid);
	build(rson,mid+1,r);
	pushup(rt);
}
void upd(int pos,int rt,int x)
{
	if(tre[rt].l==tre[rt].r)
	{
		tre[rt].x+=x;
		tre[rt].d+=x;
		return;
	}
	int mid=tre[rt].l+tre[rt].r>>1;
	if(pos<=mid)upd(pos,lson,x);
	else upd(pos,rson,x);
	pushup(rt);
}
int que(int rt,int l,int r)
{
	if(l<=tre[rt].l&&tre[rt].r<=r)return tre[rt].d;
	int mid=tre[rt].l+tre[rt].r>>1;
	int res=0;
	if(l<=mid)res=__gcd(res,que(lson,l,r));
	if(r>mid)res=__gcd(res,que(rson,l,r));
	return res;
}
int que1(int rt,int l,int r)
{
	if(l<=tre[rt].l&&tre[rt].r<=r)return tre[rt].x;
	int mid=tre[rt].l+tre[rt].r>>1;
	int res=0;
	if(l<=mid)res+=que1(lson,l,r);
	if(r>mid)res+=que1(rson,l,r);
	return res;
}
int main()
{
	#ifdef local
	freopen("1.txt","r",stdin);
	#endif
	scandd(n,m);
	int x,y,op;
	for(int i=1;i<=n;i++)scand(a[i]),cf[i]=a[i]-a[i-1];
	build(1,1,n);
	for(int i=1;i<=m;i++)
	{
		scanddd(op,x,y);
		if(x>y)swap(x,y);
		if(op==1)
		{
			int sum1=que1(1,1,x);
			int sum2=que(1,x+1,y);
			printf("%d\n",abs(__gcd(sum1,sum2)));
		}
		else
		{
			scand(op);
			upd(x,1,op);
			upd(y+1,1,-op);
		}
	}
}