cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A112468 Riordan array (1/(1-x), x/(1+x)).

This page as a plain text file.
%I A112468 #64 Aug 28 2022 08:40:31
%S A112468 1,1,1,1,0,1,1,1,-1,1,1,0,2,-2,1,1,1,-2,4,-3,1,1,0,3,-6,7,-4,1,1,1,-3,
%T A112468 9,-13,11,-5,1,1,0,4,-12,22,-24,16,-6,1,1,1,-4,16,-34,46,-40,22,-7,1,
%U A112468 1,0,5,-20,50,-80,86,-62,29,-8,1,1,1,-5,25,-70,130,-166,148,-91,37,-9,1,1,0,6,-30,95,-200,296,-314,239,-128,46,-10,1
%N A112468 Riordan array (1/(1-x), x/(1+x)).
%C A112468 Row sums are A040000. Diagonal sums are A112469. Inverse is A112467. Row sums of k-th power are 1, k+1, k+1, k+1, .... Note that C(n,k) = Sum_{j=0..n-k} C(n-j-1, n-k-j).
%C A112468 Equals row reversal of triangle A112555 up to sign, where log(A112555) = A112555 - I. Unsigned row sums equals A052953 (Jacobsthal numbers + 1). Central terms of even-indexed rows are a signed version of A072547. Sums of squared terms in rows yields A112556, which equals the first differences of the unsigned central terms. - _Paul D. Hanna_, Jan 20 2006
%C A112468 Sum_{k=0..n} T(n,k)*x^k = A000012(n), A040000(n), A005408(n), A033484(n), A048473(n), A020989(n), A057651(n), A061801(n), A238275(n), A238276(n), A138894(n), A090843(n), A199023(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 respectively (see the square array in A112739). - _Philippe Deléham_, Feb 22 2014
%H A112468 Reinhard Zumkeller, <a href="/A112468/b112468.txt">Rows n = 0..125 of triangle, flattened</a>
%H A112468 H. Belbachir and F. Bencherif, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL11/Belbachir/belbachir13.html">On some properties of bivariate Fibonacci and Lucas Polynomials</a>, JIS 11 (2008) 08.2.6.
%H A112468 Hacene Belbachir and Athmane Benmezai, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL15/Belbachir/bel22.html">Expansion of Fibonacci and Lucas Polynomials: An Answer to Prodinger's Question</a>, Journal of Integer Sequences, Vol. 15 (2012), #12.7.6.
%H A112468 Emeric Deutsch, L. Ferrari and S. Rinaldi, <a href="http://dx.doi.org/10.1016/j.aam.2004.05.002">Production Matrices</a>, Advances in Applied Mathematics, 34 (2005) pp. 101-122.
%H A112468 Kyu-Hwan Lee and Se-jin Oh, <a href="http://arxiv.org/abs/1601.06685">Catalan triangle numbers and binomial coefficients</a>, arXiv:1601.06685 [math.CO], 2016.
%H A112468 H. Prodinger, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL12/Prodinger/prodinger27.html">On the expansion of Fibonacci and Lucas Polynomials</a>, JIS 12 (2009) 09.1.6.
%F A112468 Triangle T(n,k) read by rows: T(n,0)=1, T(n,k) = T(n-1,k-1) - T(n-1,k). - _Mats Granvik_, Mar 15 2010
%F A112468 Number triangle T(n, k)= Sum_{j=0..n-k} C(n-j-1, n-k-j)*(-1)^(n-k-j).
%F A112468 G.f. of matrix power T^m: (1+(m-1)*x)*(1+m*x)/(1+m*x-x*y)/(1-x). G.f. of matrix log: x*(1-2*x*y+x^2*y)/(1-x*y)^2/(1-x). - _Paul D. Hanna_, Jan 20 2006
%F A112468 T(n, k) = R(n,n-k,-1) where R(n,k,m) = (1-m)^(-n+k)-m^(k+1)*Pochhammer(n-k,k+1)*hyper2F1([1,n+1],[k+2],m)/(k+1)!. - _Peter Luschny_, Jul 25 2014
%e A112468 Triangle starts
%e A112468   1;
%e A112468   1,  1;
%e A112468   1,  0,  1;
%e A112468   1,  1, -1,  1;
%e A112468   1,  0,  2, -2,  1;
%e A112468   1,  1, -2,  4, -3,  1;
%e A112468   1,  0,  3, -6,  7, -4,  1;
%e A112468 Matrix log begins:
%e A112468   0;
%e A112468   1,  0;
%e A112468   1,  0,  0;
%e A112468   1,  1, -1,  0;
%e A112468   1,  1,  1, -2,  0;
%e A112468   1,  1,  1,  1, -3,  0; ...
%e A112468 Production matrix begins
%e A112468   1,  1,
%e A112468   0, -1,  1,
%e A112468   0,  0, -1,  1,
%e A112468   0,  0,  0, -1,  1,
%e A112468   0,  0,  0,  0, -1,  1,
%e A112468   0,  0,  0,  0,  0, -1,  1,
%e A112468   0,  0,  0,  0,  0,  0, -1,  1.
%e A112468 - _Paul Barry_, Apr 08 2011
%p A112468 T := (n,k,m) -> (1-m)^(-n+k)-m^(k+1)*pochhammer(n-k,k+1)*hypergeom( [1,n+1],[k+2],m)/(k+1)!; A112468 := (n,k) -> T(n,n-k,-1);
%p A112468 seq(print(seq(simplify(A112468(n,k)),k=0..n)),n=0..10); # _Peter Luschny_, Jul 25 2014
%t A112468 T[n_, 0] = 1; T[n_, n_] = 1; T[n_, k_ ]:= T[n, k] = T[n-1, k-1] - T[n-1, k]; Table[T[n, k], {n, 0, 12}, {k, 0, n}]//Flatten (* _Jean-François Alcover_, Mar 06 2013 *)
%o A112468 (PARI) {T(n,k)=local(m=1,x=X+X*O(X^n),y=Y+Y*O(Y^k)); polcoeff(polcoeff((1+(m-1)*x)*(1+m*x)/(1+m*x-x*y)/(1-x),n,X),k,Y)} \\ _Paul D. Hanna_, Jan 20 2006
%o A112468 (Haskell)
%o A112468 a112468 n k = a112468_tabl !! n !! k
%o A112468 a112468_row n = a112468_tabl !! n
%o A112468 a112468_tabl = iterate (\xs -> zipWith (-) ([2] ++ xs) (xs ++ [0])) [1]
%o A112468 -- _Reinhard Zumkeller_, Jan 03 2014
%o A112468 (PARI) T(n,k) = if(k==0 || k==n, 1, T(n-1, k-1) - T(n-1, k)); \\ _G. C. Greubel_, Nov 13 2019
%o A112468 (Magma)
%o A112468 function T(n,k)
%o A112468   if k eq 0 or k eq n then return 1;
%o A112468   else return T(n-1,k-1) - T(n-1,k);
%o A112468   end if;
%o A112468   return T;
%o A112468 end function;
%o A112468 [T(n,k): k in [0..n], n in [0..12]]; // _G. C. Greubel_, Nov 13 2019
%o A112468 (Sage)@CachedFunction
%o A112468 def T(n, k):
%o A112468     if (k<0 or n<0): return 0
%o A112468     elif (k==0 or k==n): return 1
%o A112468     else: return T(n-1, k-1) - T(n-1, k)
%o A112468 [[T(n, k) for k in (0..n)] for n in (0..12)] # _G. C. Greubel_, Nov 13 2019
%o A112468 (GAP)
%o A112468 T:= function(n,k)
%o A112468     if k=0 or k=n then return 1;
%o A112468     else return T(n-1,k-1) - T(n-1,k);
%o A112468     fi;
%o A112468   end;
%o A112468 Flat(List([0..12], n-> List([0..n], k-> T(n,k) ))); # _G. C. Greubel_, Nov 13 2019
%Y A112468 Cf. A174294, A174295, A174296, A174297. - _Mats Granvik_, Mar 15 2010
%Y A112468 Cf. A072547 (central terms), A112555 (reversed rows), A112465, A052953, A112556, A112739, A119258.
%Y A112468 See A279006 for another version.
%K A112468 easy,sign,tabl
%O A112468 0,13
%A A112468 _Paul Barry_, Sep 06 2005