CCF 除法

一、试题

问题描述
　　小葱喜欢除法，所以他给了你N个数a1, a2, ⋯, aN，并且希望你执行M次操作，每次操作可能有以下两种：
　　给你三个数l, r, v，你需要将al, al+1, ⋯, ar之间所有v的倍数除以v。
　　给你两个数l, r，你需要回答al + al+1 + ⋯ + ar的值是多少。
输入格式
　　第一行两个整数N, M，代表数的个数和操作的次数。
　　接下来一行N个整数，代表N个数一开始的值。
　　接下来M行，每行代表依次操作。每行开始有一个整数opt。如果opt=1，那么接下来有三个数l, r, v，代表这次操作需要将第l个数到第r个数中v的倍数除以v；如果opt = 2，那么接下来有两个数l, r，代表你需要回答第l个数到第r个数的和。
输出格式
　　对于每一次的第二种操作，输出一行代表这次操作所询问的值。
样例输入
5 3
1 2 3 4 5
2 1 5
1 1 3 2
2 1 5
样例输出
15
14
评测用例规模与约定
　　对于30%的评测用例，1 ≤ N, M ≤ 1000；
　　对于另外20%的评测用例，第一种操作中一定有l = r；
　　对于另外20%的评测用例，第一种操作中一定有l = 1 , r = N；
　　对于100%的评测用例，1 ≤ N, M ≤ 105，0 ≤ a1, a2, ⋯, aN ≤ 106, 1 ≤ v ≤ 106, 1 ≤ l ≤ r ≤ N。

二、代码

树状数组基础知识
如果使用常规算法只能30分就超时，感谢这位老兄算法
原版

#include<iostream>#include<math.h>#include<cstring>using namespace std;typedef long long LL;int N,M;int type,l,r,v;int A[100000+10];LL TreeArr[100000+10];//求解i二进制最后一个1后面有多少个0 int lowbit(int i)  {      return i&(-i);  } void toValue(int i,int num){   //每次读入一个数时都将树状数组的对应节点更新一下，当所有数都读取完毕后，整个树状数组也就建成。     while(i<=N)      {          TreeArr[i] += num;          i += lowbit(i);      }  }  //使用树状数组求解和 LL sum(int i)  {      LL total = 0;      while(i != 0)      {          total += TreeArr[i];          i -= lowbit(i);      }      return total;  } int main(){    cin>>N>>M;//  为啥 A[10^5+10]本地可以，CCF运行出错 //  long long *out = new long long[M];//  long long *head = out;//  memset(out, 0, sizeof(long long)*M);    memset(TreeArr,0,sizeof(TreeArr));    for(int j=1; j<=N; j++){        cin>>A[j];//      逐步更新树状数组的值，每次只是更新A[j]产生的父节点         toValue(j, A[j]);    }    while(M--){        cin>>type;        if(type==1){            cin>>l>>r>>v;            if(v==1) continue;            for(int i=l; i<=r; i++){                if(A[i]>=v && A[i]%v==0){//                  v&(v-1)==0表示v是2的n次幂                     if(v>0 && v&(v-1)==0){                        int n2 = log(v)/log(2);                        toValue(i, A[i]>>n2-A[i]);                        A[i]>>n2;                    }else{    //                  当A[i]值变化时，树状数组对于受影响的节点都需要更新                         toValue(i, A[i]/v-A[i]);                        A[i] /= v;    //                  if(v==2){    //                      A[i]>>1;    //                  }else{    //                      A[i] /= v;    //                  }                       }                }            }        }else if(type==2){            cin>>l>>r;            cout<< sum(r) - sum(l-1)<<endl; //          int total=0;//          for(int i=l; i<=r; i++){//              total += A[i];//          }//          *out = total;//          out++;//          cout<<total<<endl;        }    }//  out = head;//  while(*out){//      cout<<*out<<endl;//      out++;//  }//  delete[] head;    return 0;}