bzoj5028小Z的加油店（线段树+差分）

题意：维护支持以下两种操作的序列：1 l r询问a[l...r]的gcd，2 l r x把a[l...r]全部+x

题解：一道经典题。根据gcd(a,b)=gcd(a-b,b)以及区间加可知，这题可以差分求解。然后只需要维护支持单点加减和区间查询gcd、差分和的线段树即可，复杂度O(nlog²n)，因为要求区间gcd。

#include<bits/stdc++.h>
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
using namespace std;
const int N=1e5+;
int n,m,a[N],d[N<<],s[N<<];
void pushup(int rt){d[rt]=__gcd(d[rt<<],d[rt<<|]),s[rt]=s[rt<<]+s[rt<<|];}
void build(int l,int r,int rt)
{
    if(l==r){d[rt]=s[rt]=a[l];return;}
    int mid=l+r>>;
    build(lson),build(rson);
    pushup(rt);
}
void update(int k,int v,int l,int r,int rt)
{
    if(l==r){a[l]+=v,d[rt]=s[rt]=a[l];return;}
    int mid=l+r>>;
    if(k<=mid)update(k,v,lson);else update(k,v,rson);
    pushup(rt);
}
int query_d(int L,int R,int l,int r,int rt)
{
    if(L>R)return ;
    if(L<=l&&r<=R)return d[rt];
    int mid=l+r>>,ret=;
    if(L<=mid)ret=__gcd(ret,query_d(L,R,lson));
    if(R>mid)ret=__gcd(ret,query_d(L,R,rson));
    return ret;
}
int query_s(int L,int R,int l,int r,int rt)
{
    if(L<=l&&r<=R)return s[rt];
    int mid=l+r>>,ret=;
    if(L<=mid)ret+=query_s(L,R,lson);
    if(R>mid)ret+=query_s(L,R,rson);
    return ret;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=;i<=n;i++)scanf("%d",&a[i]);
    for(int i=n;i;i--)a[i]-=a[i-];
    build(,n,);
    while(m--)
    {
        int op,l,r,x;scanf("%d%d%d",&op,&l,&r);
        if(op==)printf("%d\n",abs(__gcd(query_d(l+,r,,n,),query_s(,l,,n,))));
        else{
            scanf("%d",&x);
            update(l,x,,n,);
            if(r<n)update(r+,-x,,n,);
        }
    }
}