登录
Tribonacci阵列,T(n,k)。
21

%I#27 2022年6月7日17:14:03

%S 1,2,3,4,6,5,7,11,9,8,13,20,17,15,10,24,37,31,28,19,12,44,68,57,51,35,

%电话:22,14,81125105,94,64,41,26,16149230193173118,75,48,30,18274,

%电话:423355318217138,88,55,33,21504778653585399254162101,61,39,23

%N Tribonacci数组,T(N,k)。

%作为一个散布(和弥散),数组作为一个序列,是正整数的置换。k列由数字m组成,因此m的摩擦学表示中的最小和为T(1,k)。例如,第1列由最小和为1的数字组成。该阵列源自tribonacci表示,与Wythoff阵列A035513源自斐波那契(或Zeckendorf)表示的方式大致相同。

%C From _Abel Amene_,2012年7月29日:(开始)

%C(第1行)=A000073(偏移量=4)a(0)=0,a(1)=0、a(2)=1

%C(第2行)=A001590(偏移量=5)a(0)=0,a(1)=1,a(2)=0

%C(第3行)=A000213(偏移量=4)a(0)=1,a(1)=1、a(2)=1

%C(第4行)=A214899(偏移量=5)a(0)=2,a(1)=1,a(2)=2

%C(第5行)=A020992(偏移量=6)a(0)=0,a(1)=2,a(2)=1

%C(第6行)=A100683(偏移量=6)a(0)=-1,a(1)=2,a(2)=2

%C(第7行)=A135491(偏移量=4)a(0)=2,a(1)=4,a(2)=8

%C(第8行)=A214727(偏移量=6)a(0)=1,a(1)=1、a(2)=2

%C(第9行)=A081172(偏移量=8)a(0)=1,a(1)=1、a(2)=0

%C(第1列)=A003265

%C(第2列)=A353083

%C(结束)[由John Keith更正和扩展,2022年5月9日]

%当k>3时,F T(1,1)=1,T(1,2)=2,T。第1行是tribonacci基础;写入B(k)=T(1,k)。每行满足递归T(n,k)=T(n、k-1)+T(n和k-2)+T。T(n,1)是不在前一行中的最小数。如果T(n,1)具有摩擦学表示B(k(1))+B(k(2))+。..+B(k(m)),则T(n,2)=B(k)(2))+B(k(3))+。..+B(k(m+1))和T(n,3)=B(k)(3))+B(k(4))+。..+B(k(m+2))。(指数的持续变化也给出了第n行中的其他术语。)

%e西北角:

%电子1 2 4 7 13 24 44 81 149 274 504

%电子邮箱3 6 11 20 37 68 125 230 423 778

%e 5 9 17 31 57 105 193 355 653

%电子邮箱:8 15 28 51 94 173 318 585

%电子邮箱:10 19 35 64 118 217 399

%电子12 22 41 75 138 254

%e 14 26 48 88 162

%电子16 30 55 101

%e 18 33 61

%电子21 39

%e 23(电子23)

%p#A73中的最大指数,使得A73<=n。

%p A73楼层Idx:=进程(n)

%p局部k;

%从3do到k的p

%p如果A000073(k)=n,则

%p返回k;

%p elif A000073(k)>n则

%p返回k-1;

%p end if;

%p端do:

%p端程序:

%n的p#tribonacci膨胀系数

%p A278038:=程序(n)

%p局部k,L,nres;

%p k:=A73楼层Idx(n);

%p L:=[1];

%编号:=n-A000073(k);

%p,而k>=4 do

%p k:=k-1;

%p如果nres>=A000073(k),则

%p L:=[1,op(L)];

%编号:=编号-A000073(k);

%p其他

%p L:=[0,op(L)];

%p end if;

%p端do:

%p返回L;

%p端程序:

%p A278038inv:=程序(L)

%p加(A000073(i+2)*op(i,L),i=1..nops(L));

%p端程序:

%p A135175:=进程(n,k)

%p选项记住;

%p本地a、已知、前v、nprev、kprev、freb;

%如果n=1,则为p

%p A000073(k+2);

%p elif k>3那么

%p进程名(n,k-1)+进程名(n,k-2)+进程名称(n,k-3);

%p其他

%p如果k=1,则

%p表示从1 do开始的a

%p已知:=假;

%nprev从1到n-1的p do

%p表示kprev从1 do开始

%p如果procname(nprev,kprev)>a那么

%p断裂;

%p elif procname(nprev,kprev)=那么

%p已知:=真;

%p end if;

%p端do:

%p端do:

%p如果未知,则

%p返回a;

%p end if;

%p端do:

%p其他

%p前缀:=进程名(n,k-1);

%p freb:=A278038(上一版本);

%p返回A278038inv([0,op(freb)]);

%p end if;

%p end if;

%p端程序:

%p序列(序列(A135175(n,d-n),n=1..d-1),d=2..12);#R.J.Mathar,2022年6月7日

%Y参考A035513、A353083和A353084。

%K nonn,表

%O 1,2号机组

%A_Clark Kimberling_,2007年12月18日

%E T(3,4)修正及更多条款,作者:John Keith,2022年5月9日