登录
G.f.:A(x)=exp(总和{n>=1}A000172号(n) *x^n/n)其中Franel编号A000172号(n) =和{k=0..n}C(n,k)^3。
15

%I#10 2017年11月27日09:43:34

%序号1,2,7,301477864472266441644771044258678248444887236,

%电话:30178236120562505701417279235598667874038692948001906,

%电话:49044034999923495112433730025061782908765618070552843977113095146839953030

%N G.f.:A(x)=exp(总和{N>=1}A000172(N)*x^N/N),其中Franel数A000172,(N)=总和{k=0..N}C(N,k)^3。

%C类似于加泰罗尼亚数字的g.f.平方(A000108):

%C C(x)^2=exp(总和{n>=1}A000984(n)*x^n/n),其中中心二项式系数A000983(n)=总和{k=0..n}C(n,k)^2。

%H G.C.Greubel,n表,n=0..500时的a(n)</a>

%F A166991的自卷积。

%F a(n)~c*8^n/n^2,其中c=0.58462…-_Vaclav Kotesovec_,2017年11月27日

%通用公式:A(x)=1+2*x+7*x^2+30*x^3+147*x^4+786*x^5+4472*x^6+。。。

%对数(A(x))=2*x+10*x^2/2+56*x^3/3+346*x*^4/4+2252*x^5/5+15184*x^6/6+104960*x^7/7+…+A000172(n)*x^n/n+。。。

%t a[n_]:=和[(二项式[n,k])^3,{k,0,n}];f[x_]:=和[a[n]*x^n/n,{n,1,75}];系数列表[系列[Exp[f[x]],{x,0,50}],x](*_G.C.Greubel_,2016年5月30日*)

%o(PARI){a(n)=极系数(exp(sum(m=1,n,sum(k=0,m,二项式(m,k)^3)*x^m/m)+x*o(x^n)),n)}

%Y参见A000172(法兰编号)、A166991、A1669962、A218117、A218119。

%K nonn公司

%0、2

%A _保罗·D·汉纳,2009年11月17日