费马小定理以及快速幂应用

艾醒(AiXing-w)

已于 2022-03-12 20:30:22 修改

阅读量482

点赞数

分类专栏： C++数据结构与算法文章标签：概率论线性代数几何学

于 2022-01-15 10:22:42 首次发布

本文链接：https://blog.csdn.net/DuLNode/article/details/122505903

版权

C++数据结构与算法专栏收录该内容

22 篇文章 2 订阅

订阅专栏

费马小定理

定理
- 例子
- 代码
逆元
- 快速幂代码
- 快速幂求逆元

定理

假设p为质数，a不是p的倍数，则有
$a^{p-1}\equiv 1(mod \space p)$

例子

求
$2^{100}关于13的余数$
$2^{100}=2^{12 \times 8} \cdot2^4=16 (mod \space 13)=3 (mod \space 13)$

代码

$假设我们要求2^k \space mod \space p， k \in [1, 100]$

#include<iostream>
#include<cmath>

using namespace std;

int fun1(int a, int p, int mod){
	int ans = pow(a, p % (mod-1));
	ans %= mod;
	return ans;
}

int fun2(int a, int p, int mod){
	int ans = pow(a, p);
	ans %= mod;
	return ans;
}

int fun3(int a, int p, int mod){
	int ans = 1;
	while(p--){
		ans *= a;
		ans %= mod;
	}
	ans %= mod;
	return ans;
}

int main(){
	for(int i=1; i <= 100; i++){
		cout<<fun1(2, i, 13)<<" "<<fun1(2, i, 13)<<" "<<fun1(2, i, 13)<<endl;
	}
	return 0;
}

fun2和fun3都会受到时间或者空间的限制，而fun1直接使用了费马小定理可以迅速的求取答案

逆元

$如果a \times b = 1 (mod \space p)，则a为b互为mod \space p逆元$
根据费马小定理得
$a^{p-1} \equiv 1 (mod \space p) \\\implies a \times a^{p-2} \equiv 1 (mod \space p) \\ 所以a和a^{p-2}互为mod \space p逆元$
费马小定理一般和快速幂配合使用，下面我们介绍一下快速幂

快速幂代码

#include<iostream>
#include<cmath>

using namespace std;

int normal_pow(int a, int p){
	int ans = pow(a, p);
	return ans;
}

int qpow_1(int a, int p){
	if(p==0){
		return 1;
	}else if(p % 2 == 1){
		return qpow_1(a, p-1) * a;
	}else{
		int temp = qpow_1(a, p / 2);
		return temp * temp;
	}
}

int qpow_2(int a, int p){
	int ans = 1;
	while(p){
		if(p & 1){
			ans *= a;
		}
		a *= a;
		p >>= 1;
	}
	return ans;
}

int main(){
	for(int i=1; i <= 10; i++){
		cout<<normal_pow(2, i)<<" "<<qpow_1(2, i)<<" "<<qpow_2(2, i)<<endl;
	}
	return 0;
}