N: =1000:#将a(1)转换为a(N)
K: =天花板(N-4)/3):
Dmax:=ilog[3](天花板(K/2+1)):
A: =矢量(3*K+4):
A[1..4]:=<1,2,4,7>:
对于从0到Dmax的d do
对于从2*3^d-1到最小(4*3^d-2,k)的k,do
A[3*k+2]:=9*A[k]+3;
A[3*k+3]:=9*A[k]+4;
A[3*k+4]:=9*A[k]+5
日期:
对于k从4*3^d-1到min(2*3^(d+1)-2,k)do
A[3*k+2]:=9*A[k]+1;
A[3*k+3]:=9*A[k]+4;
A[3*k+4]:=9*A[k]+7
日期:
日期:
seq(A[i],i=1..N)#罗伯特·伊斯雷尔,2016年6月6日