（PARI）平滑（P:vec，lim）=我的（v=列表（[1]），nxt=向量（#P，i，1），indx，t）；而（1，t=vecmin（向量（#P，i，v[nxt[i]]*P[i]），&indx）；如果（t>lim，断裂）；如果（t>v[#v]，则列表输入（v，t））；nxt[indx]++）；Vec（v）；
ok（p）=p=p%24；p==1|p==5|p==7|p==11；
list（lim）=my（u=列表（[1]），v=[5]，r，t）；对于素数（p=7，如果（ok（p），如果（v[v]*p>lim，break））；v=连接（v，v[#v]*p））；v=光滑（v，lim）；对于（i=2，#v，t=2013年12月（v[i]）；如果（t>r，r=t；列表输入（u，v[i]））；Vec（u）（车辆）
