%I#19 2018年1月10日11:31:19
%S 0,0,0.0,0,52,0,11,0,0,10,51,50,0,49,10,00,48,0,9,50,049,0,048,0,0,
%电话:8,0,49,46,0,47,48,0,45,0,0,7,46,7,47,6,0,45,44,0,46,5,0,44,44,
%U 0,43,4,5,4,0,0,4,44,4,4,43,3,0,04,0,42,0,3,4,43
%N通过b_s(1)=s,b_s;a(n)是bn到达b0、b1、b3或b9中的一个项所需的步骤数,或者,如果bn从未加入这四个序列中的任何一个,则a(n)=-1。
%我们推测a(n)决不是-1。
%H Lars Blomberg,n的表,n=0..10000的a(n)</a>
%H<a href=“/index/Coi#Colombian”>哥伦比亚或自身编号和相关序列的索引条目</a>
%p读取转换;#得到数字和
%p M:=2000;
%p#f(s)返回从s开始的序列k->k+数字和(k)
%p f:=进程全局M;选项记忆;局部n,k,s1;
%p s1:=[s];k: =秒;
%n从1到M的p do k:=k+数字和(k);
%p s1:=[op(s1),k];od:结束;
%p#g(s)返回(x,p),其中x=之间的第一个公共数字
%p#f(s)和f(1)、f(3)、f。
%p#如果f(s)和f(1)、f(3)、f
%p S1:=转换(f(1),设置):
%p S3:=转换(f(3),集):
%p S9:=转换(f(9),集合):
%p g:=进程全局f,S1,S3,S9;局部t1、p、T0、t1、T2;
%p T0:=f(s):
%p T1:=转换(T0,设置);
%p如果((s mod 9)=3)或(s mod9)=6),则T2:=T1与S3相交;t1:=最小值(T2);如果(t1=无穷大),则返回(-1,-1);else成员(t1,T0,'p');返回(t1,p-1);fi;
%p elif((s mod 9)=0)则T2:=T1与S9相交;t1:=最小值(T2);如果(t1=无穷大),则返回(-1,-1);else成员(t1,T0,'p');返回(t1,p-1);fi;
%p其他T2:=T1与S1相交;t1:=最小值(T2);如果(t1=无穷大),则返回(-1,-1);else成员(t1,T0,'p');返回(t1,p-1);fi;
%p fi;
%p端;
%p[序列(g(n)[2],n=1..45)];
%Y参见A230107、A062028、A004207、A016052、A007618、A006507、A0160096。
%K nonn,基础,看
%0、6
%A _N.J.A.Sloane和_Reinhard Zumkeller,2013年10月21日
%E 2018年1月10日_Lars Blomberg_的第a(46)条及以后条款
|