上海计算机学会 2023年9月月赛乙组T4 组合数（组合数学）

信息学奥赛-王老师

于 2024-04-13 07:00:00 发布

阅读量822

点赞数 7

分类专栏：题目题解文章标签：算法 c++ 数据结构图论学习

本文链接：https://blog.csdn.net/oitutor/article/details/137497415

版权

题目题解专栏收录该内容

96 篇文章

订阅专栏

本文介绍了如何使用组合数学中的C(n,m)公式求解从n个不同数字中取出m个的方案数，特别强调了对1,000,000,007取模以及利用费马小定理计算逆元的重要性，提供了C(n,m)的编程实现

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

第四题：T4组合数

标签：组合数学
题意：求组合数 $C_n^m$ ，即从 $n$ 个不同的数字中取出 $m$ 个数字的方案数，结果对 $1, 000, 000, 007$ 取模
（ $1≤m≤n≤10^9，1≤m≤10^6$ ）
题解：直接通过通项公式 $C_n^m=\frac {n!}{m!(n-m)!}$ ，值得注意的是要进行取模，所以得用到逆元 $in v$ 。
费马小定理： $a^{p-1}≡1(mod \ p)$ ，两边同时除以 $a$ ， $a^{p-2}≡inv(a)(mod\ p)$ ，即 $inv(a)=a^{p-2}(mod \ p)$ ，得到 $a^{p-2}$ 这部分可以通过快速幂去得到。
代码：

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
const ll p = 1e9 + 7;

ll fast_power(ll a, ll b) {
    ll ans = 1;
    a %= p;
    while (b) {
        if (b & 1) ans = (ans * a) % p;
        a = (a * a) % p;
        b >>= 1;
    }
    return ans;
}

ll C(ll a, ll b) {
    if (b > a) return 0;
    if (a == b) return 1;
    ll ans1 = 1, ans2 = 1;
    for (ll i = 1; i <= b; i++) {
        ans1 = ans1 * (a - i + 1) % p;
        ans2 = ans2 * i % p;
    }
    return ans1 * fast_power(ans2, p - 2) % p;
}

int main() {
    ll n, m;
    cin >> n >> m;
    cout << C(n, m);
    return 0;
}