%I#33 2022年2月5日16:43:50
%S 1,1,3,5,15,17,63,852193251023110540955397135152184565535,
%电话:707372621433331258905231397077419430345271851623601522365525,
%电话:575218838842984526843545527296262510737418231431655756793023635726557525
%N a(N)=和{1<=k<=N,gcd(k,N)=1}2^(k-1)。
%C对于n>0,将n的约化剩余集(三角形A054431的行)解释为二进制数而形成的数。
%H Robert Israel,n表,n=1..2658的a(n)</a>
%F M*V,其中M=A054521是一个无限下三角矩阵,V=[1,2,4,8,…]是一个向量_Gary W.Adamson_,2007年1月13日
%F a(4*n)=(2^(2*n)+1)*a(2*n)[想想4n型数的约化残数集是如何形成的]。
%对于所有素数p和整数e>1,A054432(p^e)=A019320(p^e)*((2^(p^(e-1)))-1)*(2^(p^1)-1))/(2^p)-1)。
%F a(n-1)=和{k=1..n,gcd(n,k)=1}2^(k-1).-_Vladeta Jovovic_,2002年8月15日
%e对于n=6,我们有k=1和5,然后2^0+2^4=17=a(6)。
%p rrs2bincode:=进程(n)局部i,z;z:=0;对于i从1到n-1,做z:=z*2;如果(1=igcd(n,i)),则z:=z+1;fi;od;返回(z);结束;
%t f[n_]:=总和[2^k,{k,选择[Range@n,GCD[#,n]==1&]-1}];数组[f,35](*_Robert G.Wilson v_,2014年7月21日*)
%o(PARI)a(n)=总和(k=1,n,如果(gcd(k,n)==1,2^(k-1),0));\\_米歇尔·马库斯,2014年7月20日
%o(PARI)a(n)=subst(Polrev(向量(n,i,gcd(n,i)==1)),x,2);\\_米歇尔·马库斯,2014年7月21日
%Y参见A054431、A054433、A001317、A054521。
%K nonn公司
%氧1,3
%安蒂·卡图恩_
%E根据R.J.Mathar的建议,由N.J.A.Sloane编辑,2008年7月3日_
%E来自_Michel Marcus_的更多术语,2014年7月20日
|