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.

Showing 1-10 of 35 results. Next

A302764 Pascal-like triangle with A000012 as the left border and A080956 as the right border.

Original entry on oeis.org

1, 1, 1, 1, 2, 0, 1, 3, 2, -2, 1, 4, 5, 0, -5, 1, 5, 9, 5, -5, -9, 1, 6, 14, 14, 0, -14, -14, 1, 7, 20, 28, 14, -14, -28, -20, 1, 8, 27, 48, 42, 0, -42, -48, -27, 1, 9, 35, 75, 90, 42, -42, -90, -75, -35, 1, 10, 44, 110, 165, 132, 0, -132, -165, -110, -44
Offset: 1

Views

Author

Gregory Gerard Wojnar, Apr 12 2018

Keywords

Comments

Number the rows of the triangle beginning with n=0. For each row construct a degree n polynomial with regularly decreasing powers, denoting the polynomial as f_n(x); e.g., for row 2 we have f_2(x)=1x^2+2x+0. Then construct g_n(x)=x^2*f_{n-1}(x)-(n+1)x+1. It obtains that g_n(x)=(1-x)(2-(1+x)^n). These g_n(x) are the denominators of the generating functions for the following sequences: A024537 (n=2); A195350 (n=3); A301417 (n=4); A301420 (n=5); A301421 (n=6); A301424 (n=7). For these sequences the asymptotic term-to-term ratios are 1/(2^(1/n)-1). The numerators of the generating functions are 1-x(x+1)^(n-1).

Examples

			Triangle begins:
  1;
  1, 1;
  1, 2,  0;
  1, 3,  2, -2;
  1, 4,  5,  0, -5;
  1, 5,  9,  5, -5,  -9;
  1, 6, 14, 14,  0, -14, -14;
  1, 7, 20, 28, 14, -14, -28, -20;
  ...
		

Crossrefs

Programs

  • PARI
    T(n,k) = if (k==0, 1, if (k==n, (n+1)*(2-n)/2, if (k>n, 0, T(n-1,k) + T(n-1,k-1))));
    tabl(nn) = for (n=0, nn, for (k=0, n, print1(T(n, k), ", ")); print); \\ Michel Marcus, Apr 21 2018

Formula

T(n,k) = T(n-1,k) + T(n-1,k-1) with T(n, 0) = 1 and T(n, n) = (n+1)*(2-n)/2.

A005563 a(n) = n*(n+2) = (n+1)^2 - 1.

Original entry on oeis.org

0, 3, 8, 15, 24, 35, 48, 63, 80, 99, 120, 143, 168, 195, 224, 255, 288, 323, 360, 399, 440, 483, 528, 575, 624, 675, 728, 783, 840, 899, 960, 1023, 1088, 1155, 1224, 1295, 1368, 1443, 1520, 1599, 1680, 1763, 1848, 1935, 2024, 2115, 2208, 2303, 2400, 2499, 2600
Offset: 0

Views

Author

Keywords

Comments

