欧拉定理——PowMod ( HDU 5728 )(2016 Multi-University Training Contest 1 1006)

• 题目链接：
  http://acm.hdu.edu.cn/showproblem.php?pid=5728

• 分析：
  题中有 k = ∑(m,i=1) φ(i∗n) mod 1000000007，即求和函数 sum(n, m) = ∑(m,i=1) φ(i∗n) mod 1000000007. 那么考虑将质因子p的贡献单独处理，可以得到 sum(n, m) = φ(p) × sum(n/p, m) + sum( n, m/p )，求出k之后使用欧拉定理递归求值。

• 题解：

  1. sum(n, m) = φ(p) × sum(n/p, m) + sum( n, m/p )
     这个式子可以这么理解：
     欧拉函数是不完全积性函数，在i 与 p 互质（因为n是无平方因子数，所以n/p肯定和p互质）的时候有φ(i∗n) = φ(i*n/p) * φ(p)，但是当i与p不互质的时候是不成立的所以当i∈(1~m)且i与n的质因数p互质时候，可以把p提出来列为左边这项；当不互质时，即i∈(1×p， m)，令 j = i×p，那么有j ∈ (1，m/p)就得到右边这项。然后不断这样转换递归能够简化计算。

  2. 欧拉定理递归求值，即a^b mod p = a^(φ(p) + b%φ(p))，然后不断递归计算，直到φ(φ(…φ(p))) = 1，然后再递归回去。

• 参考代码：

#include <iostream>
#include <vector>
#include <cmath>
#include <algorithm>
#include <ctime>
#include <cstdio>
#include <cstring>
#include <map>

#define Max 10000001
using namespace std;
typedef long long LL;
const int mod = 1000000007;



int vis[Max], eular[Max], prime[Max], cnt, a[20];
LL sum[Max];
void init()
{
    cnt = 0;
    eular[1] = 1;
    int N = Max;
    int k;
    for (int i = 2;i<N;i++)
    {
        if (!vis[i])
            prime[cnt++] = vis[i] = i, eular[i] = i - 1;
        for (int j = 0;j<cnt && (k = prime[j] * i)<N;j++)
        {
            vis[k] = prime[j];
            if (vis[i] == prime[j])
            {
                eular[k] = eular[i] * prime[j];
                break;
            }
            else
                eular[k] = eular[i] * (prime[j] - 1);
        }
    }
}
LL pow(LL a,LL n,LL mod)
{
    LL base = a, ret = 1;
    while (n)
    {
        if (n % 2)
            ret = (ret*base) % mod;
        base = (base*base) % mod;
        n /= 2;
    }
    return ret;
}
LL cal_sum(int pos,LL n,LL m)
{
    if (n == 1)//sum(1,m)=sigma(eular(i))
        return sum[m];
    if (m == 0)//sum(n,0)=0;
        return 0;
    return ((a[pos] - 1)*cal_sum(pos - 1, n / a[pos], m) % mod + cal_sum(pos, n, m / a[pos])) % mod;
    //按照sum(n,m)=eular(p)*sum(n/p,m)+sum(n,m/p),计算，其中p为n分解出来的所有质数
}
LL cal_ans(LL k,LL mod) //欧拉定理递归计算k的k次的k次的k次的k次......（无数个k次）
{
    if (mod == 1)
        return 0;
    LL temp = cal_ans(k, eular[mod]);//b%eular(p)
    LL ans = pow(k,temp+eular[mod],mod) %mod;
    return ans;
}
int main()
{
    init();
    sum[0] = 0;
    for (int i = 1;i < Max;i++)
        sum[i] = (sum[i - 1] + eular[i]) % mod;
    int n, m, p;
    while (scanf("%d %d %d", &n, &m, &p) != EOF)
    {
        int now = n, num = 0;
        for (int i = 0;i < cnt;i++)
        {
            if (!vis[n])//n是质数直接存进a数组
            {
                a[num++] = now;
                break;
            }
            if (now%prime[i] == 0)//若n不是质数则分解质因数存进a数组
            {
                a[num++] = prime[i];
                now /= prime[i];
            }
        }
        //num -1表示n含有多少个质数
        LL temp = cal_sum(num - 1, n, m);//temp表示题目表达式中的 k 值
        LL ans = cal_ans(temp, p);
        printf("%lld\n", ans);
    }
    return 0;
}