#整除分块#洛谷 2261 余数之和

题目

给定$n$和$k$，求$$(k\mathrm{m}\mathrm{o}\mathrm{d}1)+(k\mathrm{m}\mathrm{o}\mathrm{d}2)+\dots +(k\mathrm{m}\mathrm{o}\mathrm{d}n-1)+(k\mathrm{m}\mathrm{o}\mathrm{d}n)$$

---

分析

因为$(k\mathrm{m}\mathrm{o}\mathrm{d}n)=k-\lfloor k÷n\rfloor$，所以最终答案=$$n\times k-\sum _{i=1}^n\lfloor k÷i\rfloor \times i$$
因为还是会超时，所以要用一种分块的思想，
设$$g(x)=\lfloor k÷\lfloor k÷x\rfloor \rfloor$$，那么$$g(x)\ge \lfloor k÷(k÷x)\rfloor =x$$
且$$\lfloor k÷g(x)\rfloor \ge \lfloor k÷(k÷\lfloor k÷x\rfloor )\rfloor =\lfloor k÷x\rfloor$$，所以$$\lfloor k÷g(x)\rfloor =\lfloor k÷x\rfloor$$
那么就可以用整除分块求解
然后用等差数列求答案即可

---

代码

#include <cstdio>
#include <algorithm>
long long n,k,ans;
int main(){
	scanf("%lld%lld",&n,&k); ans=n*k;
	for (register long long x=1,gx;x<=n;x=gx+1){
		gx=k/x?std::min(k/(k/x),n):n;//有约数为x
		ans-=(k/x)*(x+gx)*(gx-x+1)>>1;//等差数列
	}
	printf("%lld",ans);
	return 0;
}