线性筛法求质数(素数)表 及其原理(By OWenT)

代码(c++): (By OWenT)

 1#include<iostream>
 2using namespace std;
 3const long n=200000;
 4long prime[n]={0},num_prime=0;
 5int main()
 6{
 7    int a[n]={1,1},i,j;
 8    for(i=2;i<n;i++)
 9    {
10        if(!a[i])
11            prime[num_prime++]=i;
12        for(j=0;j<num_prime && i*prime[j]<n;j++)
13        {
14            a[i*prime[j]]=1;
15            if(!(i%prime[j]))
16                break;
17        }
18    }
19    return 0;
20}
21

线性筛法,即是筛选掉所有合数,留下质数
我们知道合数可以由一个质数数与另一个数相乘得到
而同时假设合数a=质数b×质数c×一个数d
令e=c × d,假设b ≥ e,e为合数,令f=d × b
a=f × c ,其中c
即比一个合数数大的质数和该合数的乘积可用一个更大的合数和比其小的质数相乘得到
这也是if(!( i % prime[j]))break;的含义,这也是线性筛法算质数表的关键所在