筛选（线性+埃氏）

最新推荐文章于 2022-05-13 20:22:41 发布

JWei_7

最新推荐文章于 2022-05-13 20:22:41 发布

阅读量188

点赞数 2

JWei_7

本文链接：https://blog.csdn.net/qq_54729417/article/details/115708317

版权

算法

(假设)题目：
给出一个正整数n，打印出所有从1~n的素数(即质数)，关键是要找出一个判断一个正整数n是否为素数的方法

  如果我们没有学习算法  那么我们就是用最朴素的算法
  废话不多说  上代码：

对于每一个正整数我们都做如下操作：

#include<stdio.h> //运行时间的话  O(n)
int main()
{
  int i,n;
  while(~scanf("%d",&n))
  { for(i=2;i<n;i++)
         if(n%i==0)    break;
    if(i==n)    printf("YES\n");
    else           printf("NO\n");
  }
}

那么针对上一个代码进行优化
我们可以发现循环 i<n 这里可以优化一下，打个比方 8%2=0 就发现8不是素数，所以没必要做到8%4=0去，只要i<=sqrt(n) 即可判断完；

#include<stdio.h>   //对于每个正整数的时间复杂度为O（log n）
#include<math.h>
int main()
{ int i,n,x;
  while(~scanf("%d",&n))
  { x=(int)sqrt(n);//因为sqrt返回的值为double；
    for(i=2;i<=x;i++)
         if(n%i==0)    break;
    if(i>x)    printf("YES\n");
    else           printf("NO\n");
  }

说完以上俩种其实可以发现效率并不是很高如果处理较大的数据时就会很头疼。因为都有时间限制

埃拉托斯特尼筛法（不清楚的可以点链接看一下bd怎么说）

基本思想：素数的倍数一定不是素数
实现方法：用一个长度为N+1的数组保存信息（0表示素数，1表示非素数），先假设所有的数都是素数（初始化为0），从第一个素数2开始，把2的倍数都标记为非素数（置为1），一直到大于N；然后进行下一趟，找到2后面的下一个素数3，进行同样的处理，直到最后，数组中依然为0的数即为素数。
这样时间复杂度 O（nloglogn）就降下来很多很多了

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e6+10;
int a[maxn];//标记 
int A[100000];//保存素数 
int main()
{  int n,x=0;
  scanf("%d",&n);
  memset(a,0,sizeof(a));c++函数   目的是把数组初值都复赋为0；
  for(int i=2;i<=n;i++)
  {
  	if(a[i]==0) //未被标记
	  {
	  	A[x++]=i;
	   } 
	   for(int j=i+i;j<=n;j+=i)
	       a[j]=1;//把i的倍数都筛掉； 
  }
   printf("1-%d的素数有：\n",n);
  for(int i=0;i<x;i++) //打印表
   printf("%d ",A[i]) ;
}

例如：
在这里插入图片描述
以上就是埃氏筛选法，但是自习想一想在处理比较大的数据事是不是会重复很多次，
打个比分
n=6时
i=2的时候可以筛选一次
当i=3的时候又会筛选一次这是一个较小的数据，如果n很大的时候是不是会多做很多次无用的筛选，
所以就有了

线性筛选

线性筛，复杂度为O(n)。与埃氏筛相比，不会对已经被标记过的合数再进行重复标记，故效率更高。欧拉筛将合数分解为 (最小质因数 * 一个合数) 的形式，通过最小质因数来判断当前合数是否已经被标记过。

const int maxn = 101;   // 表长
int prime[maxn], pNum = 0;    // prime记录素数，pNum记录素数个数 
bool p[maxn] = {false};        // p记录当前数是否被筛去

void sushu(int n)    // 查找记录2-n的素数
{
	for (int i = 2; i <= n; i++)
	{
		if (p[i] == false)  // 如果未被筛过，则为素数
			prime[pNum++] = i;//记录素数
		for (j = 0; j < pNum; j++)//
		{
			if (i * prime[j] > n)      // 当要标记的合数超出范围时跳出
				break;
			p[i * prime[j]] = true;     // 将已经记录的素数的倍数进行标记
			if (i % prime[j] == 0)      //关键步骤  用最小的因数去筛
				break;
		}
	}
}

if (i % prime[j] == 0)
这个是关键的一步
可以参考一下别的up主的思路
在这里插入图片描述
对于 i%prime[j] == 0 就break的解释：当 i是prime[j]的倍数时，i = kprime[j]，如果继续运算 j+1，i * prime[j+1] = prime[j] * k prime[j+1]，这里prime[j]是最小的素因子，当i = k * prime[j+1]时会重复，所以才跳出循环。
举个例子：i = 8 ，j = 1，prime[j] = 2，如果不跳出循环，
prime[j+1] = 3，8 * 3 = 2 * 4 * 3 = 2 * 12，在i = 12时会计算。
因为欧拉筛法的原理便是通过最小素因子来消除。