费马小定理

费马小定理

费马小定理(Fermat's little theorem)是数论中的一个重要定理，在1636年提出，其内容为： 假如p是质数，且gcd(a,p)=1，那么 a(p-1)≡1（mod p），例如：假如a是整数，p是质数，则a,p显然互质(即两者只有一个公约数1)，那么我们可以得到费马小定理的一个特例，即当p为质数时候， a^(p-1)≡1(mod p)

（≡   的意思是左右二边的数对p求模所得到的数是一样的；

证明:（源自数论妙趣——数学女王的盛情款待）

　任意取一个质数，比如13。考虑从1到12的一系列整数1,2,3,4,5,6,7,8,9,10,11,12，给这些数都乘上一个与13互质的数，比如3，得到3,6,9,12,15,18,21,24,27,30,33,36。对于模13来说，这些数同余于3,6,9,12,2,5,8,11,1,4,7,10。这些余数实际上就是原来的1,2,3,4,5,6,7,8,9,10,11,12，只是顺序不同而已。

　　把1,2,3,„,12统统乘起来，乘积就是12的阶乘12！。把3,6,9,„,36也统统乘起来，并且提出公因子3，乘积就是3¹²×12！。对于模13来说，这两个乘积都同余于1,2,3,„,12系列，尽管顺序不是一一对应，即3¹²×12！≡12！mod 13。两边同时除以12！得3¹²≡1 mod 13。如果用p代替13，用x代替3，就得到费马小定理

x^p－1≡1 mod p。

根据费马小定理可以求数的逆元，采用快速幂

代码：

#include<iostream>
#define LL long long
using namespace std;

LL quickm(LL a, LL b, LL p){
    LL ret = 1;
    while(b){
        if(b & 1) ret = (ret * a) % p;
        a = (a * a) % p;
        b >>= 1;
    }
    return ret;
}

LL Fermat(LL a, LL p){
        return quickm(a, p-2, p);
}

int main()
{
	LL ans,n;
	cin>>ans>>n;
	cout<<Fermat(ans,n)<<endl;
	return 0;
}