素数

素数

判断素数

判断到$\sqrt{n}$即可。

bool isPrime(int n){
		if(n<=1)
			return false;
		for(int i=2;i*i<=n;i++){
			if(n%i==0){
				return false;
			}		
		}
		return true;
}

获取素数表

埃氏筛选法----从2开始筛选掉合数

const int maxn = 101;//表长
int prime[maxn],pNum = 0;
//prime数组存放所有素数，pNum为素数个数
bool p[maxn] = {0};
//i为素数，p[i]为false；否则，p[i]为true
void find_prime(){
		for(int i=2;i<maxn;i++){
			if(p[i]==false){
				prime[pNum++]=i;//素数存入素数数组
			for(int j=i+i;j<maxn;j+=i){
				//筛去所有i的倍数
				p[j]=true;
				}
			}
		}
}

质因子分解–素数的应用

对一个正整数n来说，如果它存在[2,n]范围内的质因子，要么这些质因子全部小于等于$\sqrt{n}$，要么只存在一个大于$\sqrt{n}$的质因子，其余质因子全部小于等于$\sqrt{n}$。

Prime Factors

#include<stdio.h>
struct factor{
    int cnt;//个数
    int x;//质因子
}fac[10];
const int maxn = 100001;//表长
int prime[maxn],pNum = 0;
//prime数组存放所有素数，pNum为素数个数
bool p[maxn] = {0};
//i为素数，p[i]为false；否则，p[i]为true
bool isPrime(int n){
		if(n<=1)
			return false;
		for(int i=2;i*i<=n;i++){
			if(n%i==0){
				return false;
			}
		}
		return true;
}
void find_prime(){
		for(int i=2;i<maxn;i++){
			if(p[i]==false){
				prime[pNum++]=i;//素数存入素数数组
			for(int j=i+i;j<maxn;j+=i){
				//筛去所有i的倍数
				p[j]=true;
				}
			}
		}
}
int main(){
    int n;//输入数字 num统计质因子个数
    while(scanf("%d",&n)!=EOF){
        for(int i=0;i<10;i++){
            fac[i].cnt = 0;
            fac[i].x = 0;
        }
        isPrime(n);
        find_prime();
        int num=0;//fac数组下标
        if(n==1)
            printf("1=1\n");
        else{
            printf("%d=",n);
            //枚举根号n以内的数字
            int t_n=n;//n在不停变化，不能作为判断条件
            for(int i=0;i<pNum&&p[i]*p[i]<t_n;i++){
                if(n%prime[i]==0){
                    while(n%prime[i]==0){
                        fac[num].cnt++;
                        fac[num].x=prime[i];
                        n/=prime[i];
                    }
                    num++;
                }
                if(n==1)//若n等于1 及时退出
                    break;
            }
            if(n!=1){//若此时n不为1，则n为最后一个因子
                fac[num].cnt++;
                fac[num++].x=n;
            }
            for(int i=0;i<num-1;i++){
                printf("%d",fac[i].x);
                if(fac[i].cnt>1){
                    printf("^%d",fac[i].cnt);
                }
                printf("*");
            }
            printf("%d",fac[num-1].x);
            if(fac[num-1].cnt>1){
                printf("^%d",fac[num-1].cnt);
            }
            printf("\n");
        }
    }
    return 0;
}