%I#34 2021年12月23日06:02:57
%S 1,-1,2,-3,-20370,-407434293,-1383312,-2932533106271090,-2192834490,
%电话:32208497124,-206343936097,-7657279887698412496622532785,
%电话:12455477719752976260294034150380430,电话:2256541295745391542,电话:122593550603395508438728842979656718306780
%N a(N)=N!*[x^n]经验(n*(1-exp(x)))。
%H Alois P.Heinz,n的表格,n=0..415的a(n)</a>
%F a(n)=exp(n)*Sum_{k>=0}(-n)^k*k^n/k!.-_伊利亚·古特科夫斯基,2019年7月13日
%F a(n)=和{k=0..n}(-n)^k*斯特林2(n,k).-_Seiichi Manyama,2019年7月28日
%F a(n)=贝尔多项式(n,-n).-_Peter Luschny_,2021年12月23日
%p b:=proc(n,k)选项记忆`如果`(n=0,1,
%p-(1+加法(二项式(n-1,j-1)*b(n-j,k),j=1..n-1))*k)
%p端:
%pa:=n->b(n$2):
%p序列(a(n),n=0..30);#_Alois P.Heinz_,2017年9月25日
%t表[n!*系列系数[E^(n*(1-E^x)),{x,0,n}],{n,0,20}](*_Vaclav Kotesovec_,2017年9月25日*)
%t a[n_]:=贝尔B[n,-n];表[a[n],{n,0,20}](*_Peter Luschny_,2021年12月23日*)
%o(红宝石)
%o定义ncr(n,r)
%o如果r==0,则返回1
%o(n-r+1..n).注入(:*)/(1..r).注入
%o端
%o定义A(k,n)
%o元=[1]
%o(1..n).每个{i|ary<<k*(0..i-1).注入(0){s,j|s+ncr(i-1,j)*ary[j]}}
%o ary公司
%o端
%o定义A292866(n)
%o(0..n).map{|i|A(-i,i)[-1]}
%o端
%o p A292866(20)
%o(PARI){a(n)=总和(k=0,n,(-n)^k*stirling(n,k,2))}\\_Seiichi Manyama_,2019年7月28日
%Y A292861的主对角线。
%Y参考A242817。
%K符号
%0、3
%A _Seiichi Manyama,2017年9月25日
|