%I#31 2023年9月3日10:43:27
%S 1,3,4,3,6,12,8,3,4,18,12,24,14,24,24,3,18,30,20,38,32,36,24,48,6,42,
%电话:4,52,30,72,32,3,48,54,48,42,38,60,56,78,42,96,44,80,69,72,48,96,8,68,
%U 72,94,54,84,72108,80,90,60164,62,96,95,3,84144,68122,96144,72,66,74114,99136,96168,80158,4126,84220108132120168,90225112164128144120192,98122147,88单位
%N a(N)是N的除数d之和,它要么是1,要么是Product_{i}(p_i^e_i)形式,其中至少有一个e_i=1。
%C a(n)=n的非幂除数d之和,其中幂数是A001694(m)中m>=1的数。
%H Antti Karttunen,n的表,n=1..16385的a(n)</a>
%H<a href=“/index/Su#sums_of_divisors”>为与除数和相关的序列索引条目</a>。
%F a(n)=A000203(n)-A183099(n)=A183098(n)+1。
%F a(1)=1,a(p)=p+1,a(p*q)=(p+1)*(q+1),a(p*q*…*z)=(p+1)**(z+1),a(p^k)=p+1,对于p,q=素数,k=自然数,p*q**z=k(k>2)不同素数p,q,…,的乘积。。。,z。
%e对于n=12,这样的除数集是{1,2,3,6,12};a(12)=1+2+3+6+12=24。
%t f1[p,e]:=(p^(e+1)-1)/(p-1);f2[p,e_]:=f1[p,e]-p;a[1]=1;a[n_]:=倍@@f1@@@(f=系数整数[n])-倍@@f2@@@f+1;阵列[a,100](*_Amiram Eldar_,2023年8月29日*)
%o(PARI)A183100(n)=(1+sumdiv(n,d,d*(!ispowerful(d)));\\_Antti Karttunen,2017年10月7日
%Y参考A000203、A001694、A183098和A183099。
%K nonn,简单
%O 1,2号机组
%A _Jaroslav Krizek,2010年12月25日
|