cp's OEIS Frontend

A144331 Triangle b(n,k) for n >= 0, 0 <= k <= 2n, read by rows. See A144299 for definition and properties.

%I A144331 #25 Oct 05 2023 08:36:24
%S A144331 1,0,1,1,0,0,1,3,3,0,0,0,1,6,15,15,0,0,0,0,1,10,45,105,105,0,0,0,0,0,
%T A144331 1,15,105,420,945,945,0,0,0,0,0,0,1,21,210,1260,4725,10395,10395,0,0,
%U A144331 0,0,0,0,0,1,28,378,3150,17325,62370,135135,135135,0,0,0,0,0,0
%N A144331 Triangle b(n,k) for n >= 0, 0 <= k <= 2n, read by rows. See A144299 for definition and properties.
%C A144331 Although this entry is the last of the versions of the underlying triangle to be added to the OEIS, for some applications it is the most important.
%C A144331 Row n has 2n+1 entries.
%C A144331 A001498 has a b-file.
%H A144331 Reinhard Zumkeller, <a href="/A144331/b144331.txt">Rows n = 0..100 of triangle, flattened</a>
%H A144331 Moa Apagodu, David Applegate, N. J. A. Sloane, and Doron Zeilberger, <a href="http://arxiv.org/abs/1701.08394">Analysis of the Gift Exchange Problem</a>, arXiv:1701.08394 [math.CO], 2017.
%H A144331 David Applegate and N. J. A. Sloane, <a href="http://arxiv.org/abs/0907.0513">The Gift Exchange Problem</a> (arXiv:0907.0513, 2009)
%F A144331 E.g.f.: Sum_{n >= 0} Sum_{k = 0..2n} b(n,k) y^n * x^k/k! = exp(x*y*(1 + x/2)).
%F A144331 b(n, k) = 2^(n-k)*k!/((2*n-k)!*(k-n)!).
%F A144331 Sum_{k=0..2*n} b(n, k) = A001515(n).
%F A144331 Sum_{n >= 0} b(n, k) = A000085(k).
%F A144331 From _G. C. Greubel_, Oct 04 2023: (Start)
%F A144331 T(n, k) = 0 for 0 <= k <= n-1, otherwise T(n, k) = k!/(2^(k-n)*(k-n)!*(2*n-k)!) for n <= k <= 2*n.
%F A144331 Sum_{k=0..2*n} (-1)^k * T(n, k) = A278990(n). (End)
%e A144331 Triangle begins:
%e A144331   1
%e A144331   0 1 1
%e A144331   0 0 1 3 3
%e A144331   0 0 0 1 6 15 15
%e A144331   0 0 0 0 1 10 45 105 105
%e A144331   0 0 0 0 0  1 15 105 420  945  945
%e A144331   0 0 0 0 0  0  1  21 210 1260 4725 10395 10395
%e A144331   ...
%t A144331 Flatten[Table[PadLeft[Table[(n+k)!/(2^k*k!*(n-k)!), {k,0,n}], 2*n+1, 0], {n,0,12}]] (* _Jean-François Alcover_, Oct 14 2011 *)
%o A144331 (Haskell)
%o A144331 a144331 n k = a144331_tabf !! n !! k
%o A144331 a144331_row n = a144331_tabf !! n
%o A144331 a144331_tabf = iterate (\xs ->
%o A144331   zipWith (+) ([0] ++ xs ++ [0]) $ zipWith (*) (0:[0..]) ([0,0] ++ xs)) [1]
%o A144331 -- _Reinhard Zumkeller_, Nov 24 2014
%o A144331 (Magma)
%o A144331 A144331:= func< n,k | k le n-1 select 0 else Factorial(k)/(2^(k-n)*Factorial(k-n)*Factorial(2*n-k)) >;
%o A144331 [A144331(n,k): k in [0..2*n], n in [0..12]]; // _G. C. Greubel_, Oct 04 2023
%o A144331 (SageMath)
%o A144331 def A144331(n, k): return 0 if k<n else factorial(k)/(2^(k-n)*factorial(2*n-k)*factorial(k-n))
%o A144331 flatten([[A144331(n,k) for k in range(2*n+1)] for n in range(13)]) # _G. C. Greubel_, Oct 04 2023
%Y A144331 Cf. A144299, A278990.
%Y A144331 Row sums give A001515, column sums A000085.
%Y A144331 Other versions of this triangle are given in A001497, A001498, A111924 and A100861.
%Y A144331 See A144385 for a generalization.
%K A144331 nonn,tabf,nice
%O A144331 0,8
%A A144331 _David Applegate_ and _N. J. A. Sloane_, Dec 07 2008