指数循环节例题(模板)

1. 2019南京网络赛super_log： 求(fuck[0]^fuck[1]^fuck[2]^...^fuck[b-1] mod m )>=b  其中fuck数组=a; 数据范围:1<=a<=1e6: 0<=b<=1e6,1<=m<=1e6

上代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MAX = 1e6+100;
ll phi[MAX],fuck[MAX];
void getphi(){
	phi[1]=1;
	for (ll i = 2; i < MAX;i++){
		if(!phi[i]){
			for (ll j = i; j < MAX;j+=i){
				if(!phi[j]) phi[j]=j;
				phi[j]=phi[j]/i*(i-1);
			}
		}
	}
}
ll quick(ll a,ll b,ll c){
	ll ans=1;
    if(a>=c) a=a%c+c;
    while(b){
        if(b&1){
            ans*=a;
            if(ans>=c) ans=ans%c+c;
        }
        a*=a;
        if(a>=c) a=a%c+c;
        b>>=1;
    }
    return ans;
}
ll a,b,m;
ll zsxhj(ll id,ll mod){// 求(fuck[0]^fuck[1]^fuck[2]^...^fuck[b-1] mod m )>= b
	if(b==0) return 1%m;
    if(mod==1) return 1;
    if (id==b-1){
        if(fuck[id]>=mod) return fuck[id]%mod+mod;
        else return fuck[id];
    }
    else return quick(fuck[id],zsxhj(id+1,phi[mod]),mod);
}
int main(){
	getphi();
	int t;
	scanf("%d",&t);
	while(t--){
		scanf("%lld%lld%lld",&a,&b,&m);
		for (int i = 0; i < b;i++){
			fuck[i]=a;
		}
		ll w=zsxhj(0,m);
		printf("%lld\n",w%m);
	}
	return 0;
} 

2. FZU-1759 求a^b%mod 数据范围：(1<=a,c<=1000000000,1<=b<=10^1000000).

上代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll get_euler(ll n){
    ll res=n,a=n;
    for (ll i = 2; i*i<=a;i++){
        if(a%i==0){
            res=res/i*(i-1);
            while(a%i==0) a=a/i;
        }
    }
    if(a>1) res=res/a*(a-1);
    return res;
}
ll quick(ll a,ll b,ll c){
	ll ans=1;
    if(a>=c) a=a%c+c;
    while(b){
        if(b&1){
            ans*=a;
            if(ans>=c) ans=ans%c+c;
        }
        a*=a;
        if(a>=c) a=a%c+c;
        b>>=1;
    }
    return ans;
}
int main(){
	ll a,c;
	string b;
	while(~scanf("%lld",&a)){
		cin >> b;
		scanf("%lld",&c);
		ll p=get_euler(c);
		int len=b.size();
		ll sum=0;
		for (int i = 0; i < len;i++){
			sum=sum*10+(b[i]-'0');
			sum%=p;
		}
		printf("%lld\n",quick(a,sum,c)%c);
	}
	return 0;
} 

3. 计蒜客 Exponial  求(n^(n-1)^(n-2)^...^1 )mod m 数据范围：1<=n<=1e9,1<=m<=1e9

上代码：

#include<bits/stdc++.h>
typedef long long ll;
int get_euler(int x){
    int ans=x;
    for(int i=2;i*i<=x;i++){
        if(x%i==0){
            ans-=ans/i;
            while(x%i==0) x/=i;
        }
    }
    if(x>1) ans-=ans/x;
    return ans;
}
ll quick(ll a,ll b,ll c){
    ll ans=1;
    if(a>=c) a=a%c+c;
    while(b){
        if(b&1){
            ans*=a;
            if(ans>=c) ans=ans%c+c;
        }
        a*=a;
        if(a>=c) a=a%c+c;
        b>>=1;
    }
    return ans;
}
ll exponial(int x,int y){//求(x^(x-1)^(x-2)^...^1 )mod y 
    if(x==1) return 1%y;
    if(x%y==0) return quick(x,x-1,y);
    return quick(x,exponial(x-1,get_euler(y)),y);
}
int main(){
    int n,m; 
    while(~scanf("%d%d",&n,&m)){
        printf("%lld\n",exponial(n,m)%m);
    }
    return 0;
}