codeforces 616E Sum of Remainders 数学公式转化

Calculate the value of the sum: n mod 1 + n mod 2 + n mod 3 + … + n mod m. As the result can be very large, you should print the value modulo 109 + 7 (the remainder when divided by 109 + 7).

The modulo operator a mod b stands for the remainder after dividing a by b. For example 10 mod 3 = 1.

Input
The only line contains two integers n, m (1 ≤ n, m ≤ 1013) — the parameters of the sum.

Output
Print integer s — the value of the required sum modulo 109 + 7.

Example
Input
3 4
Output
4
Input
4 4
Output
1
Input
1 1
Output
0

题意

给你两个数n,m
问你n % 1 + n % 2 + … + n% m为几

公式化简 n%i=n-n/i*i 那么就要求 n * m - sum(n/i * i) (i从1到m)

由于 当到达i的时候 从i到 n/(n/i) 的 n/i 的值是一样的 因为如果说 n/i 得到了一个因子的下界，那么n/这个因子。则是这个因子的上界。那么就可以那个求和公式求每一段了。

求得过程有些坑点需要注意，首先是注意m的大小，m的大小可能位于n的两个因子之间
然后是，求和公式那里也需要先求下余，再算乘法，这样就需要注意因为有个除2，所以要判断用哪一项能整除2，除完后再求余

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
int main()
{
    long long n,m;
    cin>>n>>m;
    ll all = (n%mod)*(m%mod);
    all%=mod;
    long long res=0;
    ll t=min(n,m);
    for(ll i=1;i<=t;i++)
    {
        ll l=i;
        ll r=min((n/(n/i)),m);
        long long num=r-l+1;
        long long sw=(l+r);
        if(sw%2==0) sw/=2;
        else num/=2;
        num%=mod,sw%=mod;
        ll sum=(num%mod)*(sw%mod);//这里不取模居然不行。。好奇怪
        sum%=mod;
        sum*=(n/i);
        sum%=mod;
        res+=sum;
        res%=mod;
        i=r;
    }
    res%=mod;
   cout<<(all-res+mod)%mod<<endl;
}