4302 Interval GCD 线段树

描述
给定一个长度为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(x,y)=gcd(x,y-x)，可以推广到多个数，所以可以求原序列差分后的gcd，这样就可以做到单点修改。
然后还要用一个树状数组维护原序列。

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define LL long long
const int Maxn=2000010;
int n,m;
LL a[Maxn],b[Maxn],s[Maxn];
void add(int x,LL y){for(;x<=n;x+=(x&-x))s[x]+=y;}
LL get(int x){LL re=0;for(;x;x-=(x&-x))re+=s[x];return re;}
LL gcd(LL a,LL b)
{
    if(!a)return b;
    return gcd(b%a,a);
}
struct Seg{int l,r,lc,rc;LL g;}tr[Maxn<<1];
int len=0;
void build(int l,int r)
{
    int t=++len;
    tr[t].l=l;tr[t].r=r;
    if(l==r){tr[t].g=b[l];return;}
    int mid=l+r>>1;
    tr[t].lc=len+1,build(l,mid);
    tr[t].rc=len+1,build(mid+1,r);
    tr[t].g=gcd(tr[tr[t].lc].g,tr[tr[t].rc].g);
}
void modify(int x,int p,LL d)
{
    if(tr[x].l==tr[x].r){tr[x].g+=d;return;}
    int mid=tr[x].l+tr[x].r>>1,lc=tr[x].lc,rc=tr[x].rc;
    if(p<=mid)modify(lc,p,d);
    else modify(rc,p,d);
    tr[x].g=gcd(tr[lc].g,tr[rc].g);
}
LL query(int x,int l,int r)
{
    if(l>r)return 1;
    if(tr[x].l==l&&tr[x].r==r)return tr[x].g;
    int mid=tr[x].l+tr[x].r>>1,lc=tr[x].lc,rc=tr[x].rc;
    if(r<=mid)return query(lc,l,r);
    else if(l>mid)return query(rc,l,r);
    else return gcd(query(lc,l,mid),query(rc,mid+1,r));
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)scanf("%lld",&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 op[2];
        scanf("%s",op);
        int l,r;LL d;
        scanf("%d%d",&l,&r);
        if(op[0]=='Q')printf("%lld\n",abs(gcd(get(l),query(1,l,r-1))));
        else
        {
            scanf("%lld",&d);
            add(l,d);
            if(l>1)modify(1,l-1,d);
            if(r<n)add(r+1,-d),modify(1,r,-d);
        }
    }
}