乘法逆元

定义：
满足a*k≡1 (mod p)的k值就是a关于p的乘法逆元。

为什么要有乘法逆元呢？
当我们要求（a/b） mod p的值，且a很大，无法直接求得a/b的值时，我们就要用到乘法逆元。我们可以通过求b关于p的乘法逆元k，将a乘上k再模p，即(a*k) mod p=a mod p * k。其结果与(a/b) mod p等价。
eg:
40/4 mod 3 ?
4*k=1(mod 3)
k=1
=>40/4 mod 3相当于40*1%3=1

证：（其实很简单。。。）
根据b*k≡1 (mod p)有b*k=p*x+1。
k=(p*x+1)/b。
把k代入(a*k) mod p，得：
(a*(p*x+1)/b) mod p
=((a*p*x)/b+a/b) mod p
=[((a*p*x)/b) mod p +(a/b)] mod p
=[(p*(a*x)/b) mod p +(a/b)] mod p
//p*[(a*x)/b] mod p=0
所以原式等于：(a/b) mod p

证明转载：http://blog.csdn.net/liangzhaoyang1/article/details/51104333

#include<iostream>
#include<cstdio>
#include<cstring>
#define p 9973
using namespace std;
void ex_gcd(int a,int b,int &x,int &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return;
    }
    ex_gcd(b,a%b,x,y);
    int tmp =x;
    x=y;
    y=tmp-(a/b)*y;
}
int main()
{
    int t,n,b,x,y;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&b);
        ex_gcd(b,p,x,y);
        x=(x%p+p)%p;
        printf("%d\n",(n%p*x)%p);
    }

}