URAL 2062 Ambitious Experiment（树状数组）

题目地址：http://acm.timus.ru/problem.aspx?space=1&num=2062

思路：对于x位置上的数，对其有贡献的数为位置为x的因数且属于修改区间[l,r]的数。所以对于查询操作ans=a[p]+sigma(sum(j)) (j为p的约数）。区间修改时仅需修改区间端点值：c[l]+=d，c[r+1]-=d（当数x在某个无重复区间[l,r]时，每次求前x项的和时已包含该区间即已加上该数贡献，当x>r时，求前x项和时，该区间和为0，即数x贡献为0）。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn=3e5+50;
int n,q;
LL c[maxn],a[maxn];
int lowbit(int x)
{
    return x&(-x);
}
void add(int x,int w)
{
    while(x<=n)
    {
        c[x]+=w;
        x+=lowbit(x);
    }
}
LL sum(int x)
{
    LL tot=0;
    while(x)
    {
        tot+=c[x];
        x-=lowbit(x);
    }
    return tot;
}
int main()
{
    scanf("%d",&n);
    for(int i=1; i<=n; i++)
        scanf("%d",&a[i]);
    scanf("%d",&q);
    while(q--)
    {
        int id;
        scanf("%d",&id);
        if(id==2)
        {
            int l,r,w;
            scanf("%d%d%d",&l,&r,&w);
            add(l,w);
            add(r+1,-w);
        }
        else
        {
            int p;
            scanf("%d",&p);
            LL ans=a[p];
            for(int i=1; i*i<=p; i++)
            {
                if(p%i) continue;
                ans+=sum(i);
                if(i!=p/i) ans+=sum(p/i);
            }
            printf("%I64d\n",ans);
        }
    }
    return 0;
}