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