BZOJ【2982】combination

最新推荐文章于 2020-12-21 00:52:56 发布

Richard__Luan

最新推荐文章于 2020-12-21 00:52:56 发布

阅读量204

点赞数 1

分类专栏：板子数论文章标签： Lucas

本文链接：https://blog.csdn.net/Richard__Luan/article/details/82693573

版权

板子同时被 2 个专栏收录

37 篇文章

订阅专栏

数论

10 篇文章

订阅专栏

本文介绍了一种基于中国剩余定理和递归分解的方法来高效计算Lucas数，通过预处理阶乘及其逆元，避免了大数运算中的除法操作，实现了对Lucas数的快速求解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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;
}