组合数选择
从五个数12345中选择三个恰好有十种方式,分别是:
123、124、125、134、135、145、234、235、245和345
在组合数学中,我们记作:5C3 = 10。
一般来说,
nCr=n!r!(n−r)!,其中r ≤ n,n! = n×(n−1)×…×3×2×1,且0! = 1。
直到n = 23时,才出现了超出一百万的组合数:23C10 = 1144066。
若数值相等形式不同也视为不同,对于1 ≤ n ≤ 100,有多少个组合数nCr超过一百万?
代码演示
#include <iostream>
using namespace std;
#define MAX_N 1000000
int prime[MAX_N + 5] = {0};
int f[MAX_N + 5] = {0};
int cnt[MAX_N + 5] = {0};
void init() {
for (int i = 2; i <= MAX_N; i++) {
if (!prime[i]) {
prime[++prime[0]] = i;
f[i] = i + 1;
cnt[i] = i * i;
}
for (int j = 1; j <= MAX_N; j++) {
if (prime[j] * i > MAX_N) break;
prime[prime[j] * i] = 1;
if (i % prime[j] == 0) {
f[i * prime[j]] = f[i] * (cnt[i] * prime[j] - 1) / (cnt[i] - 1);
cnt[i * prime[j]] = cnt[i] * prime[j];
break;
} else {
f[i * prime[j]] = f[i] * f[prime[j]];
cnt[i * prime[j]] = prime[j] * prime[j];
}
}
}
return ;
}
int main() {
init();
int ans = 0;
for (int i = 2; i < MAX_N; i++) {
if (f[i] > ans) ans = f[i];
}
cout << ans << endl;
return 0;
}