「洛谷 4139」上帝与集合的正确用法

最新推荐文章于 2022-05-07 12:09:16 发布

forever_dreams

最新推荐文章于 2022-05-07 12:09:16 发布

阅读量199

点赞数

分类专栏： --------数学问题-------- 文章标签：欧拉定理降幂

本文链接：https://blog.csdn.net/forever_dreams/article/details/103016656

版权

--------数学问题-------- 专栏收录该内容

141 篇文章 1 订阅

订阅专栏

传送门

problem

有 $T$ 组数据，每组数据给定 $p$ ，求：

$2^{2^{2^{\cdots}}}\bmod p$

数据范围： $T\le1000$ ， $p\le10^7$ 。

solution

一开始觉得不可做，后来仔细想想还是挺简单的。

首先肯定有降幂：

$a^b\equiv \begin{cases} a^b,&b<\varphi(p)\\ a^{b\bmod \varphi(p)+\varphi(p)},&b\geq \varphi(p) \end{cases} \;\pmod p$

那么就有：

$2^{\left(2^{2^{\cdots}}\bmod \varphi(p)+\varphi(p)\right)}\bmod p$

于是就可以递归做了，递归到 $\varphi(p)=1$ 停止即可。

时间复杂度 $O(p+T\log p)$ 。

code

#include<bits/stdc++.h>
using namespace std;
const int N=1e7+5;
int T,P;
int add(int x,int y)  {return x+y>=P?x+y-P:x+y;}
int dec(int x,int y)  {return x-y< 0?x-y+P:x-y;}
int mul(int x,int y)  {return 1ll*x*y>=P?1ll*x*y%P:x*y;}
int power(int a,int b,int p){
	int ans=1;
	for(;b;b>>=1,a=1ll*a*a%p)  if(b&1)  ans=1ll*ans*a%p;
	return ans;
}
int sum,prime[N],mark[N],phi[N];
void linear_sieves(){
	phi[1]=1;
	for(int i=2;i<N;++i){
		if(!mark[i])  prime[++sum]=i,phi[i]=i-1;
		for(int j=1;j<=sum&&i*prime[j]<N;++j){
			mark[i*prime[j]]=1;
			if(i%prime[j])  phi[i*prime[j]]=phi[i]*(prime[j]-1);
			else  {phi[i*prime[j]]=phi[i]*prime[j];break;}
		}
	}
}
int solve(int p){
	if(p==1)  return 0;
	return power(2,solve(phi[p])+phi[p],p);
}
int main(){
	freopen("data.in","r",stdin);
	scanf("%d",&T);
	linear_sieves();
	while(T--){
		scanf("%d",&P);
		printf("%d\n",solve(P));
	}
	return 0;
}