c++实现筛除法求素数和分解质因数

关于代码中的maxn变量,它实际上是和输入的n是同一量级,或者说直接可以等于n;但和普通素数求法一样,sqrt(n)就可以解决问题。
知识点一:
筛除法原理:若n为素数,那他的n × \times × 2,n × \times × 3等等,就都不是素数,在之后的筛查就直接跳过。至于为何从j × \times × j开始,是因为从j × \times × 2开始,会有好多数重复计算。
知识点二:
vector变量做形参。有两种形式
1.地址形式,除在形参处加地址符外别无区别,见最终代码。
2.指针形式,形参处加指针符号*,实参处传地址&,使用时要用指针操作->

void initial(vector<int>*prime){
	bool isPrime[200];
	for(int i = 0;i<n;i++){
		isPrime[i] = true;
	}
	isPrime[0] = false;
	isPrime[1] = false;
	for(int j = 2;j<n;j++){
		if(!isPrime[j])
			continue;
		prime->push_back(j);//只能用指针操作
		for(int k = j*j;k<n;k += j){
			isPrime[k] = false;
		}
	}
}
int main(){
	vector<int>prime;
	initial(&prime);
	......
}

最终代码,形参形式

#include<iostream>
#include<cstdio>
#include<vector>
#include<cmath>

using namespace std;

const int maxn = sqrt(1e4) +1;

void initial(vector<int>&prime){
	bool isPrime[maxn];
	for(int i = 0;i<maxn;i++){
		isPrime[i] = true;
	}
	isPrime[0] = false;
	isPrime[1] = false;
	for(int j = 2;j<maxn;j++){
		if(!isPrime[j])
			continue;
		prime.push_back(j);
		for(int k = j*j;k<maxn;k += j){
			isPrime[k] = false;
		}
	}
}

int main(){
	vector<int>prime;
	int n;
	initial(prime);
	while(scanf("%d",&n) != EOF){
		int anwser = 0;
		for(int i = 0;i<prime.size() && prime[i] < n;i++)
		{
			int factor = prime[i];
			while(n%factor == 0){
				n /= factor;
				anwser++;
			}
		}
		if(n != 1)
			anwser++;

		printf("%d\n",anwser);
	}
	return 0;
}

### 回答1: #include <stdio.h> int main() { int n, i; printf("输入一个正整数: "); scanf("%d",&n); printf("分解质因数: "); for(i=2; i<=n; i++) { while(n%i == 0) { printf("%d ",i); n=n/i; } } printf("\n"); return 0; } ### 回答2: 要实现一个运行速度最快的分解质因数算法,可以使用Pollard Rho算法。这是一个基于随机性的质因数分解算法,其主要思想是通过不断迭代生成一个序列,其中的元素可能是合数,然后使用辗转相除法来判断是否找到一个质因数。 下面是使用C代码实现的简化版本: ``` #include <stdio.h> long long gcd(long long a, long long b) { if (b == 0) { return a; } else { return gcd(b, a % b); } } long long pollardRho(long long n) { long long x = 2, y = 2, d = 1; while (d == 1) { x = (x * x + 1) % n; y = (y * y + 1) % n; y = (y * y + 1) % n; d = gcd((x - y + n) % n, n); } return d; } int main() { long long n; printf("请输入一个正整数:"); scanf("%lld", &n); printf("%lld的一个质因数是:%lld\n", n, pollardRho(n)); return 0; } ``` 该代码通过求解随机迭代序列的非平凡最大公约数来找到一个质因数。与传统的试除法相比,Pollard Rho算法在较大的数上运行速度更快。在实际使用时,可以进一步改进代码以提高效率,并处理特殊情况。 ### 回答3: 要编写一个运行速度最快的分解质因数算法,可以采用试除法的优化形式,又称为Pollard's rho算法。 首先,我们需要写一个用于检测一个数是否为素数的函数。可以使用试除法,从2到sqrt(n)逐个检查是否有因数。如果找到一个因数,说明n不是素数;如果没有找到因数,说明n是素数。 接下来,我们可以使用Pollard's rho算法来分解质因数。这个算法的思想是利用了环的周期性来找到质因数。 具体步骤如下: 1. 定义一个函数f(x),表示f(x) = x^2 + 1。 2. 初始化x = 2和y = 2,计算f(x)和f(f(y))。 3. 利用欧几里得算法计算f(x)和f(f(y))的差值d。 4. 如果d是一个质数,则d即为n的一个质因数。 5. 如果d是合数,重复步骤2-4直到找到一个质因数。 通过不断重复步骤2-4,我们可以找到n的一个质因数。然后可以递归地对这个质因数进行分解,直到无法再分解为止。 通过这种方式,可以得到一个快速的质因数分解算法。当然,实现中还可以对代码进行优化,例如可以用位运算代替乘法和除法操作,以提高运行速度。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值