codeforces-914D Bash and a Tough Math Puzzle(线段树维护gcd)

链接：http://codeforces.com/problemset/problem/914/D

题意：给定一个序列，给出两种操作1 l r x表示判断l到r区间内能否通过改变一个数字使得该区间的gcd=x，2 i y把第i个数改为y

题解：通过这道题才知道gcd也是有区间性质的，那么就是一道普通的单点修改加区间查询了。可以通过判断递归到叶子判断该数能否整除x，不能的话返回1进行累加，要注意的是当一个区间的gcd已经能整除x就不需要继续递归下去了，否则会超时

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

int n, a[5*100010];
int q;

int getgcd(int a, int b){
    return b == 0? a:getgcd(b, a%b);
}

struct node{
    int l, r, gcd;
}s[4*5*100010];

void build(int l, int r, int k){
    s[k].l = l; s[k].r = r;
    if(l == r){
        s[k].gcd = a[l];
        return;
    }
    int mid = (l+r)>>1;
    build(l, mid, 2*k);
    build(mid+1, r, 2*k+1);
    s[k].gcd = getgcd(s[2*k].gcd, s[2*k+1].gcd);
}

void update(int k, int i, int c){
    if(s[k].l == s[k].r){
        s[k].gcd = c;
        return;
    }
    int mid = (s[k].l+s[k].r)>>1;
    if(i<=mid) update(2*k, i, c);
    else if(i>mid) update(2*k+1, i, c);
    s[k].gcd = getgcd(s[2*k].gcd, s[2*k+1].gcd);
}

int query(int l, int r, int k, int x){
    if(s[k].l == s[k].r){
        //cout<<"! "<<k<<" "<<s[k].gcd<<endl;
        if(s[k].gcd%x == 0) return 0;
        else return 1;
    }
    if(l>=s[k].l && r<=s[k].r && s[k].gcd%x == 0) return 0;
    //if(s[k].gcd>=2) return 2;
    int mid = (s[k].l+s[k].r)>>1;
    if(r<=mid) {
        return query(l, r, 2*k, x);
    }
    else if(l>mid) {
        return query(l, r, 2*k+1, x);
    }
    else {
        //cout<<"???"<<endl;
        int res1 = query(l, mid, 2*k, x);
        if(res1>=2) return 2;
        int res2 = query(mid+1, r, 2*k+1, x);
        return res1+res2;
    }

}

int main(){
    scanf("%d", &n);
    for(int i = 1; i<=n; i++){
        scanf("%d", &a[i]);
    }
    build(1, n, 1);
    scanf("%d", &q);
    while(q--){
        int ch, l, r, x, i, y;
        scanf("%d", &ch);
        if(ch == 1){
            scanf("%d%d%d", &l, &r, &x);
            if(query(l, r, 1, x)>=2) printf("NO\n");
            else printf("YES\n");
        }
        else if(ch == 2){
            scanf("%d%d", &i, &y);
            update(1, i, y);
            /*for(int i = 1; i<=10; i++){
                cout<<i<<" "<<" "<<s[i].l<<" "<<s[i].r<<" "<<s[i].gcd<<endl;
            }*/
        }
    }

    return 0;
}