%I#19 2022年1月24日16:53:46

%编号：1,2,8410800285768012859560000880599202560853262368358400，

%电话：1111400775560275200187327646047473280003967400888465895264384000，

%电话：103139980547138966296473600322919706180911082676956696000119851575591550917401545594880000

%N[2n]上使图像大小等于N的内函数数。

%H Alois P.Heinz，n的表格，n=0..195的a（n）</a>

%F a（n）=箍筋2（2*n，n）*n！*二项式（2*n，n）。

%F a（n）=A090657（2n，n）=A101817（2n，n）=A219859（2n、n）。

%F a（n）~n^（2*n-1/2）*2^（4*n）/（sqrt（Pi*（1-c））*c^n*（2-c）^n*exp（2*n）），其中c=-LambertW（-2*exp

%e a（1）=2:（1,1），（2,2）。

%p b:=proc（n，k）选项记忆`如果`（k=n，n！，

%p`如果`（k=0，0，n*（b（n-1，k-1）+b（n-1，k）*k/（n-k）））

%p端：

%pa:=n->b（2*n，n）：

%p序列（a（n），n=0..15）；

%t表[StirlingS2[2*n，n]*（2*n）/n！，{n，0，20}]（*_Vaclav Kotesovec_，2017年6月10日*）

%o（PARI）a（n）=斯特林（2*n，n，2）*n*二项式（2*n，n）；\\_Indranil Ghosh，2017年7月4日

%o（Python）

%o来自mpmath导入*

%o最大dps=100

%o定义a（n）：返回int（stirling2（2*n，n）*fac（n）*二项式（2*n，n））

%o打印（[a（n）代表范围（21）中的n）]#_Indranil Ghosh，2017年7月4日

%Y参见A000142、A048993、A090657、A101817、A219859。

%K nonn公司

%0、2

%A _Alois P.Heinz，2017年6月7日