Erdős conjectured that n^2 - 1 = k! has a solution if and only if n is 5, 11 or 71 (when k is 4, 5 or 7).
Second-order linear recurrences y(m) = 2y(m-1) + a(n)*y(m-2), y(0) = y(1) = 1, have closed form solutions involving only powers of integers. - Len Smiley, Dec 08 2001
Number of edges in the join of two cycle graphs, both of order n, C_n * C_n. - Roberto E. Martinez II, Jan 07 2002
Let k be a positive integer, M_n be the n X n matrix m_(i,j) = k^abs(i-j) then det(M_n) = (-1)^(n-1)*a(k-1)^(n-1). - Benoit Cloitre, May 28 2002
Also numbers k such that 4*k + 4 is a square. - Cino Hilliard, Dec 18 2003
For each term k, the function sqrt(x^2 + 1), starting with 1, produces an integer after k iterations. - Gerald McGarvey, Aug 19 2004
a(n) mod 3 = 0 if and only if n mod 3 > 0: a(A008585(n)) = 2; a(A001651(n)) = 0; a(n) mod 3 = 2*(1-A079978(n)). - Reinhard Zumkeller, Oct 16 2006
a(n) is the number of divisors of a(n+1) that are not greater than n. - Reinhard Zumkeller, Apr 09 2007
Nonnegative X values of solutions to the equation X^3 + X^2 = Y^2. To find Y values: b(n) = n(n+1)(n+2). - Mohamed Bouhamida, Nov 06 2007
Sequence allows us to find X values of the equation: X + (X + 1)^2 + (X + 2)^3 = Y^2. To prove that X = n^2 + 2n: Y^2 = X + (X + 1)^2 + (X + 2)^3 = X^3 + 7*X^2 + 15X + 9 = (X + 1)(X^2 + 6X + 9) = (X + 1)*(X + 3)^2 it means: (X + 1) must be a perfect square, so X = k^2 - 1 with k>=1. we can put: k = n + 1, which gives: X = n^2 + 2n and Y = (n + 1)(n^2 + 2n + 3). - Mohamed Bouhamida, Nov 12 2007
From R. K. Guy, Feb 01 2008: (Start)
Toads and Frogs puzzle:
This is also the number of moves that it takes n frogs to swap places with n toads on a strip of 2n + 1 squares (or positions, or lily pads) where a move is a single slide or jump, illustrated for n = 2, a(n) = 8 by
T T - F F
T - T F F
T F T - F
T F T F -
T F - F T
- F T F T
F - T F T
F F T - T
F F - T T
I was alerted to this by the Holton article, but on consulting Singmaster's sources, I find that the puzzle goes back at least to 1867.
Probably the first to publish the number of moves for n of each animal was Edouard Lucas in 1883. (End)
a(n+1) = terms of rank 0, 1, 3, 6, 10 = A000217 of A120072 (3, 8, 5, 15). - Paul Curtz, Oct 28 2008
Row 3 of array A163280, n >= 1. - Omar E. Pol, Aug 08 2009
Final digit belongs to a periodic sequence: 0, 3, 8, 5, 4, 5, 8, 3, 0, 9. - Mohamed Bouhamida, Sep 04 2009 [Comment edited by N. J. A. Sloane, Sep 24 2009]
Let f(x) be a polynomial in x. Then f(x + n*f(x)) is congruent to 0 (mod f(x)); here n belongs to N. There is nothing interesting in the quotients f(x + n*f(x))/f(x) when x belongs to Z. However, when x is irrational these quotients consist of two parts, a) rational integers and b) integer multiples of x. The present sequence represents the non-integer part when the polynomial is x^2 + x + 1 and x = sqrt(2), f(x+n*f(x))/f(x) = A056108(n) + a(n)*sqrt(2). - A.K. Devaraj, Sep 18 2009
For n >= 1, a(n) is the number for which 1/a(n) = 0.0101... (A000035) in base (n+1). - Rick L. Shepherd, Sep 27 2009
For n > 0, continued fraction [n, 1, n] = (n+1)/a(n); e.g., [6, 1, 6] = 7/48. - Gary W. Adamson, Jul 15 2010
Starting (3, 8, 15, ...) = binomial transform of [3, 5, 2, 0, 0, 0, ...]; e.g., a(3) = 15 = (1*3 + 2*5 +1*2) = (3 + 10 + 2). - Gary W. Adamson, Jul 30 2010
a(n) is essentially the case 0 of the polygonal numbers. The polygonal numbers are defined as P_k(n) = Sum_{i=1..n} ((k-2)*i-(k-3)). Thus P_0(n) = 2*n-n^2 and a(n) = -P_0(n+2). See also A067998 and for the case k=1 A080956. - Peter Luschny, Jul 08 2011
a(n) is the maximal determinant of a 2 X 2 matrix with integer elements from {1, ..., n+1}, so the maximum determinant of a 2x2 matrix with integer elements from {1, ..., 5} = 5^2 - 1 = a(4) = 24. - Aldo González Lorenzo, Oct 12 2011
Using four consecutive triangular numbers t1, t2, t3 and t4, plot the points (0, 0), (t1, t2), and (t3, t4) to create a triangle. Twice the area of this triangle are the numbers in this sequence beginning with n = 1 to give 8. - J. M. Bergot, May 03 2012
Given a particle with spin S = n/2 (always a half-integer value), the quantum-mechanical expectation value of the square of the magnitude of its spin vector evaluates to = S(S+1) = n(n+2)/4, i.e., one quarter of a(n) with n = 2S. This plays an important role in the theory of magnetism and magnetic resonance. - Stanislav Sykora, May 26 2012
Twice the harmonic mean [H(x, y) = (2*x*y)/(x + y)] of consecutive triangular numbers A000217(n) and A000217(n+1). - Raphie Frank, Sep 28 2012
Number m such that floor(sqrt(m)) = floor(m/floor(sqrt(m))) - 2 for m > 0. - Takumi Sato, Oct 10 2012
The solutions of equation 1/(i - sqrt(j)) = i + sqrt(j), when i = (n+1), j = a(n). For n = 1, 2 + sqrt(3) = 3.732050.. = A019973. For n = 2, 3 + sqrt(8) = 5.828427... = A156035. - Kival Ngaokrajang, Sep 07 2013
The integers in the closed form solution of a(n) = 2*a(n-1) + a(m-2)*a(n-2), n >= 2, a(0) = 0, a(1) = 1 mentioned by Len Smiley, Dec 08 2001, are m and -m + 2 where m >= 3 is a positive integer. - Felix P. Muga II, Mar 18 2014
Let m >= 3 be a positive integer. If a(n) = 2*a(n-1) + a(m-2) * a(n-2), n >= 2, a(0) = 0, a(1) = 1, then lim_{n->oo} a(n+1)/a(n) = m. - Felix P. Muga II, Mar 18 2014
For n >= 4 the Szeged index of the wheel graph W_n (with n + 1 vertices). In the Sarma et al. reference, Theorem 2.7 is incorrect. - Emeric Deutsch, Aug 07 2014
If P_{k}(n) is the n-th k-gonal number, then a(n) = t*P_{s}(n+2) - s*P_{t}(n+2) for s=t+1. - Bruno Berselli, Sep 04 2014
For n >= 1, a(n) is the dimension of the simple Lie algebra A_n. - Wolfdieter Lang, Oct 21 2015
Finding all positive integers (n, k) such that n^2 - 1 = k! is known as Brocard's problem, (see A085692). - David Covert, Jan 15 2016
For n > 0, a(n) mod (n+1) = a(n) / (n+1) = n. - Torlach Rush, Apr 04 2016
Conjecture: When using the Sieve of Eratosthenes and sieving (n+1..a(n)), with divisors (1..n) and n>0, there will be no more than a(n-1) composite numbers. - Fred Daniel Kline, Apr 08 2016
a(n) mod 8 is periodic with period 4 repeating (0,3,0,7), that is a(n) mod 8 = 5/2 - (5/2) cos(n*Pi) - sin(n*Pi/2) + sin(3*n*Pi/2). - Andres Cicuttin, Jun 02 2016
Also for n > 0, a(n) is the number of times that n-1 occurs among the first (n+1)! terms of A055881. - R. J. Cano, Dec 21 2016
The second diagonal of composites (the only prime is number 3) from the right on the Klauber triangle (see Kival Ngaokrajang link), which is formed by taking the positive integers and taking the first 1, the next 3, the following 5, and so on, each centered below the last. - Charles Kusniec, Jul 03 2017
Also the number of independent vertex sets in the n-barbell graph. - Eric W. Weisstein, Aug 16 2017
Interleaving of A000466 and A033996. - Bruce J. Nicholson, Nov 08 2019
a(n) is the number of degrees of freedom in a triangular cell for a Raviart-Thomas or Nédélec first kind finite element space of order n. - Matthew Scroggs, Apr 22 2020
From Muge Olucoglu, Jan 19 2021: (Start)
For n > 1, a(n-2) is the maximum number of elements in the second stage of the Quine-McCluskey algorithm whose minterms are not covered by the functions of n bits. At n=3, we have a(3-2) = a(1) = 1*(1+2) = 3 and f(A,B,C) = sigma(0,1,2,5,6,7).
.
0 1 2 5 6 7
+---------------
*(0,1)| X X
(0,2)| X X
(1,5)| X X
*(2,6)| X X
*(5,7)| X X
(6,7)| X X
.
*: represents the elements that are covered. (End)
1/a(n) is the ratio of the sum of the first k odd numbers and the sum of the next n*k odd numbers. - Melvin Peralta, Jul 15 2021
For n >= 1, the continued fraction expansion of sqrt(a(n)) is [n; {1, 2n}]. - Magus K. Chu, Sep 09 2022
Number of diagonals parallel to an edge in a regular (2*n+4)-gon (cf. A367204). - Paolo Xausa, Nov 21 2023
For n >= 1, also the number of minimum cyclic edge cuts in the (n+2)-trapezohedron graph. - Eric W. Weisstein, Nov 21 2024
For n >= 1, a(n) is the sum of the interior angles of a polygon with n+2 sides, in radians, multiplied by (n+2)/Pi. - Stuart E Anderson, Aug 06 2025

