AcWing 246. 区间最大公约数 (gcd性质 线段树)

简介: AcWing 246. 区间最大公约数 (gcd性质 线段树)

20200401134307494.png

线段树支持标记和标记的合并,值和值的合并,标记和值的合并。

维护原序列的话,不知道打了区间修改的标记后gcd的变化。

转化为差分序列后,单点修改区间查询gcd。时间复杂度nlogn

gcd(a,b)=gcd(a,b-a)

const int maxn=500010;
int n,m;
ll w[maxn];
struct node{
    int l,r;
    ll sum,d;
}tr[maxn*4];
ll gcd(ll a,ll b){
    return b==0?a:gcd(b,a%b);
}
void pushup(node &u,node &l,node &r){
    u.sum=l.sum+r.sum;
    u.d=gcd(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 tmp=w[r]-w[r-1];
        tr[u].l=l;tr[u].r=r;
        tr[u].d=tmp;tr[u].sum=tmp;
    }
    else{
        tr[u].l=l;tr[u].r=r;
        int mid=l+r>>1;
        build(u<<1,l,mid);
        build(u<<1|1,mid+1,r);
        pushup(u);
    }
}
void update(int u,int x,ll v){
    if(tr[u].l==x&&tr[u].r==x){
        ll tmp=tr[u].sum+v;
        tr[u].l=x;tr[u].r=x;
        tr[u].d=tmp;tr[u].sum=tmp;
    }
    else{
        int mid=tr[u].l+tr[u].r>>1;
        if(x<=mid) update(u<<1,x,v);
        else update(u<<1|1,x,v);
        pushup(u);
    }
}
node query(int u,int l,int r){
    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{
            node left=query(u<<1,l,r);
            node right=query(u<<1|1,l,r);
            node res;
            pushup(res,left,right);
            return res;
        }
    }
}
int main(){
    n=read();m=read();
   for(int i=1;i<=n;i++) w[i]=read();
   build(1,1,n);
   int x,y;char op[2];
    while(m--){
        cin>>op;x=read();y=read();
        if(op[0]=='Q'){
            node left=query(1,1,x),right=query(1,x+1,y);
            printf("%lld\n",abs(gcd(left.sum,right.d)));
        }
        else{
            ll d;d=read();
            update(1,x,d);
            if(y+1<=n) update(1,y+1,-d);
        }
    }
    return 0;
}
目录
相关文章
|
算法 测试技术 C++
C++二分算法习题:判断是否是完全平方数[容易]和排列箱子[容易]
C++二分算法习题:判断是否是完全平方数[容易]和排列箱子[容易]
|
算法 C语言 C++
【数论】最大公约数、约数的个数与约数之和定理
先来科普下什么是约数:当a能被b整除,我们就说b为a的约数,b的倍数为a
138 0
|
人工智能
Leetcode53/152—最大子数组和/最大子数组乘积(状态转移方程/不熟)
Leetcode53/152—最大子数组和/最大子数组乘积(状态转移方程/不熟)
114 0
|
人工智能
归并排序求逆序对 acwing例题代码
归并排序求逆序对 acwing例题代码
容斥原理 (两个例题)
容斥原理 (两个例题)
158 0
|
算法 C++
快排例题——第k个数
做道简单一点的题巩固一下
算法--递归辗转相除法求最大公约数
算法--递归辗转相除法求最大公约数
对分查找、欧几里得算法求最大公约数
对分查找、欧几里得算法求最大公约数