今天浅学一下一个简单算法
目录
普通的素数判断
通常,我们写的筛选素数的算法是这样的:
bool isPrime(int x) {
for (int i = 2; i <= sqrt(x); ++i)
{
if (x % i == 0)
{
return false;
}
}
return true;
}
如果要检查n个数时间复杂度为O(n*sqrt(n)),空间复杂度为O(1)
埃氏筛
下面来看看埃氏筛,假设要求0-100的素数个数
#include <stdio.h>
#include <stdbool.h>
#include <stdlib.h>
#define maxn 101
void sieve(bool prime[])
{
for (int i = 0; i < maxn; i++)
prime[i] = true;
prime[0] = prime[1] = false;
for (int i = 2; i < maxn; i++)
{
if (prime[i])
{
for (int j = 2 * i; j < maxn; j += i)
{
prime[j] = false;
}
}
}
}
int main()
{
bool prime[101];
sieve(prime);
int count = 0;
for (int i = 0; i < 101; i++)
{
if (prime[i])
count++;
}
printf("0-100间共有%d个素数", count);
return 0;
}
时间复杂度O(nloglogn),证明参考埃氏筛时间复杂度证明
空间复杂度O(n),因为需要开辟n个空间来存放是否是素数
比较完以后来看一下埃氏筛的原理,首先,我们要把存放素数的数组全部置为1(也可以用0表示是素数,不是素数的用下标表示,这里有时间浪费,在下面用更快的方法写),假设刚开始每个数是素数,然后因为0和1不是素数,将这两个置成0,接下来从2开始判断,2是素数,那么2的倍数就都不是素数,4,6,8,10......都不是,把他们置成0。3是素数,那么3的倍数也都不是素数,6,9,12.....置成0。接下来4不是素数,因为在2的地方已经置成0了,不进入循环。5是素数,6不是素数,以此类推,直到筛选完。
通俗易懂还好用,以后筛素数再也不用暴力循环遍历了
但是我们发现这种算法还是有漏洞,比如6这个数字,我们在i=2的时候已经判断把6置成0了,而在i=3的时候又多把它置成0一次,这就造成了时间消耗,但是这些漏洞可以用一个更为牛逼的算法解决——线性筛。
线性筛
首先,我们可以肯定的是合数是可以由很多素数相乘得到的,那么这些因子中肯定有最大的,有最大的就肯定有最小的,比如12=3*4和12=2*6,而且最大的和最小的一定同时出现,比如2是12的最小因子,6是12的最大因子。每个数都有且只有一个最大因子,那么我们把每个数按照最大因子筛选,是不是就可以避免上面埃氏筛的问题了呢,比如6,我们在最大因子为2的这一次循环判断中不筛掉,只在3的这一次循环中筛掉,12只在6中筛掉。再者,最小因子一定小于等于最大因子,可以类比判断素数算法中的sqrt(n),比方说25=5*5,最大因子和最小因子都是5,49=7*7,最大因子和最小因子都是7,也就是说在每一轮筛选中我们最多筛选到最大因子和最小因子相等时就停下。
#include <stdio.h>
#include <stdbool.h>
#include <stdlib.h>
#define maxn 101
int main()
{
bool check[101] = { 0 };//用于检查素数,0表示是素数,1表示不是素数
int prime[101] = { 0 };//用于存放素数
int count = 0;//用于记录素数个数
for (int i = 2; i < maxn; i++)//i是最大的因子
{
if (!check[i])
{
prime[count++] = i;//如果是素数,count++,同时把素数存进数组
}
for (int j = 0; j < count && i * prime[j] < maxn; j++)
{
check[i * prime[j]] = 1;//把最大因子与最小因子乘积筛掉
if (i % prime[j] == 0)//保证只有当i是最大因子的时候进行筛选
{
break;
}
}
}
printf("0-100间共有%d个素数:", count);
for (int i = 0; i < count; i++)
{
printf("%d ", prime[i]);
}
printf("\n");
if (!check[55])
printf("55是素数\n");
else
printf("55不是素数\n");
return 0;
}
当最大因子为4的时候,并没有4*3=12这一列,因为4%2可以整除,当j++时,4*3=12,但是我们想一下,这个时候4是可以拆一个2出来和3相乘在一起变成6,也就是说最大因子不是4而是6了。我们再举个例子,当最大因子为6的时候,6%2可以整除,6*3=18,但是此时6可以拆出来一个3,3*3=9,也就是说18的最大因子是9
if (i % prime[j] == 0)//保证只有当i是最大因子的时候进行筛选
{
break;
}
因此这一步就是线性筛的精华!
线性筛的时间复杂度为O(n),空间复杂度为O(n),来看一下运行后的结果:
学会这个以后筛素数又快又狠!