传送门:URAL - 2062
题意:给定n个数的初始值,和两种操作:
1.询问第i个数是多少
2.给定l,r,d,将i,2*i,3*i...位置上的数加d,l<=i<=r.
思路:采用平方分割算法,也叫块状数组,算法思想大致就是讲整个区间分成长度为sqrt(n)的块,然后每次更新如果包括完整的块就直接维护整个块,然后将剩余部分单点更新。具体讲解见挑程序设计第二版P183.
就本题来说,对第i个数有贡献的数的下标就是i的所有因数。因此询问某个点时暴力因数分解一下然后加和就好了。
代码:
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<string>
#include<math.h>
#include<algorithm>
#include<queue>
#include<stack>
#include<set>
#include<vector>
#include<map>
#define ll long long
#define pi acos(-1)
#define inf 0x3f3f3f3f
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
using namespace std;
typedef pair<int,int>P;
const int MAXN=300010;
int n,num[MAXN];
int L[3000],R[3000],id[MAXN];
ll interval[3000],own[MAXN];
void init()
{
int t=sqrt(n);
int cnt=n/t;
if(n%t)cnt++;
for(int i=1;i<=cnt;i++)
{
L[i]=(i-1)*t+1;
R[i]=(i)*t;
}
for(int i=1;i<=n;i++)
id[i]=(i-1)/t+1;
}
void update(int l,int r,int d)
{
if(id[l]==id[r])
{
for(int i=l;i<=r;i++)
own[i]+=d;
return ;
}
for(int i=l;i<=R[id[l]];i++)
own[i]+=d;
for(int i=id[l]+1;i<=id[r]-1;i++)
interval[i]+=d;
for(int i=L[id[r]];i<=r;i++)
own[i]+=d;
}
ll solve(int k)
{
ll ans=0;
for(int i=1;i*i<=k;i++)
{
if(k%i==0)
{
ans+=interval[id[i]]+own[i];
if(i*i!=k)
ans+=interval[id[k/i]]+own[k/i];
}
}
return ans;
}
int main()
{
int m,l,r,d,k;
cin>>n;
for(int i=1;i<=n;i++)
scanf("%d",num+i);
init();
cin>>m;
while(m--)
{
scanf("%d",&k);
if(k==2)
{
scanf("%d%d%d",&l,&r,&d);
update(l,r,d);
}
else
{
scanf("%d",&k);
printf("%lld\n",solve(k)+num[k]);
}
}
return 0;
}
网上有很多用树状数组解这个题的,但是基本上都是卡过去的,时间上限很紧,2700+ms。说不准什么时候T了。