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.

A140998 Triangle G(n, k), read by rows, for 0 <= k <= n, where G(n, 0) = G(n+1, n+1) = 1, G(n+2, n+1) = 2, and G(n+3, m) = G(n+1, m-1) + G(n+1, m) + G(n+2, m) for n >= 0 and m = 1..n+1.

This page as a plain text file.
%I A140998 #57 Dec 06 2020 06:27:44
%S A140998 1,1,1,1,2,1,1,4,2,1,1,7,5,2,1,1,12,11,5,2,1,1,20,23,12,5,2,1,1,33,46,
%T A140998 28,12,5,2,1,1,54,89,63,29,12,5,2,1,1,88,168,137,69,29,12,5,2,1,1,143,
%U A140998 311,289,161,70,29,12,5,2,1,1,232,567,594,367,168,70,29,12,5,2,1
%N A140998 Triangle G(n, k), read by rows, for 0 <= k <= n, where G(n, 0) = G(n+1, n+1) = 1, G(n+2, n+1) = 2, and G(n+3, m) = G(n+1, m-1) + G(n+1, m) + G(n+2, m) for n >= 0 and m = 1..n+1.
%C A140998 From _Petros Hadjicostas_, Jun 10 2019: (Start)
%C A140998 According to the attached picture, the index of asymmetry here is s = 1 and the index of obliqueness (or obliquity) is e = 0.
%C A140998 In the picture, the equation G(n, e*n) = 1 becomes G(n, 0) = 1, while the equations G(n+x+1, n-e*n+e*x-e+1) = 2^x for 0 <= x < s = 1 become G(n+1, n+1) = 1 and G(n+2, n+1) = 2.
%C A140998 Also, in the picture, the recurrence G(n+s+2, k) = G(n+1, k-e*s+e-1) + Sum_{m=1..s+1} G(n+m, k-e*s+m*e-2*e) for k = 1..n+1 becomes G(n+3, k) = G(n+1, k-1) + G(n+1, k) + G(n+2, k) for k = 1..n+1.
%C A140998 Except for a shifting of the indices by 1, this array is a mirror image of array A140993. We have G(n, k) = A140993(n+1, n-k+1) for 0 <= k <= n. Triangular array A140993 has the same index of asymmetry (i.e., s = 1) but index of obliqueness e = 1.
%C A140998 (End)
%H A140998 G. C. Greubel, <a href="/A140998/b140998.txt">Rows n = 0..100 of triangle, flatten</a>
%H A140998 Juri-Stepan Gerasimov, <a href="/A140998/a140998.jpg">Stepan's triangles and Pascal's triangle are connected by the recurrence relation ...</a>
%F A140998 From _Petros Hadjicostas_, Jun 10 2019: (Start)
%F A140998 G(n, k) = A140993(n+1, n-k+1) for 0 <= k <= n.
%F A140998 Let A(x,y) = Sum_{n,k >= 0} G(n, k)*x^n*y^k and B(x,y) = Sum_{n,k >= 1} A140993(n, k). Then A(x, y) = x^(-1) * B(x*y, y^(-1)). Thus, the g.f. of the current array is A(x, y) = (1 - x - x^2 + x^3*y)/((1 - x) * (1 - x*y) * (1 - x - x^2 - x^2*y)).
%F A140998 To find the g.f. of the k-th column (where k >= 0), we differentiate A(x, y) k times with respect to y, divide by k!, and substitute y = 0. For example, differentiating A(x, y) once w.r.t. y and setting y = 0, we get the g.f. of the k = 1 column: x/((1 - x)*(1 - x - x^2)). This is the g.f. of sequence (A000071(n+2): n >= 0) = (Fibonacci(n+2) - 1: n >= 0).
%F A140998 G.f. of column k = 2 is x^2*(1 - x + x^3)/((1 - x)*(1 - x - x^2)^2). Thus, column k = 2 is a shifted version of (A140992(n): n >= 0).
%F A140998 (End)
%e A140998 Triangle begins (with rows for n >= 0 and columns for k >= 0):
%e A140998   1;
%e A140998   1,   1;
%e A140998   1,   2,   1;
%e A140998   1,   4,   2,   1;
%e A140998   1,   7,   5,   2,   1;
%e A140998   1,  12,  11,   5,   2,   1;
%e A140998   1,  20,  23,  12,   5,   2,   1;
%e A140998   1,  33,  46,  28,  12,   5,   2,   1;
%e A140998   1,  54,  89,  63,  29,  12,   5,   2,   1;
%e A140998   1,  88, 168, 137,  69,  29,  12,   5,   2,   1;
%e A140998   1, 143, 311, 289, 161,  70,  29,  12,   5,   2,   1;
%t A140998 G[n_,k_] := G[n,k] = Which[k==0 || k==n, 1, k==n-1, 2, True, G[n-2,k-1] + G[n-2,k] + G[n-1,k]]; Table[G[n,k], {n,0,12}, {k,0,n}] (* _Jean-François Alcover_, Jun 09 2019 *)
%o A140998 (PARI) G(n,k) = if(k==0 || k==n, 1, if(k==n-1, 2, G(n-1, k) + G(n-2, k) + G(n-2, k-1)));
%o A140998 for(n=0,12, for(k=0,n, print1(G(n,k), ", "))) \\ _G. C. Greubel_, Jun 09 2019
%o A140998 (Sage)
%o A140998 def G(n,k):
%o A140998     if (k==0 or k==n): return 1
%o A140998     elif (k==n-1): return 2
%o A140998     else: return G(n-1, k) + G(n-2, k) + G(n-2, k-1)
%o A140998 [[G(n,k) for k in (0..n)] for n in (0..12)] # _G. C. Greubel_, Jun 09 2019
%Y A140998 Cf. A000071, A007318, A140992, A140993, A140994, A140995, A140996, A140997, A141020, A141021.
%K A140998 nonn,tabl
%O A140998 0,5
%A A140998 _Juri-Stepan Gerasimov_, Jul 08 2008
%E A140998 Indices in the definition corrected by _R. J. Mathar_, Aug 02 2009
%E A140998 Name edited by _Petros Hadjicostas_, Jun 10 2019