%I#17 2016年6月7日10:36:19

%S 0,1,2,3,4,6,7,8,10,12,15,16,18,20,22,24,26,28,30,31,32,36,40,42,44，

%电话：46,48,52,54,56,58,60,63,64,66,70,72,78,80,82,84,88,92,96100102104，

%电话：106108110112116120124126127130132136138140144148150156160162164168172176

%N（素数幂减1）与欧拉-托利函数值的并集。

%F A181062和A002202的接头。

%t最大值=200；

%t selNu=选择[Range[max]，PrimeNu[#]==1&]-1；

%t phiQ[m_]：=选择[Range[m+1，2*m*乘积[1/（1-1/（k*Log[k]）），{k，2，DivisorSigma[0，m]}]]，EulerPhi[#]==m&，1]！={};

%t selPhi=选择[Range[max]，phiQ]；

%t联接[{0}，联合[selNu，selPhi]]

%o（PARI）列表（lim）=我的（P=1，q，v，u=列表（[0]））；对于素数（p=2，默认（素数极限），如果（eulerphi（p*=p）>=lim，q=p；断裂））；v=vecsort（向量（P/q*lim\eulerphi（P/q），k，eulerpchi（k）），8）；v=选择（n->n<=lim，v）；对于素数（p=2，平方（lim\1+1），p=p；而（P*=P）<=lim+1，listput（u，P-1））；向量排序（concat（v，Vec（u）），，8）\\_Charles R Greathouse IV_，2013年1月8日

%Y参见A000010、A002202、A000961、A181062、A070932（乘法闭包）。

%K nonn公司

%氧1,3

%A _Jean-François Alcover，2013年1月6日

%E编辑：N.J.A.Sloane，2013年1月6日