这道题,我们没办法求出N!,因为N范围是1e6,1e6的阶乘太大了,肯定超long long,所以我们调整枚举顺序,不去枚举N!中每个数进行质数分解,我们先求出N范围内的所有质数,然后枚举每一个质数,在N阶乘中的幂,关键在于这个在N阶乘中幂怎么求出来。
#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define ll long long
#define x first
#define y second
typedef pair<int,int> PII;
const int N=1e6+10;
int n,m,k,ecnt,cnt,ans;
int vis[N],prim[N],st[N];
void get_prime(int n) {//线性筛
memset(vis,0,sizeof vis);
for(int i=2;i<=n;i++){
if(!vis[i])prim[++cnt]=i;
for(int j=1;1ll*i*prim[j]<=n;j++){//越界中断
vis[i*prim[j]]=1;
if(i%prim[j]==0)break;
}
}
}
int main() {
cin>>n;
get_prime(n);
for(int i=1;i<=cnt;i++){//枚举质数
int p=prim[i];
int s=0;
for(int j=n;j;j/=p) s+=j/p;//累加和
printf("%d %d\n",p,s);
}
return 0;
}