素數的查找

使用篩法找出4萬內的素數

因爲2,3分別是素數中最小的偶數和奇數,所以只要是2或者3的倍數的都不是素數,剩下的則是素數。


#define MAXN 40000
int SPrime(void)
{
     int i, j;
     int b_size;
     int prime[MAXN], b[MAXN];
	
     for (i = 2 ; i < MAXN; i ++ )
 	prime[i] = 1 ;
     for (i = 2 ; i < MAXN; i ++ )
     {
        if (prime[i] == 1 )
	{
             for (j = 2 ; i * j < MAXN; j ++ )
                prime[i * j] = 0 ;
	}
    }
     for (i = 0 , j = 0 ; j < MAXN; j ++ )
     {
         if (prime[j] == 1 ) 
	{
            b[i] = j;
            i ++ ;
        }
    }
    b_size = i;
    return b_size;
}

另一種方法使用位圖法

#define MAX 40000

/*
i/8 第i/8個字節裏
i%8 第i%8位
*/

size_t count_prime()
{
	char a[MAX/8 +1] ={0};
	size_t count = 0;
	for(int i = 2; i < MAX; i++)
	{
		a[i/8] |= (1<<(i%8));
	}

	for(int i = 2; i 




点赞