Lucas定理（大组合数取模）

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll f[10008];
const int p=10007;
void init()
{
	f[0]=1;
	for(int i=1;i<=10007;i++)
	f[i]=(f[i-1]*i)%p;
}
ll pow(ll a,int b)
{
	ll ans=1;
	while(b)
	{
		if(b&1)
		ans=ans*a%p;
		a=a*a%p;
		b>>=1;
	}
	return ans;
}
ll C(ll a,int b)
{
	if(b>a)
	return 0;
	return f[a]*(pow(f[b]*f[a-b],p-2))%p;
}
ll lucas(ll n,ll m)
{
	if(m==0)
	return 1;
	else return (C(n%p,m%p)*lucas(n/p,m/p))%p;
}

int main()
{
	ll n,m;
	
	while(~scanf("%lld %lld",&n,&m))
   {
		init();
	printf("%lld\n",lucas(m,n));
   }
}