%I#8 2020年5月23日05:44:07
%S 24,30,40,42,54,56,66,70,78,88,96102104114120138150174186222,
%电话:246258270282294318354360366402420426438474898534,
%电话:5405826066186306426546660678726262780786822834894906942
%N无穷好数:数k是k的一个适当的无穷除数d,使得isigma(k)-2*d=2*k,其中isigma是无穷除数函数的和(A049417)。
%等价地,等于其固有无穷除数之和的数,其中一个除数用减号表示。
%C好数(A111592)的除数是2的幂(A036537)也是无穷好数,因为它们的所有除数都是无穷大的。除数不是2的幂的项是96、150、294、360、420、486、540、630、660、726、780、960、990。。。
%H Amiram Eldar,n的表,a(n)表示n=1..10000</a>
%e 150是序列中的,因为150=1+2+3-6+25+50+75是其固有的无限除数之和,其中一个除数为6,用减号表示。
%t fun[p_,e_]:=模块[{b=整数位数[e,2],m},m=长度[b];乘积[如果[b[[j]]>0,1+p^(2^(m-j)),1],{j,1,m}]];isigma[1]=1;isigma[n_]:=倍@@fun@@FactorInteger[n];infDivQ[n_,1]=真;infDivQ[n_,d_]:=位与[IntegerExponent[n,First/@(f=FactorInteger[d])],(e=Last/@f)]==e;infAdmQ[n_]:=(ab=isigma[n]-2n)>0&&EvenQ[ab]&&ab/2<n&&Divisible[n,ab/2]&&infDivQ[n,ab/2];选择[Range[1000],infAdmQ]
%Y无限版本的A111592。
%Y A129656的后续序列。
%Y参考A036537、A049417、A077609、A328328、A334972。
%K nonn公司
%O 1,1号机组
%A阿米拉姆·埃尔达尔,2020年5月18日
|