poj1845 sumdiv 题解
Emmm...并非题解
其实是边想边写现编的
先审题：
考虑两个自然数 A 和 B。令 S 为 A\^B 的所有自然因子之和。确定 S 除以 9901 的余数.
eg. \(2^3 = 8\)。 8 的自然因子是：1、2、4、8。它们的和是 15。 15 除以 9901 的余数是 15（输出值）。
啊，好简洁的题面像我的大脑一样！
可以肯定的，\(A\)的因子一定是\(A^B\)的因子依旧废话
而且我们会发现，\(A^B\)的质因子一定与A的质因子完全相同。
很简单吧 谁问你了
思路
这时候，可以想到：(敲黑板，这是重点/)
对 A 进行 唯一分解
\[\large A = \prod _{i=1}^{k}{(P_i^{e_i})} \]
还是拆开好看：
\[\large A= P_1^{e_1} \times P_2^{e_2} \times ... \times P_k^{e_k} \]
------其中\(P_i\) 表示A的质因子，\(e_i\)表示该质因子的次数。
那\(A^B\)就是给每个\(e_i\)乘上一个B就好了。
应该没人不知道这个吧
很显然地，A的 因子 就是一个x，使得：
\[\large x= P_1^{v_1} \times P_2^{v_2} \times ... \times P_k^{v_k} \]
其中\(v_i < e_i\).
Emmm...还是不好求嘞...
求解
题目要求我们求所有x的和，
每个\(P_i\)会被算:
\[\tag {1} 1+P_i^1 + P_i^2 + P_i^3 + ... + P_i^{e_i} \]
引入一个另外的一次质因子\(P_j\), 与上面组合，答案就是：
\[\tag {2} 1 \times P_j + P_i^1 \times P_j + P_i^2 \times P_j + ... \]
可以发现，\((2)=(1)\times P_j\).
那么多次的\(P_j\)就是：
\[(1) + (1) \times P_j^1 + (1) \times P_j^2 +...(1) \times P_j^{e_j} \]
嘿，您猜怎么着， 咱把这 (1) 一提:
\[(1) \times (1+P_j^1+P_j^2+P_j^3+...+P_j^{e_j}) \]
诶呦喂，瞧瞧，我是不是在哪遇见过您？awa
这就是变了个样的 (1) 啊!
噫嘘兮，情乎喜哉！数论之易，易于切水题。
很好，我们只要求每个\(P_i\)所对应的(i),对其求积就可以了。
应该吧QAQ
拓展
没想到吧，这篇题解还有拓展awa
上面我们已经求出了这个解释，我还有什么要说的呢?
注意：
\[(I) = 1+P_i^1 + P_i^2 + P_i^3 + ... + P_i^{e_i} \]
有没有发现什么端倪?
oi,这不是个等比数列吗
知道我要说啥了吧awa
\[\tag{1} (I) = 1+P_i^1 + P_i^2 + P_i^3 + ... + P_i^{e_i} \]
\[\tag{2} P_i \times (I) = P_i^1 + P_i^2 + P_i^3 + ... + P_i^{e_i}+P_i^{e_i+1} \]
(2)-(1),得
\[(P_i^{e_i}-1) \times (I) = 1 + P_i^{e_i+1} \]
那么有
\[(I)=\frac{P_i^{e_i+1}+1}{P_i^{e_i}-1} \]
...太丑了吧
虽然丑，但它好求啊，这不比一个一个枚举快多了qwq
其实看有好多大佬还用了 乘法逆元 二次化简，我这个蒟蒻没学过就不在这里做洋相了
好不容易打了一下午的题解，浅浅要个赞不过分吧qwq
代码部分
.
.
.
.
.
.
.
.
.
自己写awa
特别鸣谢
完结撒花✿ヽ(^ ▽^)ノ✿✿