(岩浆)[1..1300]|IsZero中的n(DivisorSigma(2,n)mod Divisor西格玛(1,n))]//布鲁诺·贝塞利2013年4月10日
(哈斯克尔)
a020487 n=a020487_列表!!(n-1)
a020487_list=过滤器(\x->a001157 x`mod`a000203 x==0)[1..]
(Python)
从symy导入divisorsigma
定义ok(n):返回divisor_sigma(n,2)%divisor_sigma(n,1)==0
打印([k代表范围(11300)中的k,如果正常(k)])#迈克尔·布拉尼基2024年2月25日
(Python)#更快地生成序列的初始段
从数学导入prod
来自sympy导入因子
定义正常(n):
f=因子(n)
sigma1=prod((p**(e+1)-1)//(p-1)表示f中的p,e项())
σ2=prod((p**(2*e+2)-1)//(p**2-1)对于f中的p,e,items())
返回sigma2%sigma1==0
打印([k代表范围(11300)中的k,如果正常(k)])#迈克尔·布拉尼基,2024年2月25日
|