阶乘分解[CH3101]

欢迎大家访问我的老师的OJ———caioj.cn

题面描述

传送门

思路

算术基本定理:
$$N=p_1^{c_1}{p}_2^{c_2}......p_m^{c_m}$$

题意让我们求$\mathrm{N}!$的每一个质因子$p_i$的个数$c_i$

$\mathrm{N}!$中质因子$p$的个数实际上就是$[1,n]$中每个数包含$p$的个数之和.

显然地,$[1,n]$中至少含有一个$p$的数有$\lfloor N/p\rfloor$个

至少含有2个$p$的数有$\lfloor N/p^2\rfloor$，有$2\ast \lfloor N/p^2\rfloor$个$p$，

但由于至少含有一个$p$包含至少含有2个$p$，所以

统计第2个质因子，即累加上$\lfloor N/p^2\rfloor$,而不是$2\ast \lfloor N/p^2\rfloor$个$p$。

综上所述，$\mathrm{N}!$中质因子$p$的个数为:
$$\lfloor \frac{N}{p}\rfloor +\lfloor \frac{N}{p^2}\rfloor +\lfloor \frac{N}{p^3}\rfloor +......+\lfloor \frac{N}{p^{\lfloor {\log }_{p}N\rfloor }}\rfloor =\sum _{p^k\le N}\lfloor \frac{N}{p^k}\rfloor$$

对于每个$p$,都这样进行即可。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
#define ll long long
using namespace std;
const int N=1e6+10;
const int inf=1e6;
const int M=5e5+10;
int prime[M],m,n;bool v[N];
void g_p()
{
	m=0;memset(v,false,sizeof(v));
	for(int i=2;i<=n;i++)
	{
		if(!v[i])
		{
			prime[++m]=i;
			ll p=prime[m];int c=0;
			while(p<=n)c+=n/p,p*=prime[m];
			printf("%d %d\n",prime[m],c);
		}
		for(int j=1;j<=m&&prime[j]*i<=n;j++)
		{
			v[i*prime[j]]=1;
			if(i%prime[j]==0)break;
		}
	}
}
int main()
{
	scanf("%d",&n);
	g_p();
	return 0;
}