卢卡斯定理

最新推荐文章于 2022-05-21 15:35:21 发布

hahaAll

最新推荐文章于 2022-05-21 15:35:21 发布

阅读量1.3k

点赞数

本文链接：https://blog.csdn.net/hahaAll/article/details/49079241

版权

遇到好多次了，，，但总是写不出来。

来练几个最基础的题吧：http://acm.hdu.edu.cn/showproblem.php?pid=3037

http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1118

卢卡斯定理的基本应用：基本的组合数学，C(N+M,M)%P

卢卡斯定理模板：（hdu3037）

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <string>

using namespace std;

typedef long long LL;
#define M 1e9 + 7

typedef struct Node
{
    LL exp_mod(LL a, LL b, LL p)
    {
        LL res = 1;
        while(b != 0)
        {
            if(b&1) res = (res * a) % p;
            a = (a*a) % p;
            b >>= 1;
        }
        return res;
    }
    LL Comb(LL a, LL b, LL p)
    {
        if(a < b)   return 0;
        if(a == b)  return 1;
        if(b > a - b)   b = a - b;

        LL ans = 1, ca = 1, cb = 1;
        for(LL i = 0; i < b; ++i)
        {
            ca = (ca * (a - i))%p;
            cb = (cb * (b - i))%p;
        }
        ans = (ca*exp_mod(cb, p - 2, p)) % p;
        return ans;
    }
    LL lucas(int n, int m, int p)
    {
        LL ans = 1;
        while(n&&m&&ans)
        {
            ans = (ans*Comb(n%p, m%p, p)) % p;
            n /= p;
            m /= p;
        }
        return ans;
    }
} Lucas;
int main()
{
    int t;
    Lucas Lu;
    scanf("%d",&t);
    while(t --)
    {
        int n,m,p;
        scanf("%d%d%d",&n,&m,&p);
        printf("%lld\n",Lu.lucas(n+m,m,p));
    }
    return 0;
}