c++编译
算法笔记
注释:
(1)枚举1~sqrt(n)范围内所有质因子p,判断p是否是n的因子。如果p是n的因子,那么给fac数组中增加质因子p,并初始化其个数为0,然后只要p还是n的因子,就让n不断除以p,每次操作令p的个数加1,直到p不在是n的因子为止。
(2)如果p不是n的因子,直接跳过。如果在上面步骤结束后n仍然大于1,说明n有且仅有一个大于sqrt(n)的质因子(有可能是n本身),这时需要把这个质因子加入fac数组,并令其个数为1。
(3)题目说的是int范围内的正整数进行质因子分解,因此素数表打概开10^5大小就可以了。
(4)注意n==1需要特判输出“1=1”,否则不会输出结果。
(5)容易忘记的地方:一是在main函数开头忘记调用find_prime()函数;二是find_prime函数中把i<maxn误写成i<=maxn;三是没有处理大于sqrt(n)部分的质因子;四是在枚举质因子的过程中发生了死循环;五是没有在循环外定义变量来存储sqrt(n),而在循环条件中直接计算sqrt(n),这样当循环中使用n本身进行操作的话会导致答案错误。
代码
#include <iostream>
#include <math.h>
#include <stdio.h>
using namespace std;
const int maxn=100010;
bool is_prime(int n){
if(n==1)
return false;
int b=(int)sqrt(1.0*n);
for(int i=2;i<=b;i++){
if(n%i==0)
return false;
}
return true;
}
int prime[maxn];
int pNum=0;
void find_prime(){
for(int i=1;i<maxn;i++){
if(is_prime(i)==true){
prime[pNum++]=i;
}
}
}
struct factor{
int x,cnt;
}fac[10];
int main()
{
find_prime();
int n;
int num=0;
cin>>n;
if(n==1)
printf("1=1");
else{
printf("%d=",n);
int sqr=(int)sqrt(1.0*n);
for(int i=0;i<pNum&&prime[i]<=sqr;i++){
if(n%prime[i]==0){
fac[num].cnt=0;
fac[num].x=prime[i];
while(n%prime[i]==0){
fac[num].cnt++;
n/=prime[i];
}
num++;
}
if(n==1)break;
}
if(n!=1){
fac[num].x=n;
fac[num++].cnt=1;
}
for(int i=0;i<num;i++){
if(i>0)
cout<<"*";
cout<<fac[i].x;
if(fac[i].cnt>1)
printf("^%d",fac[i].cnt);
}
}
return 0;
}