n的阶乘的素数因子

n的阶乘的素数因子

题目：

输入正整数n（2《=n《=100），那n！=1*2*3*****n，分解成素因子相乘的形势，从小到大输出各个素数的指数，

例如825=3* 5**2 * 11应该表示成（0，1，2，0，1），分别表示有0个2，1个3，2个5，0个7，1个11。

忽略比最大素因子更大的素数（否则末尾会有无穷多个0）

样例

输入：53

输出：53！ = 49 23 12 8 4 4 3 2 2 1 1 1 1 1 1 1

解题思路：

构造prime数组存放用到的素数。创建指数因子数组p。

我们只需要判断 n！%素数因子 的情况,如果余数为零表示可以整除，那么对应的素因子的指数加一。

根据数学知识，我们可以从1到n依次判断对素数取余的情况，所以我们让m（1<=m<=n)不断地对素数因子取余（从最小的素数开始）

代码;

#include<iostream>
#include<cstdio> 
#include<cstring> 
using namespace std;
//判断素数 
int is_prime(int n)
{
	for (int i=2;i*i<=n;i++)
	{
		if(n%i ==0)
		{
			return 0;
		}
	}
	return 1;
 }
 int main()
 {
 	int n;
 	//存放素数的数组 
 	int prime[100],count=0;
 	for (int i=2;i<=100;i++)
 	{
 		if(is_prime(i))
 			prime[count++]=i;
	 }
	 
	
	while(scanf("%d",&n)==1)
	{
		printf("%d!= ",n);
		//存放指数 
		int p[100];
		memset(p,0,sizeof(p));
		int maxp=0;
		for(int i=1;i<=n;i++)
		{
			int m=i;
		    for(int j=0;j<count;j++)
			{
				while(m%prime[j]==0)//反复除以prime【j】，并累加p【j】 
				{
					m /= prime[j];
					p[j]++;
					//跟新最大素数因子的索引 
					if(j>maxp)
						maxp=j;
				}
			}
		}
		
		//输出指数 
		for(int i=0;i<=maxp;i++)
		printf(" %d",p[i]);
		
		cout<<endl;
	}
	
	return 0; 
}