%I#16 2020年3月3日08:58:33

%序号1,2,12,34,42,9017019840243445648249420464086451875207605，

%电话：8622963210924282805157051714741247772100172139653143136，

%电话：1767602945883998047182667535669053413583601577696208907421219402136256

%N将k编号为k|sigma_13（k）-phi（k）^13。

%C sigma_13（k）是k（A013961）除数的13次幂之和。

%H Amiram Eldar，<a href=“/A055707/b055707.txt”>n表，n=1..100的a（n）</a>

%t Do[If[Mod[DivisorSigma[13，n]-EulerPhi[n]^13，n]==0，打印[n]]，｛n，1，10^5｝]

%o（PARI）isok（n）=！（（sigma（n，13）-eulerphi（n）^13）%n）；\\_Michel Marcus，2014年3月2日

%Y参考A000010、A013961。

%K nonn公司

%O 1,2号机组

%A _Robert G.Wilson v_，2000年6月9日

%E 2014年3月2日修正的定义和来自Michel Marcus的更多术语