#linear_prime_sieves
#prime[]，存储n以内所有的素数，其index为pi，初值为0
#is_prime[i]，表示自然数i(i<=n)是否质数，False就不是
set is_prime[] to true
for i=2 to n
if is_prime[i]=true then prime[pi++]=i
for j=0 to pi-1
if prime[j]*i>n then exit loop_j
is_prime[prime[j]*i]=false
if i mod prime[j]=0 then exit loop_j
endif
endif