%I #7 Aug 31 2016 05:31:04
%S 0,0,0,0,1,9,60,368,2218,13496,84010,538885,3576885,24624234,
%T 176021146,1306975716,10078487587,80670305273,669699040920,
%U 5760761324172,51293043905644,472213256230868,4489889820251839
%N Bell number A000110(n) minus Bessel number A006789(n).
%H G. C. Greubel, <a href="/A153820/b153820.txt">Table of n, a(n) for n = 0..200</a>
%F a(n) = A000110(n) - A006789(n).
%t nmax = 50; m = SparseArray[{{i_, i_} :> i, Band[{1, 2}] -> 1, Band[{2, 1}] -> 1}, {nmax, nmax}]; a[n_] := MatrixPower[m, n][[1, 1]]; Table[BellB[n] - a[n], {n, 0, nmax}] (* _G. C. Greubel_, Aug 30 2016, built using _Jean-François Alcover_'s code from A006789 *)
%Y Cf. A000110, A006789.
%K nonn
%O 0,6
%A _Omar E. Pol_, Jan 02 2009
%E More terms and correction of 9th term by _Dmitry Kamenetsky_, Jan 17 2009
|