组合数模板(四种情况)

组合数的性质

从n个不同元素中，任取m(m≤n)个元素并成一组，叫做从n个不同元素中取出m个元素的一个组合;从n个不同元素中取出m(m≤n)个元素的所有组合的个数，叫做从n个不同元素中取出m个元素的组合数。C(n,m)=>C_n^m

1. C(n,m)=n! / ( (n-m)! * m! )         (m≤n)
2. C(n,m)= C(n,n-m)
3. C(n,m)=C(n-1,m-1)+C(n-1,m)
4. C(n,m)%p=C(n/p,m/p) * C(n%p,m%p)%p卢卡斯定理对较大的组合数对素数p取模

求组合数 1

题目链接
给定n组询问，每组询问给定两个整数a，b，请你输出C(a,b) mod (10⁹+7)的值。

输入格式

第一行包含整数n。

接下来n行，每行包含一组a和b。

输出格式

共n行，每行输出一个询问的解。

数据范围

1≤n≤10000,
1≤b≤a≤2000

输入样例：

3
3 1
5 3
2 2

输出样例：

3
10
1

算法分析

数目比较少，而且a,b的值也比较小我们可以用递推的方法,利用性质3. C(a,b)=C(a-1,b-1)+C(a-1,b).

代码实现

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
int n;
const int maxn=2e3+5;
int C[maxn][maxn];
const int mod=1e9+7;
void init()
{
    for(int i=0;i<maxn;i++)
    {
        for(int j=0;j<=i;j++)
        {
            if(j==0)
                C[i][j]=1;
            else
                C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
        }
    }
}
int main()
{
    int t;
    cin>>t;
    init();
    int a,b;
    while(t--)
    {
        cin>>a>>b;
        cout<<C[a][b]<<endl;
    }
}

求组合数 2

题目链接
给定n组询问，每组询问给定两个整数a，b，请你输出C(a,b) mod (10⁹+7)的值。

输入格式

第一行包含整数n。
接下来n行，每行包含一组a和b。

输出格式

共n行，每行输出一个询问的解。

数据范围

1≤n≤10000,
1≤b≤a≤10⁵

输入样例：

3
3 1
5 3
2 2

输出样例：

3
10
1

算法分析

数目比较少，但是a和b比较大,但是一个阶乘就超过int范围了,所以我们要用到long long还有逆元的定理,边乘边除,利用1. C(n,m)=n! / ( (n-m)! * m! )

代码实现

#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=1e5+5;
const int mod=1e9+7;
typedef long long ll;
ll fac[maxn],infac[maxn];
//一个存乘以的阶乘,一个存除以的阶乘
ll ksm(ll a,ll b,ll p)
{
    ll res=1;
    while(b)
    {
        if(b&1) res=res*a%p;
        a=a*a%p;
        b>>=1;
    }
    return res;
}
void init()
{
    fac[0]=infac[0]=1;
    for(int i=1;i<maxn;i++)
    {
        fac[i]=(ll)fac[i-1]*i%mod;
        infac[i]=(ll)infac[i-1]*ksm(i,mod-2,mod)%mod;
    }
}
int main()
{
    init();
    ll t;
    cin>>t;
    while(t--)
    {
        ll a,b;
        cin>>a>>b;
        cout<<fac[a]*infac[b]%mod*infac[a-b]%mod<<endl;//性质1
    }
}

求组合数 3

题目链接
给定n组询问，每组询问给定三个整数a,b,p，其中p是质数，请你输出C ^b_a mod p的值。

输入格式

第一行包含整数n。
接下来n行，每行包含一组a,b,p。

输出格式

共n行，每行输出一个询问的解。

数据范围

1≤n≤20,
1≤b≤a≤10¹⁸,
1≤p≤10⁵,

输入样例：

3
5 3 7
3 1 5
6 4 13

输出样例：

3
3
2

算法分析

组合数比较大,还有要取余,需要利用组合数的卢卡斯定理
C(n,m)%p=C(n/p,m/p) * C(n%p,m%p)%p卢卡斯定理对较大的组合数对素数p取模

代码实现

#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
ll n,a,b,p;
ll ksm(ll a,ll b,ll p)
{
    ll res=1;
    while(b)
    {
        if(b&1) res=res*a%p;
        a=a*a%p;
        b>>=1;
    }
    return res;
}
ll C(ll a,ll b,ll p)
{
    if(b>a) return 0;
    ll x=1,y=1;
    for(int i=1,j=a;i<=b;i++,j--)
    {
        x=x*j%p;
        y=y*i%p;
    }
    return x*ksm(y,p-2,p)%p;
}
ll lucas(ll a,ll b,ll p)
{
    if(a<p&&b<p)    return C(a,b,p);
    else
        return C(a%p,b%p,p)*lucas(a/p,b/p,p)%p;
}
int main()
{
    cin>>n;
    while(n--)
    {
        cin>>a>>b>>p;
        cout<<lucas(a,b,p)<<endl;
    }
    return 0;
}

求组合数 4

题目链接
输入a,b，求C^b_a的值。
注意结果可能很大，需要使用高精度计算。

输入格式

共一行，包含两个整数a和b。

输出格式

共一行，输出C^b_a的值。

数据范围

1≤b≤a≤5000

输入样例：

5 3

输出样例：

10

算法分析

单纯的考察高精度的算法,我们可以用高精度的原理是,我们组合数求解中除法的上下等式会包含相同的质因子,可以将其变为最简然后计算

代码实现

#include<iostream>
#include<algorithm>
#include<vector>
using namespace std;
const int maxn=5e3+5;
int pri[maxn],cnt;
bool nopri[maxn];
int sum[maxn];
void getpri(int n)//分解质因子,存入pri中
{
    for(int i=2;i<=n;i++)
    {
        if(!nopri[i])   pri[cnt++]=i;
        for(int j=0;pri[j]<=n/i;j++)
        {
            nopri[pri[j]*i]=true;
            if(i%pri[j]==0) break;
        }
    }
}
int get(int a,int p)//获得a数中有几个p
{
    int res=0;
    while(a)
    {
        res+=a/p;
        a/=p;
    }
    return res;
}
vector<int> mul(vector<int> a,int b)//高精度乘法
{
    vector<int> ans;
    int t=0;
    for(int i=0;i<a.size();i++)
    {
        t+=a[i]*b;
        ans.push_back(t%10);
        t/=10;
    }
    while(t)
    {
        ans.push_back(t%10);
        t/=10;
    }
    return ans;
}
int main()
{
    int a,b;
    cin>>a>>b;
    getpri(a);
    for(int i=0;i<cnt;i++)
    {
        int p=pri[i];
        sum[i]=get(a,p)-get(b,p)-get(a-b,p);//最简后有几个质因子p;
    }
    vector<int> res;
    res.push_back(1);
    for(int i=0;i<cnt;i++)
        for(int j=0;j<sum[i];j++)
            res=mul(res,pri[i]);
    for(int i=res.size()-1;i>=0;i--)    cout<<res[i];
    return 0;
}