现在才看Lucas,真是菜QAQ。
Lucas首先是解决大组合数取模问题的,C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)%p,循环求解这个公式就是Lucas。
#include <iostream>
#include <cmath>
#include <vector>
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;
typedef long long int ll;
const int maxn=1e6+5;
const int mod=1000003;
int t;
ll n,m;
ll f[maxn];
void fac(ll x)//处理出0-mod内的阶乘
{
f[0]=1;//一定要处理0,0的阶乘是1,容易忘
for(ll i=1;i<=x;i++)
f[i]=(f[i-1]*i)%x;
}
ll quick_pow(ll a,ll b)//快速幂
{
ll res=1;
while(b>0)
{
if(b%2) res=(res*a)%mod;
a=(a*a)%mod;
b/=2;
}
return res;
}
ll Lucas(ll n,ll k)//Lucas
{
ll res=1;
while(n && k)
{
ll nn=n%mod,kk=k%mod;
if(nn<kk) return 0;
res=(res*f[nn]*quick_pow(f[kk]*f[nn-kk]%mod,mod-2))%mod;//注意最后取模
n/=mod;
k/=mod;
}
return res;
}
int main()
{
fac(mod);
scanf("%d",&t);
for(int i=1;i<=t;i++)
{
scanf("%lld%lld",&n,&m);
ll ans=Lucas(n,m,mod);
printf("Case %d: %lld\n",i,ans);
}
return 0;
}
C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)%p