C(5, 2) 是多少
运行下面这段程序:
本节模型:comb(n, k, p) 分子分母各自连乘,分母乘费马逆元。
#include <array>
#include <iostream>
#include <string>
#include <vector>
using namespace std;
long long inv(long long a, long long p) {
// 费马小定理求逆元:p 是质数时 a^(p-2) 就是 a 的逆元
long long r = 1;
a %= p;
long long e = p - 2;
while (e > 0) {
if (e & 1) r = r * a % p;
a = a * a % p;
e >>= 1;
}
return r;
}
long long comb(long long n, long long k, long long p) {
// 组合数 C(n, k) mod p(p 是质数):分子分母各自连乘,分母乘费马逆元
if (k < 0 || k > n) return 0;
long long num = 1;
for (long long i = n - k + 1; i <= n; i++) num = num * i % p;
long long den = 1;
for (long long i = 1; i <= k; i++) den = den * i % p;
return num * inv(den, p) % p;
}
int main() {
cout << comb(5, 2, 1000000007) << endl;
}
(本题用 g++ -std=c++17 -O0 编译。)
全部评论