求逆元的方法

逆元简介

如果存在ｘ使得a*x ≡ 1 mod ( p ),　我们就说x 是 a 的一个逆元。

拓展欧几里得算法求逆元

首先拓展欧几里得可以解决的问题是 ax + py = 1，这样通过拓展欧几里得算法求解出来的ｘ既为ａ的逆元。因为对左右同时取模p，就变回了ａｘ≡ 1 mod ( p )。得到的结果通过对ｘ取模p就能的到最小的正逆元，如果需要当然也可以求出来负数逆元。
代码如下：

void exgcd(int a, int b, int &x, int &y)
{
    if (!b)
    {
        x = -1, y = 0;
    }
    else
    {
        exgcd(b, a%b, y, x);
        y -= (a/b)*x;
    }
}

费马小定理解决求逆元

费马小定理(Fermat’s little theorem)是数论中的一个重要定理，在1636年提出。如果p是一个质数，而整数a不是p的倍数，则有a^（p-1）≡1（mod p）。

通过对a^（p-1）≡1（mod p）进行变换的到，a*a^p-2≡1（mod p），所以ａ的逆元是a^p-2，通过快速幂进行求解。但别忘记的时，ｐ必须是一个素数，在做题的时候给出的取模数字有时候可能不是素数，就不能用这种方法进行计算了。

#include <iostream>

using namespace std;
typedef long long ll;
const int mod = prime;

ll qul(ll a, ll b)
{
    ll sum = 1;
    
    while (b)
    {
        if (b&1)
            sum = sum*a%mod;
        
        a = a*a%mod;
        b>>=1;
    }
}

int main()
{
    ll ni = qul(ni, mod-2);
    
    return 0;
}

打表求逆元

首先给出证明。
　
代码如下：

#include <iostream>

using namespace std;

typedef long long ll;
const int MOD = 1e7+9;
ll inv[MOD];

void init()
{
	int[1] = 1;
	
    for (int i=2; i<MOD; i++)
    {
        inv[i] = (MOD - MOD/i)%MOD*ni[MOD%i]%MOD)%MOD;
    }
}

int main()
{
    init();
    
    return 0;
}