mathemtica 使用素数筛求素数

算小于N的质数列表1，创建从2到n的列表

2，令p=2，就是第一个指数

3，把大于P的P的倍数从列表里删除

4，p变成表中剩下的大于p的数

5， 重复3，4

 

?

n = 1000;

a = Range[2, n];

p = 2;

While[p < Last@a,

 a = Complement[a, Select[a, Mod[#, p] == 0 && # > p &]];

 p = First@Select[a, # > p &]]