Examples

			G.f. = 3*x + 8*x^2 + 15*x^3 + 24*x^4 + 35*x^5 + 48*x^6 + 63*x^7 + 80*x^8 + ...
		

References

  • E. R. Berlekamp, J. H. Conway and R. K. Guy, Winning Ways, Academic Press, NY, 2 vols., 1982, see index under Toads and Frogs Puzzle.
  • Martin Gardner, Perplexing Puzzles and Tantalizing Teasers, p. 21 (for "The Dime and Penny Switcheroo").
  • R. K. Guy, Unsolved Problems in Theory of Numbers, Section D25.
  • Derek Holton, Math in School, 37 #1 (Jan 2008) 20-22.
  • Edouard Lucas, Récréations Mathématiques, Gauthier-Villars, Vol. 2 (1883) 141-143.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

Formula

G.f.: x*(3-x)/(1-x)^3. - Simon Plouffe in his 1992 dissertation
a(n) = A000290(n+1) - 1.
A002378(a(n)) = A002378(n)*A002378(n+1); e.g., A002378(15)=240=12*20. - Charlie Marion, Dec 29 2003
a(n) = A067725(n)/3. - Zerinvary Lajos, Mar 06 2007
a(n) = Sum_{k=1..n} A144396(k). - Zerinvary Lajos, May 11 2007
a(n) = A134582(n+1)/4. - Zerinvary Lajos, Feb 01 2008
A143053(a(n)) = A000290(n+1), for n > 0. - Reinhard Zumkeller, Jul 20 2008
a(n) = Real((n+1+i)^2). - Gerald Hillier, Oct 12 2008
A053186(a(n)) = 2*n. - Reinhard Zumkeller, May 20 2009
a(n) = (n! + (n+1)!)/(n-1)!, n > 0. - Gary Detlefs, Aug 10 2009
a(n) = floor(n^5/(n^3+1)) with offset 1 (a(1)=0). - Gary Detlefs, Feb 11 2010
a(n) = a(n-1) + 2*n + 1 (with a(0)=0). - Vincenzo Librandi, Nov 18 2010
Sum_{n>=1} 1/a(n) = 3/4. - Mohammad K. Azarian, Dec 29 2010
a(n) = 2/(Integral_{x=0..Pi/2} (sin(x))^(n-1)*(cos(x))^3), for n > 0. - Francesco Daddi, Aug 02 2011
a(n) = A002378(n) + floor(sqrt(A002378(n))); pronic number + its root. - Fred Daniel Kline, Sep 16 2011
a(n-1) = A008833(n) * A068310(n) for n > 1. - Reinhard Zumkeller, Nov 26 2011
G.f.: U(0) where U(k) = -1 + (k+1)^2/(1 - x/(x + (k+1)^2/U(k+1))); (continued fraction, 3-step). - Sergei N. Gladkovskii, Oct 19 2012
a(n) = 15*C(n+4,3)*C(n+4,5)/(C(n+4,2)*C(n+4,4)). - Gary Detlefs, Aug 05 2013
a(n) = (n+2)!/((n-1)! + n!), n > 0. - Ivan N. Ianakiev, Nov 11 2013
a(n) = 3*C(n+1,2) - C(n,2) for n >= 0. - Felix P. Muga II, Mar 11 2014
a(n) = (A016742(n+1) - 4)/4 for n >= 0. - Felix P. Muga II, Mar 11 2014
a(-2 - n) = a(n) for all n in Z. - Michael Somos, Aug 07 2014
A253607(a(n)) = 1. - Reinhard Zumkeller, Jan 05 2015
E.g.f.: x*(x + 3)*exp(x). - Ilya Gutkovskiy, Jun 03 2016
For n >= 1, a(n^2 + n - 2) = a(n-1) * a(n). - Miko Labalan, Oct 15 2017
Sum_{n>=1} (-1)^(n+1)/a(n) = 1/4. - Amiram Eldar, Nov 04 2020
From Amiram Eldar, Feb 17 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = 2.
Product_{n>=1} (1 - 1/a(n)) = -sqrt(2)*sin(sqrt(2)*Pi)/Pi. (End)
a(n) = A000290(n+2) - n*2. See Bounded Squares illustration. - Leo Tavares, Oct 05 2021
From Leo Tavares, Oct 10 2021: (Start)
a(n) = A008585(n) + 2*A000217(n-1). See Trapezoids illustration.
2*A005563 = A054000(n+1). See Trapagons illustration.
a(n) = 2*A000217(n) + n. (End)
a(n) = (n+2)!!/(n-2)!! for n > 1. - Jacob Szlachetka, Jan 02 2022

Extensions

Partially edited by Joerg Arndt, Mar 11 2010
More terms from N. J. A. Sloane, Aug 01 2010

A139600 Square array T(n,k) = n*(k-1)*k/2+k, of nonnegative numbers together with polygonal numbers, read by antidiagonals upwards.

Original entry on oeis.org

0, 0, 1, 0, 1, 2, 0, 1, 3, 3, 0, 1, 4, 6, 4, 0, 1, 5, 9, 10, 5, 0, 1, 6, 12, 16, 15, 6, 0, 1, 7, 15, 22, 25, 21, 7, 0, 1, 8, 18, 28, 35, 36, 28, 8, 0, 1, 9, 21, 34, 45, 51, 49, 36, 9, 0, 1, 10, 24, 40, 55, 66, 70, 64, 45, 10, 0, 1, 11, 27, 46, 65, 81, 91, 92, 81, 55, 11
Offset: 0

