Python和C++求10000以内的所有素数（Eratosthenes筛选 求素数）

Eratosthenes筛选求素数

“除了自身之外，无法被其它整数整除的数”称之为质数，要求质数很简单，但如何快速的求出质数则一直是程序设计人员与数学家努力的课题，在这边介绍一个著名的 Eratosthenes求质数方法。

解法

首先知道这个问题可以使用循环来求解，将一个指定的数除以所有小于它的数，若可以整除就不是质数，然而如何减少循环的检查次数？如何求出小于N的所有质数？

首先假设要检查的数是N好了，则事实上只要检查至N的开根号就可以了，道理很简单，假设A*B = N，如果A大于N的开根号，则事实上在小于A之前的检查就可以先检查到B这个数可以整除N。不过在程序中使用开根号会精确度的问题，所以可以使用 i*i <= N进行检查，且执行更快。

再来假设有一个筛子存放1～N，例如：

2　3　4　5　6　7　8　9　10　11　12　13　14　15　16　17　18　19　20　21 ........ N

先将2的倍数筛去：

2　3　5　7　9　11　13　15　17　19　21 ........ N

再将3的倍数筛去：

2　3　5　7　11　13　17　19　........ N

再来将5的倍数筛去，再来将7的质数筛去，再来将11的倍数筛去........，如此进行到最后留下的数就都是质数，这就是Eratosthenes筛选方法（Eratosthenes Sieve Method）。

检查的次数还可以再减少，事实上，只要检查6n+1与6n+5就可以了

Python代码：

def findAllprime(n):
     pt=[True]*n
     prime=[]
     for p in range(2,n):
          if not pt[n]:continue
          prime.append(p)
          for i in range(p*p,n,p)
              pt[i]=False
     return prime , pt

prime,pt=findAllprime(10000)

C++代码：

const int N = 10000;  
bool arrPrimeTable[N + 1];  
int nCount = 0;  
   
void InitPrimeTable()  
{  
    //初始化  
    for (int i = 2;i <= N;i++)  
    {  
        arrPrimeTable[i] = true;  
    }  
    //进行剔除  
    for ( i = 2;i * i <= N;i++)  
    {  
        if (arrPrimeTable[i])  
        {  
            //剔除i的倍数：t = x * i  
            for (int t = i + i;t <= N;t += i)  
            {  
                arrPrimeTable[t] = false;  
            }  
        }  
    }  
    //输出  
    for ( i = 2;i <= N;i++)  
    {  
        if (arrPrimeTable[i])  
        {  
            cout<<i<<" ";  
            nCount++;  
            if(nCount%10==0)
                cout<<endl;
        }  
    }  
   // cout<<endl<<"总计: "<<nCount<<" 个"<<endl;  
}  

或者可以用这种方法，质数只能在6的左右，

即：num % 6 == 1 或者num % 6 == 5

证明，反证法：如果num%6可以等于2,3,4

6x+2,6x+3,6x+4可以提取2或3

即2（x+1）不是质数

bool isPrime(int num)  
{  
    if (num == 2 || num == 3)  
    {  
        return true;  
    }  
    if (num % 6 != 1 && num % 6 != 5)  
    {  
        return false;  
    }  
    for (int i = 5; i*i <= num; i += 6)  
    {  
        if (num % i == 0 || num % (i+2) == 0)  
        {  
            return false;  
        }  
    }  
    return true;  
}