%I#32 2020年9月15日09:52:31

%S 1,3,4,5,6,12,8,9,13,18,12,20,14,24,24,17,18,39,20,30,32,36,24,31，

%电话：42、40、40、30、72、32、33、48、54、48、65、38、60、56、54、42、96、44、60、78、72、48、68，

%U 57,93,72,70,54120,72,72,80,90,60120,62,96104,65,84144,68,90,96单位

%N N的除数d之和，使得d或N/d是奇数。

%C可能被称为UnitaryOrdinarySigma（n）：如果n=乘积p_i^r_i，那么UOSiga（n_Kohmoto Yasutoshi，2005年6月11日

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

%H<a href=“/index/Su#sums_of_divisors”>与除数和相关的序列的索引项</a>

%对于奇素数p，F与a（2^e）=2^e+1和a（p^e）=（p^（e+1）-1）/（p-1）相乘。

%计算公式：和{m>0}m*x^m*（1+x^m+x^（2*m）-x^。

%F Dirichlet g.F.：zeta（s）*zeta（s-1）*（2^（2-3s）-2-（1-2s）-2^（1-s）+1）/（1-2^_R.J.Mathar，2011年6月2日

%F和{k=1..n}a（k）~7*Pi^2*n^2/96.-_Vaclav Kotesovec_，2019年2月8日

%e UOSigma（2^4*7^2）=单位sigma（2^4）*sigma（7^2”）=17*57=969。

%p A069184:=进程（n）局部a，f，p，e；a:=1；对于ifactors（n）[2]中的f，做p:=op（1，f）；e：=op（2，f）；如果p=2，则a:=a*（2^e+1）；否则a:=a*（p^（e+1）-1）/（p-1）；结束条件：；结束do；a；结束程序：#R.J.Mathar，2011年6月2日

%t表[Sum[d*Boole[OddQ[d]||OddQ[n/d]]，{d，Divisors[n]}]，{n，1，69}]（*_Jean-François Alcover_，2013年3月26日*）

%tf[2，e_]：=2^e+1；f[p_，e_]：=（p^（e+1）-1）/（p-1）；a[1]=1；a[n_]：=倍@@f@@FactorInteger[n]；阵列[a，100]（*_Amiram Eldar_，2020年9月15日*）

%o（PARI）a（n）=汇总（n，d，d*（（d%2）||（n/d）%2））；\\_米歇尔·马库斯（Michel Marcus），2014年4月10日

%o（PARI）a（n）=我的（e=估价（n，2））；西格玛（n>>e）*if（e，2^e+1,1）\\_Charles R Greathouse IV_，2014年4月10日

%Y参见A069733、A107749、A092356。

%K mult，nonn（非n）

%O 1,2号机组

%A _弗拉德塔·乔沃维奇，2002年4月10日

%E根据R.J.Mathar的建议，由N.J.A.Sloane于2008年8月29日编辑_