素数筛选法

简单素数筛选法

我们知道,判断一个数是否是素数,只需判断在其算术平方根之前是否有数能 被其整除,因此,简单的素数筛选法便是根据这一思想产生了,即对每一个数的算术平方根开始之前的数都进行判断。

代码如下:
#include <stdio.h>
#include <math.h>
#define MAX 100
int prime[MAX+1];
int main()
{
    int j,k=0;
    for(int i=2;i<=MAX;i++)
    {
        for(j=2;j<=sqrt(i);j++)
            if(i%j==0 )
                break;
       if(j>sqrt(i))
            prime[k++]=i;
    }
    for(int i=0;i<k;i++)
        printf("%d ",prime[i]);
    printf("\n");
    return 0;
}

快速素数筛选法

然而简单的素数筛选法的耗时比较长,尤其是求的素数范围比较大的时候,在oj做题时,基本上有简单素数筛选法的代码,都超时了。因此我就去学习一种快速筛选法。这种筛选法主要思想就是假如一个数是素数,那么它的2倍,3倍....n倍必定是合数。而且素数除2之外全部都是奇数,我们建立一个数组,令数组奇数下标等于0,数组偶数下标等于1,同时,在确定一个数是素数的时候,将其倍数下标的数也赋值为1,这样一来,只有数组下标等于0的数才是奇数,数组下标等于1的,不用考虑,因此省去了不少时间

代码如下:
#include <stdio.h>
#define MAX 100
int a[MAX+1];
int prime[MAX+1];
int main()
{
    int k=0;
    prime[k]=2;
    for(int i=3;i<=MAX;i++)
    {
        if(i%2==0)
            a[i]=1;
        else
            a[i]=0;
    }
    for(int i=3;i<=MAX;i+=2)
    {
        if(a[i]==0)
            prime[++k]=i;
        for(int j=i*2;j<=MAX;j+=i)
            a[j]=1;
    }
    for(int i=0;i<=k;i++)
        printf("%d ",prime[i]);
    printf("\n");
    return 0;
}

代码优化

快速筛选法虽然比较好,但是其仍然有一些计算是重复了,比如偶数下标的数不需要计算;但是我们仍然给其赋值;而且有的数甚至重复赋值,比如在计算3的倍数的时候,将15赋值为1,在计算5的倍数的时候,又将15赋值为1,这样就花费了不少时间。因此我们可以将下标为偶数的不予计算,只计算奇数下标的即可。另外,我们发现,当只计算奇数的时候,a数组中下标为0对应的则是3,下标为1时 对应的是5下标为2时,对应的是7......,即i*2+3的关系,另外,因为我们不对偶数进行计算,因此当筛选的倍数只能是3倍,5倍....

代码如下:

#include <stdio.h>
#include <math.h>
#define MAX 100
int a[MAX+1];
int prime[MAX+1];
int main()
{
    int k=0;
    prime[k]=2;
    for(int i=0;i*2+3<=MAX;i++)
        a[i]=0;
    for(int i=0;i*2+3<=MAX;i++)
    {
        if(a[i]==0)
            prime[++k]=i*2+3;
        for(int j=3*(i*2+3);j<=MAX;j+=2*(i*2+3))
            a[(j-3)/2]=1;
    }
    for(int i=0;i<=k;i++)
        printf("%d ",prime[i]);
    printf("\n");
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值