%I#21 2022年9月8日08:46:23

%S 0,1,1,4,1,0,11,5,11,9,13,13,26,10,10,1,8,17,1,16,7,19,18,0,1，

%T 28,1,57,19,23,22,34,1,25,23,24,41,1,65,45,29,1,57,9,68,25,75,1,39，

%U 25,25,29,35,1,88,1,37,74120,29,37,1,91,31,24,1103单位

%N a（N）=和{d|N}（σ（N）mod d）。

%H Antti Karttunen，n表，n=1..65537的a（n）</a>

%H Carlos Rivera，<a href=“https://www.primepuzzles.net/puzzles/puzz_1065.htm“>拼图1065。大于45的整数，因此</a> ，主要困惑和问题联系。

%F a（A007691（n））=0。

%F a（A000040（n））=1。

%F a（A008578（n））=τ（n）-1。

%对于数字4、45、6048、14421…，F a（n）=n。。。

%e当n=4时；a（4）=（7模1）+（7模2）+（7mod 4）=0+1+3=4。

%t a[n_]：=块[{s=DivisorSigma[1，n]}，DivisorSum[n，Mod[s，#]&]]；阵列[a，72]（*Giovanni Resta_，2018年9月7日*）

%o（岩浆）[&+[SumOfDivisors（n）mod d:d in Divisor（n）]：n in[1..1000]]

%o（PARI）a（n）=我的（sn=σ（n））；sumdiv（n，d，sn%d）；\\_米歇尔·马库斯，2018年9月7日

%o（Python）

%o来自sympy导入除数

%o定义a（n）：divs=除数（n）；s=总和（divs）；返回总和（divs中d的s%d）

%o打印（[a（n）代表范围（1,73）中的n）]#_Michael S.Branicky_，2021年11月27日

%Y参考A000005、A000040、A000203、A007691、A008578、A300657。

%K nonn公司

%O 1,4型

%2018年9月7日，雅罗斯拉夫·克里泽克