给你一个小于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;
}