苦恼的小明

版权声明：本文为博主原创文章，转载请注明出处：http://blog.csdn.net/cggwz https://blog.csdn.net/cggwz/article/details/83316723

cgg传送门

这题的解法是费马小定理+欧拉函数+快速幂

整理：费马小定理：
a^(φ(m))≡1(mod m) ((a,m)=1) （phi为欧拉函数）

所以得：
(a^b)mod m=(a^(b%phi(m)))mod m

欧拉函数线筛打表就好了

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int phi[10020];
int n;
int num[1234570],phm[1234570];
void init(){
	for(int i=0;i<=10010;i++){
		phi[i]=i;
	}
	phi[1]=1;
	for(int i=2;i<=10010;i++){
		if(phi[i]==i){
			for(int j=i;j<=10010;j+=i){
				phi[j]=phi[j]/i*(i-1);
			}
		}
	}
}
int fpow(int a,int n,int p){
	int ans=1;
	while(n>0){
		if(n&1){
			ans=ans*a%p;
		} 
		a=a*a%p;
		n>>=1;
	}
	return ans;
}
int main(){
	init();
	scanf("%d",&n);
	phm[1]=10007;
	for(int i=1;i<=n;i++){
		scanf("%d",&num[i]);
		phm[i+1]=phi[phm[i]];
	}
	int cur=num[n]%phm[n];
	for(int i=n-1;i>=1;i--){
		cur=fpow(num[i],cur,phm[i]);
	}
	printf("%d",cur);
	return 0;
}