素數篩法

我也蠢蠢欲動，寫個素數篩法。有一些問題必須先儲存一個素數表，然後查表。

這時就要用到素數篩。具體的篩法見代碼，具體原理這兩篇文章有很好的介紹：

http://blog.aikilis.com/1299.html

http://blog.aikilis.com/1287.html

My Code主要是用bit做合數標記，在素數範圍較大時能夠有效地節省記憶體。

#include <stdio.h>#include <math.h>#define  MAXN  1000000int  f[MAXN/32+1];int  prime[87000];  // 1.2 * MAXN / ln( MAXN )int  cnt = 0;#define  TEST( n )  ( f[n>>5] & (1<<(n&31) ) )#define  SET( n )  ( f[n>>5] |= (1<<(n&31)) )void  sieve( int x ) {        int  rt = sqrt( (double)x ) + 1;        prime[cnt++] = 2;        for( int i = 3; i <= x; i += 2 ) {                if( !TEST(i) ) {                        prime[cnt++] = i;                        if( i < rt ) for( int j = i * i; j <= x; j += i ) SET( j );                }        }}int main() {        sieve( 1000000 );        printf( "cnt = %d\n", cnt );        return 0;}