Lucas定理

40 篇文章 0 订阅

37. Lucas定理

1. Definition

​ Let n = s p + q , m = t p + r . ( 0 ≤ , r ≤ p − 1 ) n=sp+q,m=tp+r.(0\leq,r\le p-1) n=sp+q,m=tp+r.(0,rp1)

Then : C n m ≡ C s t C q r ( m o d   p ) C_n^m\equiv C_s^tC_q^r(mod\ p) CnmCstCqr(mod p)

2. deduction

百度百科

3.Example

Link

Solution:

#include<bits/stdc++.h>
using namespace std;
#define maxn 100005
#define ll long long
ll nn; 
ll qpow(ll a,ll b,ll p)
{
    if(b==1) return a;
    ll s=1;
    while(b!=0)
    {
        if(b&1) 
        s=s*a%p;
        a=(a*a)%p;
        b>>=1;
    }
    return s%p;
}
ll pow(ll y,int z,int p){
    y%=p;ll ans=1;
    for(int i=z;i;i>>=1,y=y*y%p)if(i&1)ans=ans*y%p;
    return ans;
}
ll C(ll n,ll m,ll p)
{
    if(n<m) return 0;
    if(m>(n-m)) m=n-m;
    ll a=1;ll b=1;
    for(ll i=1;i<=m;i++)
    {
        a=(a*(i))%p;
        b=(b*(n-i+1))%p;
    }
    return (b*pow(a,p-2,p))%p;
}
ll C1(ll n,ll m,ll p){
	if(n<m) return 0;
	if(m>n-m) m=n-m;
	ll a=1,b=1;
	for(int i=0;i<m;i++){
		a=(a*(n-i))%p;
		b=(b*(i+1))%p;
	}
	return a*pow(b,p-2,p)%p;
}
ll lucas(ll n,ll m,ll p)
{
    if(m==0) return 1;
    return lucas(n/p,m/p,p)*C(n%p,m%p,p)%p;
}

int main()
{
    cin>>nn;
    //cout<<qpow(2,10,maxn);
    for(int i=1;i<=nn;i++)
    {
        int x,y,p;
        cin>>x>>y>>p;
        cout<<lucas(x+y,x,p)%p<<endl;
    }
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值