%I#31 2022年4月9日23:01:57

%S 1,3,4,5,6,7,8,9,10,12,13,14,15,17,18,20,21,24,26,28,30,31,32,33,36，

%电话：38、39、40、42、44、48、50、54、56、57、60、62、63、65、68、72、73、74、78、80、82、84、85，

%电话：90,91,93,96,98102104108110112114120121122

%N数N，使N=σ_k（m），对于某些k>=1。

%C Sigma_k（n）=总和[d|n，d^k]。

%C Sigma_0（n）可以是任何正整数，因此在此序列中被忽略。

%C该序列的渐近密度为0（Niven，1951，Rao和Murty，1979）_Amiram Eldar，2020年7月23日

%H Giovanni Resta，n的表，n=1..10000的a（n）</a>

%H Ivan Niven，<a href=“https://doi.org/10.1090/S0002-9904-1951-09543-9“>序列的渐近密度，美国数学学会，第57卷（1951年），第420-434页。

%H R.Sita Rama Chandra Rao和G.Sri Rama Chandri Murty，<a href=“https://doi.org/10.4153/CBM-1979-018-5“>关于Niven的一个定理，《加拿大数学公报》，第22卷，第1期（1979年），第113-115页。

%H Eric W.Weisstein，<a href=“http://mathworld.wolfram.com/DivisorFunction.html“>MathWorld:除数函数</a>

%e Sigma_2（4）=1+4+16=21，因此21在序列中。

%t高达[n_]：=选择[扁平接头[{1，除数Sigma[最大范围[1,地板@原木[#，n]]，#]&&@范围[2，n]}]，#<=n&&]；截至[122]（*Giovanni Resta_，2013年2月5日*）

%o（PARI）列表（lim）=如果（lim<3，返回（如果（lim<1，[]，[1]））；my（v=列表（[1]））；对于（k=1，logint（（lim=1）-1,2），forfactored（m=2，sqrtnint（lim-1，k），my（t=sigma（m，k））；如果（t<=lim，列表输入（v，t））；集（v）\\_Charles R Greathouse IV_，2022年4月9日

%Y参考A0000203、A001157、A001158、A001159、A001160。

%Y参见A013954、A013955、A013956、A013957、A013958、A013959、A013960、A013961、A013962、A013963、A013964、A013965、A013966、A013967、A013968、A013969、A013970、A013971、A013972。

%Y参考A000005。

%K nonn公司

%O 1,2号机组

%A _乔恩·佩里（Jon Perry），2013年2月5日