今日学习在线编程题：阶乘的质因子

题目来源：码蹄集

https://matiji.net/exam/brushquestion/14/3181/1DC60EA6DF83A333301CFFE1407FBA59

时间限制：1000ms
内存限制：65535kb

题目描述：小码哥有一天想把20！表示出来，但是小码哥意识到这个数字算出来将非常大，于是你给小码哥出主意，用算术基本定理的形式表示，于是这个任务理所应当的被小码哥推给了你，正好你想干点简单的事来换换脑子，就接受了。你要把N！按照算术基本定理的形式表示。

基本算数定理：一个数一定可以用一下形式表示：

x = p[1]^{c[1]} * p[2]^{c[2]} * ······ * p[i]^{c[i]} 

p[i]为质数数组，c[i]为每一个质数的幂的次数。

对于不同的x来说，有不同的p和c数组。

输入格式：输入一个整数N。​​​​​​​
输出格式：N! 分解质因数后的结果，共若干行，每行两个数p,c,表示含有 p^c 项。按照p[i]从小到大的顺序输出
 

输入样例：5
输出样例：2 3
                  3 1
                  5 1

参考程序：

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <queue>
using namespace std;
#define ll long long 
#define N 1100000
ll b[N], z[N] ,num = 0;
void zhi(){
    for (ll i = 2; i <= N ;i++){
        if (b[i] == 0){
            z[++num] = i;
            for (ll j = 2; j * i <= N; j++){
                b[i * j] = 1;
            }
        }
    }
}
int main(){
    ll n;
    cin >> n;
    zhi();
    if (n == 1) cout << 2 <<" " << 0 << endl;
    else 
    for (ll i = 1; z[i] <= n; i++){
        ll k = z[i] ;
        ll ans = 0;
        for (ll j = 1; k <= n; j++){
            ans = ans + n / k;
            k = k * z[i] ;
        }
        cout << z[i] <<" " << ans << endl;
    }
    return 0;
}