数学体操 模运算

---

概述

求下列式子的值：$$\sum _{i=1}^{n}nmodi.$$

其中，$1\le n\le 10^9.$

----------------------------------------
#思路
我们可以将$nmodi$写成$n-\lfloor \frac{n}{i}\rfloor \times i$，那么原式就变成了：$$Ans=\sum _{i=1}^{n}nmodi=\sum _{i=1}^{n}n-\lfloor \frac{n}{i}\rfloor \times i=(\sum _{i=1}^{n}n)-(\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor \times i).$$

整理一下：$$Ans=n^2-\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor \times i.$$

那么现在的关键就是求后面那个$\sum$的值了。

我们可以观察一下$n=11$时，每个$i$的$\lfloor \frac{n}{i}\rfloor$值：$$i:1234567891011$$$$\lfloor \frac{n}{i}\rfloor :115322111111$$

观察细致的读者也许发现了，同一个$\lfloor \frac{n}{i}\rfloor$值总是连续出现的，而且出现的次数似乎也有规律？

是的，令$i$的贡献$\lfloor \frac{n}{i}\rfloor =x$，那么$x$出现的次数为$\lfloor \frac{n}{x}\rfloor -\lfloor \frac{n}{x+1}\rfloor$次。

这个结论可以这么理解：贡献大于等于$x$的$\lfloor \frac{n}{i}\rfloor$一共有$\lfloor \frac{n}{x}\rfloor$个，贡献大于$x$的$\lfloor \frac{n}{i}\rfloor$一共有$\lfloor \frac{n}{x+1}\rfloor$个，两者相减即是贡献等于$x$的个数了。

所以，我们可以分块计算贡献。对于$\lfloor \frac{n}{i}\rfloor$相等的一个块，我们计算$\sum \lfloor \frac{n}{i}\rfloor \times i$可以提出公因数$\lfloor \frac{n}{i}\rfloor$，剩下$\sum i$是个等差数列，可以用求和公式直接计算。总复杂度均摊$O(\sqrt{n}).$

----------------------------------------

代码

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define ll long long
#define For(i,j,k) for(register int i=j; i<=(int)k; ++i)
#define Forr(i,j,k) for(register int i=j; i>=(int)k; --i)
#define INF 0x3f3f3f3f
using namespace std;

ll n, top, v, len, Ans;

int main(){
	scanf("%lld", &n);
	
	top = n;
	Ans = n*n;
	while(top >= 1){
		v = n/top;//计算当前位置的贡献.
		len = n/v - n/(v+1);//计算当前位置贡献出现的次数.
		Ans -= v*(2*top-len+1)*len/2;//更新答案.
		top -= len;//更新当前位置.
	}
	
	printf("%lld", Ans);
	return 0;
}

----------------------------------------

小结

本题重点在于将模运算转化成常见的四则运算，再观察贡献表（姑且这么叫吧）找出每个数的贡献所具有的规律，然后就可以快速计算了。

----------------------------------------

$——wrotebymiraclejzd.$