(PARI)a(n)=如果(n<0,0,n=a(n-1));直到(0==(σ(n,0)*n)%σ(n,1),n++);n)/*迈克尔·索莫斯2004年2月6日*/
(哈斯克尔)
导入数据。比率(分母)
导入数据。列表(genericLength)
a001599 n=a001599_列表!!(n-1)
a001599_list=过滤器((==1)。分母。hm)[1..]其中
hm x=genericLength ds*recip(总和$map(recip.fromIntegral)ds)
其中ds=a027750_低x
(GAP)级联([1],过滤([2,4..2000000],n->IsInt(n*Tau(n)/Sigma(n))); #穆尼鲁·A·阿西鲁2018年11月26日
(Python)
从symy导入divisorsigma到sigma
定义ok(n):返回(n*σ(n,0))%sigma(n,1)==0
打印([n代表范围(1,10**4)中的n,如果正常(n)])#迈克尔·布拉尼基2021年1月6日
(Python)
从itertools导入计数,islice
从functools导入reduce
从数学导入prod
来自sympy导入因子
定义A001599号_gen(startvalue=1):#术语生成器>=startvalue
对于计数中的n(max(startvalue,1)):
f=因子(n)
s=产品((p**(e+1)-1)//(p-1)对于f项中的p,e())
如果不减少(λx,y:x*y%s,(e+1代表f.values()中的e),1)*n%s:
产量n