对筛质数三种方法的理解[普通筛,诶氏筛,线性筛]

最新推荐文章于 2025-04-11 22:28:29 发布

泥烟

最新推荐文章于 2025-04-11 22:28:29 发布

阅读量507

点赞数

分类专栏：数据结构与算法文章标签：算法

本文链接：https://blog.csdn.net/qq_39391544/article/details/123670350

版权

数据结构与算法专栏收录该内容

41 篇文章

订阅专栏

这篇博客探讨了三种不同的质数筛选算法：普通筛法、埃拉托斯特尼筛法和线性筛法。通过代码实现和时间复杂度分析，展示了线性筛法在效率上的显著优势，特别是在处理大规模数据时，其时间复杂度为O(n)，远低于普通筛法的O(nlogn)和埃拉托斯特尼筛法的O(nloglogn)。对于10^7和10^8规模的数据，线性筛法分别只需要0.098s和0.948s，显著优于其他两种方法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

//primes数组用来存放质数
int primes[N], cnt;

//i如果是质数则false,否则为true
bool st[N];

1.普通筛法 O(nlogn)

void get_primes(){
    for(int i=2;i<=n;i++)
    {
        if(!st[i]) primes[cnt++] = i; //把素数存起来
        for(int j = i; j <= n; j += i)//不管是合数还是质数，都用来筛掉后面它的倍数
            st[j] = true;
    }
}

2.诶氏筛法 O(nloglogn)

void get_primes(){
    for(int i = 2; i <= n; i ++)
    {
        if(!st[i]){
            primes[cnt++] = i; //把素数存起来
            for(int j = i; j <= n; j += i)//仅用质数就可筛掉所有合数
                st[j] = true;
        }   
    }
}

3.线性筛法 O(n)

void get_primes(){
	for(int i = 2; i <= n; i ++)
	{
		if(!st[i]) primes[cnt++] = i;
		for(int j = 0; primes[j] <= n/i; j ++)
		{
			st[primes[j]*i] = true; //用最小质因子筛合数(已经挑出来的质因子的i倍)
			if(i % primes[j] == 0) break; // 此时primes[j]是i的i的最小质因子
		}
	}
}

if(i % primes[j] == 0) break;

对于👆的理解:

对于一个数c=ab(b为c的最小质因数），当通过该算法的循环循环至cb时，易得此时c%b==0,如果此时继续循环至b后面的一个素数d，则有：cd=abd=(ad)b，因为d>b,所以ad>c。当循环从c继续查找到ad时我们发现当ad再次与素数b想乘时，就又对c*d进行了一次操作，出现了冗余，所以在if(n%prime[j]==0)成立时要将该层循环break掉；

举个例子，对于一个数9，92=18将18标记为合数，循环继续；93=27将27标记为合数，此时发现9%3=0，循环退出。如果将循环继续下去会出现筛除95=45的情况，而45=153，在15时会被在筛去一次，故不可行