区间线段公约数(线段树+差分)

题目连接

#include <bits/stdc++.h>

using namespace std;
#define ll long long
const int maxn=5e5+10;
int m,n,ans,p;
ll w[maxn];
struct node
{
    int l,r;
    ll sum,d;

} tr[maxn*4];

ll gcb(ll x,ll y)
{
    return y?gcb(y,x%y):x;
}

void pushup(node &u,node &l,node&r)
{
    u.sum=l.sum+r.sum;
    u.d=gcb(l.d,r.d);
}

void pushup(int u)
{
    pushup(tr[u],tr[u<<1],tr[u<<1|1]);
}

void build(int u,int l,int r)
{
    if(l==r) //叶子节点
    {
        ll b=w[r]-w[r-1]; //求差分
        tr[u]= {l,r,b,b};
    }
    else
    {
        int mid=l+r>>1;
        tr[u].l=l,tr[u].r=r;
        build(u<<1,l,mid),build(u<<1|1,mid+1,r); //向左和向右递归建树
        pushup(u);//更新当前结点的前缀和,和最大公约数
    }
}

void modify(int u,int x,ll v)
{
    if(tr[u].l==x&&tr[u].r==x)
    {
        ll b=tr[u].sum+v;
        tr[u]= {x,x,b,b};
    }
    else
    {
        int mid=tr[u].l+tr[u].r>>1;
        if(x<=mid)
            modify(u<<1,x,v);
        else modify(u<<1|1,x,v);
        pushup(u);
    }
}

node query(int u,int l,int r)
{
    if(l>r)return {0};
    if(tr[u].l>=l&&tr[u].r<=r)return tr[u];
    else
    {
        int mid=tr[u].l+tr[u].r>>1;
        if(r<=mid)
        {
            return query(u<<1,l,r);
        }
        else if(l>mid)
        {
            return query(u<<1|1,l,r);
        }
        else
        {
            auto left=query(u<<1,l,r);
            auto right=query(u<<1|1,l,r);
            node res;
            pushup(res,left,right);
            return res;
        }
    }
}

int main()
{

    scanf("%d%d",&n,&m);
    for(int i=1; i<=n; i++)scanf("%lld",&w[i]);
    build(1,1,n);
    int l,r;
    ll d;
    char op[2];
    while(m--)
    {
        scanf("%s%d%d",op,&l,&r);
        if(*op=='Q')
        {
            auto left=query(1,1,l),right=query(1,l+1,r);
            printf("%lld\n",abs(gcb(left.sum,right.d)));
        }
        else
        {
            scanf("%lld",&d);
            //指定结点进行修改
            modify(1,l,d);
            if(r+1<=n)modify(1,r+1,-d);
        }
    }



    return 0;
}