洛谷 P3935 Calculating 题解 & 因数个数定理

题目传送门

题目大意： 若 $x$ 可以分解成 $p_1^{k_1}{p}_2^{k_2}{p}_3^{k_3}...p_s^{k_s}$，那么 $f(x)=(k_1+1)(k_2+1)...(k_s+1)$，求出 ${\sum }_{i=l}^{r}f(i)$

前置知识——因数个数定理

描述： 将一个数 $x$ 质因数分解后，可以表示为 $x=p_1^{c_1}{p}_2^{c_2}...p_k^{c_k}$，那么 $x$ 的因数个数就是 $\prod _{i=1}^k(c_i+1)$。

证明： 由于 $x$ 的每个因子也可以表示上面的形式，只不过 $c$ 的值不同，设因子的质因数分解后为 $p_1^{c_1^{\prime }}{p}_2^{c_2^{\prime }}...p_k^{c_k^{\prime }}$，对于 $c_i^{\prime }$，它的取值范围是 $0$ ~ $c_i$，即有 $(c_i+1)$ 种取值方案，所以总因子个数就是 $\prod _{i=1}^k(c_i+1)$。

题解

发现题目定义的 $f(i)$ 其实就是 $i$ 的因数个数。

那么考虑换一个角度，统计每一个因数 $x$ 会被多少个 $i$ 计算到。

显然，$x$ 的每个倍数都会将 $x$ 计算一次，将询问差分一下，变成 ${\sum }_{i=1}^{r}f(i)-{\sum }_{i=1}^{l-1}f(i)$，那么问题就变成了求 $1$ ~ $n$ 范围内有多少个 $x$ 的倍数，即

求 ${\sum }_{i=1}^n\lfloor \frac{n}{i}\rfloor$

然而这是个除法分块的模板。

代码如下：

#include <cstdio>
#include <cstring>
#define mod 998244353ll
#define ll long long

ll a,b;
ll solve(ll x)
{
	ll l=1,r,re=0;
	while(l<=x)
	{
		r=x/(x/l);
		re=(re+(r-l+1)*(x/l)%mod)%mod;
		l=r+1;
	}
	return re;
}

int main()
{
	scanf("%lld %lld",&a,&b);
	printf("%lld",(solve(b)-solve(a-1)+mod)%mod);
}