关于组合数与Lucas定理

这是属于数论的一部分，找了代码及题目。刷了几道题目之后发现其实大部分也都是在模版上做修改。

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

正常用组合数做必然爆炸。

因此

先给出两种组合数的写法。

第一种中规中矩的算，变算变除。减少了溢出的可能性，但数据一大还是要溢出的。

long long Combination(int n,int m)//组合数
{
long long nr = 1;
if( n-m < m ) m = n-m; // C(n, m) = C(n, n-m)
    for (int j = 0; j < m; j++)
    {
        nr *= (n - j);
        nr/=(j+1);
    }
    return nr;
}

第二种也是类似的，不过使用了dp写法。可用作计算二项式系数。

long long binomial_coefficient(int n, int m)//计算二项式系数，dp写法
{
int i, j;
if( n-m < m ) m = n-m; // C(n, m) = C(n, n-m)
    for (i = 0; i <= n; i++)
        bc[i][0] = 1;

    for (j = 0; j <= n; j++)
        bc[j][j] = 1;

    for (i = 1; i <= n; i++)
        for (j = 1; j < i; j++)
            bc[i][j] = bc[i - 1][j - 1] + bc[i - 1][j];

    return bc[n][m];
}

但是当n，m再大，再怎么去算都要爆炸。

这就用到了Lucas定理。简单来说就是一条公式

C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)

然后代码实现。

还有一个最关键的问题是涉及到模运算时，组合数的下方不应该是直接除去某些数字，因此需要乘以这个数的逆元

#include<bits/stdc++.h>
#define ll long long
using namespace std;
int n,m;
const int p=10007;
int qpow(int a, int b)
{
    int ans;
    for(ans = 1; b; b >>= 1, a = a * a % p)
        if(b & 1)
            ans = ans * a % p;
    return ans;
}
int getc(int n, int m)
{
    if(n < m)
        return 0;
    if(m > n - m)
        m = n - m;
    ll s1 = 1, s2 = 1;
    for(int i = 0; i < m; i++)
    {
        s1 = s1 * (n - i) % p;
        s2 = s2 * (i + 1) % p;
    }
    return s1 * qpow(s2, p - 2) % p;
}
int lucas(int n, int m)
{
    if(m == 0)
        return 1;
    return getc(n % p, m % p) * lucas(n / p, m / p) % p;
}
int main() 
{
    scanf("%d%d", &n, &m);
    printf("%d\n", getc(n, m));
    printf("%d\n", lucas(n, m));
    return 0;
}

再附上另一个模版

大同小异

#include <iostream>
#include <cstdio>
#include <cstring>
usingnamespace std;

typedef longlong llg;

constint N =150000;

llg n, m, p, fac[N];

void init()
{
    int i;
    fac[0] =1;
    for(i =1; i <= p; i++)
        fac[i] = fac[i-1]*i % p;
}

llg pow(llg a, llg b)
{
    llg tmp = a % p, ans =1;
    while(b)
    {
        if(b &1)  ans = ans * tmp % p;
        tmp = tmp*tmp % p;
        b >>=1;
    }
    return  ans;
}

llg C(llg n, llg m)
{
    if(m > n)  return0;
    return  fac[n]*pow(fac[m]*fac[n-m], p-2) % p;
}

llg Lucas(llg n, llg m)
{
    if(m ==0)  return1;
    elsereturn  (C(n%p, m%p)*Lucas(n/p, m/p))%p;
}

int main()
{
    int t;
    scanf("%d", &t);
    while(t--)
    {
        scanf("%I64d%I64d%I64d", &n, &m, &p);
        init();
        printf("%I64d\n", Lucas(n+m, n));
    }
    return0;
}

另外还有Lucas定理的延伸，迄今为止只做过一道关于这个的题目。

用作求二项式系数中奇数的个数

int main(){
    while(scanf("%d",&n) != EOF){
        sum = 0;
        while(n){
            sum += (n & 1);
            n >>= 1;
        }
        printf("%d\n", (1 << sum));
    }
	return 0;
}