UVa 10692 Huge Mods (指数循环节)

UVa 10692 Huge Mods

题目大意:

给出模数$m$和正整数$a_1,a_2...a_n,$求出$a_1^{a_2^{...^a_n}}\ mod\ m$的值.
(注意指数运算的顺序:$2^{3^4}=2^{(3^4)}=2^{81}$)
($2\leq m\leq 10000,1\leq n \leq 10,1\leq a_i \leq 1000$)

题目分析:

主要是要运用到欧拉定理的推广——指数循环节.

$A^B\ mod\ C=A^{B\ mod\ \varphi (C)+\varphi(C)}\ mod\ C$

但是需要注意的是,指数运算自上而下的模数也会改变,如

$a^{b^{c^d}}\ mod\ p$

设$x=b^{c^d}$

$a^x=a^{x\ mod\ \varphi(p)+\varphi(p)}\ mod\ p$

那么对于$b^{c^d}$而言,其模数应为$\varphi(p)$

按照顺序,递归求解即可.

代码:

#include<cstdio>
#include<iostream>
#include<algorithm>

using namespace std;

const int maxn=10+5;
const int maxm=10000;

int phi[maxm+1];

void init(int n)
{
    phi[1]=1;
    for(int i=2;i<=n;i++) if(!phi[i])
        for(int j=i;j<=n;j+=i) {
            if(!phi[j]) phi[j]=j;
            phi[j]=phi[j]/i*(i-1);
        }
}

int pow_mod(int x,int y,int mod)
{
    int ret=1;
    while(y>0) {
        if(y&1) ret=ret*x%mod;
        x=x*x%mod;y>>=1;
    }
    return ret;
}

int e[maxn],n;

void dfs(int pos,int mod)
{
    if(pos<n-2) dfs(pos+1,phi[mod]);//递归边界a^b 
    e[pos]=pow_mod(e[pos],e[pos+1]%phi[mod]+phi[mod],mod);
}

int main()
{
    init(maxm);
    int mod,kase=0;
    while(scanf("%d",&mod)==1) {
        scanf("%d",&n);
        for(int i=0;i<n;i++) scanf("%d",&e[i]);
        if(n>1) dfs(0,mod);
        else e[0]%=mod;//n=1时,特殊处理 
        printf("Case #%d: %d\n",++kase,e[0]);
    }
    return 0;
}