素数筛法

我也蠢蠢欲动，写个素数筛法。有一些问题必须先存储一个素数表，然后查表。

这时就要用到素数筛。具体的筛法见代码，具体原理这两篇文章有很好的介绍：

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

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

我的代码主要是用bit做合数标记，在素数范围较大时能够有效地节省内存。

#include <stdio.h>
#include <math.h>

#define  MAXN  1000000
int  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;
}