pat 甲级 1059 Prime Factors

1059 Prime Factors

分数 25

全屏浏览题目

切换布局

作者 HE, Qinming

单位 浙江大学

Given any positive integer N, you are supposed to find all of its prime factors, and write them in the format N = p1​k1​×p2​k2​×⋯×pm​km​.

Input Specification:

Each input file contains one test case which gives a positive integer N in the range of long int.

Output Specification:

Factor N in the format N = p1​^k1​*p2​^k2​*…*pm​^km​, where pi​'s are prime factors of N in increasing order, and the exponent ki​ is the number of pi​ -- hence when there is only one pi​, ki​ is 1 and must NOT be printed out.

Sample Input:

97532468

Sample Output:

97532468=2^2*11*17*101*1291

思路：打出素数表，用n除以这个素数并计数，一直到无法除尽，即可得到其一项素因子的个数和该数因子的种类，储存到结构题里面，即一个结构体储存素因子的个数和该素因子的种类；

考点：打印素数表，枚举法；

          一个人数的质因子只有两种情况，都小于sqr，第二种是有一个数大于sqr，其余数都小于sqr；

代码：

#include<iostream>
#include<cmath>
using namespace std;
const int maxn=100010;
int primetable[maxn]={0};
bool isprime(int num){//判断是不是质数； 
	if(num==1)return false;
	int sqr=(int)sqrt(1.0*num);
	for(int i=2;i<=sqr;i++){
		if(num%i==0)return false;
	}
	return true;
}
int primenum=0;
void find_prime(){//寻找所有可能的质数； 
	for(int i=1;i<maxn;i++){
		if(isprime(i)==true)primetable[primenum++]=i;
	}
}
struct factor{//质因子，质因子的个数，从小到大存储； 
	int cnt;
	int x;
}f[10];
int main(){
	int n,num=0;
	find_prime();
	cin>>n;
	if(n==1){//1是特例; 
		cout<<"1=1";
	}else {
		cout<<n<<"=";
		int sqr=(int )sqrt(1.0*n);
		//枚举法; 
		for(int i=0;i<primenum&&primetable[i]<=sqr;i++){//第一种情况，如果质因子都是小于sqr,sqr是向下取的,所以sqr可以取值,如15开方得3.8,int向下取是3，3是15得质因子； 
			if(n%primetable[i]==0){
				f[num].x=primetable[i];
				f[num].cnt=0;
				while(n%primetable[i]==0){
					f[num].cnt++;
					n/=primetable[i];
				}
				num++; 
				if(n==1){
					break;
				}
			}
		}
		if(n!=1){//如果到了sqr还没被除得只剩1，说明有一个大于sqr的因子，该因子就是此时的n（n>sqr），15=3*5，15/3==5，5>sqr； 
			f[num].x=n;
			f[num++].cnt=1;
		}
	}
	for(int i=0;i<num;i++){
		if(i>0)cout<<"*";
		cout<<f[i].x;
		if(f[i].cnt>1){
			cout<<"^"<<f[i].cnt;
		}
	}
	return 0;
}

纯净版：

#include<iostream>
#include<cmath>
using namespace std;
const int maxn=100010;
int primetable[maxn]={0};
bool isprime(int num){
	if(num==1)return false;
	int sqr=(int )sqrt(1.0*num);
	for(int i=2;i<=sqr;i++){
		if(num%i==0)return false;
	}
	return true;
}
int primenum=0;
void find_prime(){
	for(int i=1;i<maxn;i++){
		if(isprime(i)==true){
			primetable[primenum++]=i;
		} 
	}
}
struct factor{
	int cnt;
	int x;
}f[10];
int main(){
	int n,num=0;
	cin>>n;
	find_prime();
	if(n==1){
		cout<<"1=1";
	}else {
		cout<<n<<"=";
		int sqr=(int )sqrt(1.0*n);
		for(int i=0;i<primenum&&primetable[i]<=sqr;i++){
			if(n%primetable[i]==0){
				f[num].x=primetable[i];
				f[num].cnt=0;
				while(n%primetable[i]==0){
				f[num].cnt++;
				n/=primetable[i];	
				}
				num++; 
			}
			if(n==1)break;
		}
		if(n!=1){
			f[num].cnt=1;
			f[num++].x=n;
		}
	}
	for(int i=0;i<num;i++){
		if(i>0)cout<<"*";
		cout<<f[i].x;
		if(f[i].cnt>1)cout<<"^"<<f[i].cnt;
	}
    return 0;
}

埃式筛法求素数表：

#include<iostream>
#include<cmath>
using namespace std;
const int maxn=100010;
int prime[maxn]={0};
bool primetable[maxn]={false};
int primenum=0;
void find_prime(){//筛法求素数表; 
	for(int i=2;i<maxn;i++){
		if(primetable[i]==false){//如果是素数; 
			prime[primenum++]=i;
			for(int j=i+i;j<maxn;j+=i){
				primetable[j]=true;//标记为不是素数; 
			}
		}
	}
}
struct factor{
	int cnt=0;
	int x;
}f[10];
int main(){
	int n,num=0;
	find_prime();//没写括号也是可以过编译的，但是相当于没用调用函数，有时候发现函数没有调用，可能不是函数的问题，而是忘记用括号了;
	cin>>n;
	if(n==1){
		cout<<"1=1";
	}else {
		cout<<n<<"=";
		int sqr=(int )sqrt(1.0*n);
		for(int i=0;i<primenum&&prime[i]<=sqr;i++){
			if(n%prime[i]==0){
				while(n%prime[i]==0){
				f[num].cnt++;
				f[num].x=prime[i];
				n/=prime[i];	
				}
				num++;
			}
			if(n==1)break;//因为是在sqr内讨论,一直找质因数相除，最后肯定是得1，除非没有质因数，即一个数有质因子，一定可以被整除为1，否则因子只有它和它本身； 
		} 
		if(n!=1){
			f[num].cnt=1;
			f[num++].x=n;
		}
	}
	for(int i=0;i<num;i++){
		if(i!=0)cout<<"*";
		cout<<f[i].x;
		if(f[i].cnt>1)cout<<"^"<<f[i].cnt;
	}
	return 0; 
}