codeforce_914D(GCD+线段树)

题意:给出一个数组,查询某一段的GCD是否猜测正确,假如这段最多有一个数不是所有数GCD的倍数,那么就算猜对了,否则就是猜错。
思路:用线段树维护每一段的GCD,单点更新不用说,那么在查询的时候,需要稍微做点手脚,向下查询的时候,如果说左子树的GCD不能整除x,就说明至少在左子树这个区间有一个数不是x的倍数,这些数可能在要找的区间[a,b]里,也可能不在,而是在左子树-[a,b]的集合里,如果是后者,线段树在寻找GCD的过程中就会将这些外面的部分自动撇弃,如果是前者,线段树找到最后l==r的时候,就说明这个数就是那个毒瘤,现在ans++,如果ans>1,那么不用再查询,直接输出“NO”

#include<bits/stdc++.h>
#define maxn 500005*4
#define ls rt<<1
#define rs rt<<1|1
using namespace std;
int v[maxn],a[maxn],ans;
int gcd(int a,int b)
{
    if(b==0) return a;
    else return gcd(b,a%b);
}
void up(int rt){v[rt]=gcd(v[ls],v[rs]);}

void bulid(int l,int r,int rt)
{
    if(l==r) {v[rt]=a[l];return;}
    int mid=(l+r)>>1;
    bulid(l,mid,ls);
    bulid(mid+1,r,rs);
    up(rt);
}

void modfiy(int l,int r,int rt,int pos,int val)
{
    if(l==r) {v[rt]=val;return;} int mid=(l+r)>>1;
    if(pos<=mid) { modfiy(l,mid,ls,pos,val);}
    else {modfiy(mid+1,r,rs,pos,val);}
    up(rt);
}
void query(int l,int r,int s,int t,int rt,int x)
{
    if(l==r) {ans++;return ;}
    if(ans>1) return;
    int mid=(l+r)>>1;
    if(s<=mid&&v[ls]%x!=0){query(l,mid,s,t,ls,x);}
    if(t>mid&&v[rs]%x!=0) {query(mid+1,r,s,t,rs,x);}
}
int main()
{
    int n,m,op,l,r,x;
    scanf("%d",&n);
    for(int i=1;i<=n;i++){ scanf("%d",&a[i]);}
    memset(v,0,sizeof(0));
    bulid(1,n,1);
    scanf("%d",&m);
    for(int i=0;i<m;i++){
        scanf("%d",&op);
        if(op==1){
            ans=0;
            scanf("%d%d%d",&l,&r,&x);
            query(1,n,l,r,1,x);
            if(ans<=1) cout<<"YES"<<endl;else cout<<"NO"<<endl;
        }
        else{
            int pos;
            scanf("%d%d",&pos,&x);
            modfiy(1,n,1,pos,x);
        }
    }
    return 0;
}
  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值