判断+筛选 素数

举例:筛选从1—n的所有素数
  1. 第一种

暴力判断,无优化。直接根据素数的定义判断即可

#include<stdio.h>
int main()
{
    int n=1000,flag=1;   
    int prime[25]= {0};
    for(int i=2; i<=n; i++)
    {
        flag=1;
        for(int j=2; j<i; j++)
            if(i%j==0)
                flag=0;
        if(flag)
            prime[i]++;
    }
    for(int i=1; i<=n; i++)
        if(prime[i])
            printf("%d\n",i);
}

  1. 第二种

素数的倍数一定不是素数,对其进行标记。最后遍历一遍数组,可得未被标记的即为筛选出的素数。

  • 比第一种快了一些,但仍有可以改进的地方,因为在标记的时候有重复标记,因此这个地方也可以再优化一下
#include<stdio.h>
int prime[1100]={0};  //初始化0,假设全部是素数,标记为0;
int main()
{
    int n=1000,cnt=0;    //cnt计数,表示下标
    for(int i=2;i<=n;i++)
    {
        if(prime[cnt++]==0)//这里如果用a[i]=i的话,就会将原有的标记值改变,所以需要cnt来改变;
            prime[cnt++]=i;
        for(int j=i*2;j<=n;j+=i)   //j从i的两倍开始,每次j+=i,就是以此类推,如i的3倍,4倍……
        {
            prime[j]=1;    //对倍数进行标记
        }
    }
    for(int i=1;i<=n;i++)
        if(prime[i]==0)
        printf("%d\n",i);
}


  1. 第三种

线性筛法:不选重复过的数,大大加快了代码运行时间

注意这个分解:任意一个合数都可分解为素数相乘,如4=2x2;12=2x2x3;15=3x5 ……

  • 对2 ~ n的数进行遍历,那么就有两种情况。
  • 1.当i 为合数时,那么就可以将i = p1 * p2 * p3……,(p1<p2<p3<……),此时可知p1是最小质因数。我们就可以筛出i* p1 的数。
  • 2.当i 为素数时,因数只有1和它本身,无法进一步分解。此时代码中的cnt(下图代码) 不仅可以将数依次存入prime数组,也可以记录素数的个数,那么就有cnt个素数,此时筛出i *cnt1,cnt2,cnt3……

#include<stdio.h>
int main()
{
    int check[1000]= {0};
    int prime[1000]= {0};
    int n=1000,cnt=0;    //cnt同上,计数,表示下标
    for(int i=2; i<=n; i++)  //对2-n的数进行遍历
    {
        if(!check[i])   //如果是素数,就记录进prime数组里
            prime[cnt++]=i; //同时,cnt也记录了素数的个数
        for(int j=0; j<=cnt&&i*prime[j]<=n; j++) //筛选合数
        {
            check[i*prime[j]]=1;   
            if(i%prime[j]==0)  //对于合数,筛出它与最小质因数的乘积即可
                break;
        }
    }
    for(int i=0; i<=n; i++)
        if(prime[i])
            printf("%d\n",prime[i]);
}

** 表格会更加直观一点 **

在这里插入图片描述

_ 第三种方法最快 _

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值