%I#16 2022年6月26日09:01:39

%S 1,121451600170321791321883117199293902133321012316793121，

%电话25577181324287421068697329039439709738393883291996，

%电话：45675345280069155405974390088616853048954734169786421860831500723011109867095322250145563654356205885737

%N在N的所有集合分区中，大小>=10的块数。

%H Alois P.Heinz，n的表，n的a（n）=10..576</a>

%H维基百科，<a href=“https://en.wikipedia.org/wiki/Partition_of_a_set（英文）“>集合的分区</a>

%F a（n）=贝尔（n+1）-和{j=0..9}二项式（n，j）*贝尔（n-j）。

%F a（n）=和{j=0..n-10}二项式（n，j）*Bell（j）。

%例如：（exp（x）-和{k=0..9}x^k/k！）*经验（经验（x）-1）_伊利亚·古特科夫斯基，2022年6月26日

%p b:=proc（n）选项记忆`如果`（n=0，1，添加(

%p b（n-j）*二项式（n-1，j-1），j=1..n）

%p端：

%p g：=proc（n，k）选项记忆`如果`（n<k，0，

%pg（n，k+1）+二项式（n，k）*b（n-k））

%p端：

%pa:=n->g（n，10）：

%p序列（a（n），n=10..30）；

%t表[总和[二项式[n，j]BellB[j]，{j，0，n-10}]，{n，10，30}]（*_Indranil Ghosh_，2017年7月6日*）

%o（Python）

%o来自sympy输入bell，二项式

%o定义a（n）：返回和（[二项式（n，j）*bell（j）for j in range（n-9）]）

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

%A283424的Y列k=10。

%Y参考A000110。

%K nonn公司

%O 10,2号机组

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