登录
n个X n个二元矩阵的唯一X射线数,精确到n个。

%I#25 2018年5月6日11:59:43

%S 1,1,2,5,8,13,21,31,45,65,9212717523731856173591241,

%电话1597204726073305417452476569819710189126211558819189,

%电话:2355128829351894284152033630397619791903110603132831159215190463227416

%N N个二元矩阵的唯一X射线数,正好有N个X射线。

%矩阵的X射线定义为反对角线和的序列。

%C独特的X射线可以重建二进制矩阵。

%H Alois P.Heinz,n的表格,n=0..10000的a(n)</a>

%H C.Bebeacua、T.Mansour、A.Postnikov和S.Severini,<A href=“https://arxiv.org/abs/math/0506334“>关于排列的X射线,arXiv:math/0506334[math.CO],2005。

%H<a href=“/index/Mat#binmat”>与二进制矩阵相关的序列的索引项</a>

%F a(n)~exp(Pi*sqrt(2*n/3))/(2^(9/4)*3^(1/4)*n^(3/4))_瓦茨拉夫·科特索维奇,2018年5月6日

%e a(3)=5:00021、00300、02001、10020、12000。

%e a(4)=8:0000301、0004000、0030001、0200020、1000021、1000300、1030000、1200001。

%p b:=proc(n,i)选项记忆;(m->`if`(n>m,0,

%p`如果`(n=m或n=0,1,加上(b(n-i*j,min(n-i*j,i-1))*

%p`如果`(j=1,2,1),j=0..分钟(2,n/i)))(i*(i+1))

%p端:

%p a:=n->`如果`(n=0,1,1+b(n,n-1)):

%p序列(a(n),n=0..60);

%tb[n_,i_]:=b[n,i]=函数[m,如果[n>m,0,如果[n==m|n==0,1,和[b[n-i*j,Min[n-i**j,i-1]]*如果[j==1,2,1],{j,0,Min[2,n/i]}]]][i*(i+1)];

%t a[n_]:=如果[n==0,1,1+b[n,n-1]];

%t表[a[n],{n,0,60}](*_Jean-François Alcover_,2017年11月7日,在_Alois P.Heinz_*之后)

%Y参考A290052、A290134。

%K nonn公司

%0、3

%A _Alois P.Heinz,2017年7月20日