Views

Author

Omar E. Pol, Apr 27 2008

Keywords

Comments

A general formula for polygonal numbers is P(n,k) = (n-2)*(k-1)*k/2 + k, where P(n,k) is the k-th n-gonal number.
The triangle sums, see A180662 for their definitions, link this square array read by antidiagonals with twelve different sequences, see the crossrefs. Most triangle sums are linear sums of shifted combinations of a sequence, see e.g. A189374. - Johannes W. Meijer, Apr 29 2011

Examples

			The square array of nonnegatives together with polygonal numbers begins:
=========================================================
....................... A   A   .   .   A    A    A    A
....................... 0   0   .   .   0    0    1    1
....................... 0   0   .   .   1    1    3    3
....................... 0   0   .   .   6    7    9    9
....................... 0   0   .   .   9    3    6    6
....................... 0   1   .   .   5    2    0    0
....................... 4   2   .   .   7    9    6    7
=========================================================
Nonnegatives . A001477: 0,  1,  2,  3,  4,   5,   6,   7, ...
Triangulars .. A000217: 0,  1,  3,  6, 10,  15,  21,  28, ...
Squares ...... A000290: 0,  1,  4,  9, 16,  25,  36,  49, ...
Pentagonals .. A000326: 0,  1,  5, 12, 22,  35,  51,  70, ...
Hexagonals ... A000384: 0,  1,  6, 15, 28,  45,  66,  91, ...
Heptagonals .. A000566: 0,  1,  7, 18, 34,  55,  81, 112, ...
Octagonals ... A000567: 0,  1,  8, 21, 40,  65,  96, 133, ...
9-gonals ..... A001106: 0,  1,  9, 24, 46,  75, 111, 154, ...
10-gonals .... A001107: 0,  1, 10, 27, 52,  85, 126, 175, ...
11-gonals .... A051682: 0,  1, 11, 30, 58,  95, 141, 196, ...
12-gonals .... A051624: 0,  1, 12, 33, 64, 105, 156, 217, ...
...
=========================================================
The column with the numbers 2, 3, 4, 5, 6, ... is formed by the numbers > 1 of A000027. The column with the numbers 3, 6, 9, 12, 15, ... is formed by the positive members of A008585.
		

Crossrefs

A formal extension negative n is in A326728.
Triangle sums (see the comments): A055795 (Row1), A080956 (Row2; terms doubled), A096338 (Kn11, Kn12, Kn13, Fi1, Ze1), A002624 (Kn21, Kn22, Kn23, Fi2, Ze2), A000332 (Kn3, Ca3, Gi3), A134393 (Kn4), A189374 (Ca1, Ze3), A011779 (Ca2, Ze4), A101357 (Ca4), A189375 (Gi1), A189376 (Gi2), A006484 (Gi4). - Johannes W. Meijer, Apr 29 2011
Sequences of m-gonal numbers: A000217 (m=3), A000290 (m=4), A000326 (m=5), A000384 (m=6), A000566 (m=7), A000567 (m=8), A001106 (m=9), A001107 (m=10), A051682 (m=11), A051624 (m=12), A051865 (m=13), A051866 (m=14), A051867 (m=15), A051868 (m=16), A051869 (m=17), A051870 (m=18), A051871 (m=19), A051872 (m=20), A051873 (m=21), A051874 (m=22), A051875 (m=23), A051876 (m=24), A255184 (m=25), A255185 (m=26), A255186 (m=27), A161935 (m=28), A255187 (m=29), A254474 (m=30).

