Решето Ератосфена

#include <iostream>
#include <ctime>
 
using namespace std;
 
template<int N>
 
class bitset
{
private:
   char* bits;
public:
    bitset() : bits(new char[(N - 1) / 8 + 1]){}
    bool test(int n) 
    { 
        return (bits[n >> 3] & (1 << (n & 7))) != 0; 
    }
    void set(int n) 
    { 
        bits[n >> 3] |= 1 << (n & 7); 
    }
    void reset(int n) 
    { 
        bits[n >> 3] &= ~(1 << (n & 7));
    }
};
 
int main()
{ 
   const long N = 10000000;
   clock_t cstart = clock();
   bitset<N + 1> b;
   int count = 0;
   int i;
   for (i = 2; i <= N; i++)
      b.set(i);
   i = 2;
   while (i * i <= N)
   { 
      if (b.test(i))
      { 
         count++;
         int k = 2 * i;
         while (k <= N)
         { 
             b.reset(k);
             k += i;
         }
      }
      i++;
   }
   
   while (i <= N)
   {
      if (b.test(i))
         count++;
      i++;
   }
 
   clock_t cend = clock();
   double millis = 1000.0 * (cend - cstart) / CLOCKS_PER_SEC;
 
   cout << count << " primes\n" << millis << " milliseconds\n";
   system("pause");
   return 0;
}

bool test(int n) 
{ 
 int tmp1=n/8;
 int tmp2=n%8;
 unsigned char bt1=bits[tmp1];
 unsigned char bt2=0x01<<tmp2;
 unsigned char bt3=bt1&bt2;
 bool bl=bt3!=0;
 return bl;
}