51nod1225 余数之和（分块）

1225 余数之和

1.0 秒 131,072.0 KB 80 分 5级题
F(n) = (n % 1) + (n % 2) + (n % 3) + … (n % n)。其中%表示Mod，也就是余数。
例如F(6) = 6 % 1 + 6 % 2 + 6 % 3 + 6 % 4 + 6 % 5 + 6 % 6 = 0 + 0 + 0 + 2 + 1 + 0 = 3。
给出n，计算F(n), 由于结果很大，输出Mod 1000000007的结果即可。

输入

输入1个数N(2 <= N <= 10^12)。

输出

输出F(n) Mod 1000000007的结果。

输入样例

6

输出样例

3

代码

#include<bits/stdc++.h>

using namespace std;
typedef long long ll;
const ll mod = 1e9 + 7;
ll N;

ll quick_pow(ll a, ll n, ll Mod) {
    ll ans = 1;
    while (n) {
        if (n & 1) {
            ans = (ans * a) % Mod;
        }
        a = (a * a) % Mod;
        n >>= 1;
    }
    return ans;
}

ll inv = quick_pow(2, mod - 2, mod);

int main() {
    while (~scanf("%lld", &N)) {
        ll ans = (N % mod) * (N % mod) % mod;
        for (ll i = 1; i <= N;) {
            ll tmp = N / i;
            ll t = N / tmp;
            // i 首项 t 末项 t-i+1 项数
            tmp = (i + t) % mod * ((t - i + 1) % mod) % mod * inv % mod * tmp % mod;
            ans = (ans % mod - tmp % mod + mod) % mod;
//            cout << tmp << " " << t << " " << i + t << " " << t - i + 1 << " " << endl;
            i = t + 1;
        }
        printf("%lld\n", ans);
    }
    return 0;
}