OEIS哀悼西蒙斯感谢西蒙斯基金会支持包括OEIS在内的许多科学分支的研究。
登录
OEIS由OEIS基金会的许多慷慨捐赠者.

 

标志
提示
(来自的问候整数序列在线百科全书!)
A352687型 按行读取的三角形,一个与Narayana相关的三角形,其行是加泰罗尼亚数字的两倍的细化(对于n>=2)。
1, 0, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 4, 4, 1, 0, 1, 7, 12, 7, 1, 0, 1, 11, 30, 30, 11, 1, 0, 1, 16, 65, 100, 65, 16, 1, 0, 1, 22, 126, 280, 280, 126, 22, 1, 0, 1, 29, 224, 686, 980, 686, 224, 29, 1, 0, 1, 37, 372, 1512, 2940, 2940, 1512, 372, 37, 1 (列表;桌子;图表;参考;;历史;文本;内部格式)
抵消
0,9
评论
这是纳拉亚纳三角形序列中的第二个三角形。第一个是A090181号其第n行是加泰罗尼亚语(n)的细化,而这里T的第n行则是2*加泰罗尼亚语(n-1)的细化。我们可以证明T(n,k)<=A090181号(n,k)表示所有n,k。这个序列中的第三个三角形是A353279型,其中也给出了一般情况下的递归。
这里,我们给出了行多项式的一个递推公式,它对应于Sulanke组合证明的经典Narayana多项式的递推公式(参见链接)。
多项式只有实数零,形成Sturm序列。这是根据Chen等人的论文中给出的路线重复得出的。
一些有趣的序列是多项式序列在一个固定点上的求值(参见交叉引用),例如雅各布斯塔尔数的反转A001045号本质上是-(-2)^n*P(n,-1/2)。
多项式也可以表示为广义Narayana多项式之间的差异,请参阅公式部分。
链接
Xi Chen、Arthur Li Bo Yang、James Jing Yu Zhao、,Callan推广Narayana多项式的递归性,J.系统。科学。综合体(2021年)。
彼得·卢什尼,多项式的图解.
罗伯特·苏兰克,Narayana分布,J.Statist。计划。推理,2002,101:311-326,公式2。
公式
显式公式(加法形式):
T(n,n)=1,T(n>0,0)=0,否则T(n、k)=二项式(n,k)*二项式。
具有相同边界条件的乘法公式:
T(n,k)=二项式(n,k)^2*(k*(2*k^2+(n+1)*(n-2*k))。
二元生成函数:
T(n,k)=[x^n][y^k](1-x+(1+y)*(1-x*(y-1)-sqrt((x*y+x-1)^2-4*x^2*y))/2)。
基于多项式的递归:
T(n,k)=[x^k](((2*n-3)*(x+1)*P(n-1,x)-(n-3。
基于行的递归(请参阅第二个Python程序):
T(n,k)=(((B(k)+B(k-1))*(2*n-3)-(A(k)-2*A(k-1)+A(k-2))*(n-3))/n),其中A(k)=T(n-2,k)和B(k)=T(n-1,k),对于n>=3。
超几何表示:
对于n>=2,T(n,k)=[x^k]x*(x+1)*超几何([1-n,2-n],[2],x)。
行总和:
和{k=0..n}T(n,k)=(2/n)*二项式(2*(n-1),n-1)=A068875号(n-1)对于n>=2。
Narayana多项式的推广如下所示
N{N,k}(x)=和{j=0..N-2*k}。
N{N,0}(x)是经典的Narayana多项式A001263号而N{N,1}(x)是A145596号基于(3,2)。我们的多项式是n>=1的差分P(n,x)=n{n,0}(x)-n{n,1}(x)。
设RS(T,n)表示T的第n行的行和,则RS(T、n)-RS(A090181号,n)=-4*二项式(2*n-3,n-3)/(n+1)=A115143号(n+1),对于n>=3。
例子
三角形开始:
[0] 1;
[1] 0, 1;
[2] 0, 1, 1;
[3] 0, 1, 2, 1;
[4] 0, 1, 4, 4, 1;
[5] 0, 1, 7, 12, 7, 1;
[6] 0, 1, 11, 30, 30, 11, 1;
[7] 0, 1, 16, 65, 100, 65, 16, 1;
[8] 0, 1, 22, 126, 280, 280, 126, 22, 1;
[9] 0, 1, 29, 224, 686, 980, 686, 224, 29, 1;
MAPLE公司
T:=(n,k)->如果n=k,则1 elif k=0,否则0
二项式(n,k)^2*(k*(2*k^2+(n+1)*(n-2*k))
seq(seq(T(n,k),k=0..n),n=0..10);
#备选方案:
gf:=1-x+(1+y)*(1-x*(y-1)-sqrt((x*y+x-1)^2-4*x^2*y))/2:
serx:=展开(系列(gf,x,16)):系数:=n->系数(serx,x,n):
seq(seq(coeff(n),y,k),k=0..n),n=0..10);
#使用多项式递归:
P:=proc(n,x)选项记住;如果n<3,则[1,x,x+x^2][n+1]else
((2*n-3)*(x+1)*P(n-1,x)-(n-3)x(x-1)^2*P(n-2,x))/n结束:
Trow:=n->seq(系数(P(n,x),x,k),k=0..n):seq(Trow(n),n=0..10);
#用广义Narayana多项式表示:
N:=(N,k,x)->加法((k+1)/(N-k))*二项式(N-k,j-1)*二项式(N-k,j+k)*x^(j+k),j=0..N-2*k):seq(打印(如果其他(N=0,1,展开(N(N,0,x)-N(N,1,x)),N=0..7);
数学
H[0,_]:=1;H[1,x_]:=x;
H[n_,x_]:=x*(x+1)*超几何2F1[1-n,2-n,2,x];
Hrow[n_]:=系数列表[H[n,x],x];表[Hrow[n],{n,0,9}]//表格
黄体脂酮素
(Python)
从二项式的数学导入梳
定义T(n,k):
如果k==n:返回1
如果k==0:返回0
返回((二项式(n,k)**2*(k*(2*k**2+(n+1)*(n-2*k)))
//(n**2*(n-1)*(n-k+1))
def Trow(n):返回[T(n,k)for k in range(n+1)]
对于范围(10)中的n:打印(Trow(n))
(Python)#带缓存的递归速度更快:
从functools导入缓存
@高速缓存
定义行(n)(_R):
如果n<3:返回([1],[0,1],[0,1,1])[n]
A=T_行(n-2)+[0,0]
B=T_行(n-1)+[1]
对于范围(n-1,1,-1)中的k:
B[k]=((B[k]+B[k-1])*(2*n-3)
-(A[k]-2*A[k-1]+A[k-2])*(n-3))//n)
返回B
对于范围(10)中的n:打印(T_row(n))
交叉参考
囊性纤维变性。A090181号A001263号(纳拉亚纳),A353279型(情况3),A000108美元(加泰罗尼亚语),A145596号A172392号(中心术语),A000124号(子对角线,第2列),15143年.
基本上是加泰罗尼亚数字的两倍:A284016型(同时A068875号A002420型).
多项式序列的值:A068875号(行总和):P(1),A154955号:P(-1),2013年2月:P(2)/2,A125695号(同时A152681号):P(-2),A054872号:P(3)/2,P(3A234939型A336729型:P(-3)/6,A082298号:P(4)/5,A238113型:2^n*P(1/2),A154825号A091593号:2^n*P(-1/2)。
关键词
非n
作者
彼得·卢什尼2022年4月26日
状态
经核准的

查找|欢迎光临|维基|注册|音乐|地块2|演示|索引|浏览|更多|网络摄像头
贡献新序列。或评论|格式|样式表|变换|超级搜索|最近
OEIS社区|维护人OEIS基金会。

许可协议、使用条款、隐私政策。.

上次修改时间:美国东部夏令时2024年5月28日16:48。包含372916个序列。(在oeis4上运行。)