|
（*列表minDiv[n，b]中的第d个元素是最小的数字k<=n，正好有d<=b个除数，否则为零。只要任一不等式失败，计算就会停止。*）
最小刻度[n_，b_]：=
模块[{list=Array[0&，b]，k=1，d}，
而[k<=n，d=DivisorSigma[0，k]；
如果[d<=b&&list[[d]]==0，list[[d]]=k]；
如果[d<=b，k++，k=n+2]]；列表]
a256259【n，b】：=
映射[DivisorSigma[1，#]&，排序[Select[minDiv[n，b]，#！=0 &]]]
a256259[100000，300]（*计算序列的前60个元素*）
| 