（算法竞赛进阶指南）余数之和（约数）

AC__dream

已于 2022-03-03 16:23:27 修改

阅读量324

点赞数

分类专栏：数学知识文章标签：算法优化数学问题连续整数取模运算区间求和

于 2022-03-03 16:22:39 首次发布

本文链接：https://blog.csdn.net/AC__dream/article/details/123257119

版权

数学知识专栏收录该内容

67 篇文章 2 订阅

订阅专栏

分析：由k%i=k-k/i*i(k/i是下取整)可得 $\sum_{i=1}^{n}$ (k%i)= $\sum_{i=1}^{n}$ （k-k/i*i）=n*k- $\sum_{i=1}^{n}$ (k/i*i)

而减号右边的部分可以用简便方法来进行计算，我们发现，k/i的值可能对于多个i都是相同的，而且这样的i是连续的，我们知道了这个如果能求出来满足k/i=j的区间[l,r]那么就可以直接求出这一块区间的k/i*i的和就是(l+r)*(r-l+1)*j了，现在问题就是求满足k/i=j的左右边界，或者换句话说就是知道左边界如何求取右边界，现在不妨假设左边界为l，那么右边界就是k/(k/l)，找几个样例试一下就不难发现这个式子，但是需要注意的是k和n的大小关系不确定，所以我们在计算右边界的时候需要注意当前求出的右边界与n的关系，还需要注意k/l=的情况，防止除0操作。

下面是代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
using namespace std;
int main()
{
	long long n,k;
	cin>>n>>k;
	long long ans=n*k;
	for(int l=1,r;l<=n;l=r+1)
	{
		if(k/l==0) break;
		r=min(k/(k/l),n);
		ans-=k/l*(l+r)*(r-l+1)/2;
	}
	printf("%lld",ans);
	return 0;
}