Codeforces616 E. Sum of Remainders（数论分块）

题意：

$给定n和m，要求计算{\sum }_{i=1}^{m}n\%i$

$数据范围：n,m<=1e13$

解法：

$$\begin{gathered} \sum _{i=1}^{m}n\%i \\ =\sum _{i=1}^{m}n-\frac{n}{i}\ast i \\ =n\ast m-\sum _{i=1}^m\frac{n}{i}\ast i \\ 发现有\frac{n}{i},那么可以用整除分块优化到O(sq) \end{gathered}$$

code：

#include <bits/stdc++.h>
#define int long long
using namespace std;
const int mod=1e9+7;
int n,m;
int inv2;
int ppow(int a,int b,int mod){
    int ans=1%mod;a%=mod;
    while(b){
        if(b&1)ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}
int cal(int x){
    return x%mod*((x+1)%mod)%mod*inv2%mod;
}
signed main(){
    ios::sync_with_stdio(0);
    inv2=ppow(2,mod-2,mod);
    cin>>n>>m;
    int ans=n%mod*(m%mod)%mod;
    for(int l=1,r;l<=m;l=r+1){
        if(n/l)r=min(m,n/(n/l));
        else break;//l>n的时候,n/l=0,因此提前退出.
        int temp=(n/l)%mod*((cal(r)-cal(l-1))%mod)%mod;
        ans=(ans-temp)%mod;
    }
    ans=(ans%mod+mod)%mod;
    cout<<ans<<endl;
    return 0;
}

---