一个性质:n%k=n-k*(n/k);(利用整除的性质)

注意:sigma(i)不需要考虑除法自动向下取整的问题(本身就全都是整数,通项公式本身是绝对正确的,那么在/2之前必然已经是2的倍数了)


 

#include<cstdio>
#include<iostream>
#define ll long long
using namespace std;

int main(){
	ll n,k;//不开long long见祖宗
	scanf("%lld%lld",&n,&k);
	
	ll sum=n*k;
	n=min(n,k);//特判n<k 
	for(ll l=1,r;l<=n;l=r+1){
		 r=k/(k/l);
		 if(r>=n)r=n;
		 ll tmp=k/l;
		 sum-=tmp*(r-l+1)*(l+r)/2;
	} 
	printf("%lld\n",sum);
	return 0;
}