cp's OEIS Frontend

A217391 Partial sums of the squares of the ordered Bell numbers (number of preferential arrangements) A000670.

%I A217391 #27 Jun 24 2022 23:37:47
%S A217391 1,2,11,180,5805,298486,22228975,2258856824,300194704049,
%T A217391 50529463186170,10505093602625139,2643441560563225468,
%U A217391 791779611505017309493,278371498870260182630654,113516551713466910954246903,53143864598655784249290736512,28309328562668956145157858372537
%N A217391 Partial sums of the squares of the ordered Bell numbers (number of preferential arrangements) A000670.
%H A217391 Vincenzo Librandi, <a href="/A217391/b217391.txt">Table of n, a(n) for n = 0..200</a>
%F A217391 a(n) = Sum_{k=0..n} t(k)^2 where t = A000670 (ordered Bell numbers).
%F A217391 a(n) ~ (n!)^2 / (4 * (log(2))^(2*n+2)). - _Vaclav Kotesovec_, Nov 08 2014
%t A217391 t[n_] := Sum[StirlingS2[n, k]k!, {k, 0, n}]; Table[Sum[t[k]^2, {k, 0, n}], {n, 0, 100}]
%o A217391 (Maxima)
%o A217391 t(n):=sum(stirling2(n,k)*k!,k,0,n);
%o A217391 makelist(sum(t(k)^2,k,0,n),n,0,40);
%o A217391 (Magma)
%o A217391 A000670:=func<n | &+[StirlingSecond(n,i)*Factorial(i): i in [0..n]]>;
%o A217391 [&+[A000670(k)^2: k in [0..n]]: n in [0..14]]; // _Bruno Berselli_, Oct 03 2012
%o A217391 (PARI) for(n=0,30, print1(sum(k=0,n, (sum(j=0,k, j!*stirling(k,j,2)))^2), ", ")) \\ _G. C. Greubel_, Feb 07 2018
%Y A217391 Cf. A000670, A006957, A005649, A217388, A217389, A217392.
%Y A217391 Partial sums of A122725.
%K A217391 nonn
%O A217391 0,2
%A A217391 _Emanuele Munarini_, Oct 02 2012