习题 47：因子分解阶乘版★

题目描述：
给你一个小于10000的非负整数n，然后对n!的结果进行因子分解


输入：
多组测试数据，每行就是数值n


输出：
输出分解结果，详见样例


样例输入：
1
5
7


样例输出：
1! = 1
5! = 2^3 * 3 * 5
7! = 2^4 * 3^2 * 5 * 7


难度：very easy

site:http://www.yzfy.org/dis/listpost.php?tid=164&extra=page%3D1

TIPS：
素数：
n被2，3，4，，，，sqrt（n）整除就不是素数，不能除尽就是素数，

给定两个数m,n
求m!分解质因数后因子n的个数。
这道题涉及到了大数问题，如果相乘直接求的话会超出数据类型的范围。
下面给出一种效率比较高的算法，我们一步一步来。
m!=1*2*3*……*(m-2)*(m-1)*m
可以表示成所有和n倍数有关的乘积再乘以其他和n没有关系的
??? =(n*2n*3n*......*kn)*ohter???? other是不含n因子的数的乘积???因为 kn<=m?而k肯定是最大值? 所以k=m/n
??? =n^k*(1*2*......*k)*other?? 
??? =n^k*k!*other　　　　　
从这个表达式中可以提取出k个n，然后按照相同的方法循环下去可以求出k!中因子n的个数。
每次求出n的个数的和就是m!中因子n的总个数
然后说N的阶乘：
例如:20!1.先求出20以内的素数,(2,3,5,7,11,13,17,19)2.再求各个素数的阶数
e(2)=[20/2]+[20/4]+[20/8]+[20/16]=18;
e(3)=[20/3]+[20/9]=8;
e(5)=[20/5]=4;
...
e(19)=[20/19]=1;
所以20!=2^18*3^8*5^4*...*19^1

// 47.cpp : Defines the entry point for the console application.
//
#define PB_47
#define CP_VC6

#include "stdafx.h"
#include "stdio.h"
#include "math.h"

int prime[100];
int pos;

void findprime(int n, int *a)
{
	int i, j, flag;
	for (i=2; i<=n; i++)
	{
		flag = 0;
		for(j=2; j<=sqrt(i); j++)
		{
			if(i%j == 0) flag = 1; break;
		}
		if( flag == 0) { a[pos++] = i;}

	}
}


int main(int argc, char* argv[])
{
	int n, i, tmp, result, input;
	result = 0;
	pos = 0;
	while( scanf("%d",&n) != EOF )
	{
		if(n == 0) 
		{
			printf("0! = 1\n");
			continue;
		}
		else if(n == 1)
		{
			printf("1! = 1\n");
			continue;
		}
		findprime(n,prime);
		printf("%d! = ",n);
		i = 0;
		while(--pos >= 0)
		{
			printf("%d", prime[i]);
			tmp = prime[i++];
			input = n;
			while(input)
			{
				result += input/tmp;
				input  /= tmp;
			}
			if(result != 1) printf("^%d",result);
			if(pos != 0)printf(" * ");
			result = 0;
		}
		pos = 0;
		printf("\n");
	}
	return 0;
}