CH 3101 阶乘分解【质因数分解】

题目链接

题目描述：

给定整数 N(1 ≤ N ≤ $10^6$)，试把阶乘 N! 分解质因数，按照算术基本定理的形式输出分解结果中的 $p_i$ 和 $c_i$ 即可。

---

关于质因数分解：

算术基本原理：

任何一个大于 1 的正整数都能唯一分解成有限个质数的乘积，可写做：
$$N=p_1^{c1}{p}_2^{c2}...p_m^{cm}$$
其中 $c_i$ 都是正整数，$p_i$都是质数，且满足$p_1<p_2<...<p_m$。

试除法：

结合质数判定的“试除法”和质数筛选的“Eratosthenes筛法”，我们可以扫描 2 ~ $\lfloor N^{1/2}\rfloor$的每个数d，若d能整除N，则从N中除掉所有的因子d，同时累计除去d的个数。
因为一个合数的因子一定在扫描到这个合数之前就从N中被除掉了，所以在上述过程中能整除N的一定是质数，最终就到了质因数分解的结果，易知时间复杂度为$O(N^{1/2})$。

---

题解：

若把1 ~ N每个数都分解质因数，再把结果合并，时间复杂度过高，为$O(N{N}^{1/2})$，显然，N !的每个质因子都不会超过N，我们可以先筛选出1~N的每个质数p，然后考虑阶乘N !中一共包含多少个质因子p。
N !中质因子p的个数就等于1 ~ N每个数包含质因子p的个数，在1 ~ N中，p的倍数，即至少包含 1 个质因子 p 的显然有 $\lfloor N/p\rfloor$个，而$p^2$的倍数，即至少包含两个质因子p的有$\lfloor N/p^2\rfloor$个。不过其中的一个质因子已经在$\lfloor N/p\rfloor$里统计过了，所以只需要在统计第二个质因子，即累加上$\lfloor N/p^2\rfloor$，而不是2 * $\lfloor N/p^2\rfloor$。
综上所述：N !中质因子p的个数为：
$$\lfloor \frac{N}{p}\rfloor +\lfloor \frac{N}{p^2}\rfloor +\lfloor \frac{N}{p^3}\rfloor +...+\lfloor \frac{N}{p^{\lfloor lo{g}_{p}N\rfloor }}\rfloor =\sum _{p^k\le N}\lfloor \frac{N}{p^k}\rfloor$$
对于每个p，只需要$O(lo{g}_{2}N)$计算上述和式，故对于整个 N ! 分解质因数复杂度为$O(Nlo{g}_{2}N)$

---

AC代码

#include <iostream>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <unordered_set>
#include <unordered_map>
#include <deque>
#include <list>
#include <iomanip>
#include <algorithm>
#include <fstream>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <cstdlib>
//#pragma GCC optimize(2)
using namespace std;
typedef long long ll;
//cout << fixed << setprecision(2);
//cout << setw(2);
const int N = 1e6 + 6, M = 1e9 + 7;

bool vis[N];
int prime[200000], c[200000];

int main() {
    //freopen("/Users/xumingfei/Desktop/ACM/test.txt", "r", stdin);
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int n, cnt = 0;
    cin >> n;
    for (int i = 2; i <= n; i++) {
        if (vis[i] == 0) prime[cnt++] = i;
        for (int j = 0; j < cnt && prime[j] * i <= n; j++) {
            vis[prime[j] * i] = 1;
            if (i % prime[j] == 0) break;
        }
    }
    for (int i = 0; i < cnt; i++) {
        for (long long j = prime[i]; j <= n; j *= prime[i])
            c[i] += n / j;
    }
    for (int i = 0; i < cnt; i++) {
        cout << prime[i] << " " << c[i] << '\n';
    }
    return 0;
}