%I#33 2022年9月8日08:45:07

%S 2,3,5,7,9,21,27,65,69,70,81,84,90110123123130153154189222，

%电话：228243252259264327329333340343350365372381402434，

%电话：450516528580588618621650684729731731738740741756765774

%N数N，使得以9为底的数字之和是素数除数之和的除数（A008472）。

%序列是无限的，因为对于m=9^k，k>=0，digsum（m_9）=1_Marius A.Burtea_，2019年7月10日

%H Marius A.Burtea，n的表，n的A（n）=1..5189</a>

%e 21=23_9，digsum（23_9）=5，素数除数（21）={3，7}，sopf（21）=3+7=10=5*2。

%p A076387：=proc（n）局部i，j，t，t1，sod，sopd；t:=空；对于i从2到n，做t1:=i；超氧化物歧化酶：=0；t1<>0时，sod:=sod+（t1mod9）；t1:=地板（t1/9）；od；sopd:=0；j:=1；当ithprime（j）<=i do时，如果i mod ithprime；fi；j:=j+1；od；如果sopd-mod-sod=0，则t:=t，i；fi；od；t；结束；

%o（PARI）{对于（ixp=2783，

%o casi=ixp；cvst=0；dsu=0；M=系数（ixp）；smt=0；

%o表示（i=1，矩阵大小（M）[1]，smt=smt+M[i，1]）；

%o while（casi！=0，

%o cvd=casi%9；dsu=dsu+cvd；casi=（casi-cvd）/9）；

%o if（smt%dsu==0，print1（ixp，“，”））｝\\_Douglas Latimer_，2012年5月8日

%o（Magma）[1..800]中的n：n |&+PrimeDivisors（n）mod&+Intseq（n，9）eq 0]；//_Marius A.Burtea，2019年7月10日

%Y参见A075657、A076380-A076387。

%K nonn，基础

%O 1,1号机组

%2002年10月8日，A层面包车拉蒙