HDU 5728 PowMod 欧拉函数 递归

---

感觉智商被掏空…

---

原题见HDU 5728

定义$k=\sum_{i=1}^{m} \varphi (i·n)\ mod\ 1000000007$
$n$是无质因子平方项的数.
求$ans=k^{k^{k^{k^{...^k}}}}(\ mod \ p)$，其中k有无穷多个
数据范围：$1 \leq n, m, p \leq 10^{7}$

---

各种引理

1. 若素数p满足$p^2\mid n$，则$\varphi(n)=p· \varphi({n\over p})$ 见性质3
2. 若素数p满足$p\mid n, ~p^2\nmid n$，则$\varphi(n)=(p-1)·\varphi({n\over p})$
3. $a^b \equiv a^{\varphi (p) + b\%\varphi (p)} (\mod p)$
   当$gcd(a,p)=1$，有欧拉定理$a^{\varphi(p)}\equiv 1(\mod p)$，等式成立。
   当$gcd(a,p)=g>1$，只要证$a^{\varphi(p)}\equiv a^{2\varphi(p)}(\mod p)$，即可得到$a^{\varphi(p)}\equiv a^{\varphi(p)的任意多倍}$
   只要证$a^{\varphi(p)}(a^{\varphi(p)}-1)=p的倍数(*)$
   令$a=gx,p=gy,则gcd(x,y)=1$，有$g\mid a^{\varphi(p)}$
   又$\varphi(y)\mid \varphi(p)，则y\mid a^{\varphi(y)}-1\mid a^{\varphi(p)}-1$即证(*)成立。
   $a^{\varphi(p)}\equiv 0或1(\mod p)$

---

分析

题意很清晰啦。数据范围这么大，肯定要先化简。
$n$是无质因子平方项的数.这个条件很关键，说明任取n中的因子g，满足$gcd(g,{n\over g})=1(*)$
取i，满足$gcd(i,n)=g$，则$gcd({i\over g},{n\over g})=1$，综合$(*)可得gcd(g^2\times{i\over g},{n\over g})=1$
则由$\varphi (i·n)=\varphi(i·g·{n\over g})=\varphi(i·g)\varphi({n\over g})$
又$g\mid i$,则$\varphi(i·g)=g·\varphi(i)$

$$\begin{align*} \sum_{i=1}^{m} \varphi (i·n)=&\sum_{g\mid n}\sum_{i=1}^{m} g·\varphi(i)\varphi({n\over g})[gcd(i,n)=g]\\ =&\sum_{g\mid n}g·\varphi({n\over g})\sum_{i=1}^{m} \varphi(i)[gcd(i,n)=g]\\ =&\sum_{g\mid n}g·\varphi({n\over g})\sum_{i=1}^{m/g} \varphi(i·g)[gcd(i,n)=1] \end{align*}$$
然后…莫比乌斯反演展开…并没有什么卵用

所以看了题解，不按套路啊…

任取n中的质数p，若$p\mid i$，则$p^2\mid i·n$，则由引理1：$\varphi (i·n)=p·\varphi (i·{n\over p})$
若$p\nmid i$，则$p^2\nmid i·n$，则由引理2：$\varphi (i·n)=(p-1)\varphi (i·{n\over p})$

$$\begin{align*} \sum_{i=1}^{m} \varphi (i·n)=&\sum_{p\mid i,i=1}^{m} p·\varphi (i·{n\over p})+\sum_{p\nmid i,i=1}^{m} (p-1)\varphi (i·{n\over p})\\ =&\sum_{p\mid i,i=1}^{m} (p-1+1)·\varphi (i·{n\over p})+\sum_{p\nmid i,i=1}^{m} (p-1)\varphi (i·{n\over p})\\ ^{前面的p-1项合并到后面}=&\sum_{p\mid i,i=1}^{m} \varphi (i·{n\over p})+(p-1)\sum_{i=1}^{m} \varphi (i·{n\over p})\\ =&\sum_{i=1}^{m/p} \varphi (i·n)+(p-1)\sum_{i=1}^{m} \varphi (i·{n\over p}) \end{align*}$$
令 $S(n,m)=\sum\limits_{i=1}^{m} \varphi (i·n)$，则 $$S(n,m)=S(n,{m\over p})+(p-1)S({n\over p},m)$$
处理一下欧拉函数前缀和，递归处理即可。
已经询问过的S(n,m)可以保存在 map<pair<n,m>, S>里，不用反复查询。

然后看看$ans=k^{k^{k^{k^{...^k}}}}(\ mod \ p)$
两种错误的想法，样例都不能过，不想清楚还是很抓狂…
认为$ans=k^{ans}(\ mod \ p)$或者$ans=ans^k(\ mod \ p)$的，就不要玩了。
一般情况下，$a^x\not\equiv a^{x\% p}(\mod p)$

令$A(k, p) = k^{k^{k^{k^{...^k}}}}(\ mod \ p)$
由引理3，$k^b \equiv k^{\varphi (p) + b\%\varphi (p)} (\mod p)$，$b=k^{k^{k^{k^{...^k}}}}$，可得

$$A(k,p)=k^{\varphi(p)+A(k,\varphi(p))}$$

---

代码

#include <bits/stdc++.h>
using namespace std;
#define N 10000010
#define LL long long
int MOD = 1000000007;
int num[N], prim[N], phi[N], sum[N];
int pn = 0;
void init(){
    memset(num, -1, sizeof(num));
    phi[0] = phi[1] = 1;
    for(int i = 2;i < N;i++){
        if(num[i]){
            prim[pn++] = i;
            phi[i] = i-1;
        }
        for(int j = 0;j < pn && 1LL*i*prim[j] < N;j++){
            if(i % prim[j] == 0){
                phi[i*prim[j]] = phi[i] * prim[j];
                num[i*prim[j]] = 0;
                break;
            }
            phi[i*prim[j]] = phi[i] * (prim[j] - 1);
            num[i*prim[j]] = 0;
        }
    }
    sum[0] = 0;
    for(int i = 1;i < N;i++)
        sum[i] = 1LL*(sum[i-1] + phi[i]) % MOD;
}
int s(int n, int m){
    if(n == 1)
        return sum[m];
    if(m == 0)
        return 0;
    for(int i = 0;i < pn && prim[i] <= n;i++){
        if(n % prim[i] == 0){
            int p = prim[i];
            return (1LL*(p-1)*s(n/p, m)%MOD + s(n, m/p)%MOD)%MOD;
        }
    }
}
LL poo(LL a, int k, int m){
    LL res = 1;
    while(k){
        if(k & 1LL) res = res*a%m;
        k >>= 1;
        a = a*a%m;
    }
    return res;
}
int A(int k, int p){
    if(p == 1) return 0;
    if(k == 1) return 1;
    return poo(1LL*k%p, phi[p]+A(k, phi[p]), p);
}
int main(){
    init();
    int n, m, p;
    while(~scanf("%d%d%d", &n, &m, &p)){
        int k = s(n, m);
        int ans = A(k, p);
        printf("%d\n", ans);
    }
    return 0;
}