【数学扩展欧拉定理】luogu_4139 上帝与集合的正确用法

最新推荐文章于 2024-07-29 23:24:37 发布

原创最新推荐文章于 2024-07-29 23:24:37 发布 · 302 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#数学 #扩展欧拉定理

收录于

数学专栏收录该内容

53 篇文章

订阅专栏

本文探讨了在解决特定数学问题时，如何利用欧拉定理及其推论进行模幂运算的优化。通过定义函数f(p)，文章详细阐述了如何递归地计算2的指数次方模p的值，特别关注于当指数非常大时的计算方法。代码部分展示了使用C++实现这一算法的过程，包括幂运算、欧拉函数计算以及主函数调用。

题意

求 $2^{2^{2...}}\ mod\ p$

思路

设 $f(p)=2^{2^{2...}}\ mod\ p$
根据欧拉定理的推论， $p，(b>φ(p))a^b\ mod\ p=a^{b\% \varphi(p)+\varphi(p)}\ mod\ p，(b>\varphi(p))$
所以 $p2^{2^{2..}}\ mod\ p=2^{f(\varphi(p))+\varphi(p)}\ mod \ p$
递归即可。

代码

#include<cstdio>

int t, p;
int phi[10000001];

int power(int a, int b, int m) {
	long long res = 1;
	for (; b; b >>= 1) {
		if (b & 1) res = res * a % m; 
		a = (long long)a * a % m;
	}
	return res;
}

int solve(int p) {
	if (p == 1) return 0;
	return power(2, solve(phi[p]) + phi[p], p);
}

int main() {
	phi[1] = 1;
	for (int i = 2; i <= 10000000; i++)
		if (!phi[i]) 
			for (int j = i; j <= 10000000; j += i) {
				if (!phi[j]) phi[j] = j;
				phi[j] = phi[j] / i * (i - 1);
			}
	scanf("%d", &t);
	for (; t; t--) {
		scanf("%d", &p);
		printf("%d\n", solve(p));
	}
}