【题解】LightOJ1067 逆元

题目链接
费小马定理求逆元
a/b=1mod( M );
只要 M 是一个素数，而且 b 不是 M 的倍数，就可以用一个逆元整数 b’，通过 a/b=a*b’*(mod M)来以乘换除
费马小定理说，对于素数 M 任意不是 M 的倍数的 b，都有：b^(M-1)=1 (mod) M;
于是可以拆成：b*b^(M-2)=1(mod)M;
所以：a/b=a/b*(b*b^(M-2))=a*(b^(M-2))(mod M)
也就是说我们要求的逆元就是b^(M-2)(mod M)

#include<cstdio>
const int N=1e6+5;
const int mod=1000003;
typedef long long ll;
ll fact[N];
void init()
{
    fact[1]=fact[0]=1;
    for(ll i=2;i<N;i++)
    fact[i]=(i*fact[i-1])%mod;
}
ll niyuan(ll a,ll p)//快速幂 
{
    ll res=1%mod;
    ll t=a%mod;
    while(p)
    {
        if(p%2)
        res=res*t%mod;
        t=t*t%mod;
        p/=2;
    }
    return res;
}
ll work(int n,int k)
{
    ll fm=(fact[k]*fact[n-k])%mod;
    ll ans1=niyuan(fm,mod-2);
    return (ans1*fact[n])%mod;
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int t,n,k,ca=0;
    init();
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&k);
        if(k*2>n)
        k=n-k;//组合数对称性，减少计算
        printf("Case %d: %lld\n",++ca,work(n,k)); 
    }
}