%I#21 2019年12月4日18:18:35
%编号:6,7061612402090893217980197802032026980295124318049742,
%电话:516885805879000100130116870128570175370176715201376208280,
%电话:221536275770280670282680302176373065427924435470764483616618772642124
%N编号N,使sigma(N)=6*phi(N)。
%如果p>2&2^p-1是素数(梅森素数),那么序列中是5*2^(p-2)*(2^p-1)。所以5*2^(A000043-2)*(2^A000043-1)是这个序列的子序列。
%H Amiram Eldar,n的表,n=1..100000的a(n)(使用Jud McCranie的数据计算,Donovan Johnson的术语1..1000)
%H Kevin A.Broughan和Daniel Delbourgo,<A href=“https://cs.uwaterloo.ca/journals/JIS/VOL16/Broughan/broughan26.html“>关于除数和与欧拉总函数I之比,整数序列杂志,第16卷(2013年),第13.8.8条。
%H Kevin A.Broughan和Qizhi Zhou,<A href=“https://cs.uwaterloo.ca/journals/JIS/VOL17/Broughan/bro32.html“>关于除数和与欧拉总函数之比II,整数序列杂志,第17卷(2014年),第14.9.2条。
%ep>2,q=2^p-1(q是素数);m=5*2^(p-2)*q所以sigma(m)=6*(2^,p-1)-1)*2^p=6*phi(m),因此m在序列中。
%eσ(79000)=187200=6*31200=6*phi(7900),因此79000在序列中,但79000不是5*2^(p-2)*(2^p-1)的形式。
%t做[If[DivisorSigma[1,m]==6*EulerPhi[m],Print[m]],{m,1000000}]
%o(PARI)is(n)=σ(n)==6*eulerphi(n
%o(PARI)v=列表();对于因子(n=6,10^6,如果(sigma(n)==6*eulerphi(n),listput(v,n[1]));Vec(v)\\_Charles R Greathouse IV_,2017年5月9日
%Y参见A000043、A062699、A068390、A104901。
%K nonn公司
%O 1,1号机组
%A_Farideh Firoozbakht,2005年4月1日
|