10亿素数算法

1. //计算10亿需要60MB内存,比普通的bit位算法节省一半内存
2.  unsigned int PrimeCounter(unsigned int MaxNumber)
3. {
4.     const unsigned MachineLength = sizeof(unsigned int) << 3;
5.     const unsigned int MaxSize = MaxNumber / (MachineLength<<1) + 1;
6.     unsigned int* Number = new unsigned int[MaxSize];
7.     memset(Number,0xFF,MaxSize*sizeof(unsigned int));
8.     for(unsigned long long num = 3;num<MaxNumber;num+=2)
9.     if(Number[num/(MachineLength<<1)] & (0x1L<<((num>>1) % MachineLength))) {
10.         for(unsigned long long NotPrime = num * num;NotPrime < MaxNumber ;NotPrime += num ) 
11.             if(NotPrime % 2)
12.                 Number[NotPrime / (MachineLength<<1)] &= ~(0x1L << ((NotPrime>>1) % (MachineLength)));               
13.         }
14.     unsigned int counter = 1;
15.     for(unsigned int num = 3;num < MaxNumber;num +=2)
16.         if( Number[num /(MachineLength<<1)] & (0x1L<<(num>>1 % (MachineLength)))){
17.             //cout << num << "/t";
18.             counter ++;
19.         }
20.     delete []Number;
21.     return counter;
22. }
23. int main()
24. {
25.     clock_t pretime = clock();
26.     std::cout << "/nPrimeCounter(10000000)= "<<PrimeCounter(10000000) << std::endl;
27.     std::cout << clock() - pretime <<"ms has used.";
28.     std::cout <<std::endl;
29.     system("pause");
30.     return 0;
31. }