ACM：因子和阶层

题目：输入正整数n（2≤n≤100），把阶乘n!=1×2×3ׄ×n分解成素因子相乘的形式，从 小到大输出各个素数（2、3、5、„）的指数。例如825=3×52 ×11应表示成(0,1,2,0,1)，表示分别有0、1、2、0、1个2、3、5、7、11。你的程序应忽略比最大素因子更大的素数（否则末尾会有无穷多个0）。 

 样例输入：  5

                       53 

样例输出： 5!=3 1 1 

                      53!=49 23 12 8 4 4 3 2 2 1 1 1 1 1 1 1

 

分析：打印素数表。

 

#include <iostream>
#include <string>
#include <assert.h>
#include <math.h>
using namespace std;

const int MAXN = 1000;
int n, x[MAXN], y[MAXN], G[MAXN][MAXN], d[MAXN];

int is_prime(int x) {               //判断一个数是否为素数。
	assert(x >= 0);
	if(x == 1) return 0;
	int m = floor(sqrt(x) + 0.5);       //用这种方式避免溢出！如果用 i * i < n 的话，i太大的时候会超过int的最大值！
	for(int i = 2; i <= m; ++i) {
		if(x % i == 0) return 0;
	}
	return 1;
}

int prime[100];   //素数表！ 
int count = 0;

int main() {
	for(int i = 2; i <= 100; ++i) {             //构造素数表！
		if(is_prime(i))  prime[count++] = i;
	}
	int n, p[100];       //p[100]记录各素数的指数！
	while(cin >> n) {
		cout << n << "! " << "=";
		memset(p, 0, sizeof(p));
		int maxp = 0;
		for(int i = 1; i <= n; ++i) {         //将n！一个一个分解出来！
			int m = i;
			for(int j = 0; j < count; ++j) {
				while(m % prime[j] == 0) {
					m /= prime[j];
					p[j]++;
					if(j > maxp) maxp = j;
				}
			}
		}
		for(int i = 0; i <= maxp; ++i) {
			cout << " " << p[i];
		}
		cout << endl;
	}
	return 0;
}