%I#110 2022年1月23日12:10:37

%S 2,2,2,1,2,2,2,2,2,2,3,2,2,2,1,2,2,2,3,2,2,2,3,1,2,2,3,2,2,12,2,2，

%温度2,1,2,2,2,2,2,2,2,2,1,3,2,1,2,2,2,2,1,1,2,2,1,2,1,2,2，2,2,2,2，

%U 2,2,2,1,2,2,2,2,1,2,2,2,2,2,2,2,2,2,22,2,2,2，2,2,2,3,2,2,2

%N a（N）是正k的数量（可以大于N），因此A000005（N）/N=A000005。

%C第一个1出现在a（4）处。相应的k值为4。

%C前2出现在a（1）处。相应的k值为1、2。

%C前3出现在a（9）处。相应的k值为9、18、24。

%C前4出现在a（243）处。相应的k值为24340548610。

%C每个数字都出现在这个序列中吗？

%C前5个位于a（3189375），k值{3189375,5467500,6378750,6561000,8748000}，而前6个位于3176523，k值为{3176523,4084101,6353046,6806835,8168202,13613670}。如果出现7，则超过10^8。-_Charles R Greathouse IV，2021年10月5日

%H David A.Corneth，n的表，A（n）表示n=1..10000</a>

%H Charles R Greathouse IV，PARI/gp脚本</a>

%e对于n=9，A000005（9）/9=1/3=A000005k（k）/k的k为9、18和24。因此，a（9）=3。

%t数组[Function[r，Count[Range[Ciling[4/r^2]]，_？（DivisorSigma[0，#]/#==r&）]][Divisor西格玛[0，#]/#]&，105]（*或*）

%t块[{nn=7，m，s}，m=2^（2 nn）；s=KeySort@PositionIndex[Array[DivisorSigma[0，#]/#&，m]]；s=反向@KeyDrop[s，TakeWhile[Keys@s，4/#^2>m&]]；长度/@Array[Lookup[s，Divisor Sigma[0，#]/#]&，2^nn]]（*_Michael De Vlieger_，2021年10月4日*）

%o（PARI）a（n）=｛my（q=numdiv（n）/n）；sum（i=1，4/q^2，numdiv（i）/i==q）；｝\\_Michel Marcus_，2021年10月4日

%o（PARI）a（n）=my（q=numdiv（n）/n，s=分母（q），res=0）；对于步骤（i=s，4/q^2，s，如果（numdiv（i）==q*i，res++））；res\\_David A.Corneth，2021年10月7日

%o（巴黎）查看Greathouse链接

%Y参考A000005、A090387、A090395、A005382、A005833、A348184。

%K nonn公司

%O 1，1

%2021年10月4日，A_Tejo Vrush_