今天看到一個素數的算法，利用素數集中在6x+1/6x-1的原理
整理了一下 算法如下 打印max以下num個素數：
http://lmlf001.blog.sohu.com/
void print_prime(int max,int num)
{
if(max<5)return;
unsigned long t,m,n,p;
unsigned long x;
int i,j,a,b,k;
x=max;
int z=0;
x=x/6; //x>=5
for(t=x;t>=0,z<num;t--)
{
i=1;j=1; k=t%10;
m=6*t;
/**i,j的值 是是否進行驗證的標志也是對應的6t-1和6t+1的素性標志**/
if(((k-4)==0)||((k-9)==0)||((m+1)%3==0))j=0;
/*此處是簡單驗證6*t-1,6*t+1 是不是素數，借以提高素數純度**/
if(((k-6)==0)||((m-1)%3==0))i=0;
/***先通過初步判斷去除末尾是5，及被3整除的數***/
for(p=1;p*6<=sqrt(m+1)+2;p++ )
{
n=p*6;
/**將6*p-1和6*p+1看作偽素數來試除*****/
k=p%10;
a=1;b=1;
/**同樣此處a,b的值也是用來判斷除數是否為素數提高除數的素數純度**/
if(((k-4)==0)||((k-9)==0))a=0;
if(((k-6)==0))b=0;
if(i){
/*如果i非零就對m-1即所謂6*t-1進行驗證，當然還要看除數n+1,n-1,素性純度*/
if(a){if((m-1)%(n+1)==0)i=0;} /***一旦被整除就說明不是素數故素性為零即將i 賦值為零***/
if(b){if((m-1)%(n-1)==0)i=0;}
}
if(j){
/**如果j非零就對m+1即所謂6*t+1進行驗證，當然還要看除數n+1,n-1,素性純度*/
if(a){if((m+1)%(n+1)==0)j=0;} /***一旦被整除就說明不是素數故素性為零即將j 賦值為零***/
if(b){if((m+1)%(n-1)==0)j=0;}
}
if((i+j)==0)break;
/**如果已經知道6*t-1,6*t+1都不是素數了那就結束試除循環***/
}
if(j){cout<<m+1<<endl;z++;if(z>=num)break;}
if(i){cout<<m-1<<endl;z++;if(z>=num)break;}
}
}