2019南京网络赛-B-super_log（欧拉降幂）

题目链接

题意

求幂塔函数，b个a mod m .

扩展欧拉定理

欧拉定理概念

代码参考

理解：这两种情况必须要分，但是特例很少，所以很多代码没有考虑到第二种情况（B>=phi(C)）但是仍然可以过，应该是数据水，没有特例的情况，也可能不会出现特例？？

比如 6^1 % 2560 = 6 ， 6^ (1+phi(2560))  %2560   =>  6^1025 =  4096  ，2560的欧拉值是1024，小于1025，所以这就是第二种情况的一个特例。但是很多其他的都符合，不过考虑加上phi(C). 用指数循环节理解：在开头一部分，比较特殊，在之后才有循环，但是开头的这个并不一定是按照循环来的，所以开头的和后面循环的不一定相等（但是大多数情况都是相等的，只有少数特例），所以才有上述公式的两种情况，表示循环节前的和循环节。

求欧拉函数模板

打表，用于处理多个数据，数据范围最大大概为5亿（一维数组大小）

ll phi[N],prime[N];
ll tot;
void init(){
    phi[1]=1;
    for(ll i=2;i<N;i++){
        if(!phi[i]){
            phi[i]=i-1;
            prime[tot++]=i;
        }
        for(ll j=0;j<tot && i*prime[j]<N;j++){
            if(i%prime[j]) phi[i*prime[j]] = phi[i]*(prime[j]-1);
            else{
                phi[i*prime[j]] = phi[i]*prime[j];
                break;
            }
        }
    }
}

单次计算，如果多次计算，每次都要重复计算，但是可以处理很大到long long 的数据

ll phi(ll n)
{
    ll k = (ll)sqrt(n+0.5);
    ll ans = n;
    for(int i=2;i<=k;i++) {
        if(n%i==0){
            ans = ans/i*(i-1);
            while(n%i==0) n /= i;
        }
    }
    if(n>1) ans = ans/n*(n - 1);
    return ans;
}

欧拉降幂

有一种黑科技可以不用考虑B和phi(C)的比较  ，用一个自定义的MOD函数可以回避这种比较讨论，黑科技证明

ll MOD(ll x,ll m)
{
	return x>m?x%m+m:x;
}
ll Pow(ll x,ll y,ll m)
{
	ll ans=1;
	while(y){
		if(y&1) ans = MOD(ans*x,m);
		x = MOD(x*x,m);
		y >>= 1;
	}
	return ans;
}

用了自定义的MOD之后，需要注意的是当取模的数p，if(p==1)return a;  不是返回0。有很多其他的降幂都是写的返回0，但是在其他的地方分类比较了，但感觉更多的是因为数据水？欧拉定理对不互质的两种情况的公式，很多博客的代码都是回避的，也就是直接用快速幂，没有考虑取模的是大于phi还是小于phi，但是都能AC.

---

#include <bits/stdc++.h>
using namespace std;
#define N 1000010

typedef long long ll;

ll phi[N],prime[N];
ll tot;
void init(){
    phi[1]=1;
    for(ll i=2;i<N;i++){
        if(!phi[i]){
            phi[i]=i-1;
            prime[tot++]=i;
        }
        for(ll j=0;j<tot && i*prime[j]<N;j++){
            if(i%prime[j]) phi[i*prime[j]] = phi[i]*(prime[j]-1);
            else{
                phi[i*prime[j]] = phi[i]*prime[j];
                break;
            }
        }
    }
}

ll MOD(ll x,ll m)
{
	return x>m ? x%m+m : x;
}
ll Pow(ll x,ll y,ll m)
{
	ll ans=1;
	while(y){
		if(y&1) ans = MOD(ans*x,m);
		x = MOD(x*x,m);
		y >>= 1;
	}
	return ans;
}

ll solve(ll a,ll b,ll p){
	if(p==1)return a;
	if(b==0)return 1;
	ll P=solve(a,b-1,phi[p]);
	return Pow(a,P,p);
}

int main()
{
	init();
	int t;scanf("%d",&t);
	while(t--){
		ll a,b,mod;
		scanf("%lld%lld%lld",&a,&b,&mod);
		ll ans = solve(a,b,mod);
		printf("%lld\n",ans%mod);
	}
	return 0;
}