A007953号:=proc(n)add(d,d=convert(n,base,10));结束进程:
A007954号:=程序(n)mul(d,d=转换(n,基数,10));结束进程:
A061762号:=进程(n)A007953号(n)+A007954号(n) ;结束进程:
具有是A000079 := 程序(n个) 如果 n个 < 1 然后 假; 否则如果 n个 = 1 然后 真实的; 其他的 如果 类型(n个,'即使') 然后 是( 无操作系统(数字理论):对于 n个 从 1 到 2000 做:我:=长度[因子集](n):编号0:=n个:秒:=0:第页:) = 1) ; 其他的 假; 结束 如果; 结束 如果; 结束 程序:对于 米
从1到l do:q:=n0:u:=irem(q,10):v:=iquo(q,100):n0:=v:s:=s+u:p:=p*u:od:
x: =s+p:z:=0:对于从0到50的a,而(z=0)do:b:=2^a:如果x=b,
`,n):z:=1:否则fi:od:od:
isA185400:=进程(n)为A000079(A061762美元(n) );结束进程:
对于从1到1300的n,如果是A185400(n),则执行printf(“%a,”,n);结束条件:;结束日期:#R.J.Mathar,2011年2月8日