%I#13 2020年1月19日12:48:11

%S 4,8,9,10,16,18,20,21,22,24,25,26,27,28,32,33,34,36,38,39,40,44,45,46，

%电话：48,49,50,51,52,54,55,57,58,60,62,63,64,65,66,68,69,72,74,75,76,77,80，

%U 81,82,84,85,86,87,88,90,91,92,93,94,95,96,98,99100102104106108单位

%N不平衡复合数。

%e 46=2*23和sigma（46）/phi（46）=72/22不是整数，所以46在序列中。

%t Do[s=DivisorSigma[1，n]/EulerPhi[n]；如果[！IntegerQ[s]&&！PrimeQ[n]，打印[n]]，{n，1，256}]

%t选择[Range[150]，CompositeQ[#]&&！IntegerQ[DivisorSigma[1，#]/EulerPhi[#]]&]（*需要Mathematica版本10或更高版本*）（*Harvey P.Dale_，2020年1月19日*）

%o（PARI）lista（nn）=复合（n=1，nn，if（分母（sigma（n）/eulerphi（n））！=1，打印1（n，“，”）；\\_米歇尔·马库斯，2018年7月11日

%Y参见A020492、A078538、A078529、A078540。

%K nonn公司

%O 1,1号机组

%A _Labos Elemer，2002年12月4日