這種求組合數問題是要把最終答案算出來，而不是把最終答案模上一個數

就是直接用定義求組合數，分解質因數，然後用上高精度乘法

還是要吐槽一句python自帶高精，python大法好，之後還要轉python，酸了酸了

如何求a！中有多少個質因子p呢

一直乘到p的若干次方比a大了為止，這步操作時間複雜度是logp

思路：
1：先把1 ~ 5000內的質數篩出來

2：求每個質數的次數

3：用高精度乘法把所有質因子乘出來，算出結果

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 const int N = 5010;
 4 int primes[N], cnt;
 
 5 bool st[N];
 6 int sum[N];
 7 void get_primes(int n) {
 8     for (int i = 2; i <= n; i++) {
 9         if (!st[i]) {
10             primes[cnt++] = i;
11         }
12         for (int j = 0; primes[j] <= n / i; j++) {
13             st[primes[j] * i] = true;
14             if (i % primes[j] == 0
 
) {
15                 break;
16             }
17         }
18     }
19 }
20 int get(int n, int p) { //求n!包含的p的次數
21     int res = 0;
22     while (n) {
23         res += n / p;
24         n /= p;
25     }
26     return res;
27 }
28 vector<int> mul(vector<int> &a, int b) {
29     vector<int
 
> c;
30     int t = 0;
31     for (int i = 0; i < a.size(); i++) {
32         t += a[i] * b;
33         c.push_back(t % 10);
34         t /= 10;
35     }
36     while (t) {
37         c.push_back(t % 10);
38         t /= 10;
39     }
40     return c;
41 }
42 int main() {
43     int a, b;
44     cin >> a >> b;
45     get_primes(a); //1.篩選出所有質因子
46     for (int i = 0; i < cnt; i++) { //2.求出每一個質數的次數
47         int p = primes[i];
48         sum[i] = get(a, p) - get(b, p) - get(a - b, p);
49     }
50     //3.高精度乘法
51     vector<int> res;
52     res.push_back(1); //最開始是1
53     for (int i = 0; i < cnt; i++) { //列舉所有質因數
54         for (int j = 0; j < sum[i]; j++) { //這個質因數的次數
55             res = mul(res, primes[i]);
56         }
57     }
58     for (int i = res.size() - 1; i >= 0; i--) { //輸出
59         cout << res[i];
60     }
61     cout << endl;
62     return 0;
63 }