%I#23 2019年9月26日09:14:22
%S 0,10,11,21,2,12,31,3,13,41,4,14,51,5,15,61,6,16,71,7,17,81,8,18,91,
%电话:9,1910111111131141151161171191201,22,32,23,42,24,52,
%U 25,62,26,72,27,82,28,92,29102211122122131232241142
%N a(1)=0。对于n>1,a(n)是序列中尚未出现的最小正整数k,因此k的最低有效位等于a(n-1)的最高有效位。
%C除了0和10之外,序列中没有10的倍数(A008592项),即序列中除0或10之外的任何项都是A067251项。
%除了10的倍数之外,序列中是否还有其他数字没有出现?换句话说,这是A067251 UNION[0,10]的置换吗?
%F A010879(a(n))=A000030(a(n-1))。
%t L={0};Do[k=整数位数[Last@L][[1];而[MemberQ[L,k],k+=10];附加到[L,k],{80}];L(*_Giovanni Resta_,2019年9月24日*)
%o(PARI)isinv(vec,k)=对于(t=1,#vec,if(vec[t]==k,return(1)));0
%o有效(x,y)=我的(d=数字(x),e=数字(y));d[#d]==e[1]
%o术语(n)=我的(v=[0,10]);while(1,如果(#v>=n,返回(v));for(k=1,oo,if(isvalid(k,v[#v]))&&!isinv(v,k),v=连接(v,[k]);断裂))
%o术语(100)\\打印最初的100个术语
%o(Python)
%o n,a,msdc=0,0,[1,0,0,0-0,0,1,0,0,0,0]
%o,当n≤62时:
%o打印(n,a)
%o s=str(a)
%o msd=int(s[0])
%o n,a=n+1,msdc[msd]*10+msd
%o msdc[msd]=msdc[msd]+1#_A.H.M.Smeets_,2019年9月25日
%Y参考A000030、A008592、A010879、A067251。
%K nonn,基础,简单
%O 1,2号机组
%A_Felix Fröhlich,2019年9月24日
|