%I#13 2023年7月25日20:24:15
%S 1,3,7,4,15,13,12,6,31,40,39,31,28,24,18,8,6312120156,91124,93,
%电话:57,60,78,72,48,42,32,24,12127364363781280624468400195403,
%电话:372342217228171133124240234248168192144,96,90104,96,72,56,48,36,1425510923906847312432801600
%N a(N)=西格玛(A163511(N))。
%H Antti Karttunen,n的表格,n=0..16383的a(n)</a>
%H Antti Karttunen,<a href=“/A324184/A324184.txt”>数据补充:n,a(n)为n=0..65537计算</a>
%F a(n)=A000203(A163511(n))。
%F对于n>=1,a(n)=A324054(A054429(n))。
%o(PARI)A324184(n)=如果(!n,1,my(p=2,mp=p*p,m=1);while(n>1,if(n%2,p=nextprime(1+p));mp=p*p,如果(2==n)||!(n%4),mp*=p,m*=(mp-1)/(p-1));n>>=1);(m*(mp-1)/(p-1)));
%o(PARI)
%o A163511(n)=如果(!n,1,my(p=2,t=1);while(n>1,if(!(n%2),(t*=p),p=下一素数(p+1));n>>=1);(t*p));
%o A324184(n)=西格玛(A163511(n));
%o(Python)
%o从sympy导入nextprime
%o定义A324184(n):
%o如果n:
%o c,p=1,1
%o而n:
%o c*=((p:=nextprime(p))**(s:=(~n&n-1).bit_length()+1)-1)//(p-1)
%o n>>=秒
%o返回c*(p**(s+1)-1)//(p**s-1)
%o 2023年7月25日返回1#Chai Wah Wu_
%Y参见A000203、A054429、A163511、A324054、A324183、A324185、A326186、A324.187、A32.4188、A324 189。
%K nonn公司
%0、2
%A _Antti Karttunen,2019年2月17日