%I#27 2021年10月28日20:15:03

%S 2,4,8,16,24,32，6412825632051210240483712409671688192，

%电话：105121619216384327683302437884172849280515525452857280，

%电话：6259265536664326796868832679281600843528769691968920169388947201212812988131072

%N满足tau（m）+omega（m）=#（{phi（x）=m}）的整数m。

%C A000079的所有偶数项都包含在此序列中。

%C a（5）=24是第一项而不是A000079的项，a（10）=320是第二项。

%H Robert Israel，n的表，n=1..1000的a（n）</a>

%H Max Alekseyev，<a href=“http://home.gwu.edu/~maxal/gpscripts/“>PARI脚本用于各种问题</a>（请参阅此处的invphi.gp）。

%Fτ（m）+Ω（m）=#（{φ（x）=m}）。

%F整数m，使A163523（m）=A014197（m）。

%e2是一个项，因为tau（2）=2，omega（2）=1，并且#（{phi（x）=2}）=3。

%e24是一个术语，因为tau（24）=8，omega（24）=2，并且#（{phi（x）=24}）=10。

%p过滤器：=proc（n）使用numtheory；tau（n）+nops（因子集（n））=nops（invphi（n）

%p选择（过滤器，[seq（i，i=2..10^5,2）]）；#_罗伯特·伊斯雷尔，2021年10月28日

%t块[{nn=10^5，s}，s=函数[s，函数[t，Take[#，nn]&@ReplacePart[t，Map[#->Length@Lookup[s，#]&，Keys]]@ConstantArray[0，Max@Keys]]@KeySort@PositionIndex@Array[EulerPhi，Floor[nn^（3/2）]+10]；选择[Range@nn，DivisorSigma[0，#]+PrimeNu[#]==s[[#]]&]]（*_Michael De Vlieger_，2018年7月21日*）

%o（PARI）isok（m）=numdiv（m）+omega（m）==#invphi（m）；\\_米歇尔·马库斯，2018年6月8日

%Y参考A000005、A001222、A163523。

%Y参考A000010、A058277、A014197。

%Y参考A000079。

%K nonn公司

%O 1，1

%A _Torlach Rush，2018年6月7日

%E来自米歇尔·马库斯的更多条款，2018年6月8日