$$
g(n)=\sum_{i=1}^n i\frac{(\lfloor\frac{n}{i}\rfloor+1)\times (\lfloor\frac{n}{i}\rfloor)}{2}
$$
同样也可以用整除分块实现。
---
代码
```cpp
#include <cstdio>
using namespace std;
const int mod=10007;
const int i4=(mod+1)>>2;
const int i6=(mod+1)/6;
int l,r;
int mo(int x,int y){return x<y?x-y+mod:x-y;}
void Mo(int &x,int y){x=x+y>=mod?x+y-mod:x+y;}
int sf(int n){
int _n=n%mod;
return 1ll*i6*_n*(_n+1)*(_n<<1|1)%mod;
}
//1^2+2^2+...+n^2=n*(n+1)*(2n+1)/6
int query0(int n){//f(n)
int ans=0;
for (int l=1,r;l<=n;l=r+1)
r=n/(n/l),Mo(ans,(r-l+1ll)*(n/l)%mod);
return ans;
}
int query1(int n){//g(n)
int ans=0;
for (int l=1,r;l<=n;l=r+1)
r=n/(n/l),Mo(ans,i4*(1ll*(l+r)%mod*(r-l+1)%mod)*(1ll*(n/l)*(n/l+1)%mod)%mod);
return ans;
}
int query(int n){
int ans=0;
for (int l=1,r;l<=n;l=r+1){
r=n/(n/l),Mo(ans,12*mo(sf(r),sf(l-1))%mod*query0(n/l)%mod);//左边这一部分
Mo(ans,24*(r-l+1ll)%mod*query1(n/l)%mod);//右边这一部分
}
return ans;
}
int main(){
scanf("%d%d",&l,&r);
printf("%d",mo(query(r),query(l-1)));
return 0;
}
```