%I#15 2019年1月21日20:06:41

%S 1,2,4,8,10,20,22,40,80100101103107110111113117130131133，

%电话137170171173177200206220260301305310311317331，

%电话：350371400404440350353060262070170971071171371773177179080088809079701000110031010101216

%N数字的性质是数字和加上数字的乘积是2的幂。

%e 371位于序列中，因为（3+7+1）+（3*7*1）=11+21=32=2^5。

%e 116291在序列中，因为（1+1+6+2+9+1）+（1*1*6*2*9*1）=20+108=128=2^7。

%p A007953:=过程（n）加（d，d=转换（n，基数，10））；结束进程：

%p A007954：=程序（n）mul（d，d=转换（n，基数，10））；结束进程：

%p A061762：=程序（n）A007953（n）+A007954（n）；结束进程：

%p is A000079：=proc（n），如果n<1，则为false；elif n=1，则为真；否则，如果类型（n，‘even’）则为（nops（numtheory[factorset]（n））=1）；否则为假；结束条件：；结束条件：；结束进程：

%p为A185400：=进程（n）为A000079（A061762（n））；结束进程：

%p表示n从1到1300，do如果是A185400（n），则打印f（“%a，”，n）；结束条件：；结束日期：#R.J.Mathar，2011年2月8日

%t pwrs2Q[n_]：=模块[{idn=整数位数[n]，x，y}，x=总数[idn]+次数@@idn；y=圆形[Log[x]/Log[2]；2^y==x]

%t选择[范围[1100]，pwrs2Q]（*哈维·P·戴尔，2011年2月16日*）

%Y抄送A061762。

%K nonn，基础

%O 1,2号机组

%2011年2月3日，拉格瑙市