Sumdiv(约数和问题)

题目地址

看到这题的题解，大佬都说是小学奥数，~~蔡得我不敢鸡声~~。

求 \(a^b\) 所有的约数之和 mod \(9901\) \((1<=a,b<=5*10^7)\)

做这道题，我还赶紧去看了一下唯一分解定理

我们先把 \(a\) 分解质因数

\[a=p_1^{c_1}*p_2^{c_2}*...*p_n^{c_n}\]

比如说 \(12\) 可以分成 \(2^2+3^1\) 啦

因为 同指数幂相乘，指数不变，底数相乘 ，所以就有：

\[a^b=p_1^{c_1*b}*p_2^{c_2*b}*...*p_n^{c_n*b}\]

根据唯一分解定理，\(a^b\) 的约数和就是

\[(1+p_1+p_1^2+...+p_1^{c_1*b})*(1+p_2+p_2^2+...+p_2^{c_2*b})*...*(1+p_3+p_3^2+...+p_3^{c_3*b})\]

大佬看出了是等比数列，而我这个蒟蒻没有看出来

因为等比数列的求和公式要用除法，除法不满足 \(\text{mod}\) 的分配律

所以我们就迎来了这个题目的重点——分治

设 \(\text{sum}(p,c)\)，为 \((1+p+p^2+...+p^{c})\)

\[\text{sum}(p,c)=(1+p+...+p^{\frac{c-1}{2}})+(p^{\frac{c+1}{2}}+...+p^c)\]

\[=1*(1+p+...+p^{\frac{c-1}{2}})+\frac{c+1}{2}*(1+p+...+p^{\frac{c-1}{2}})\]

\[=(1+\frac{c+1}{2})*\text{sum}(p,\frac{c-1}{2})\]

\[\text{sum}(p,c)=(1+\frac{p}{2})*\text{sum}(p,\frac{p}{2}-1)*p^c\]

结合快速幂，时间复杂度上可以过得去

讲了这么多（虽然是看书），我忘了告诉你这个题目我是口胡的。