(鼠尾草)
(圣人)定义A057949号_列表(绑定):
.... numterms=(bound-1)//4+1
.... M=[1]*数字项
.... 对于范围(1,数字项)中的k:
........ 如果M[k]==1:
............ k功率=k
............ 而kpower<数字项:
................ 步长=4*kpower+1
................ 对于范围内的j(kpower、numterms、step):
.................... M[j]*=4*k+1
................ kpower=4*kpower*k+kpower+k
....##现在M[k]包含p^e项的乘积,其中p是S-素数
....##e是最大的,因此p^e除以4*k+1
.... 如果M[k]>4*k+1],则返回[4*k+1表示范围内的k(数字项)
结束# _# _埃里克·施密特(Eric M.Schmidt),2016年12月11日
|