%I#30 2024年1月26日13:31:16

%S 2,3,4,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83，

%电话89,97101103107109113127131137139149151157163167173，

%电话1791811911931971992112232272233239241251257269271277

%N个素数和4个素数。

%C sopf（n）是不同素数除以n（A008472）的总和。由于sopf（n）=n，如果n是素数，那么这个序列就是n，因此n^sopf。

%素因子之和为n.-Arkadiusz-Wesolowski的C数n，2012年1月17日

%C数n，使2n正好有四个除数_韦斯利·伊万·赫特，2013年7月1日

%C数n，使n^2不除n！.-_Charles R Greathouse IV_，2013年11月4日

%F a（n）=A046022（n+1）_Omar E.Pol_，2012年11月27日

%p with（numtheory）：digits:=200:nn:=200:对于从1到nn的a do:t1:=ifactors（a）[2]：t2:=sum（t1[i][1]，i=1..nops（t1））：如果a^t2=t2^a，则打印f（`%d，`，a）：else fi:od:

%t插入[Prime[Range[60]]，4,3]（*哈维·P·戴尔，2024年1月26日*）

%o（PARI）a（n）=如果（n>3，素数（n-1），n+1）\\查尔斯·格里特豪斯IV，2011年8月26日

%Y参考A008472、A046022。

%K nonn，简单

%O 1,1号机组

%2010年9月4日，拉格瑙市

%E切换注释和名称_Charles R Greathouse IV_，2013年11月4日