%I#14 2022年6月16日21:18:56
%S 1,1,1,1,1,1,1,1,1,1,2,1,1,8,1,1,11,1,3,2,1,1,7,1,1,18,1,8,1,1,20,
%T 1,12,3,1,1,2,28,1,2,1,30,32,1,1,2,1,1,1,1,32,3,1,1,4,6,2,1,8,1,16,1,
%U 48,20,1,3,44,12,7,1,1,2,54,20,2,1,4,1,1,1,9,32,1,56,35,1,32,72,66,2,1,72,31,1,20单位
%N a(N)=A047994(N)/gcd(A04799四(N),A344005(N))。
%C分数A047994(n)/A344005(n)的分子。
%H Antti Karttunen,n的表,n=1..20000的a(n)</a>
%H Antti Karttunen,<a href=“/A354986/A354986.txt”>数据补充:n,a(n)计算n=1.65537</a>
%F a(n)=A047994(n)/A354985(n)=A047995(n)/gcd(A047994n),A344005(n))。
%ts[n_]:=模[{m=1},而[!可除[m*(m+1),n],m++];m] ;f[p_,e_]:=p^e-1;uphi[n_]:=倍@@f@@FactorInteger[n];uphi[1]=1;a[n_]:=分子[uphi[n]/s[n]];阵列[a,100](*_Amiram Eldar_,2022年6月16日*)
%o(PARI)
%o A047994(n)={my(f=因子(n)~);prod(i=1,#f,(f[1,i]^f[2,i])-1);};
%o A344005(n)=对于(m=1,oo,如果(m*(m+1))%n==0,返回(m));\\来自A344005
%o A354986(n)={my(u=A047994(n));(u/gcd(u,A344005(n);
%Y参见A047994、A344005、A346607、A346408(术语位置>1)、A354928(1的位置)、A354/985、A354987(分母)。
%K nonn,压裂
%O 1,12号
%A _Antti Karttunen,2022年6月16日