%I#8 2016年3月17日04:22:41

%S 3,4,5,71425816127854912643169995170357121868338919592029，

%电话：103039104647112093137317149851218269261883266923323233，

%电话：3362734498815058915243716101716175699078739996431119253113422717287190052045171

%N个数，其算术导数等于其数字的某些固定幂之和。

%H Paolo P.Lava，<a href=“/A269719/A269719_1.txt”>序列术语及其固定功率</a>

%e3^0=1和3'=1；

%e4^1=4和4'=4；

%e 1 ^3+4 ^3+2 ^3=73和143’=73。

%p with（numtheory）：p:=进程（q）局部a，b，c，d，j，k，n，ok；对于从3到q的n，做a:=[]；b： =n；ok：=0；

%pd:=n*加（op（2，p）/op（1，p），p=ifactors（n）[2]）；a： =[]；b： =n；确定：=0；

%p表示k从1到ilog10（n）+1 do，如果（b mod 10）>1，则ok：=1；fi；a： =[（b mod 10），op（a）]；b： =trunc（b/10）；od；b： =-1；c： =0；

%p如果ok=1，则当c<d时，执行b:=b+1；

%p如果b>0，则c:=加（a[k]^b，k=1..nops（a））；否则，对于k从1到nops（a），如果a[k]=0，那么c:=0；断裂；

%p其他c：=c+1；fi；od；fi；od；如果c＝d，则lprint（n，b）；fi；fi；od；结束：P（10^9）；

%t f[n_]：=如果[Abs@n<2，0，n总计[#2/#1&@@@FactorInteger@Abs@n]；Select[Range[3，10^5]，Function[k，IntegerQ@SelectFirst[Range[0，10]，Fuction[d，If[MemberQ[d，0]&&#==0，Total@Power[d/.0->Nothing，#]==f@k，Total@Power[d、#]==f@k]]@IntegerDigits@k&]]（*_Michael De Vlieger_，2016年3月4日，第10版，f（n）after _Michale Somos_ at A003415*）

%Y参考A003415。

%K nonn，简单

%O 1,1号机组

%A_Paolo P.Lava，2016年3月4日