Programs

  • Magma
    T:= func< n,k | k*(n*(k-1)+2)/2 >;
    A139600:= func< n,k | T(n-k, k) >;
    [A139600(n,k): k in  [0..n], n in [0..12]]; // G. C. Greubel, Jul 12 2024
    
  • Maple
    T:= (n, k)-> n*(k-1)*k/2+k:
    seq(seq(T(d-k, k), k=0..d), d=0..14);  # Alois P. Heinz, Oct 14 2018
  • Mathematica
    T[n_, k_] := (n + 1)*(k - 1)*k/2 + k; Table[T[n - k - 1, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Robert G. Wilson v, Jul 12 2009 *)
  • Python
    def A139600Row(n):
        x, y = 1, 1
        yield 0
        while True:
            yield x
            x, y = x + y + n, y + n
    for n in range(8):
        R = A139600Row(n)
        print([next(R) for  in range(11)]) # _Peter Luschny, Aug 04 2019
    
  • SageMath
    def T(n,k): return k*(n*(k-1)+2)/2
    def A139600(n,k): return T(n-k, k)
    flatten([[A139600(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Jul 12 2024

Formula

T(n,k) = n*(k-1)*k/2+k.
T(n,k) = A057145(n+2,k). - R. J. Mathar, Jul 28 2016
From Stefano Spezia, Apr 12 2024: (Start)
G.f.: y*(1 - x - y + 2*x*y)/((1 - x)^2*(1 - y)^3).
E.g.f.: exp(x+y)*y*(2 + x*y)/2. (End)

Extensions

Edited by Omar E. Pol, Jan 05 2009

A214292 Triangle read by rows: T(n,k) = T(n-1,k-1) + T(n-1,k), 0 < k < n with T(n,0) = n and T(n,n) = -n.

Original entry on oeis.org

0, 1, -1, 2, 0, -2, 3, 2, -2, -3, 4, 5, 0, -5, -4, 5, 9, 5, -5, -9, -5, 6, 14, 14, 0, -14, -14, -6, 7, 20, 28, 14, -14, -28, -20, -7, 8, 27, 48, 42, 0, -42, -48, -27, -8, 9, 35, 75, 90, 42, -42, -90, -75, -35, -9, 10, 44, 110, 165, 132, 0, -132, -165, -110, -44, -10
Offset: 0

Views

Author

Reinhard Zumkeller, Jul 12 2012

Keywords

Examples

			The triangle begins:
    0:                              0
    1:                            1   -1
    2:                          2   0   -2
    3:                       3    2   -2   -3
    4:                     4    5   0   -5   -4
    5:                  5    9    5   -5   -9   -5
    6:                6   14   14   0  -14  -14   -6
    7:             7   20   28   14  -14  -28  -20   -7
    8:           8   27   48   42   0  -42  -48  -27   -8
    9:        9   35   75   90   42  -42  -90  -75  -35   -9
   10:     10   44  110  165  132   0 -132 -165 -110  -44  -10
   11:  11   54  154  275  297  132 -132 -297 -275 -154  -54  -11  .
		

Crossrefs

Programs

  • Haskell
    a214292 n k = a214292_tabl !! n !! k
    a214292_row n = a214292_tabl !! n
    a214292_tabl = map diff $ tail a007318_tabl
       where diff row = zipWith (-) (tail row) row
  • Mathematica
    row[n_] := Table[Binomial[n, k], {k, 0, n}] // Differences;
    T[n_, k_] := row[n + 1][[k + 1]];
    Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jul 31 2018 *)

Formula

T(n,k) = A007318(n+1,k+1) - A007318(n+1,k), 0<=k<=n, i.e. first differences of rows in Pascal's triangle;
T(n,k) = -T(n,k);
row sums and central terms equal 0, cf. A000004;
sum of positive elements of n-th row = A014495(n+1);
T(n,0) = n;
T(n,1) = A000096(n-2) for n > 1; T(n,1) = - A080956(n) for n > 0;
T(n,2) = A005586(n-4) for n > 3; T(n,2) = A129936(n-2);
T(n,3) = A005587(n-6) for n > 5;
T(n,4) = A005557(n-9) for n > 8;
T(n,5) = A064059(n-11) for n > 10;
T(n,6) = A064061(n-13) for n > 12;
T(n,7) = A124087(n) for n > 14;
T(n,8) = A124088(n) for n > 16;
T(2*n+1,n) = T(2*n+2,n) = A000108(n+1), Catalan numbers;
T(2*n+3,n) = A000245(n+2);
T(2*n+4,n) = A002057(n+1);
T(2*n+5,n) = A000344(n+3);
T(2*n+6,n) = A003517(n+3);
T(2*n+7,n) = A000588(n+4);
T(2*n+8,n) = A003518(n+4);
T(2*n+9,n) = A001392(n+5);
T(2*n+10,n) = A003519(n+5);
T(2*n+11,n) = A000589(n+6);
T(2*n+12,n) = A090749(n+6);
T(2*n+13,n) = A000590(n+7).

A118801 Triangle T that satisfies the matrix products: C*[T^-1]*C = T and T*[C^-1]*T = C, where C is Pascal's triangle.

Original entry on oeis.org

1, 1, -1, 1, -3, 1, 1, -7, 5, -1, 1, -15, 17, -7, 1, 1, -31, 49, -31, 9, -1, 1, -63, 129, -111, 49, -11, 1, 1, -127, 321, -351, 209, -71, 13, -1, 1, -255, 769, -1023, 769, -351, 97, -15, 1, 1, -511, 1793, -2815, 2561, -1471, 545, -127, 17, -1, 1, -1023, 4097, -7423, 7937, -5503, 2561, -799, 161, -19, 1
Offset: 0

Views

Author

Paul D. Hanna, May 02 2006

Keywords

Comments

Matrix inverse is triangle A118800. Row sums are: (1-n). Unsigned row sums equal A007051(n) = (3^n + 1)/2. Row squared sums equal A118802. Antidiagonal sums equal A080956(n) = (n+1)(2-n)/2. Unsigned antidiagonal sums form A024537 (with offset).
T = C^2*D^-1 where matrix product D = C^-1*T*C = T^-1*C^2 has only 2 nonzero diagonals: D(n,n)=-D(n+1,n)=(-1)^n, with zeros elsewhere. Also, [B^-1]*T*[B^-1] = B*[T^-1]*B forms a self-inverse matrix, where B^2 = C and B(n,k) = C(n,k)/2^(n-k). - Paul D. Hanna, May 04 2006
Riordan array ( 1/(1 - x), -x/(1 - 2*x) ) The matrix square is the Riordan array ( (1 - 2*x)/(1 - x)^2, x ), which belongs to the Appell subgroup of the Riordan group. See the Example section below. - Peter Bala, Jul 17 2013

Examples

			Formulas for initial columns are, for n>=0:
T(n+1,1) = 1 - 2^(n+1);
T(n+2,2) = 1 + 2^(n+1)*n;
T(n+3,3) = 1 - 2^(n+1)*(n*(n+1)/2 + 1);
T(n+4,4) = 1 + 2^(n+1)*(n*(n+1)*(n+2)/6 + n);
T(n+5,5) = 1 - 2^(n+1)*(n*(n+1)*(n+2)*(n+3)/24 + n*(n+1)/2 + 1).
Triangle begins:
1;
1,-1;
1,-3,1;
1,-7,5,-1;
1,-15,17,-7,1;
1,-31,49,-31,9,-1;
1,-63,129,-111,49,-11,1;
1,-127,321,-351,209,-71,13,-1;
1,-255,769,-1023,769,-351,97,-15,1;
1,-511,1793,-2815,2561,-1471,545,-127,17,-1;
1,-1023,4097,-7423,7937,-5503,2561,-799,161,-19,1; ...
The matrix square, T^2, starts:
1;
0,1;
-1,0,1;
-2,-1,0,1;
-3,-2,-1,0,1;
-4,-3,-2,-1,0,1; ...
where all columns are the same.
The matrix product C^-1*T*C = T^-1*C^2 is:
1;
-1,-1;
0, 1, 1;
0, 0,-1,-1;
0, 0, 0, 1, 1; ...
where C(n,k) = n!/(n-k)!/k!.
		

Crossrefs

Cf. A118800 (inverse), A007051 (unsigned row sums), A118802 (Row squared sums), A080956 (antidiagonal sums), A024537 (unsigned antidiagonal sums).
A145661, A119258 and A118801 are all essentially the same (see the Shattuck and Waldhauser paper). - Tamas Waldhauser, Jul 25 2011

Programs

  • Mathematica
    Table[(1 + (-1)^k*2^(n - k + 1)*Sum[ Binomial[n - 2 j - 2, k - 2 j - 1], {j, 0, Floor[k/2]}]) - 4 Boole[And[n == 1, k == 0]], {n, 0, 10}, {k, 0, n}] // Flatten (* Michael De Vlieger, Nov 24 2016 *)
  • PARI
    {T(n,k)=if(n==0&k==0,1,1+(-1)^k*2^(n-k+1)*sum(j=0,k\2,binomial(n-2*j-2,k-2*j-1)))}

Formula

T(n,k) = 1 + (-1)^k*2^(n-k+1)*Sum_{j=0..[k/2]} C(n-2j-2,k-2j-1) for n>=k>=0 with T(0,0) = 1.
For k>0, T(n,k) = -T(n-1,k-1) + 2*T(n-1,k). - Gerald McGarvey, Aug 05 2006
O.g.f.: (1 - 2*t)/(1 - t) * 1/(1 + t*(x - 2)) = 1 + (1 - x)*t + (1 - 3*x + x^2)*t^2 + (1 - 7*x + 5*x^2 - x^3)*t^3 + .... - Peter Bala, Jul 17 2013
From Tom Copeland, Nov 17 2016: (Start)
Let M = A200139^(-1) = (unsigned A118800)^(-1) and NpdP be the signed padded Pascal matrix defined in A097805. Then T(n,k) = (-1)^n* M(n,k) and T = P*NpdP = (A239473)^(-1)*P^(-1) = P*A167374*P^(-1) = A156644*P^(-1), where P is the Pascal matrix A007318 with inverse A130595. Cf. A112857.
Signed P^2 = signed A032807 = T*A167374. (End)

A156644 Mirror image of triangle A080233.

Original entry on oeis.org

1, 0, 1, -1, 1, 1, -2, 0, 2, 1, -3, -2, 2, 3, 1, -4, -5, 0, 5, 4, 1, -5, -9, -5, 5, 9, 5, 1, -6, -14, -14, 0, 14, 14, 6, 1, -7, -20, -28, -14, 14, 28, 20, 7, 1, -8, -27, -48, -42, 0, 42, 48, 27, 8, 1, -9, -35, -75, -90, -42, 42, 90, 75, 35, 9, 1
Offset: 0

Views

Author

Philippe Deléham, Feb 12 2009

Keywords

Comments

Inverse of A239473. Equals A007318*A167374. - Tom Copeland, Nov 14 2016

Examples

			Triangle begins as:
   1;
   0,  1;
  -1,  1, 1;
  -2,  0, 2, 1;
  -3, -2, 2, 3, 1;
  -4, -5, 0, 5, 4, 1; ...
		

Crossrefs

Programs

  • Magma
    A156644:= func< n,k | ((2*k-n+1)/(k+1))*Binomial(n,k) >;
    [A156644(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Feb 28 2021
  • Mathematica
    Table[Binomial[n, k] -Binomial[n, k+1], {n,0,10}, {k,0,n}]//Flatten (* Michael De Vlieger, Nov 24 2016 *)
  • Sage
    def A156644(n,k): return ((2*k-n+1)/(k+1))*binomial(n,k)
    flatten([[A156644(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Feb 28 2021
    

Formula

T(n,k) = A080233(n,n-k) = (-1)^(n-k)*A097808(n,k).
T(n,k) = ((2*k-n+1)/(k+1))*binomial(n,k).
T(n,k) = T(n-1,k-1) + T(n-1,k), k>0, with T(n,0) = 1-n = A024000(n), T(n,n) = 1.
T(n,k) = binomial(n,k) - binomial(n,k+1) = Sum_{i=-k-1..k+1} (-1)^(i+1) * binomial(n,k+1+i) * binomial(n+2,k+1-i). - Mircea Merca, Apr 28 2012
Sum_{k=0..n} T(n, k) = A000012(n) = 1^n. - G. C. Greubel, Feb 28 2021

A269953 Triangle read by rows: T(n, k) = Sum_{j=0..n} binomial(-j-1, -n-1)*S1(j, k) where S1 are the Stirling cycle numbers A132393.

Original entry on oeis.org

1, -1, 1, 1, -1, 1, -1, 2, 0, 1, 1, 0, 5, 2, 1, -1, 9, 15, 15, 5, 1, 1, 35, 94, 85, 40, 9, 1, -1, 230, 595, 609, 315, 91, 14, 1, 1, 1624, 4458, 4844, 2779, 924, 182, 20, 1, -1, 13209, 37590, 43238, 26817, 9975, 2310, 330, 27, 1
Offset: 0

Views

Author

Peter Luschny, Apr 12 2016

Keywords

Comments

Replacing the Stirling cycle numbers in the definition by the Stirling set numbers leads to A105794.
From Wolfdieter Lang, Jun 19 2017: (Start)
The triangle t(n, k) = (-1)^(n-k)*T(n, k) is the matrix product of P = A007318 (Pascal) and s1 = A048994 (signed Stirling1). This is Sheffer (exp(t), log(1+t)).
The present triangle T is therefore the Sheffer triangle (exp(-t), -log(1-t)). Note that P is Sheffer (exp(t), t) (of the Appell type). (End)
The triangle T(n,k) is a representative of the parametric family of triangles T(m,n,k), whose columns are the coefficients of the standard expansion of the function f(x) = (-log(1-x))^(k)*exp(-m*x)/k! for the case m=1. See A381082. - Igor Victorovich Statsenko, Feb 14 2025

Examples

			Triangle starts:
   1;
  -1,  1;
   1, -1,  1;
  -1,  2,  0,  1;
   1,  0,  5,  2,  1;
  -1,  9, 15, 15,  5,  1;
   1, 35, 94, 85, 40,  9,  1.
		

Crossrefs

Columns k=0..4 give A033999, A002741, A381064, A381065, A381066.
Cf. A000166 (row sums), A080956 (diag n,n-1).
KummerU(-n,1-n-x,z): this sequence (z=-1), A094816 (z=1), |A137346| (z=2), A327997 (z=3).

Programs

  • Maple
    A269953 := (n,k) -> add(binomial(-j-1,-n-1)*abs(Stirling1(j,k)), j=0..n):
    seq(print(seq(A269953(n, k), k=0..n)), n=0..9);
    # Alternative:
    egf := exp(-t)*(1-t)^(-x): ser := series(egf, t, 12): p := n -> coeff(ser, t, n):
    seq(n!*seq(coeff(p(n), x, k), k=0..n), n=0..9); # Peter Luschny, Oct 28 2019
  • Mathematica
    Flatten[Table[Sum[Binomial[-j-1,-n-1] Abs[StirlingS1[j,k]], {j,0,n}], {n,0,9},{k,0,n}]]
    (* Or: *)
    p [n_] := HypergeometricU[-n, 1 - n - x, -1];
    Table[CoefficientList[p[n], x], {n, 0, 9}] (* Peter Luschny, Oct 28 2019 *)

Formula

From Wolfdieter Lang, Jun 19 2017: (Start)
E.g.f. of row polynomials R(n, x) = Sum_{k=0..n} T(n,k)*x^k: exp(-t)/(1 - t)^x.
E.g.f. of column k sequence: exp(-x)*(-log(1-x))^k/k!, k >= 0. (End)
From Peter Bala, Oct 26 2019: (Start)
Let R(n, x) = (-1)^n*Sum_{k >= 0} binomial(n,k)*k!* binomial(-x,k) the n-th row polynomial of this triangle.
R(n, x) = c_n(-x;-1), where c_n(x;a) denotes the n-th Poisson Charlier polynomial.
The series representation e = Sum_{k >= 0} 1/k! is the case n = 0 of the more general result e = n!*Sum_{k >= 0} 1/(k!*R(n,k)*R(n,k+1)), n = 0,2,3,4,.... (End)
R(n, x) = KummerU(-n, 1-n-x, -1). - Peter Luschny, Oct 28 2019

A100218 Riordan array ((1-2*x)/(1-x), (1-x)).

Original entry on oeis.org

1, -1, 1, -1, -2, 1, -1, 0, -3, 1, -1, 0, 2, -4, 1, -1, 0, 0, 5, -5, 1, -1, 0, 0, -2, 9, -6, 1, -1, 0, 0, 0, -7, 14, -7, 1, -1, 0, 0, 0, 2, -16, 20, -8, 1, -1, 0, 0, 0, 0, 9, -30, 27, -9, 1, -1, 0, 0, 0, 0, -2, 25, -50, 35, -10, 1, -1, 0, 0, 0, 0, 0, -11, 55, -77, 44, -11, 1, -1, 0, 0, 0, 0, 0, 2, -36, 105, -112, 54, -12, 1
Offset: 0

Views

Author

Paul Barry, Nov 08 2004

Keywords

Examples

			Triangle begins as:
   1;
  -1,  1;
  -1, -2,  1;
  -1,  0, -3,  1;
  -1,  0,  2, -4,  1;
  -1,  0,  0,  5, -5,   1;
  -1,  0,  0, -2,  9,  -6,   1;
  -1,  0,  0,  0, -7,  14,  -7,  1;
  -1,  0,  0,  0,  2, -16,  20, -8,  1;
  -1,  0,  0,  0,  0,   9, -30, 27, -9,  1;
		

Crossrefs

Row sums are A100219.
Matrix inverse of A100100.
Apart from signs, same as A098599.
Very similar to triangle A111125.

Programs

  • Magma
    A100218:= func< n,k | n eq 0 select 1 else (-1)^(n+k)*(Binomial(k,n-k) + Binomial(k-1,n-k-1)) >;
    [A100218(n,k): k in [0..n], n in [0..13]]; // G. C. Greubel, Mar 28 2024
    
  • Mathematica
    T[0,0]:= 1; T[1,1]:= 1; T[1,0]:= -1; T[n_, k_]:= T[n, k]= If[k<0 || k>n, 0, T[n- 1,k] +T[n-1,k-1] -2*T[n-2,k-1] +T[n-3,k-1]]; Table[T[n,k], {n,0,14}, {k,0,n} ]//Flatten (* G. C. Greubel, Mar 13 2017 *)
  • SageMath
    def A100218(n,k): return 1 if n==0 else (-1)^(n+k)*(binomial(k,n-k) + binomial(k-1,n-k-1))
    flatten([[A100218(n,k) for k in range(n+1)] for n in range(14)]) # G. C. Greubel, Mar 28 2024

Formula

Sum_{k=0..n} T(n, k) = A100219(n) (row sums).
Number triangle T(n, k) = (-1)^(n-k)*(binomial(k, n-k) + binomial(k-1, n-k-1)), with T(0, 0) = 1. - Paul Barry, Nov 09 2004
T(n,k) = T(n-1,k) + T(n-1,k-1) - 2*T(n-2,k-1) + T(n-3,k-1), T(0,0)=1, T(1,0)=-1, T(1,1)=1, T(n,k)=0 if k < 0 or if k > n. - Philippe Deléham, Jan 09 2014
From G. C. Greubel, Mar 28 2024: (Start)
T(n, n-1) = A000027(n), n >= 1.
T(n, n-2) = -A080956(n-1), n >= 2.
T(2*n, n) = A280560(n).
T(2*n-1, n) = A157142(n-1), n >= 1.
T(2*n+1, n) = -A000007(n) = A154955(n+2).
T(3*n, n) = T(4*n, n) = A000007(n).
Sum_{k=0..n} (-1)^k*T(n, k) = A355021(n).
Sum_{k=0..floor(n/2)} T(n-k, k) = (-1)^n*A098601(n).
Sum_{k=0..floor(n/2)} (-1)^k*T(n-k, k) = -1 + 2*A077961(n) + A077961(n-2). (End)
From Peter Bala, Apr 28 2024: (Start)
This Riordan array has the form ( x*h'(x)/h(x), h(x) ) with h(x) = x*(1 - x) and hence belongs to the hitting time subgroup of the Riordan group (see Peart and Woan for properties of this subgroup).
T(n,k) = [x^(n-k)] (1/c(x))^n, where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the g.f. of the Catalan numbers A000108. In general the (n, k)-th entry of the hitting time array ( x*h'(x)/h(x), h(x) ) has the form [x^(n-k)] f(x)^n, where f(x) = x/( series reversion of h(x) ). (End)

A262221 a(n) = 25*n*(n + 1)/2 + 1.

Original entry on oeis.org

1, 26, 76, 151, 251, 376, 526, 701, 901, 1126, 1376, 1651, 1951, 2276, 2626, 3001, 3401, 3826, 4276, 4751, 5251, 5776, 6326, 6901, 7501, 8126, 8776, 9451, 10151, 10876, 11626, 12401, 13201, 14026, 14876, 15751, 16651, 17576, 18526, 19501, 20501, 21526, 22576, 23651
Offset: 0

Views

Author

Bruno Berselli, Sep 15 2015

Keywords

Comments

Also centered 25-gonal (or icosipentagonal) numbers.
This is the case k=25 of the formula (k*n*(n+1) - (-1)^k + 1)/2. See table in Links section for similar sequences.
For k=2*n, the formula shown above gives A011379.
Primes in sequence: 151, 251, 701, 1951, 3001, 4751, 10151, 12401, ...

References

  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 51 (23rd row of the table).

Crossrefs

Cf. centered polygonal numbers listed in A069190.
Similar sequences of the form (k*n*(n+1) - (-1)^k + 1)/2 with -1 <= k <= 26: A000004, A000124, A002378, A005448, A005891, A028896, A033996, A035008, A046092, A049598, A060544, A064200, A069099, A069125, A069126, A069128, A069130, A069132, A069174, A069178, A080956, A124080, A163756, A163758, A163761, A164136, A173307.

Programs

  • Magma
    [25*n*(n+1)/2+1: n in [0..50]];
  • Mathematica
    Table[25 n (n + 1)/2 + 1, {n, 0, 50}]
    25*Accumulate[Range[0,50]]+1 (* or *) LinearRecurrence[{3,-3,1},{1,26,76},50] (* Harvey P. Dale, Jan 29 2023 *)
  • PARI
    vector(50, n, n--; 25*n*(n+1)/2+1)
    
  • Sage
    [25*n*(n+1)/2+1 for n in (0..50)]
    

Formula

G.f.: (1 + 23*x + x^2)/(1 - x)^3.
a(n) = a(-n-1) = 3*a(n-1) - 3*a(n-2) + a(n-3).
a(n) = A123296(n) + 1.
a(n) = A000217(5*n+2) - 2.
a(n) = A034856(5*n+1).
a(n) = A186349(10*n+1).
a(n) = A054254(5*n+2) with n>0, a(0)=1.
a(n) = A000217(n+1) + 23*A000217(n) + A000217(n-1) with A000217(-1)=0.
Sum_{i>=0} 1/a(i) = 1.078209111... = 2*Pi*tan(Pi*sqrt(17)/10)/(5*sqrt(17)).
From Amiram Eldar, Jun 21 2020: (Start)
Sum_{n>=0} a(n)/n! = 77*e/2.
Sum_{n>=0} (-1)^(n+1) * a(n)/n! = 23/(2*e). (End)
E.g.f.: exp(x)*(2 + 50*x + 25*x^2)/2. - Elmo R. Oliveira, Dec 24 2024

A176145 a(n) = n*(n-3)*(n^2-7*n+14)/8.

Original entry on oeis.org

0, 1, 5, 18, 49, 110, 216, 385, 638, 999, 1495, 2156, 3015, 4108, 5474, 7155, 9196, 11645, 14553, 17974, 21965, 26586, 31900, 37973, 44874, 52675, 61451, 71280, 82243, 94424, 107910, 122791, 139160, 157113, 176749, 198170, 221481, 246790, 274208, 303849
Offset: 3

Views

Author

Michel Lagneau, Apr 10 2010

Keywords

Comments

Number of points of intersection of diagonals of a general convex n-polygon. (both inside and outside the polygon).
n(n-3)/2 (n >= 3) is the number of diagonals of an n-gon (A080956). The number of points (inside or outside), distinct of tops, where these diagonals intersect is : (1/2)( n(n-3)/2)(n(n-3)/2 - 2(n-4) -1) = n(n-3)(n^2 - 7n + 14) / 8.

Examples

			For n=3, a(3) = 0 (no diagonals in triangle),
For n=4, a(4) = 1 (2 diagonals => 1 point of intersection),
For n=5, a(5) = 5 (5 diagonals => 5 points of intersection),
For n=6, a(6) = 18 (9 diagonals => 18 points of intersection).
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 797.

Crossrefs

Programs

  • Magma
    [n*(n-3)*(n^2 - 7*n + 14) / 8: n in [3..60]]; // Vincenzo Librandi, May 21 2011
    
  • Maple
    for n from 3 to 50 do: x:=n*(n-3)*(n^2 - 7*n +14)/8 : print(x):od:
  • Mathematica
    Table[n*(n - 3)*(n^2 - 7*n + 14)/8, {n, 3, 42}] (* Bobby Milazzo, Jun 24 2013 *)
    Drop[CoefficientList[Series[x^4(1+3x^2-x^3)/(1-x)^5,{x,0,50}],x],3] (* or *) LinearRecurrence[{5,-10,10,-5,1},{0,1,5,18,49},50] (* Harvey P. Dale, Mar 14 2022 *)
  • PARI
    vector(100,n,(n+2)*(n-1)*(n^2-3*n+4)/8) \\ Derek Orr, Jan 21 2015

Formula

G.f.: x^4*(1+3*x^2-x^3)/(1-x)^5. [Colin Barker, Jan 31 2012]
a(n) = 5*a(n-1) -10*a(n-2) +10*a(n-3) -5*a(n-4) + a(n-5), with a(3)= 0, a(4)= 1, a(5)=5, a(6)= 18, a(7) = 49. [Bobby Milazzo, Jun 24 2013]
a(n) = Sum_{k=(n-3)..(n-2)*(n-3)/2} k. - J. M. Bergot, Jan 21 2015

Extensions

Edited by N. J. A. Sloane, Apr 19 2010
Showing 1-10 of 35 results. Next