素数--超强筛选法

 

一般筛选法求素数：用一个bool 数组Prime[size]，筛掉所有的合数，剩下的即为素数。

在这里向大家介绍一种优化后的“筛选法”。这是我从浙江师范大学ACM网站上学来的，跟大家分享一下。http://acm.zjnu.cn/show.asp?tab=arithmetic&id=34

优化：

bool 数组只存奇数(3,5,7,9,11,13.....2*i+3...)，这样筛选量就相对减小了，但这还不是超强之处！筛选过程才强！

易知，每个奇数跟其下标的关系为2*i+3（如3=2*0+3，5=2*1+3）；当我们选定3是素数时，它的倍数9，15，21....就要赋为false。

如何在只知道下标的情况下找到它的倍数呢！？！9的下标是3，15的下标是6...都是相差3；5的倍数的下标之间相差5....           15既是3的倍数也是5的倍数，在找到3时已经赋为false了，我们就不必再给它赋值了，但是如何跳过15（下标i==6）直接到25（i==11）去呢！？！......

外围循环为i时，内层循环的筛法是从 i+(2*i+3)*(i+1)     即i*(2*i+6)+3开始筛的。11=1*(2*1+6)+3,由5(i=1)跳到25(i=11)。

#include<iostream> #include<cmath> using namespace std; const int size=100000000; const int size2=50000001; bool Prime[size2]; void getPrime(){ int sq=(int)sqrt(size),i,j; for(i=0;i<size2;i++)Prime[i]=true; //只存奇数 Prime[0] =3; //a[i]=2*i+3 int temp; for(i=0;i<=sq;i++){ if(Prime[i]){ temp=2*i+3; for(j=i*(i*2+6)+3;j<size2;j+=temp) //起点 j=i+temp*(i+1) Prime[j]=false; } } }