题目
思路
其实完成这道题只需要2个思路:
- 高精度除法,举个栗子:11 mod 5=11-5*(11/5)(没错,是取下整)=1
- 分块
来,看看柿子:
∑
i
=
1
n
k
m
o
d
i
=
n
∗
k
−
∑
i
−
1
n
k
i
\sum_{i=1}^{n} k mod i=n*k-\sum_{i-1}^{n}\frac{k}{i}
∑i=1nkmodi=n∗k−∑i−1nik
接下来对于k/i这个部分分块,AC
code:
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
long long n,k,s,sx,mx;
int main()
{
cin>>n>>k;
s=n*k;
for (sx=1;sx<=n;sx=mx+1)
{
if (k<sx) mx=n;
else mx=min(n,k/(k/sx));
s-=(mx-sx+1)*(sx+mx)/2*(k/sx);//式子+等差数列求和公式,这部分抄题解的人一定要抄对,不抄的当我没说
}
cout<<s;
return 0;
}