Interval GCD CH4302

##题目大意：
给定一个长度为N的数列A，以及M条指令 (N≤5*10^5, M<=10^5)，每条指令可能是以下两种之一：
“C l r d”，表示把 A[l],A[l+1],…,A[r] 都加上 d。
“Q l r”，表示询问 A[l],A[l+1],…,A[r] 的最大公约数(GCD)

##解题思路：
线段树+GCD = 正解
GCD 更相减损术 | (扩展)欧几里得算法
都可以
欧几里得算法能过我就用了欧几里得求gcd
t[now].g=gcd(t[leftson].g,t[rightson].g);
询问的答案就是
gcd(a[l],Ask(1,l+1,r));

我发现很多人的写法和我不同，于是我用了一下别人的写法。接着我心里默默地想：再也不用别人的写法了，本来想改，但是我懒我没多少时间去改，反正能AC的就是好代码

##Accepted code：

#include<cstdio>
#include<iostream>
#include<algorithm>
#define ls (now<<1)
#define rs ((now<<1)+1)
#define ll long long
using namespace std;
const int N=5e5+5;
const int M=1e5+5;
ll a[N],b[N],num[N<<2],s[N<<2];
int l[N<<2],r[N<<2],lc[N<<2],rc[N<<2],n,m;

//以下装逼部分 
inline void read(ll &f) {
	char c=getchar();f=0;int ag=1;
	while(!isdigit(c)) {if (c=='-') ag=-1; c=getchar();}
	while(isdigit(c)){f=(f<<3)+(f<<1)+c-48;c=getchar();}
	f*=ag; return;
}
void write(ll x) {
	if(x>9) write(x/10);putchar(x%10+48);return;
}
void writeln(ll x) {
	if (x<0) putchar('-'),x*=-1; write(x); putchar('\n'); return;
}
//以上装逼部分

//以下树状数组 
void add(int x,ll y) {
	for(;x<=n;x+=(x&-x)) s[x]+=y;
}
ll get(int x) {
	ll res=0;
	for(;x;x-=(x&-x)) res+=s[x];
	return res;
}
//以上树状数组 

//以下GCD 
ll gcd(ll a,ll b)
{
    return !a?b:gcd(b%a,a);
}
//以上GCD 

//以下线段树 
//建树来袭 
int len=0;
void build(int L,int R)
{
    int t=++len;
    l[t]=L;r[t]=R;
    if(L==R){num[t]=b[L];return;}
    int mid=(L+R)>>1;
    lc[t]=len+1,build(L,mid);
    rc[t]=len+1,build(mid+1,R);
    num[t]=gcd(num[lc[t]],num[rc[t]]);
}
//修改来袭 
void change(int now,int x,ll dat) {
	if (l[now]==r[now]) {num[now]+=dat;return;}
	int mid=(l[now]+r[now])>>1;
	if (x<=mid) change(lc[now],x,dat);
	else change(rc[now],x,dat);
	num[now]=gcd(num[lc[now]],num[rc[now]]);
}
//询问来袭 
ll Ask(int now,int x,int y) {
	if (x>y) return 1;
	if (l[now]==x&&r[now]==y) return num[now];
	int mid=(l[now]+r[now])>>1;
	if (y<=mid) return Ask(lc[now],x,y);
	else if (x>mid) return Ask(rc[now],x,y);
	else
	return gcd(Ask(lc[now],x,mid)
			,Ask(rc[now],mid+1,y));
}
//以上线段树
 
int main() {
//	freopen("input","r",stdin);
//	freopen("1.txt","w",stdout);
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++) read(a[i]);
	for (int i=1;i<=n;i++) b[i]=a[i+1]-a[i];
	add(1,a[1]);
	for (int i=2;i<=n;i++) add(i,b[i-1]);
    build(1,n-1);
	while(m--)
    {
        char c[2];
        scanf("%s",c);
        int l,r;
        scanf("%d%d",&l,&r);
        if(c[0]=='Q')
			writeln(abs(gcd(get(l),Ask(1,l,r-1))));
        else
        {
        	ll dat; read(dat);
            add(l,dat);
            if(l>1)change(1,l-1,dat);
            if(r<n)add(r+1,-dat),change(1,r,-dat);
        }
    }
}