题意
求出n所能分解的素因子。
思路
这个应该是套路题,就是打印素数表,然后利用fac数组,枚举1~sqrt(n)范围内的所有素因子p,判断p是否为n的因子。
还要注意如果在枚举完之后,n仍然>1,那么必定有一个大于sqrt(n)的素因子(很有可能就是n本身),这部分也需要处理。
#include<cstdio>
#include<cmath>
using namespace std;
const int maxn=100010;
struct factor{
int x,cnt; //x为质因子,cnt为质因子的个数
}fac[10];
bool is_prime(int n){ //判断是否为素数
if(n==1) return false;
int sqr=(int)sqrt(1.0*n);
for(int i=2;i<=sqr;i++){
if(n%i==0) return false;
}
return true;
}
int prime[maxn],pNum=0;
void Find_Prime(){ //打印素数表
for(int i=1;i<maxn;i++){
if(is_prime(i)==true){
prime[pNum++]=i; //prime[]就是素数表
}
}
}
int main(){
Find_Prime(); //必须先打印出素数表
int n,num=0; //num为n的不同素因子的个数
scanf("%d",&n);
if(n==1) printf("1=1"); //特别判定n=1的情况
else {
printf("%d=",n);
int sqr=(int)sqrt(1.0*n); //枚举<sqrt(n)的素因子
for(int i=0;i<pNum&&prime[i]<=sqr;i++){
if(n%prime[i]==0){ //prime[i]是n的素因子
fac[num].x=prime[i]; //记录该因子
fac[num].cnt=0; //初始化个数
while(n%prime[i]==0){ //计算出该因子的个数
fac[num].cnt++;
n/=prime[i];
}
num++; //不同素因子的个数+1
}
if(n==1) break; //及时退出循环
}
if(n!=1){ //如果无法被根号n以内的素因子除尽
fac[num].x=n; //那么一定有一个大于根号n的素因子
fac[num++].cnt=1;
}
//按格式输出结果
for(int i=0;i<num;i++){
if(i>0) printf("*"); //第二个素因子开始要用乘号*
printf("%d",fac[i].x); //素因子
if(fac[i].cnt>1){ //如果个数>1,则需要把个数也输出
printf("^%d",fac[i].cnt);
}
}
}
return 0;
}
注意点:
需要在循环外定义sqr,要是在循环条件中直接计算sqrt(n)会导致答案出错。