卢卡斯定理(模板)

最新推荐文章于 2022-10-20 19:24:36 发布

痼

最新推荐文章于 2022-10-20 19:24:36 发布

阅读量371

点赞数

CC 4.0 BY-SA版权

分类专栏：数论文章标签：卢卡斯定理

本文链接：https://blog.youkuaiyun.com/qq_43628761/article/details/98229478

数论专栏收录该内容

5 篇文章

订阅专栏

卢卡斯定理是一个在组合数学中用于计算特定组合数值的公式，尤其适用于大数情况。本文将探讨该定理及其在解决编程模板题如Luogu P3807中的应用。

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

卢卡斯定理

卢卡斯定理主要用来求大数组合数。
模板题：https://www.luogu.org/problem/P3807

#include <stdio.h>
#include <time.h>
#include <string.h>
#include <algorithm>
#include <stack>
#include <vector>
#include <set>
#include <iostream>
#include <queue>
#include <string>
#include <math.h>
typedef long long LL;
using namespace std;
const int inf = 0x3f3f3f3f;
const int mod = 1e9 + 7;
#define mid(l,r) (( l + r ) / 2)
#define lowbit(x) (( x & ( - x )))
#define lc(root) ((root * 2))
#define rc(root) ((root * 2 + 1))
#define me(array,x) (memset( array , x , sizeof( array ) ) )
const int maxn = 200004;

LL POW(LL a,LL b,int p)
{
    LL ans = 1;
    while(b){
        if(b&1)
            ans = (ans%p*a%p) % p;
        a = (a%p*a%p) % p;
        b >>= 1;
    }
    return ans%p;
}

LL C(int n,int m,int p)
{
    if(m > n)return 0;
    if(m == n)return 1;
    if(m > n - m)m = n - m;
    LL a = 1, b = 1 ;
    for(LL 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(int n,int m,int p){
    if(m == 0)
        return 1;
    return C(n%p,m%p,p)*lucas(n/p,m/p,p)%p;
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--){
        int n,m,p;
        scanf("%d%d%d",&n,&m,&p);
        cout<<lucas(n+m,m,p)<<"\n";
    }
    return 0;
}