%I#14 2020年4月13日09:32:48
%S 1,2,4,1,7,3,12,7,3,1,20,15,8,3,1,33,30,19,9,3,154,58,42,23,10,3,1,
%电话88109,89,55,27,11,3143201182125,69,31,12,3,12365363273,
%U 166,84,35,13,3,1376655709579383212100,39,14,3,1609
%与A207611联合生成的多项式u(N,x)的系数的N三角形;请参阅“公式”部分。
%F u(n,x)=u(n-1,x)+v。
%e前五排:
%第1页
%第2页
%e 4…1
%e 7…3…1
%e 12…7…3…1
%t u[1,x_]:=1;v[1,x_]:=1;z=16;
%t u[n,x_]:=u[n-1,x]+v[n-1、x]
%tv[n,x_]:=u[n-1,x]+x*v[n-1、x]+1
%t表[系数[u[n,x]],{n,1,z}]
%t表[Factor[v[n,x]],{n,1,z}]
%t cu=表[系数列表[u[n,x],x]、{n,1,z}];
%t表格[cu]
%t压扁[%](*A207610*)
%t表格[展开[v[n,x]],{n,1,z}]
%t cv=表[系数列表[v[n,x],{n,1,z}];
%t表格[cv]
%t压扁[%](*A207611*)
%o(Python)
%o从sympy导入Poly
%o从sympy.abc导入x
%o定义u(n,x):如果n==1,则返回1,否则u(n-1,x)+v(n-1、x)
%o定义v(n,x):如果n==1,则返回1,否则u(n-1,x)+x*v(n-1、x)+1
%o定义a(n):返回多边形(u(n,x),x).all_coeffs()[::-1]
%o对于范围(1,13)中的n:打印(a(n))#_Indranil Ghosh,2017年5月28日
%Y参考A207611。
%Y参见A000071(第1列)、A023610(第2列)。
%K nonn,标签
%O 1,2号机组
%A_Clark Kimberling_,2012年2月19日