PAT A1059 Prime Factors (25point(s))

题目链接

#include<cstdio>
#include<math.h>
int prime[100001],p[100001],pnum=0;
struct factor{
    long long x,cnt;//x是质因子，cnt是个数
}fac[10];
void find_prime(){//筛法打印素数表
    for(int i=2;i<10001;i++){
        if(p[i]==0){
            prime[pnum++]=i;
            for(int j=i+i;j<10001;j+=i){
                p[j]=1;
            }
        }
    }
}
int main(){
    find_prime();//别忘了调用
    long long n,num=0;//num是质因子的个数
    scanf("%lld",&n);
    if(n==1) {//需特判
    printf("1=1");
    return 0;
    }
    printf("%lld=",n);
    int sqr=(int)sqrt(n*1.0);
    for(int i=0;i<pnum&&prime[i]<=sqr;i++){
        //枚举根号n以内的质因子
        //别忘了加上i<pnum的限制条件，否则会出现除以零的错误
        if(n%prime[i]==0){//该素数是质因子
         fac[num].x=prime[i];
         fac[num++].cnt=0;
        }
        while(n%prime[i]==0){//计算其个数
            fac[num-1].cnt++;
            n/=prime[i];
        }

    }
    if(n!=1){//如果无法被根号n以内的质因子除尽
        fac[num].x=n;//一定有且仅有1个大于根号n的质因子
        fac[num++].cnt=1;
    }
     for(int i=0;i<num;i++){//按格式输出
       printf("%lld",fac[i].x);
       if(fac[i].cnt>1) printf("^%lld",fac[i].cnt);
       //质因子的个数为1则不需要输出
       if(i!=num-1) printf("*");
   }
   return 0;
}