阶乘分解质因数[经典题组合数学枚举质因子]

引入问题：

---

给定整数$N$，试把阶乘 $N!$ 分解质因数，按照算术基本定理的形式输出分解结果中的 $p_i$ 和 $c_i$即可。
$N!$分解质因数后的结果，共若干行，每行一对$p_i$,$c_i$，表示含有$p_i^{ci}$项。按照$p_i$从小到大的顺序输出。
$N\in [1,1e6]$

---

解题思路：

---

错误：
1)很明显乘出来再分解是很不现实的。
2)如果一个从$[1,n]$去枚举一个个去分解最坏下数量级$n\sqrt{n}\approx 1e9[n=1e6]$妥妥超时
正解：
1)我们可以枚举质因子$p$看$[1,n]$里面有多少个是p的倍数$\lfloor \frac{n}{p}\rfloor$
2)假如说有多个$p$呢？那么我们就继续看看有多少个数是$p^2$的倍数$\lfloor \frac{n}{p^2}\rfloor$
3)同理一直直到某个$x$使得$p^x>n$
4)那么$c=\lfloor \frac{n}{p}\rfloor +\lfloor \frac{n}{p^2}\rfloor +....+\lfloor \frac{n}{p^x}\rfloor$

---

代码

---

题目链接qwq戳我

#include <bits/stdc++.h>

using namespace std;
const int N = 1e6 + 10;
bool vis[N];
int prime[N], cnt;
int n;
void init(int n)
{
    for(int i = 2; i <= n; ++ i)
    {
        if(!vis[i]) prime[cnt ++] = i;
        for(int j = 0; j < cnt && i * prime[j] <= n; ++ j)
        {
            vis[i * prime[j]] = true;
            if(i % prime[j] == 0) break;
        }
    }
}
int main()
{
    init(N);
    cin >> n;
    for(int i = 0; i < cnt; ++ i)
    {
        int p = prime[i];
        int s = 0;
        for(int j = n; j; j /= p) s += j / p;
        if(s) cout << p << " " << s << endl;
    }
    return 0;
}