%I#37 2024年2月19日01:55:39
%S 1,2,3,4,5,6,7,8,9,15,18,24,45,4813514442883784767568641575,
%电话:3936669984139968
%N对m进行编号,使m除以(m的数字乘积)*(m的位数之和)。
%C下一项(如果存在)大于4*10^7。-_米歇尔·腾·沃德_
%C序列是有限的,上面有10^84的界,因为如果10^k<=n<10^(k+1)(n的位数的乘积)*(n的位数的和)<=k*9^(k+2),对于k>=84,这小于10^k。-_Henry Bottomley,2000年5月18日
%不允许使用零位数的C数字_Harvey P.Dale_,2011年7月16日
%C对于2.5*10^9没有其他条款_Robert G.Wilson v_,2011年7月17日
%C序列已完成。-_Giovanni Resta_,2013年3月20日
%C如果只对非零数字进行数字乘积,那么1088也在序列中_Giovanni Resta_,2013年3月22日
%H Giovanni Resta,用于计算完整序列的方法</a>
%e 139968在序列中,因为它除以(1*3*9*9*6*8)*(1+3+9+9+6+8)。-_Giovanni Resta_,2013年3月20日
%t okQ[n_]:=模块[{idn=IntegerDigits[n]},!成员Q[idn,0]&&Divisible[(总计[idn]*次@@idn),n]](*_哈维·P·戴尔,2011年7月16日*)
%t(*完整序列*)dig[nD_]:=块[{ric,sol={},check},check[mu_,minN_]:=Block[{di=数字计数@minN,k=1,r},While[(r=mu/k)>=minN,If[IntegerQ[r]&DigitCount[r]==di,AppendTo[sol,r]];k++]];ric[n_,prod_,sum_,lastd_,cnt_]:=块[{t},如果[cnt==nD,检查[prod*sum,n],Do[t=nD-cnt-1;如果[n*10^(t+1)<=d*prod*9^t*(sum+d+9*t),ric[10*n+d,d*prod,d+sum,d,cnt+1],Break[]],{d,9,lastd,-1}]];里克[0,1,0,1,0];打印[“nDig=”,nD,“sol=”,sol=排序@sol]; 溶胶];压扁[dig/@Range[84]](*_Giovanni Resta_,2013年3月20日*)
%Y参见A038369、A049102、A049105、A049106。
%K nonn、base、fini、full
%O 1,2号机组
%A _利维尔·杰拉德_
|