欧拉定理及扩展

欧拉定理:对于任意互素的a和n，有$a^{\phi (n)}\equiv 1(modn)$
证明：
记小于n且与n互素的数为 R = {$x_1,x_2...x_{\phi }(n)$}
令S = { a x 1 % n , a x 2 % n , ⋯   , a x φ ( n ) % n } \{ax_1\%n,ax_2\%n,\cdots,ax_{\varphi(n)}\%n\} {ax1​%n,ax2​%n,⋯,axφ(n)​%n}

$\forall i\in [1,\phi (n)]$
$\because (a,n)=1,(x_i,n)=1$
$\therefore (a{x}_i,n)=1$
$由最大公约数的性质可得(a{x}_i\%n,n)=1$
$所以S中所有元素都与n互质，且都小于n。$

$\because x_1\%n=x_1,x_2\%n=x_2,\cdots ,且a\%n为定值$
$\therefore S中无重复元素$
$\therefore S=R$

$\therefore {\prod }_{i=1}^{\phi (n)}a{x}_i\%n={\prod }_{i=1}^{\phi (n)}{x}_i$
$\therefore {\prod }_{i=1}^{\phi (n)}a{x}_i\equiv {\prod }_{i=1}^{\phi (n)}{x}_i(\mathrm{m}\mathrm{o}\mathrm{d}n)⟹$ $a^{\phi (n)}{\prod }_{i=1}^{\phi (n)}{x}_i\equiv {\prod }_{i=1}^{\phi (n)}{x}_i(\mathrm{m}\mathrm{o}\mathrm{d}n)$
$又({\prod }_{i=1}^{\phi (n)}{x}_i,n)=1$
$故a^{\phi (n)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$

扩展欧拉定理(也是常说的欧拉降幂)
对于任意的a,b,m，$若b>=\phi (m)$ ：$a^b\equiv a^{b\mathrm{m}\mathrm{o}\mathrm{d}\phi (m)+\phi (m)}(modm)$

求a ^ a ^ a ^ a ^ … ^ a % c，一共b个a

#include <cstdio>
typedef long long ll;

ll euler(ll n)    //计算模数的欧拉函数 
{
	ll res = n,a = n;
	for (int i = 2; i * i <= n; i++)
	{
		if( a % i == 0 )
		{
			res -= res / i;
			while( a % i == 0 ) a /= i;
		}
	}
	if( a > 1 ) res -= res / a;
	return res;
}

ll q_pow(ll a,ll b,ll mod)    //快速幂 
{
	ll ans = 1;
	int flag = 0;
	while( b )
	{
		if( b & 1 ) 
		{
			ans *= a;
			if( ans > mod ) flag = 1;
			ans %= mod;
		}
		a *= a;
		if( a > mod ) flag = 1;
		a %= mod;
		b >>= 1;
	}
	return ans + flag * mod;   
	//不能简单的返回%mod的值,大于mod时要加mod（欧拉定理的约束） 
}

ll cal(ll a,ll l,ll r,ll mod)
{
	if( mod == 1 ) return 1;   //模数为1时直接返回1回溯,因为后面等于多少到这一层%1都等于0 
	if( l == r )      //递归到最后,回溯 
	{
		if( a > mod ) a = a % mod + mod;
		return a;
	}
	ll b = cal(a,l+1,r,euler(mod));  //计算下一层的指数(注意模数变为当前值的欧拉值) 
	return q_pow(a,b,mod);    //快速幂计算 
}

int main()
{
	int t;
	scanf("%d",&t);
	while( t-- )
	{
		ll a,b,c;
		ll ans;
		scanf("%lld%lld%lld",&a,&b,&c);
		if( b == 0 ) ans = 1;
		else if( b == 1 ) ans = a;
		else ans = cal(a,1,b,c);
		ans %= c;
		printf("%lld\n",ans);
	}
	return 0;
}