举例:筛选从1—n的所有素数
- 第一种
暴力判断,无优化。直接根据素数的定义判断即可
#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);
}
- 第二种
素数的倍数一定不是素数,对其进行标记。最后遍历一遍数组,可得未被标记的即为筛选出的素数。
- 比第一种快了一些,但仍有可以改进的地方,因为在标记的时候有重复标记,因此这个地方也可以再优化一下
#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);
}
- 第三种
线性筛法:不选重复过的数,大大加快了代码运行时间
注意这个分解:任意一个合数都可分解为素数相乘,如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]);
}
** 表格会更加直观一点 **
_ 第三种方法最快 _