%I#15 2019年6月26日08:11:38
%S 6,30,42,7010516821042057071474484012541848209024362730,
%电话29703135372050166270687880924010098107881186812180,
%电话:12192125401356614630158341618820790213182402488225080259082604026796321304389048360
%非限制平衡数:形式为m*N的平衡数,其中m,N>1都是平衡的。
%如果sigma(n)是phi(n)的倍数,则正整数n是平衡数(A020492)。由于phi和sigma是乘法函数,如果m和n是平衡数并且gcd(m,n)=1,那么mn也是平衡数。这个序列只包含这些非原罪术语。
%H Amiram Eldar,n的表,n=1..5000的a(n)</a>
%e2和3是平衡数,gcd(2,3)=1,所以6是非本原平衡数。
%t balQ[n_]:=可除[DivisorSigma[1,n],EulerPhi[n]];nonprimQ[n_]:=balQ[n]&&模[{d=Divisors[n],ans=False},Do[If[GCD[d[[k]],n/d[k]]==1&&balQ[d[[k]]]&&barQ[n/d[[k]],ans=真;Break[]],{k,2,Floor[Length[d]/2]}];ans];选择[Range[50000],nonprimQ](*_Amiram Eldar_,2019年6月26日*)
%Y参考A020492、A291565。
%K nonn公司
%O 1,1号机组
%A _Jud McCranie_,2017年8月26日
|