PAT A1059 Prime Factors

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;
}