%I#36 2023年8月26日03:51:56

%S 1,3,4,7,6,12,8,15,10,18,12,28,14,24,24,31,18,30,20,42,32,36,24,60,26，

%电话：42,28,56,30,72,32,63,48,54,48,70,38,60,56,90,42,96,44,84,60,72,48，

%U 124,50,78,72,98,54,84,72120,80,90,60168,62,96,80127,84144,68126,96

%N普通单位西格玛（N）：如果N=乘积p_i^r_i，则OUSigma（N）=西格玛。

%C n的除数d之和，使得gcd（d，n/d）是2的幂（A000079）_Amiram Eldar，2023年8月26日

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

%F a（n）=A000203（p2）*A034448（n/p2），其中p2=A006519（n）_R.J.Mathar_，2008年6月15日

%F与a（2^e）=2^（e+1）-1相乘，对于p>2，e>0，a（p^e）=p^e+1。

%F和{k=1..n}a（k）~c*n^2，其中c=（4/7）*zeta（2）/zeta（3）=（4/7）*A306633=0.781961….-_Amiram Eldar，2022年11月1日

%F Dirichlet g.F.：（4^s/（4^s-2））*泽塔*泽塔（s-1）/泽塔（2*s-1）_Amiram Eldar，2023年8月26日

%e OUSigma（2^4*7^2）=西格玛（2^4）*单位西格玛。

%p A107749:=进程（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）-1）；否则a:=a*（p^e+1）；结束条件：；结束do；a；结束程序：#R.J.Mathar，2011年6月2日

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

%o（PARI）a（n）=本地（fm）；fm=系数（n）；prod（k=1，matsize（fm）[1]，if（fm[k，1]==2,2^（fm[k，2]+1）-1，fm[k，1]^fm[k,2]+1））

%Y参见A000079、A069184、A091321、A000203、A034448、A006519、A306633。

%K nonn，简单，多

%O 1,2号机组

%A _Yasutoshi Kohmoto，2005年6月11日，2007年2月24日

%E根据Andrew S.Plewe的建议，由N.J.A.Sloane编辑，2007年6月8日

%E更多条款来自R.J.Mathar_，2008年6月15日

%E姓名由Franklin T.Adams-Waters更正，2013年8月24日