(PARI){a(n)=和(k=0,n,二项式(n,k)*stirling(n,k,1)^2)};
(岩浆)[(&+[二项式(n,k)*StirlingFirst(n,k)^2:k in[0..n]]):n in[0..20]]//G.C.格鲁贝尔,2019年8月7日
(Sage)[sum(二项式(n,k)*stirling_number1(n,k)^2 for k in(0..n))for n in(0..20)]#G.C.格鲁贝尔2019年8月7日
(GAP)列表([0..20],n->总和([0..n],k->二项式(n,k)*Stirling1(n,k)^2)#G.C.格鲁贝尔2019年8月7日