%I#40 2021年9月5日18:20:49
%S 1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,3,3,1,1,11,1,5,6,5,11,1,1,8,13,13,
%电话:8,1,1,1,13,28,40,28,13,1,1,1,21,60117117,60,21,1,1,11,34129,
%电话:348472348129,34,1,1,1,55277102916161029277,55,1,1
%N使用集成方形瓷砖的k X N矩形瓷砖的数量A(N,k);方阵A(n,k),n>=0,k>=0。
%C对于A(1,1)、A(2,2)、……的图纸。。。,A(5,5)见A224239。
%H Alois P.Heinz,反对角线n=0..30,扁平</a>
%H Steve Butler、Jason Ekstrand、Steven Osborne,<a href=“https://doi.org/10.1007/978-3-030-37853-0_5“>通过在图中行走计算瓷砖数,基于项目的数学本科生研究指南,Birkhäuser,Cham(2020),见第169页。
%e A(3,3)=6,因为一个3X3矩形有6个瓷砖,使用的是整数正方形:
%e._____.._____.._____.._____.._____._____.._____。
%e ||||_||_| |_|| | _|_|_| | |_||_|_|_|
%e |||___|_||_|___||_| | | | |_||_|_|_|
%e | _____ | | _ | _ | | _ |_ | | _| | _|_ | ___ | | ___ |_||_|_|_|
%e方阵A(n,k)开始:
%e 1,1,1,1,1,1,1,1。。。
%e 1,1,1,1,1,1,1,1。。。
%e 1、1、2、3、5、8、13、21。。。
%e 1、1、3、6、13、28、60、129。。。
%e 1、1、5、13、40、117、348、1029。。。
%e 1、1、8、28、117、472、1916、7765。。。
%e 1、1、13、60、348、1916、10668、59257。。。
%e 1、1、21、129、1029、7765、59257、450924。。。
%p b:=proc(n,l)选项记忆;局部i,k,s,t;
%p如果max(l[])>n,则0 elif n=0或l=[],则1
%p elif min(l[])>0,则t:=min(l[]);b(n-t,映射(h->h-t,l))
%如果l[k]=0,则p表示k do,然后打破fiod;s: =0;
%p表示i从k到nops(l),而l[i]=0表示s:=s+
%p b(n,[l[j]$j=1..k-1,1+i-k$j=k.i,l[j]$j=i+1..nops(l)])
%p od;秒
%功率因数
%p端:
%p A:=(n,k)->`如果`(n>=k,b(n,[0$k]),b(k,[0$n])):
%p序列(序列(A(n,d-n),n=0..d),d=0..14);
%p#下面是列出实际解剖的程序的第二个版本。它为每个解剖生成一个成对的列表:
%p b:=进程(n,l,ll)局部i,k,s,t;
%p如果max(l[])>n,则0 elif n=0或l=[],则lprint(ll);1
%p elif min(l[])>0,则t:=min(l[]);b(n-t,映射(h->h-t,l),ll)
%如果l[k]=0,则p表示k do,然后打破fiod;s: =0;
%p表示i从k到nops(l),而l[i]=0表示s:=s+
%p b(n,[l[j]$j=1..k-1,1+i-k$j=k.i,l[j]$j=i+1..nops(l)],
%p[ll[],[k,1+i-k]])
%p od;秒
%功率因数
%p端:
%p A:=(n,k)->b(k,[0$n],[]):
%p A(5,5);
%p#在每个列表中[a,b]表示放置边长为b的正方形
%p最左侧可能的位置,上角位于行a中。例如
%p[[1,3],[4,2],[4,1],[1,2]
%第_____.___页。
%第||页|
%第页||___|
%第页|_____|_|_|
%第||页|_|
%第页|___|___|_|
%tb[n_,l_List]:=b[n,l]=模[{i,k,s,t},其中[Max[l]>n,0,n==0|l=={},1,Min[l]>0,t=Min[l];b[n-t,l-t],True,k=位置[l,0,1][[1,1]];s=0;对于[i=k,i<=长度[l]&l[i]]==0,i++,s=s+b[n,连接[l[[1;;k-1]],表[1+i-k,{j,k,i}],l[[i+1;;-1]]]];s] ];a[n_,k_]:=如果[n>=k,b[n,数组[0&,k]],b[k,数组[0-,n]]];表[表[a[n,d-n],{n,0,d}],{d,0,14}]//Flatten(*_Jean-François Alcover_,2013年12月13日,翻译自Maple第一个程序*)
%Y列(或行)k=0+1,2-10给出:A000012,A000045(n+1),A002478,A054856,A05485%,A219925,A219926,A219927,A219928,A219929。
%Y主对角线给出A045846。
%Y参考A113881,A226545。
%K nonn,表
%O 0,13号
%A _Alois P.Heinz,2012年12月1日