组合数、容斥原理和博弈论

瑆鈅掖

于 2024-05-26 21:56:23 发布

阅读量1.1k

点赞数 9

CC 4.0 BY-SA版权

分类专栏：算法文章标签：算法 c++ 笔记

本文链接：https://blog.csdn.net/2302_77099705/article/details/139212967

求组合数I

问题描述

给定 n 组询问，每组询问给定两个整数 a，b，请你输出C^a_bmod(109+7)的值。

输入格式

第一行包含整数 n。

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

输出格式

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

数据范围

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

输入样例

输出样例

3
10
1

问题分析

在这里插入图片描述

代码

#include<bits/stdc++.h>
using namespace std;
const int N=2010,MOD=1e9+7;
int c[N][N];
void init()
{
    for(int i=0;i<N;i++)
        for(int j=0;j<=i;j++)
        {
            if(!j) c[i][j]=1;
            else c[i][j]=(c[i-1][j]+c[i-1][j-1])%MOD;
        }
}
int main()
{
    int n;
    cin>>n;
    init();
    while(n--)
    {
        int a,b;
        cin>>a>>b;
        cout<<c[a][b]<<endl;
    }
    return 0;
}

求组合数 II

问题描述

给定 n 组询问，每组询问给定两个整数 a，b，请你输出C^a_bmod(109+7)的值。

输入格式

第一行包含整数 n。

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

输出格式

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

数据范围

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

输入样例

输出样例

3
10
1

问题分析

在这里插入图片描述

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=100010,mod=1e9+7;
int fact[N],infact[N];
int qmi(int a,int b,int p)
{
    int res=1;
    while(b)
    {
        if(b&1) res=(ll)res*a%p;
        a=(ll)a*a%p;
        b>>=1;
    }
    return res;
}
int main()
{
    fact[0]=infact[0]=1;
    for(int i=1;i<=N;i++)
    {
        fact[i]=(ll)fact[i-1]*i%mod;
        infact[i]=(ll)infact[i-1]*qmi(i,mod-2,mod)%mod;
    }
    int n;
    cin>>n;
    while(n--)
    {
        int a,b;
        cin>>a>>b;
        cout<<(ll)fact[a]*infact[b]%mod*infact[a-b]%mod<<endl;
    }
    return 0;
}