%I#33 2022年9月8日08:46:07
%第10,12,24,32,36,48,54,60,72,76,84,96,98109112123125127129131页,
%电话:132133135137139141145147149151153155157159161163167,
%电话:169171172173175177179181183185191193195197199201209211213215217219224233235239241245
%N数字N,这样正好有两个数字m,其中m+(m的某个数字)=N。
%C数字12、112、1112。。。,111...112, ... 是序列的项_Marius A.Burtea,2020年2月18日
%D Eric Angelini,发布至序列粉丝邮件列表,2014年4月20日。
%H David A.Corneth,n的表,A(n)表示n=1..13111(术语<=10^6)
%e 12=6+6=11+1。
%e 32=26+6=31+1。
%e 112=106+6=111+1。
%p M:=2000;
%p M2:=M+10;
%p A:=阵列[0..M2];
%对于从0到M2的n,p做A[n]:=0;日期:
%p代表从0到M do的n
%p t1:=换算(n,基数,10);
%p t2:=转换(t1,设置);t3:=转换(t2,列表);
%对于i从1到nops(t3)的p,执行A[n+t3[i]]:=A[n+t3[i]]+1;日期:
%日期:
%p ans:=[];
%p表示n从0到M do,如果A[n]=2,则ans:=[op(ans),n];fi;日期:
%p[seq(ans[i],i=1..nops(ans))];
%t A241177[n_]:=模块[{m,c=0},
%t做[c=c+Count[m+Union[IntegerDigits[m]],n],{m,0,n}];c] ;
%t选择[范围[0,245],A241177[#]==2&](*_Robert Price_,2019年3月20日*)
%o(PARI)小于等于(n)={my(v=矢量(n+9));对于(i=1,n,d=集合(数字(i));对(j=1,#d,v[i+d[j]]++));针对(i=n+1,n+9,v[i]=0);选择(x->x==2,v,1)}\\_David A.Corneth_,2019年3月20日
%o(Magma)f:=func<n,m|存在(c){a:a在Intseq(m)|m+a方程n}>中;[1..250]|#[1..k]|f(k,m)]eq 2]中的[k:k_Marius A.Burtea,2020年2月18日
%Y相关序列:A241173、A241174、A241 175、A24 1176、A241177、A241178、A241179、A24.1180、A241181、A2410182、A2141183。
%K nonn,基础
%O 1,1号机组
%A _N.J.A.Sloane,2014年4月23日
|