%I#41 2024年4月25日09:12:06

%S 2,4,8,9,10,12,16,18,24,25,26,28,32,34,36,40,48,50,52,58,63,64,72,74，

%电话：75,76,80,81,82,84,88,90,98100104106108112117120121122124，

%电话：128130136144146152156160162170171172175

%对k进行编号，使k的除数数不除以k的除法和。

%D József Sándor、Dragoslav S.Mitrinovic和Borislav Crstic，《数论手册I》，Springer科学与商业媒体，2005年，第三章，第119页，第III.51节。

%H Reinhard Zumkeller，n的表，n=1..10000的a（n）</a>

%H维基百科，<a href=“http://en.wikipedia.org/wiki/Arithmetic_number“>算术数字</a>。

%F A054025（a（n））>0.-_Reinhard Zumkeller，2012年1月6日

%F A245656（a（n））=0.-_Reinhard Zumkeller，2014年7月28日

%p为A049642:=进程（n）

%p如果modp（numtheory[sigma]（n），numtheori[tau]（n

%p假；

%p其他

%p为真；

%p end if；

%p端程序：

%p A049642:=进程（n）

%p选项记忆；

%如果n=1，则为p

%第2页；

%p其他

%p表示来自procname（n-1）+1 do的a

%p如果是A049642（a），则

%p返回a；

%p end if；

%p端do：

%p end if；

%p结束程序：#_R.J.Mathar_，2015年10月26日

%t选择[Range[175]，Mod[DivisorSigma[1，#]，DivisorSigma[0，#]]>0&]（*Jayanta Basu_，2013年3月28日*）

%o（哈斯克尔）

%o a049642 n=a049642_列表！！（n-1）

%o a049642_list=过滤器（==0）。a245656）[1..]

%o——Reinhard Zumkeller，2012年1月6日

%o（GAP）a:=过滤（[1..180]，n->Sigma（n）mod Tau（n）>0）；；打印（a）；#_Muniru A Asiru_，2019年1月25日

%o（PARI）为（n）={my（f=系数（n））；σ（f）%numdiv（f）>0；}\\阿米拉姆·埃尔达尔，2024年4月25日

%A003601的Y互补。

%Y参考A000005、A000203。

%Y参考A054025、A245656。

%K nonn，简单

%O 1,1号机组

%A _N.J.A.斯隆_