C++学习之【使用位操作符求素数分析】

#include <iostream>
using namespace std;

void getPrime_1()
{
const int MAXN = 100;
bool flag[MAXN];
int primes[MAXN / 3 + 1], pi=0;
primes[pi++]=2; //2 是一个素数，先记录下来
int i, j;
for(i=0;i<MAXN;++i)flag[i]=false;//全部置假，没访问一个，相应位置真
for (i = 3; i < MAXN; i+=2){//大于2的偶数一定不是素数，所以只要判断奇数即可
if (!flag[i])//如果是素数
{
primes[pi++] = i;
for (j = i; j < MAXN; j += i)//i的倍数一定都不是素数
flag[j] = true;
}
}
for(i=0;i<pi;++i)
cout<<primes[i]<<" ";
cout<<endl;
}
这里使用了素数表，每一个bool型占用1个字节，共8位二进制位。而且这里除了多用了很多的无用bool变量在flag数组里，可以看到，我们使用的flag数组只用到2号位之后的所有奇数位。因此，这里可以进行压缩改进一下：

将flag数组减少一半
使用位操作符使空间占用减少为原来的八分之一
这里使用到对指定位置置1的操作：对于一个整数X可以通过将1左移n位后，与X进行或操作，使X的第n位置1。

int X=0;
int n=10;
X |= 1<< n; // 将 X 的第 n 位置 1
所以根据上面两条，优化后的代码如下：

void getPrime_2()
{
const int MAXN = 200;
const int BitN=(MAXN/2)/32+1;//
int flag[BitN];
int primes[MAXN / 3 + 1], pi=0;
primes[pi++]=2;  //2 是一个素数，先记录下来
int i, j;
for(i=0;i<BitN;++i)flag[i]=0;//全置0，每访问过一个，相应位置1
for (i = 3; i < MAXN; i+=2){ //大于2的偶数一定不是素数，所以只要判断奇数即可
if (!((flag[(i/2) / 32] >> ((i/2) % 32)) & 1))
{
primes[pi++] = i;
//i的倍数一定都不是素数，其中，j加上一个i后为偶数，上一级已经不考虑了，所以还要加上一个i
for (j = i; j < MAXN; j =j+i+i)
flag[(j/2) / 32] |= (1 << ((j/2) % 32));
}
}
for(i=0;i<pi;++i)
cout<<primes[i]<<" ";
cout<<endl;
}
首先，根据最大数，判断需要32的整型多少个：