华华开始学信息学

题目链接：https://ac.nowcoder.com/acm/problem/23054

题意：

初始数组为长度为n的全0数组，m次操作
操作1：使每个（i-0）mod d==0 的下标为i的数组元素+k
操作2：查询l~r范围的数组元素和

思路：

分块思想，
使lazy[i]表示原数组下标能被i整除的每个元素增加量
显然当d较小时操作数量陡增
既然如此，当d<=sqrt（n）时用lazy数组O(1)存起来，当d>sqrt（n）直接用树状数组change方法
查询时答案分为两个部分：树状数组的query部分+lazy数组下标从1~n元素*（r/i-（l-1）/i）（意思是该范围内能被i整除的元素个数）

#include<bits/stdc++.h>
#define int long long
#define rep(i,a,n) for(int i=a;i<=n;i++)
#define pb push_back
#define endl "\n"
#pragma GCC optimize(3)
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int inf=0x3f3f3f3f;
const ll llmax=LLONG_MAX;
const int maxn=1e5+5;
int n,m;
int s[maxn];
int lowbit(int x){
	return x&-x;
}
void change(int x,int k){
	while(x<=n){
		s[x]+=k;
		x+=lowbit(x);
	}
}
int query(int x){
	int t=0;
	while(x){
		t+=s[x];
		x-=lowbit(x);
	}
	return t;
}
int lazy[maxn];
signed main()
{
	ios::sync_with_stdio(false),cin.tie(0);
	cin>>n>>m;
	while(m--){
		int op;cin>>op;
		if(op==1){
			int d,k;cin>>d>>k;
			if(d>sqrt(n)){
				for(int i=d;i<=n;i+=d){
					change(i,k);
				}
			}else{
				lazy[d]+=k;
			}
		}else{
			int ans=0;
			int l,r;cin>>l>>r;
			ans+=(query(r)-query(l-1));
			for(int i=1;i<=sqrt(n);i++){
				ans+=((r/i)*lazy[i]-(l-1)/i*lazy[i]);
			}
			cout<<ans<<endl;
		}
	}
	return 0;
}