%I#41 2024年2月19日01:56:50

%S 1,4,5,15,42,44,4712133634034793025486937694751322379097，

%电话：37913137913328012052069834526409056264197152941920152942012，

%电话：8350344420617011639545591337939545913379831123930105069433687695336033687693812

%N对m进行编号，使m除以和{k=1..m}A000005（k）。

%D Julian Havil，“伽马：探索欧拉常数”，普林斯顿大学出版社，普林斯顿和牛津，第112-113页，2003年。

%H Donovan Johnson，n的表，n=1..39的a（n）（商<=40）

%如果和{i=1..m}d（i）=m*k，k是整数，其中d（i。

%e当k=15时，总和为1+2+2+3+2+4+2+4+3+4+4+2+6+2+4+4=45，可被15整除。

%t s=0；Do[s=s+除数西格玛[0，n]；如果[Mod[s，n]==0，打印[n]]，{n，1，2*10^9}]

%t k=10^6；a[1]=1；a[n]：=a[n]=除数Sigma[0，n]+a[n-1]；nd=a/@范围@k; 选择[范围@k，可分割[nd[[#]]，#]&]（*_Ivan N.Ianakiev_，2016年4月30日*）

%t模块[{nn=400000}，选择[Thread[{Range[nn]，Accumulate[DivisorSigma[0，Range[n]]}]，Divisible[#[2]]，#[[1]]&]][[All，1]]（*程序生成序列的前19个项。要生成更多项，请增加nn常数。*）（*_Harvey P.Dale_，2022年7月3日*）

%o（PARI）列表a（nn）={my（s=0）；对于（n=1，nn，s+=numdiv（n）；如果（！（s%n），打印1（n，“，”）；）；}\\_Michel Marcus_，2015年12月14日

%o（鼠尾草）

%o定义A050226_列表（长度）：

%o a，L=0，[]

%o表示（1..len）中的n：

%o a+=σ（n，0）

%o如果n.divides（a）：L.append（n）

%o返回L

%o A050226_list（10000）#_Peter Luschny_，2015年12月18日

%Y参考A000005、A006218、A057494、A085567、A085829。

%K nonn很好

%O 1,2号机组

%A _Labos Elemer，1999年12月20日

%E更多条款摘自2000年9月21日_Robert G.Wilson v_

%E Nomoto的进一步条款，2001年8月3日

%E a（26）-a（30）摘自Donovan Johnson_，2008年12月21日