%I#29 2019年11月23日04:08:16

%S 3,9,27,35,81175243245729875104512251715218743752256125，

%电话：65618575114951200519683198552187524871261252902930625，

%电话：42875550065574755843559049600256428584035876859275109375126445130625137885140335153125

%N编号N，使psi（N）=2*phi（N）。

%C无平方项为3、35、1045、24871、29029、50065、58435、64285。。。该序列的无平方项在A062699中。请注意，A062699还具有非平方项：2011009、3189625、3722875。。。

%C如果n在序列中，那么所有与n具有相同质数因子集的数字也都在序列中。——罗贝尔·伊斯雷尔，2017年9月15日

%C所有术语都是奇数。可被3整除的词是3的幂_罗伯特·伊斯雷尔，2017年9月18日

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

%e3^k是所有k>0的术语，因为psi（3^k）=4*3^（k-1）=2*phi（3^ k）。

%p pp:=n->mul（（p+1）/（p-1），p=numtheory:-系数集（n））：

%p选择（pp=2，[seq（i，i=1..200000,2）]）；#_罗伯特·伊斯雷尔，2017年9月15日

%t psi[n_]：=n*总和[MoebiusMu[d]^2/d，{d，除数@n}]; 选择[Range@200000，2EulerPhi[#]==psi[#]&]

%o（PARI）a001615（n）=我的（f=系数（n））；prod（i=1，#f~，f[i，1]^f[i、2]+f[i（1]^（f[i）-1））；

%o isok（n）=a001615（n）==2*eulerphi（n）；\\A001615查尔斯·格里特豪斯四世之后

%Y参见A000010、A001615、A062699、A291051、A291932。

%K nonn，简单

%O 1,1号机组

%A _Robert G.Wilson v _，_Amiram Eldar _和_Altug Alkan _，2017年9月15日