%I#17 2021年2月15日22:31:10

%S 1,2,4,7,12,13,14,21,22,23,37,40,43,64,67,7011112113120121122，

%电话：129130131192193194201202203210212334337340361364，

%电话：367388391394577580583604607610631

%以3为底表示和{i=0..m}d（i）*3^i的N个数，对于m-i奇数，d（i）=1。

%H Robert Israel，n的表，n的a（n）=1..10000</a>

%F From _Robert Israel_，2016年6月6日：（开始）

%Fα（3n+3）=9a（n）+4。

%F如果A110654（n）在A132141中，则a（3n+2）=9a（n）+3和a

%F否则a（3n+2）=9a（n）+1，a（3n+4）=9 a（n）+7。

%F G.F.满足G（x）=9（x^2+x^3+x^4）G（x^3）+（x+2x^2+4x^3+6x^4-x^5）/（1-x^3。

%F（结束）

%p N:=1000:#得到a（1）到a（N）

%p K:=细胞（N-4）/3）：

%p Dmax:=ilog[3]（天花板（K/2+1））：

%p A:=矢量（3*K+4）：

%p A[1..4]：=<1,2,4,7>：

%p代表d从0到Dmax do

%p代表k从2*3^d-1到min（4*3^d-2，k）do

%pA[3*k+2]：=9*A[k]+3；

%pA[3*k+3]：=9*A[k]+4；

%p A[3*k+4]：=9*A[k]+5

%日期：

%p代表k从4*3^d-1到min（2*3^（d+1）-2，k）do

%pA[3*k+2]：=9*A[k]+1；

%pA[3*k+3]：=9*A[k]+4；

%pA[3*k+4]：=9*A[k]+7

%日期：

%日期：

%p序列（A[i]，i=1..N）；#_罗伯特·伊斯雷尔（Robert Israel），2016年6月6日

%K nonn，基础

%O 1,2号机组

%百灵鸟金伯利_

%E姓名由Robert Israel_更正，2016年6月6日