26.C++中的Bitset埃拉托斯特尼_筛选法（查找质数或素数）及改进

Bitset是C++语言的一个类库，用来方便地管理一系列的bit位而不用程序员自己来写代码，

在 bitset 头文件中，它是一种类似数组的结构，它的每一个元素只能是０或１，每个元素仅用１bit空间。

什么是埃拉托斯特尼筛选法？
埃拉托色尼筛选法(the Sieve of Eratosthenes)简称埃氏筛法，是古希腊数学家埃拉托色尼提出的一种筛选法。 是针对自然数列中的自然数而实施的，用于求一定范围内的质数。
埃氏筛法步骤
先把N个自然数按照次序排列起来，
（1）然后把1删除（现今数学界1既不是质数也不是合数）
（2）读取队列中当前最小的数2，然后把2的倍数删去
（3）读取队列中当前最小的数3，然后把3的倍数删去
（4）读取队列中当前最小的数5，然后把5的倍数删去
（5）读取队列中当前最小的数7，然后把7的倍数删去
一直这样做下去，就会把不超过N的全部合数筛选掉，留下的就是不超过N的全部质数。

代码实现：

#include <iostream>
#include <bitset>
#include <time.h>

using namespace std;

int main()
{
	int const max_number(1000);
	bitset<max_number + 1> numbers;//初始化101个0！

	clock_t startTime, endTime;
	startTime = clock();//计时开始

	numbers.set();//set使101个0变为101个1
	numbers[1] = 0; //1既不是质数也不是合数，设为0
	for (int i = 1; i!=max_number; ++i)
	{  //筛选倍数，并赋值为0
		if (numbers[i])
		{
			for (int j = i*2; j < max_number+1; j+=i)
			{
				numbers[j] = 0;
			}
		}
	}
	cout << "the number of primes less than " << max_number
		<< " is " << numbers.count() << endl;  //计算bitset里面1的个数
	for (int i = 1; i != max_number; ++i)
	{
		if (numbers[i]) cout << i << ",";
	}

	endTime = clock();//计时结束
	cout << "The run time is: " << (double)(endTime - startTime) / CLOCKS_PER_SEC << "s" << endl;
	return 0;
}

修改后：
上面代码的效率不高，然后针对上面的效率进行修改:

#include <iostream>
#include <bitset>
#include <cmath>
#include <time.h>

using namespace std;

int main()
{
	int const max_number(1000);
	int const max_test((int)sqrt((double)max_number));
	clock_t startTime, endTime;

	startTime = clock();//计时开始

	bitset<max_number + 1> numbers;//初始化max_number+1个0
	numbers.set();//set使max_number+1变为1
	numbers[1] = 0; //1既不是质数也不是合数，设为0
	for (int i = 1; i != max_test; ++i)//例如，100以内的质数筛选到10就可以了，因为任何一个非质数，至少有一个因数小于他的平方根
	{  //筛选倍数，并赋值为0
		if (numbers[i])
		{
			for (int j = i * i; j < max_number + 1; j += i)   
			{
				numbers[j] = 0;
			}
		}
	}
	cout << "the number of primes less than " << max_number
		<< " is " << numbers.count() << endl;  //计算bitset里面1的个数
	for (int i = 1; i != max_number; ++i)
	{
		if (numbers[i]) cout << i << ",";
	}

	endTime = clock();//计时结束
	cout << "The run time is: " << (double)(endTime - startTime) / CLOCKS_PER_SEC << "s" << endl;
	return 0;
}