逆元

咱们首先来看个线性同余方程：

若是对于方程 ax = b（a不为0），因为a存在倒数，所以很容易求解。若是在mod m的运算下，也有知足这样a的倒数同样的数存在的话，方程就有解了。而这个解x就叫作a关于m的逆元，记作或是inv(a)。若是能求出逆元，那么就有x = inv(a) * ax = inv(a) * b, 就能够求出x了。

 那么咱们怎么求出inv(a)呢？

其实就是解

咱们设ax = mt + 1;

移项得ax - mt = 1；

咱们设y = -m，得ax + my = 1;

咦，这方程怎么那么眼熟，好像可用扩展欧几里得算法求得。 

固然对于方程知足的条件必须是gcd(a, m) = 1，不然的话逆元是不存在的。

附上伪代码：

int gcd(int a, int b){
    return !b ? gcd(b, a % b) : a;
}

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

int inv(int a, int m){
    int x, y;
    int d = extgcd(a, m, x, y); 
    if(gcd(a, m) == 1)return (m + x % m) % m;
    else return -1;//-1表示不存在逆元
}

固然，逆元还有其余的求法。在这以前咱们须要知道一些姿式。

费马小定理：

　　在p是素数的状况下，对任意整数x都有

　　

其中若是x不能被p整除则有：

继续变形：

咦，咱们发现就是x关于p的逆元。即：

所以就能够用矩阵快速幂运算来求出逆元。

在不是素数的状况下， 咱们也能够经过欧拉定理来求解。

欧拉定理：若a, n均为正整数，且a, n互质，则

而费马小定理仅仅是其一个特例而已。

固然, 咱们还能够用另外的方法。

咱们知道p % b =  p - (p / b) * b(这里的/表示整数除法ex : 7 / 2 = 3)

设x = p % b, y = p / b;

因而就有x + by = p;

咱们两边同时取余p得（x + by) % p = 0;

x % p = (-y) * b % p;

x * inv(b) % p = (-y) % p;

inv(b) = (-y) * inv(x) % p;

inv(b) = (p - y) * inv(x) % p;

将x, y代入得：

inv(b) = (p - p / b) * inv(p % b) % p;

附上伪代码：

const int MOD = (int)1e9 + 7;//按题目要求的取余数
const int N = 1000000 + 5;

int inv[N + 5];

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