Algorithm,Number Theory,Prime

/***************************************************/
// 测素数,根号阶
bool is_prime(int u)
{
     if(u == 0 || u == 1) return false;
     if(u == 2)      return true;
     if(u%2 == 0)      return false;
     for(int i=3; i <= sqrt(u) ;i+=2)
          if(u%i==0)     return false;
     return true;
}

/***************************************************/
// 线性筛素数
const int M = 1000; // M : size
bool mark[M]; // true : prime number

void sieve_prime()
{
     memset(mark, true, sizeof(mark));
     mark[0] = mark[1] = false;// 0 not used, 1 is not prime
     for(int i=2; i <= sqrt(M) ;i++) {// 2-> 4, 6, 8, ... //3->9, 12, 15...// 5->25, 30, ...
          if(mark[i]) {
               for(int j=i*i; j < M ;j+=i)
                    mark[j] = false;
          }
     }
}