BZOJ【2982】combination

版权声明:转载记得声明噢! https://blog.csdn.net/Richard__Luan/article/details/82693573

Lucas 就当板子了
沙雕n==0
Bzoj过了loj无故re。。
Code

#include<cstdio>
#include<iostream>
using namespace std;
const long long mod=10007;
int t;
long long jc[20010];
int n,m;
long long njc[20010];
long long exgcd(long long a,long long b,long long &x,long long &y)
{
	if(b==0)
	{
		x=1,y=0;
		return a;
	}
	long long t=exgcd(b,a%b,y,x);
	y-=a/b*x;
	return t;
}
long long C(int n,int m)
{
	if(m==0||n==m)return 1;
	return jc[n]*njc[m]%mod*njc[n-m]%mod;
}
long long L(int n,int m)
{
	if(n<mod&&m<mod)return C(n,m);
	else return L(n/mod,m/mod)*L(n%mod,m%mod)%mod;
}
int main()
{
	jc[0]=1;
	for(int i=1;i<=10006;i++)
	{
		jc[i]=i*jc[i-1]%mod;
	}
	long long y;
	exgcd(jc[10006],mod,njc[10006],y);
	njc[10006]=(njc[10006]%mod+mod)%mod;
	for(int i=10005;i>=1;i--)
	{
		njc[i]=njc[i+1]*(i+1)%mod;
	}
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d%d",&n,&m);
		printf("%lld\n",L(n,m)%mod);
	}
	return 0;
}

猜你喜欢

转载自blog.csdn.net/Richard__Luan/article/details/82693573