【模板】Lucas定理

最新推荐文章于 2024-08-01 16:31:59 发布

QAQQQQQQQQQQQ

最新推荐文章于 2024-08-01 16:31:59 发布

阅读量161

点赞数

分类专栏： ——数论、数学—— 文章标签： Lucas定理

本文链接：https://blog.csdn.net/sizeof_you/article/details/80977821

版权

——数论、数学—— 专栏收录该内容

17 篇文章 0 订阅

订阅专栏

Lucas定理是用来求 c(n,m) mod p，p为素数的值。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define maxn 100005
#define ll long long
using namespace std;
int n,m,p,t;

ll qpow(ll x,int k){
  ll ret=1;
  while(k){
    if(k&1) ret=(ret*x)%p;
    x=(x*x)%p;
    k>>=1;
  }
  return ret;
}

ll C(int x,int y){
  ll ret=1;
  if(x<y) return 0;//这玩意坑害了我好久···一定要加QAQ
  for(int i=y+1;i<=x;i++) ret=(ret*1ll*i)%p;
  for(int i=2;i<=x-y;i++) ret=(ret*qpow(i,p-2))%p;
  // for(int i=2;i<=y;i++) ret=(ret*qpow(i,p-2))%p;
  return ret;
}

ll Lucas(int n,int m){
  if(n<p && m<p) return C(n,m)%p;
  return (C(n%p,m%p)*Lucas(n/p,m/p))%p;
}

int main(){
  scanf("%d",&t);
  while(t--){
    scanf("%d%d%d",&n,&m,&p);
    n=n+m;
    printf("%lld\n",Lucas(n,m));
  }
  return 0;
}