POJ 1845 Sumdiv（高中数学，推公式，分治）

转自bin神：http://www.cnblogs.com/kuangbin/archive/2012/08/10/2631225.html

bin神语录（来自群acfun）：
　　１．年轻人，多刷题少水群
　　２．弃蓝翔杯保平安
　　３．acm大法好

题目大意：求Ａ的Ｂ次方的所有约数的和，结果模9901．Ａ和Ｂ的范围都是50000000。
解题思路：bin神讲的很清楚，可能我的代码和bin神的思路有些不太一样，反正无关紧要了。关键的地方是分治下去，究竟公式是怎么样的（我也没有仔细看bin神的公式），都差不多的，但是最好还是自己推导一遍，只要能得到正确的答案就行了。
Show me the code!

#include <iostream>
#define ll long long//注意５千万乘５千万爆int
using namespace std;
const int mod = 9901;
ll pow(ll a, ll b) {
    ll ans = 1;
    while (b) {
        if (b & 1) ans = ans * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return ans;
}
ll solve(ll p, ll k) {
    if (k == 0) return 1;
    if (k == 1) return 1 + p;
    if (k & 1) {
        ll t = solve(p, k / 2);
        return (t + t * pow(p, k / 2 + 1) % mod) % mod;
    }
    else {
        int t = solve(p, k / 2 - 1);
        return (t + t * pow(p, k / 2 + 1) % mod + pow(p, k / 2)) % mod;
    }
}
int main() {
    ll a, b;
    cin >> a >> b;
    ll ans = 1;
    for (int i = 2; i * i <= a; ++i) {
        ll t = 0;
        while (a % i == 0) {
            t++;
            a /= i;
        }
        ans = ans * solve(i, t * b) % mod;
    }
    if (a != 1) ans = ans * solve(a, b) % mod;//for循环中停止条件是i*i <= a，所以ａ最后可能不会变成１
    if (a == 0) ans = 0;//测试数据中可能会出现ａ==0的情况
    cout << ans << endl;
    return 0;
}