计算整数N 以内所有素数——埃拉托色尼筛法

问题：给定正整数N(N>1)，求N以内所有素数

埃拉托色尼筛法利用合数质因数分解原理，首先将所有数假设为素数，标为1，然后，将已知素数的倍数(非素数)标记为0，最后，所有标为1的数即为素数。代码如下：

#include <stdio.h>
#include <stdlib.h>

int main(int argc, char **argv)
{
	int n = 0;
	scanf("%d",&n);
	int *num = malloc(sizeof(int) * n);
	if (num == NULL)
	{
		printf("memory error\n");
		return 0;
	}
	
	int i = 0, j = 0;
	for (i = 0; i < n; i++)
	{
		num[i] = 1;
	}
	
	for (i = 2; i < n; i++)
	{
		if (num[i])//素数的整数倍不是素数
		{
			for (j = i; i * j < n; j++)
			{
				num[i*j] = 0;
			}			
		}
		
	}
	for (i = 2; i < n; i++)
	{
		if (num[i])
		{
			printf("%d ",i);
		}
		
	}
	
	
	return 0;
}

改进方案：

采用int型数组，会占用太多空间，可以改为位数组。

#include <stdio.h>
#include <stdlib.h>
#include <limits.h>
#define N 10000010


#define BITMASK(b) (1 << ((b) % CHAR_BIT))
#define BITSLOT(b) ((b) / CHAR_BIT)
#define BITSET(a, b) ((a)[BITSLOT(b)] |= BITMASK(b))
#define BITCLEAR(a, b) ((a)[BITSLOT(b)] &= ~BITMASK(b))
#define BITTEST(a, b) ((a)[BITSLOT(b)] & BITMASK(b))
#define BITNSLOTS(nb) ((nb + CHAR_BIT - 1)/CHAR_BIT)

char num[BITNSLOTS(N)];

int main(int argc, char **argv)
{
	int n = 0;
	scanf("%d",&n);
	
	int i = 0, j = 0;
	for (i = 0; i < n; i++)//初始化
	{
		BITSET(num,i);
	}
	
	for (i = 2; i < n; i++)
	{
		if (BITTEST(num, i))//素数的整数倍不是素数
		{
			for (j = 2; i * j < n && i * j < N; j++)
			{
				BITCLEAR(num, i*j);
			}			
		}
		
	}
	
	for (i = 2; i < n; i++)
	{
		if (BITTEST(num, i))
		{
			printf("%d ",i);
		}
		
	}
	
	return 0;
}