914D Bash and a Tough Math Puzzle

传送门

分析

用线段树维护区间gcd，每次查询找到第一个不是x倍数的点，如果这之后还有gcd不能被x整除的区间则这个区间不合法

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cctype>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<ctime>
#include<vector>
#include<set>
#include<map>
#include<stack>
using namespace std;
int d[],cnt,a[];
inline int gcd(int x,int y){return y==?x:gcd(y,x%y);}
inline void build(int le,int ri,int wh){
      if(le==ri){
          d[wh]=a[le];
          return;
      }
      int mid=(le+ri)>>;
      build(le,mid,wh<<);
      build(mid+,ri,wh<<|);
      d[wh]=gcd(d[wh<<],d[wh<<|]);
      return;
}
inline void update(int le,int ri,int wh,int pl,int k){
      if(le==ri){
          d[wh]=k;
          return;
      }
      int mid=(le+ri)>>;
      if(mid>=pl)update(le,mid,wh<<,pl,k);
        else update(mid+,ri,wh<<|,pl,k);
      d[wh]=gcd(d[wh<<],d[wh<<|]);
      return;
}
inline void q(int le,int ri,int wh,int x,int y,int k){
      if(le>=x&&ri<=y){
          if(d[wh]%k==)return;
          if(cnt>=){
            cnt++;
            return;
          }
      }
      if(le==ri){
          cnt++;
          return;
      }
      int mid=(le+ri)>>;
      if(mid>=x)q(le,mid,wh<<,x,y,k);
      if(mid<y)q(mid+,ri,wh<<|,x,y,k);
      return;
}
int main(){
      int n,m,i,j,k,x,y,z;
      scanf("%d",&n);
      for(i=;i<=n;i++)
        scanf("%d\n",&a[i]);
      build(,n,);
      scanf("%d",&m);
      for(i=;i<=m;i++){
          scanf("%d",&k);
          if(k==){
            scanf("%d%d%d",&x,&y,&z);
            cnt=;
            q(,n,,x,y,z);
            if(cnt>)puts("NO");
              else puts("YES");
          }else {
            scanf("%d%d",&x,&y);
            update(,n,,x,y);
          }
      }
      return ;
}