(Thomas Baruchel的计划项目);(定义素数'(2 3 5 7…999983))(计算n)返回A088028型(n) (如果素数列表太短,则为f)计算n<=16需要合理的时间(对于较小的n值,比“野蛮”方法慢,但很快就会变得更快)。结果被证明是最小的。
(define(compute*n mmax prod offset)(do((i offset(+i 1)))(l(length primes))((>=(*prod(do)((j 0(+j 1))(p 1))s(不精确->精确(楼层(sqrt(+p 1)))))))(如果(=(*s)(+p1))(设置!mmax p))))
(define(compute n)(let*(p(reverse(cdr primes)))(mmax(apply*(cons(car p)(list-tail p(-(length p)(-n 1))