题目链接
费小马定理求逆元
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));
}
}