短小精悍的线性时间素数筛法

输入n，求n以内的所有素数

算法用两个数组存储数据:
一个是prime[]，存储n以内所有的素数，其index为pi，初值为0
一个是is_prime[i]，表示自然数i(i<=n)是不是质数。
算法如下：
linear_prime_sieves
1:  set is_prime[] to true
2:  for i=2 to n
3:    if is_prime[i]=true then prime[pi++]=i
4:    for j=0 to pi-1
5:      if prime[j]*i>n then exit loop_j
6:      is_prime[prime[j]*i]=false
7:      if i mod prime[j]=0 then exit loop_j
8:    endif
9:  endif

    这个算法有两层循环，第一层遍历2到n之间的所有自然数i，看看它是不是质数，如果是，则把i放进prime数组中。第二层循环是对所有未来的数进行筛选。对于当前正在处理的i，显然它乘以任何一个已经找到的素数的结果肯定是和数，它们将会被剔除。整个算法最核心的一句是第7行：当i可以被某个已经找到的质数整除的时候，循环退出，不再进行剔除工作。这样做的原因是：当prime[j]是i的因子的时候，设i=prime[j]*k，首先，我们可以肯定的说，prime[j]是i的最小质因数，这是因为第二层循环是从小到大遍历素数的；其次，我们可以肯定的说，i已经无需再去剔除prime[j']*i (j'>