简单素数筛

 1 void find_prime()
 2 {
 3     mem(is_prime, 1), cnt = 0;//初始化为全是素数 
 4     for(int i = 2; i <= 10000; i ++)
 5     {
 6         if(is_prime[i])//注意不能写成is_prime[i] == 1 因为mem()是十六进制的1 
 7         {
 8             prime[++ cnt] = i;
 9             for(int j = i; j <= 10000; j += i)//所有倍数都不是素数 
10                 is_prime[j] = 0;  //筛为不是素数 
11         }
12     }
13 }