素数筛法

来源：互联网发布：量化派人工智能编辑：程序博客网时间：2024/05/09 04:54

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

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

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

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

我的代码主要是用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;}