URAL 2062 Ambitious Experiment（分块）

【题目链接】 http://acm.timus.ru/problem.aspx?space=1&num=2062

【题目大意】

　　给出两个操作，操作一给出区间[l,r]，对l到r中的每一个下标i，
　　i，2*i，3*i……位置都增加x，操作二要求查询一个位置的当前值

【题解】

　　在修改的时候，我们将增量只标识在i上，查询的时候，查询下标的因子和即可。
　　考虑到这种查询方法需要sqrt(n)枚举判断因子，用二叉数据结构可能会超时，
　　因此我们用分块nsqrt(n)修改，O(1)查询特定位置。

【代码】

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
using namespace std;
typedef long long LL;
const int MAX_SIZE=550,MAX_N=300010;
LL block[MAX_SIZE],u[MAX_SIZE][MAX_SIZE];
int n,m,a[MAX_N],op,x,l,r,size;
void add(int l,int r,int v){
    int L=l/size,R=r/size;
    if(L==R)for(int i=l%size;i<=r%size;i++)u[L][i]+=v;
    else{
        for(int i=l%size;i<size;i++)u[L][i]+=v;
        for(int i=L+1;i<R;i++)block[i]+=v;
        for(int i=0;i<=r%size;i++)u[R][i]+=v;
    }
}
LL query(int x){
    LL ans=0;
    for(int i=1;i*i<=x;i++){
        if(x%i==0){
            int t=i-1;
            ans+=block[t/size]+u[t/size][t%size];
            if(i*i!=x){
                int t=x/i-1;
                ans+=block[t/size]+u[t/size][t%size];
            }
        }
    }return ans;
}
int main(){
    while(~scanf("%d",&n)){
        memset(block,0,sizeof(block));
        memset(u,0,sizeof(u));
        for(int i=1;i<=n;i++)scanf("%d",&a[i]);
        scanf("%d",&m);
        size=sqrt(n);
        while(m--){
            scanf("%d",&op);
            if(op==1){
                scanf("%d",&x);
                printf("%lld\n",query(x)+a[x]);
            }else{
                scanf("%d%d%d",&l,&r,&x);
                add(--l,--r,x);
            }
        }
    }return 0;
}