欧拉筛的优化

最新推荐文章于 2024-05-20 22:00:16 发布

GuGu_boom

最新推荐文章于 2024-05-20 22:00:16 发布

阅读量317

点赞数

分类专栏：编程杂疑文章标签：算法数据结构

本文链接：https://blog.csdn.net/X720819/article/details/129166251

版权

编程杂疑专栏收录该内容

5 篇文章 0 订阅

订阅专栏

如何判断一个素数：

1. 素数的定义：

2.代码实现：

int JudgePrime(int number)
{
    if(number==0||number==1)return 0;
    for(int i = 2;i<number;i++)
    {
        if(number%i==0)return 0;
    }
    return 1;
}

从头到尾遍历，时间复杂度O(n);

3.优化：

一个非素数必可以分解为两个因数之积，而两个因数有必然存在大小之分（或是相等），也就是说一个小因数必然对应一个大因数，把所有情况列举出来，越靠近中间两因数数值越接近，直至相等停止（这里这么说是因为不仅考虑了整数）。

如图：

代码实现：

#include <math.h> //提供sqrt函数
int JudgePrime(int number)
{
    if(number==0||number==1)return 0;
    for(int i = 2;i<=sqrt(number);i++)
    {
        if(number%i==0)return 0;
    }
    return 1;
}
/*
for循环判断条件也可以这样写：i*i<=number
但不建议，因为易造成数据溢出
*/

遍历到根号n停止，时间复杂度O(√n)

数组极大的情况下如何判断一个素数：

判断某个数是不是素数，时间很短，但若是判断一个极大数组里面有哪些是素数，该素数是第几个素数的情况下，挨个数判断的方法所用的时间就过长了，就像我们如果有不认识的字，翻字典是一种可行的方法，并且在使用二分查找方法的情况下缩短了一定的时间（时间复杂度O（logn）），但我们更倾向于用浏览器，我们输入某个字，按下搜索立马就可以得到答案（时间复杂度O(1)）。

原理：我们输入一个问题，立马就能得到答案，类似输入array[i]，立刻能得到该位置上存储的值，这是因为一个元素映射另一元素，（一个键key对应一个值value）

我们可以创建一个数组用来存储素数，其下标对应素数的顺序（第几个素数），原数组下标表示对应数值，其内存储元素0或1表示是否是素数。

但目前还存在一个问题，我们仍需要遍历原数组，挨个判断一个数是不是素数；

解决方法：埃氏筛

原理：非素数必存在除1和本身之外的因数，一个素数✖除1外的任意数会得到一个数值，而该素数和另一个任意数一定是这一数值的因数，也就是说两数相乘（素数×非1任意数）之积必定为非素数。

实现：以素数为基数，从2开始遍历所有整数，直到两数之积越界停止，整个过程得到的所有积，将其开关（数组内存储元素0，1实现开关的开闭）设置为闭合

代码：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define MAXSIZE 10005
int nums[MAXSIZE]={1, 1};//对应0，1两个非素数，赋值为1表示开关闭合，其他元素默认为0
int storage[MAXSIZE];
int main()
{
    for(int i = 2,j = 0;  i < MAXSIZE; i++)
    {
        if(!nums[i])//判断开关开闭，即是否为素数，如果是进入循环
        {
            storage[j++]=i;//存储素数并自增以记录素数个数
            for(int k = i * 2; k < MAXSIZE; k += i)//以i为基数，从2i开始遍历，每次递增一个i值
            {
                nums[k] = 1;//闭合开关
             }
        }
    }
    int n;
    while(~scanf("%d",&n))//查询素数
        printf("%d\n",storage[n-1]);
    return 0;
}

时间复杂度O（nloglogn）

这里时间复杂度的计算涉及公式

欧拉筛的原理：

因为埃氏筛涉及到重复筛除，所以引入线性筛

原理：

每个合数一定能分解为若干质数之积

代码：

#include <stdio.h>
#include <string.h>
#define MAX 100001
_Bool mark[MAX];//默认初始化为0
int prime[MAX]={2};//第一个素数是2
int main()
{
    for(int i=2,k=0;i<MAX;i++)//2初始化时已经赋值，无需考虑
	{
		if (mark[i]==0)//在100001里面找到质数并且标记
            prime[k++]=i;	//如果没有被标记为1，就是质数
		for(int j=0;j<k;j++)// //j小于当前所有的质数的个数
		{
			if(i*prime[j]>MAX)break;// 如果超出给出的范围，那么就退出循环
			mark[i*prime[j]]=1;//用质数数依次×i，结果标记为合数（也就是标记为1）。
			if(i%prime[j]==0)break;//关键:只标记一次
		}
	}
    int number;
    while(~scanf("%d",&number))//查询
        printf("\t%d\n",prime[number-1]);
    return 0;
}

时间复杂度O(n)