使用篩法找出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