算法1

(根据商特征分块)

分析：
根据 k mod i = k - floor(k / i ) i,可将j(n,k)转化成如下式子：
f(n,k) = k mod 1+ k mod 2 + …k mod n
　　　 = k - floor(k/1) * 1 + k - floor(k/2) * 2+… k - floor(k /n) n
　　　 = n * k - sum(floor(k/i) * i) 其中 i = 1~ n 。
　　接下来我们拿n= 9 ,k = 4 为例分析下k/i：
　　4/1　4/2　4/3　4/4　4/5　4/6　4/7　4/8　4/9
　　4　　2　　1　　1　　0　　0　　0　　0　　0
其中sum(floor(k/i) * i) = 4 * 1 + 2 * 2 + 1 (3+ 4）+ 0 (5+6+7+8+9）,括号里的都是等差数列，可以由高斯求和的方法处理。那么我们就需要知道哪几个连续的商是相同的呢？那么在当前的i数时，找到与k/i想等的最大的gx,使得k / gx = k / i. 因此gx = k / floor(k/i)，为了避免算出来的gx超出n,因此准确来说gx = min(k / floor(k/i) , n);这样便可以将i~ gx之间的数用高斯求和处理了。然后继续从 gx+ 1为起点寻找连续块。

C++ 代码

#include<bits/stdc++.h>
using namespace std;
int main(){
    long long n,k,ans;
    cin >> n >> k;
    ans = n * k;
    int gx;
    for(int i = 1;i<= n; i= gx + 1 ){
        gx = k/i!=0 ?min( k / (k / i), n ) : n; // 找到与k/i相等的 k/x中最大的x,因为k/i ~ k/x相等的商，使得这段里是等差序列。 
        ans = ans - k/i * (i + gx)* (gx - i + 1) / 2;   
    }
    cout << ans << endl;
    return 0;
}