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.

Previous Showing 11-20 of 205 results. Next

A123125 Triangle of Eulerian numbers T(n,k), 0 <= k <= n, read by rows.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 4, 1, 0, 1, 11, 11, 1, 0, 1, 26, 66, 26, 1, 0, 1, 57, 302, 302, 57, 1, 0, 1, 120, 1191, 2416, 1191, 120, 1, 0, 1, 247, 4293, 15619, 15619, 4293, 247, 1, 0, 1, 502, 14608, 88234, 156190, 88234, 14608, 502, 1, 0, 1, 1013, 47840, 455192, 1310354, 1310354, 455192, 47840, 1013, 1
Offset: 0

Views

Author

Philippe Deléham, Sep 30 2006

Keywords

Comments

The beginning of this sequence does not quite agree with the usual version, which is A173018. - N. J. A. Sloane, Nov 21 2010
Each row of A123125 is the reverse of the corresponding row in A173018. - Michael Somos, Mar 17 2011
A008292 (subtriangle for k>=1 and n>=1) is the main entry for these numbers.
Triangle T(n,k), 0 <= k <= n, read by rows given by [0,1,0,2,0,3,0,4,0,5,0,...] DELTA [1,0,2,0,3,0,4,0,5,0,6,...] where DELTA is the operator defined in A084938.
Row sums are the factorials. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008
If the initial zero column is deleted, the result is A008292. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008
This result gives an alternative method of calculating the Eulerian numbers by an Umbral Calculus expansion from Comtet. - Roger L. Bagula, Nov 21 2009
This function seems to be equivalent to the PolyLog expansion. - Roger L. Bagula, Nov 21 2009
A raising operator formed from the e.g.f. of this entry is the generator of a sequence of polynomials p(n,x;t) defined in A046802 that specialize to those for A119879 as p(n,x;-1), A007318 as p(n,x;0), A073107 as p(n,x;1), and A046802 as p(n,0;t). See Copeland link for more associations. - Tom Copeland, Oct 20 2015
The Eulerian numbers in this setup count the permutation trees of power n and width k (see the Luschny link). For the associated combinatorial statistic over permutations see the Sage program below and the example section. - Peter Luschny, Dec 09 2015 [See Elder et al. link. Peter Luschny, Jul 13 2022]
From Wolfdieter Lang, Apr 03 2017: (Start)
The row polynomials R(n, x) = Sum_{k=0..n} T(n, k)*x^k are the numerator polynomials of the o.g.f. G(n, x) of n-powers {m^n}_{m>=0} (with 0^0 = 1): G(n, x) = R(n, x)/(1-x)^(n+1). See the Aug 14 2008 formula, where f(x,n) = R(n, x). The e.g.f. of R(n, t) is given in Copeland's Oct 14 2015 formula below.
The first nine column sequences are A000007, A000012, A000295, A000460, A000498, A000505, A000514, A001243, A001244. (End)
With all offsets 0, let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of this entry, A123125. Then the row polynomials of A046802 (the h-polynomials of the stellahedra) are given by h_n(x) = A_n(x;1); the row polynomials of A248727 (the face polynomials of the stellahedra), by f_n(x) = A_n(1 + x;1); the Swiss-knife polynomials of A119879, by Sw_n(x) = A_n(-1;1 + x); and the row polynomials of the Worpitsky triangle (A130850), by w_n(x) = A(1 + x;0). Other specializations of A_n(x;y) give A090582 (the f-polynomials of the permutohedra, cf. also A019538) and A028246 (another version of the Worpitsky triangle). - Tom Copeland, Jan 24 2020
Let b(n) = (1/(n+1))*Sum_{k=0..n-1} (-1)^(n-k+1)*T(n, k+1) / binomial(n, k+1). Then b(n) = Bernoulli(n, 1) = -n*Zeta(1 - n) = Integral_{x=0..1} F_n(x) for n >= 1. Here F_n(x) are the signed Fubini polynomials (A278075). (See also Rzadkowski and Urlinska, example 1.) - Peter Luschny, Feb 15 2021
Patrick J. Burchell (see link) describes the following method: To get the k-th row of the triangle write the nonnegative integers with a fixed exponent k as a sequence, 0^k, 1^k, 2^k, ..., and then apply the first differences to them k + 1 times. - Peter Luschny, Apr 02 2023

Examples

			The triangle T(n, k) begins:
  n\k 0 1    2     3      4       5       6      7     8    9 10...
  0:  1
  1:  0 1
  2:  0 1    1
  3:  0 1    4     1
  4:  0 1   11    11      1
  5:  0 1   26    66     26       1
  6:  0 1   57   302    302      57       1
  7:  0 1  120  1191   2416    1191     120      1
  8:  0 1  247  4293  15619   15619    4293    247     1
  9:  0 1  502 14608  88234  156190   88234  14608   502    1
 10:  0 1 1013 47840 455192 1310354 1310354 455192 47840 1013  1
...  Reformatted. - _Wolfdieter Lang_, Feb 14 2015
------------------------------------------------------------------
The width statistic over permutations, n=4.
  [1, 2, 3, 4] => 3; [1, 2, 4, 3] => 2; [1, 3, 2, 4] => 2; [1, 3, 4, 2] => 2;
  [1, 4, 2, 3] => 2; [1, 4, 3, 2] => 1; [2, 1, 3, 4] => 3; [2, 1, 4, 3] => 2;
  [2, 3, 1, 4] => 2; [2, 3, 4, 1] => 3; [2, 4, 1, 3] => 2; [2, 4, 3, 1] => 2;
  [3, 1, 2, 4] => 3; [3, 1, 4, 2] => 3; [3, 2, 1, 4] => 2; [3, 2, 4, 1] => 3;
  [3, 4, 1, 2] => 3; [3, 4, 2, 1] => 2; [4, 1, 2, 3] => 4; [4, 1, 3, 2] => 3;
  [4, 2, 1, 3] => 3; [4, 2, 3, 1] => 3; [4, 3, 1, 2] => 3; [4, 3, 2, 1] => 2;
Gives row(4) = [0, 1, 11, 11, 1]. - _Peter Luschny_, Dec 09 2015
------------------------------------------------------------------
From _Wolfdieter Lang_, Apr 03 2017: (Start)
Recurrence: T(5, 3) = (6-3)*T(4, 2) + 3*T(4, 3) = 3*11 + 3*11 = 66.
O.g.f. column k=2: (x/(1 - 2*x))*E_x*(x/(1-x)) = (x/(1-x))^2/(1-2*x).
E.g.f. column k=2: A(2, x) = x*A(1, x) + x*E(1, x) = x*1 + x*(exp(x)-1) = x*exp(x), hence E(2, x) = (1 + int(x*exp(-x),x ))*exp(2*x) = exp(x)*(exp(x) - (1+x)). See A000295. (End)
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, Holland, 1978, page 245. [Roger L. Bagula, Nov 21 2009]
  • Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics, 2nd ed.; Addison-Wesley, 1994, p. 268, Row reversed table 268. - Wolfdieter Lang, Apr 03 2017
  • Douglas C. Montgomery and Lynwood A. Johnson, Forecasting and Time Series Analysis, MaGraw-Hill, New York, 1976, page 91. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008

Crossrefs

See A008292 (subtriangle for k>=1 and n>=1), which is the main entry for these numbers. Another version has the zeros at the ends of the rows, as in Concrete Mathematics: see A173018.
T(2n,n) gives A180056.

Programs

  • Haskell
    a123125 n k = a123125_tabl !! n !! k
    a123125_row n = a123125_tabl !! n
    a123125_tabl = [1] : zipWith (:) [0, 0 ..] a008292_tabl
    -- Reinhard Zumkeller, Nov 06 2013
    
  • Maple
    gf := 1/(1 - t*exp(x)): ser := series(gf, x, 12):
    cx := n -> (-1)^(n + 1)*factor(n!*coeff(ser, x, n)*(t - 1)^(n + 1)):
    seq(print(seq(coeff(cx(n), t, k), k = 0..n)), n = 0..9); # Peter Luschny, Feb 11 2021
    A123125 := proc(n, k) option remember; if k = n then 1 elif k <= 0 or k > n then 0 else k*procname(n-1, k) + (n-k+1)*procname(n-1, k-1) fi end:
    seq(print(seq(A123125(n, k), k=0..n)), n=0..10); # Peter Luschny, Mar 28 2021
    # Alternative (Patrick J. Burchell):
    t := a -> Statistics:-Difference([0, a]): Trow := k -> (t@@(k+1))([seq(n^k, n = 0..k)]):
    seq(print(Trow(n)), n = 0..6); # Peter Luschny, Apr 02 2023
  • Mathematica
    f[x_, n_] := f[x, n] = (1 - x)^(n + 1)*Sum[k^n*x^k, {k, 0, Infinity}];
    Table[CoefficientList[f[x, n], x], {n,0,9}] // Flatten (* Roger L. Bagula, Aug 14 2008 *)
    t[n_ /; n >= 0, 0] = 1; t[n_, k_] /; k<0 || k>n = 0; t[n_, k_] := t[n, k] = (n-k) t[n-1, k-1] + (k+1) t[n-1, k]; T[n_, k_] := t[n, n-k];
    Table[T[n, k], {n,0,10}, {k, 0, n}] // Flatten (* Jean-François Alcover, May 26 2019 *)
    A123125[n_, k_] := Sum[(-1)^j*(n - j - k + 1)^n * Binomial[n + 1, j], {j, 0, n - k}];
    Table[A123125[n, k], {n, 0, 9}, {k, 0, n}] // TableForm  (* Peter Luschny, Aug 12 2022 *)
  • Python
    from math import isqrt, comb
    def A123125(n):
        a = (m:=isqrt(k:=n+1<<1))+(k>m*(m+1))
        b = comb(a+1,2)-n
        return sum(-(b-j)**(a-1)*comb(a,j) if j&1 else (b-j)**(a-1)*comb(a,j) for j in range(b)) # Chai Wah Wu, Nov 13 2024
  • Sage
    def statistic_eulerian(pi):
        if not pi: return 0
        h, i, branch, next = 0, len(pi), [0], pi[0]
        while True:
            while next < branch[len(branch)-1]:
                del(branch[len(branch)-1])
            current = 0
            h += 1
            while next > current:
                i -= 1
                if i == 0: return h
                branch.append(next)
                current, next = next, pi[i]
    def A123125_row(n):
        L = [0]*(n+1)
        for p in Permutations(n):
            L[statistic_eulerian(p)] += 1
        return L
    [A123125_row(n) for n in range(7)] # Peter Luschny, Dec 09 2015
    

Formula

Sum_{k=0..n} T(n,k) = n! = A000142(n).
Sum_{k=0..n} 2^k*T(n,k) = A000629(n).
Sum_{k=0..n} 3^k*T(n,k) = abs(A009362(n+1)).
Sum_{k=0..n} 2^(n-k)*T(n,k) = A000670(n).
Sum_{k=0..n} T(n,k)*3^(n-k) = A122704(n). - Philippe Deléham, Nov 07 2007
G.f.: f(x,n) = (1 - x)^(n + 1)*Sum_{k>=0} k^n*x^k. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008. f is not the g.f. of the triangle, it is the polynomial of row n. See an Apr 03 2017 comment above - Wolfdieter Lang, Apr 03 2017
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A000142(n), A000629(n), A123227(n), A201355(n), A201368(n) for x = 0, 1, 2, 3, 4, 5 respectively. - Philippe Deléham, Dec 01 2011
E.g.f. (1-t)/(1-t*exp((1-t)x)). A123125 * A007318 = A130850 = unsigned A075263, related to reversed A028246. A007318 * A123125 = A046802. Evaluating the row polynomials at -1, giving the alternating-sign row sum, generates A009006. - Tom Copeland, Oct 14 2015
From Wolfdieter Lang, Apr 03 2017: (Start)
T(n, k) = A173018(n, n-k), 0 <= k <= n. Row reversed Euler's triangle. See Graham et al., p. 268.
Recurrence (from A173018): T(n, 0) = 1 if n=0 else 0; T(n, k) = 0 if n < k and T(n, k) = (n+1-k)*T(n-1, k-1) + k*T(n-1, k) else.
T(n, k) = Sum_{j=0..k} (-1)^(k-j)*binomial(n-j, k-j)*S2(n, j)*j!, 0 <= k <= n, else 0. For S2(n, k)*k! see A131689.
The recurrence for the o.g.f. of the sequence of column k is
G(k, x) = (x/(1 - k*x))*(E_x - (k-2))*G(k-1, x), with the Euler operator E_x = x*d_x, for k >= 1, with G(0, x) = 1. (Proof from the recurrence of T(n, k)).
The e.g.f of the sequence of column k is found from E(k, x) = (1 + int(A(k, x),x)*exp(-k*x))*exp(k*x), k >= 1, with the recurrence
A(k, x) = x*A(k-1, x) +(1 + (1-k)*(1-x))*E(k-1, x) for k >= 1, with A(0,x)= 0. (Proof from the recurrence of T(n, k)). (End)
T(n, k) = Sum_{j=0..n-k} (-1)^j*(n-j-k+1)^n*binomial(n + 1, j). - Peter Luschny, Aug 12 2022
G.f.: Sum_{m >= 0} x^m/(1/(1-x)-m*t). - Mamuka Jibladze, Mar 12 2025

A065941 T(n,k) = binomial(n-floor((k+1)/2), floor(k/2)). Triangle read by rows, for 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 1, 4, 3, 3, 1, 1, 1, 5, 4, 6, 3, 1, 1, 1, 6, 5, 10, 6, 4, 1, 1, 1, 7, 6, 15, 10, 10, 4, 1, 1, 1, 8, 7, 21, 15, 20, 10, 5, 1, 1, 1, 9, 8, 28, 21, 35, 20, 15, 5, 1, 1, 1, 10, 9, 36, 28, 56, 35, 35, 15, 6, 1, 1, 1, 11, 10, 45, 36, 84, 56, 70, 35, 21, 6, 1
Offset: 0

Views

Author

Len Smiley, Nov 29 2001

Keywords

Comments

Also the q-Stirling2 numbers at q = -1. - Peter Luschny, Mar 09 2020
Row sums give the Fibonacci sequence. So do the alternating row sums.
Triangle of coefficients of polynomials defined by p(-1,x) = p(0,x) = 1, p(n, x) = x*p(n-1, x) + p(n-2, x), for n >= 1. - Benoit Cloitre, May 08 2005 [rewritten with correct offset. - Wolfdieter Lang, Feb 18 2020]
Another version of triangle in A103631. - Philippe Deléham, Jan 01 2009
The T(n,k) coefficients appear in appendix 2 of Parks's remarkable article "A new proof of the Routh-Hurwitz stability criterion using the second method of Liapunov" if we assume that the b(n) coefficients are all equal to 1 and ignore the first column. The complete version of this triangle including the first column is A103631. - Johannes W. Meijer, Aug 11 2011
Signed ++--++..., the roots are chaotic using f(x) --> x^2 - 2 with cycle lengths shown in A003558 by n-th rows. Example: given row 3, x^3 + x^2 - 2x - 1; the roots are (a = 1.24697, ...; b = -0.445041, ...; c = -1.802937, ...). Then (say using seed b with x^2 - 2) we obtain the trajectory -0.445041, ... -> -1.80193, ... -> 1.24697, ...; matching the entry "3" in A003558(3). - Gary W. Adamson, Sep 06 2011
From Gary W. Adamson, Aug 25 2019: (Start)
Roots to the polynomials and terms in A003558 can all be obtained from the numbers below using a doubling series mod N procedure as follows: (more than one row may result). Any row ends when the trajectory produces a term already used. Then try the next higher odd term not used as the leftmost term, then repeat.
For example, for N = 11, we get: (1, 2, 4, 3, 5), showing that when confronted with two choices after the 4: (8 and -3), pick the smaller (abs) term, = 3. Then for the next row pick 7 (not used) and repeat the algorithm; succeeding only if the trajectory produces new terms. But 7 is also (-4) mod 11 and 4 was used. Therefore what I call the "r-t table" (for roots trajectory) has only one row: (1, 2, 4, 3, 5). Conjecture: The numbers of terms in the first row is equal to A003558 corresponding to N, i.e., 5 in this case with period 2.
Now for the roots to the polynomials. Pick N = 7. The polynomial is x^3 - x^2 - 2x + 1 = 0, with roots 1.8019..., -1.2469... and 0.445... corresponding to 2*cos(j*Pi/N), N = 7, and j = (1, 2, and 3). The terms (1, 2, 3) are the r-t terms for N = 7. For 11, the r-t terms are (1, 2, 4, 3, 5). This implies that given any roots of the corresponding polynomial, they are cyclic using f(x) --> x^2 - 2 with cycle lengths shown in A003558. The terms thus generated are 2*cos(j*Pi), with j = (1, 2, 4, 3, 5). Check: Begin with 2*j*Pi/N, with j = 1 (1.9189...). The other trajectory terms are: --> 1.6825..., --> 0.83083..., -1.3097...; 545...; (a 5 period and cyclic since we can begin with any of the constants). The r-t table for odd N begins as follows:
3...............1
5...............1, 2
7...............1, 2, 3
9...............1, 2, 4
...............3 (singleton terms reduce to "1") (9 has two rows)
11...............1, 2, 4, 3, 5
13...............1, 2, 4, 5, 3, 6
15...............1, 2, 4, 7
................3, 6 (dividing through by the gcd gives (1, 2))
................5. (singleton terms reduce to "1")
The result is that 15 has 3 factors (since 3 rows), and the values of those factors are the previous terms "N", corresponding to the r-t terms in each row. Thus, the first row is new, the second (1, 2), corresponds to N = 5, and the "1" in row 3 corresponds to N = 3. The factors are those values apart from 15 and 1. Note that all of the unreduced r-t terms in all rows for N form a complete set of the terms 1 through (N-1)/2 without duplication. (End)
From Gary W. Adamson, Sep 30 2019: (Start)
The 3 factors of the 7th degree polynomial for 15: (x^7 - x^6 - 6x^5 + 5x^4 + 10x^3 - 6x^2 - 4x + 1) can be determined by getting the roots for 2*cos(j*Pi/1), j = (1, 2, 4, 7) and finding the corresponding polynomial, which is x^4 + x^3 - 4x^2 - 4x + 1. This is the minimal polynomial for N = 15 as shown in Table 2, p. 46 of (Lang). The degree of this polynomial is 4, corresponding to the entry in A003558 for 15, = 4. The trajectories (3, 6) and (5) are j values for 2*cos(j*Pi/15) which are roots to x^2 - x - 1 (relating to the pentagon), and (x - 1), relating to the triangle. (End)
From Gary W. Adamson, Aug 21 2019: (Start)
Matrices M of the form: (1's in the main diagonal, -1's in the subdiagonal, and the rest zeros) are chaotic if we replace (f(x) --> x^2 - 2) with f(x) --> M^2 - 2I, where I is the Identity matrix [1, 0, 0; 0, 1, 0; 0, 0, 1]. For example, with the 3 X 3 matrix M: [0, 0, 1; 0, 1, -1; 1, -1, 0]; the f(x) trajectory is:
....M^2 - 2I: [-1, -1, 0; -1, 0, -1; 0, -1, 0], then for the latter,
....M^2 - 2I: [0, 1, 1; 1, 0, 0; 1, 0, -1]. The cycle ends with period 3 since the next matrix is (-1) * the seed matrix. As in the case with f(x) --> x^2 - 2, the eigenvalues of the 3 chaotic matrices are (abs) 1.24697, 0.44504... and 1.80193, ... Also, the characteristic equations of the 3 matrices are the same as or variants of row 4 of the triangle below: (x^3 + x - 2x - 1) with different signs. (End)
Received from Herb Conn, Jan 2004: (Start)
Let x = 2*cos(2A) (A = Angle); then
sin(A)/sin A = 1
sin(3A)/sin A = x + 1
sin(5A)/sin A = x^2 + x - 1
sin(7A)/sin A = x^3 + x - 2x - 1
sin(9A)/sin A = x^4 + x^3 - 3x^2 - 2x + 1
... (signed ++--++...). (End)
Or Pascal's triangle (A007318) with duplicated diagonals. Also triangle of coefficients of polynomials defined by P_0(x) = 1 and for n>=1, P_n(x) = F_n(x) + F_(n+1)(x), where F_n(x) is Fibonacci polynomial (cf. A049310): F_n(x) = Sum_{i=0..floor((n-1)/2)} C(n-i-1,i)*x^(n-2*i-1). - Vladimir Shevelev, Apr 12 2012
The matrix inverse is given by
1;
1, 1;
0, -1, 1;
0, 1, -2, 1;
0, 0, 1, -2, 1;
0, 0, -1, 3, -3, 1;
0, 0, 0, -1, 3, -3, 1;
0, 0, 0, 1, -4, 6, -4, 1;
0, 0, 0, 0, 1, -4, 6, -4, 1;
... apart from signs the same as A124645. - R. J. Mathar, Mar 12 2013

Examples

			Triangle T(n, k) begins:
n\k 0  1  2  3   4   5  6   7  8  9 ...
---------------------------------------
[0] 1,
[1] 1, 1,
[2] 1, 1, 1,
[3] 1, 1, 2, 1,
[4] 1, 1, 3, 2,  1,
[5] 1, 1, 4, 3,  3,  1,
[6] 1, 1, 5, 4,  6,  3,  1,
[7] 1, 1, 6, 5, 10,  6,  4,  1,
[8] 1, 1, 7, 6, 15, 10, 10,  4,  1,
[9] 1, 1, 8, 7, 21, 15, 20, 10,  5, 1,
---------------------------------------
From _Gary W. Adamson_, Oct 23 2019: (Start)
Consider the roots of the polynomials corresponding to odd N such that for N=7 the polynomial is (x^3 + x^2 - 2x - 1) and the roots (a, b, c) are (-1.8019377..., 1.247697..., and -0.445041...). The discriminant of a polynomial derived from the roots is the square of the product of successive differences: ((a-b), (b-c), (c-a))^2 in this case, resulting in 49, matching the method derived from the coefficients of a cubic. For our purposes we use the product of the differences, not the square, resulting in (3.048...) * (1.69202...) * (1.35689...) = 7.0. Conjecture: for all polynomials in the set, the product of the differences of the roots = the corresponding N. For N = 7, we get x^3 - 7x + 7. It appears that for all prime N's, these resulting companion polynomials are monic (left coefficient is 1), and all other coefficients are N or multiples thereof, with the rightmost term = N. The companion polynomials for the first few primes are:
  N =  5:  x^2 - 5;
  N =  7:  x^3 - 7x + 7;
  N = 11:  x^5 - 11x^3 + 11x^2 + 11x - 11;
  N = 13:  x^6 - 13x^4 + 13x^3 + 26x^2 - 39x + 13;
  N = 17:  x^8 - 17x^6 + 17x^5 + 68x^4 - 119x^3 + 17x^2 + 51x - 17;
  N = 19:  x^9 - 19x^7 + 19x^6 + 95x^5 - 171x^4 - 19x^3 + 190x^2 - 114x + 19. (End)
		

Crossrefs

Cf. A065942 (central stalk sequence), A000045 (row sums), A108299.
Reflected version of A046854.
Some triangle sums (see A180662): A000045 (Fi1), A016116 (Kn21), A000295 (Kn23), A094967 (Fi2), A000931 (Ca2), A001519 (Gi3), A000930 (Ze3).

Programs

  • Haskell
    a065941 n k = a065941_tabl !! n !! k
    a065941_row n = a065941_tabl !! n
    a065941_tabl = iterate (\row ->
       zipWith (+) ([0] ++ row) (zipWith (*) (row ++ [0]) a059841_list)) [1]
    -- Reinhard Zumkeller, May 07 2012
    
  • Magma
    [Binomial(n - Floor((k+1)/2), Floor(k/2)): k in [0..n], n in [0..15]]; // G. C. Greubel, Jul 10 2019
    
  • Maple
    A065941 := proc(n,k): binomial(n-floor((k+1)/2),floor(k/2)) end: seq(seq(A065941(n,k), k=0..n), n=0..15); # Johannes W. Meijer, Aug 11 2011
    A065941 := proc(n,k) option remember: local j: if k=0 then 1 elif k=1 then 1: elif k>=2 then add(procname(j,k-2), j=k-2..n-2) fi: end: seq(seq(A065941(n,k), k=0..n), n=0..15);  # Johannes W. Meijer, Aug 11 2011
    # The function qStirling2 is defined in A333143.
    seq(print(seq(qStirling2(n, k, -1), k=0..n)), n=0..9);
    # Peter Luschny, Mar 09 2020
  • Mathematica
    Flatten[Table[Binomial[n-Floor[(k+1)/2],Floor[k/2]],{n,0,15},{k,0,n}]] (* Harvey P. Dale, Dec 11 2011 *)
  • PARI
    T065941(n, k) = binomial(n-(k+1)\2, k\2); \\ Michel Marcus, Apr 28 2014
    
  • Sage
    [[binomial(n - floor((k+1)/2), floor(k/2)) for k in (0..n)] for n in (0..15)] # G. C. Greubel, Jul 10 2019

Formula

T(n, k) = binomial(n-floor((k+1)/2), floor(k/2)).
As a square array read by antidiagonals, this is given by T1(n, k) = binomial(floor(n/2) + k, k). - Paul Barry, Mar 11 2003
Triangle is a reflection of that in A066170 (absolute values). - Gary W. Adamson, Feb 16 2004
Recurrences: T(k, 0) = 1, T(k, n) = T(k-1, n) + T(k-2, n-2), or T(k, n) = T(k-1, n) + T(k-1, n-1) if n even, T(k-1, n-1) if n odd. - Ralf Stephan, May 17 2004
G.f.: sum[n, sum[k, T(k, n)x^ky^n]] = (1+xy)/(1-y-x^2y^2). sum[n>=0, T(k, n)y^n] = y^k/(1-y)^[k/2]. - Ralf Stephan, May 17 2004
T(n, k) = A108299(n, k)*A087960(k) = abs(A108299(n, k)). - Reinhard Zumkeller, Jun 01 2005
From Johannes W. Meijer, Aug 11 2011: (Start)
T(n,k) = A046854(n, n-k) = abs(A066170(n, n-k)).
T(n+k, n-k) = A109223(n,k).
T(n, k) = sum(T(j, k-2), j=k-2..n-2), 2 <= k <= n, n>=2;
T(n, 0) =1, T(n+1, 1) = 1, n >= 0. (End)
For n > 1: T(n, k) = T(n-2, k) + T(n-1, k), 1 < k < n. - Reinhard Zumkeller, Apr 24 2013

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

Original entry on oeis.org

0, 1, 5, 17, 49, 129, 321, 769, 1793, 4097, 9217, 20481, 45057, 98305, 212993, 458753, 983041, 2097153, 4456449, 9437185, 19922945, 41943041, 88080385, 184549377, 385875969, 805306369, 1677721601, 3489660929, 7247757313, 15032385537, 31138512897, 64424509441
Offset: 0

Views

Author

Keywords

Comments

a(n) also gives number of 0's in binary numbers 1 to 111..1 (n+1 bits). - Stephen G Penrice, Oct 01 2000
Numerator of m(n) = (m(n-1)+n)/2, m(0)=0. Denominator is A000079. - Reinhard Zumkeller, Feb 23 2002
a(n) is the number of directed column-convex polyominoes of area n+2 having along the lower contour exactly one vertical step that is followed by a horizontal step (a reentrant corner). - Emeric Deutsch, May 21 2003
a(n) is the number of bits in binary numbers from 1 to 111...1 (n bits). Partial sums of A001787. - Emeric Deutsch, May 24 2003
Genus of graph of n-cube = a(n-3) = 1+(n-4)*2^(n-3), n>1.
Sum of ordered partitions of n where each element is summed via T(e-1). See A066185 for more information. - Jon Perry, Dec 12 2003
a(n-2) is the number of Dyck n-paths with exactly one peak at height >= 3. For example, there are 5 such paths with n=4: UUUUDDDD, UUDUUDDD, UUUDDUDD, UDUUUDDD, UUUDDDUD. - David Callan, Mar 23 2004
Permutations in S_{n+2} avoiding 12-3 that contain the pattern 13-2 exactly once.
a(n) is prime for n = 2, 3, 7, 27, 51, 55, 81. a(n) is semiprime for n = 4, 5, 6, 8, 9, 10, 11, 13, 15, 19, 28, 32, 39, 57, 63, 66, 75, 97. - Jonathan Vos Post, Jul 18 2005
A member of the family of sequences defined by a(n) = Sum_{i=1..n} i*[c(1)*...*c(r)]^(i-1). This sequence has c(1)=2, A014915 has c(1)=3. - Ctibor O. Zizka, Feb 23 2008
Starting with 1 = row sums of A023758 as a triangle by rows: [1; 2,3; 4,6,7; 8,12,14,15; ...]. - Gary W. Adamson, Jul 18 2008
Equivalent formula given in Brehm: for each q >= 3 there exists a polyhedral map M_q of type {4, q} with [number of vertices] f_0 = 2^q and [genus] g = (2^(q-3))*(q-4) + 1 such that M_q and its dual have polyhedral embeddings in R^3 [McMullen et al.]. - Jonathan Vos Post, Jul 25 2009
Sums of rows of the triangle in A173787. - Reinhard Zumkeller, Feb 28 2010
This sequence is related to A000079 by a(n) = n*A000079(n)-Sum_{i=0..n-1} A000079(i). - Bruno Berselli, Mar 06 2012
(1 + 5*x + 17*x^2 + 49*x^3 + ...) = (1 + 2*x + 4*x^2 + 8*x^3 + ...) * (1 + 3*x + 7*x^2 + 15*x^3 + ...). - Gary W. Adamson, Mar 14 2012
The first barycentric coordinate of the centroid of Pascal triangles, assuming that numbers are weights, is A000295(n+1)/A000337(n), no matter what the triangle sides are. See attached figure. - César Eliud Lozada, Nov 14 2014
a(n) is the n-th number that is a sum of n positive n-th powers for n >= 1. a(4) = 49 = A003338(4). - Alois P. Heinz, Aug 01 2020
a(n) is the sum of the largest elements of all subsets of {1,2,..,n}. For example, a(3)=17; the subsets of {1,2,3} are {1}, {2}, {3}, {1,2}, {1,3}, {2,3}, {1,2,3}, and the sum of the largest elements is 17. - Enrique Navarrete, Aug 20 2020
a(n-1) is the sum of the second largest elements of the subsets of {1,2,..,n} that contain n. For example, for n = 4, a(3)=17; the subsets of {1,2,3,4} that contain 4 are {4}, {1,4}, {2,4}, {3,4}, {1,2,4}, {1,3,4}, {2,3,4}, {1,2,3,4}, and the sum of the second largest elements is 17. - Enrique Navarrete, Aug 24 2020
a(n-1) is also the sum of diameters of all subsets of {1,2,...,n} that contain n. For example, for n = 4, a(3)=17; the subsets of {1,2,3,4} that contain 4 are {4}, {1,4}, {2,4}, {3,4}, {1,2,4}, {1,3,4}, {2,3,4}, {1,2,3,4}; the diameters of these sets are 0,3,2,1,3,3,2,3 and the sum is 17. - Enrique Navarrete, Sep 07 2020
a(n-1) is also the number of additions required to compute the permanent of general n X n matrices using trellis methods (see Theorems 5 and 6, pp. 10-11 in Kiah et al.). - Stefano Spezia, Nov 02 2021

References

  • F. Harary, Topological concepts in graph theory, pp. 13-17 of F. Harary and L. Beineke, editors, A seminar on Graph Theory, Holt, Rinehart and Winston, New York, 1967.
  • V. G. Gutierrez and S. L. de Medrano, Surfaces as complete intersections, in Riemann and Klein Surfaces, Automorphisms, Symmetries and Moduli Spaces, edited by Milagros Izquierdo, S. Allen Broughton, Antonio F. Costa, Contemp. Math. vol. 629, 2014, pp. 171-.
  • F. Harary, Graph Theory. Addison-Wesley, Reading, MA, 1969, p. 119.
  • G. H. Hardy, A Theorem Concerning the Infinite Cardinal Numbers, Quart. J. Math., 35 (1904), p. 90 = Collected Papers of G. H. Hardy, Vol. VII, p. 430.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n) = T(3, n), array T given by A048472. A036799/2.
Cf. A003338.
Main diagonal of A336725.

Programs

  • GAP
    List([0..30],n->(n-1)*2^n+1); # Muniru A Asiru, Oct 24 2018
  • Magma
    [(n-1)*2^n + 1: n in [0..40]]; // Vincenzo Librandi, Nov 21 2014
    
  • Maple
    A000337 := proc(n) 1+(n-1)*2^n ; end proc: # R. J. Mathar, Oct 10 2011
  • Mathematica
    Table[Sum[(-1)^(n - k) k (-1)^(n - k) Binomial[n + 1, k + 1], {k, 0, n}], {n, 0, 28}] (* Zerinvary Lajos, Jul 08 2009 *)
    Table[(n - 1) 2^n + 1, {n, 0, 40}] (* Harvey P. Dale, Jun 21 2011 *)
    LinearRecurrence[{5, -8, 4}, {0, 1, 5}, 40] (* Harvey P. Dale, Jun 21 2011 *)
    CoefficientList[Series[x / ((1 - x) (1 - 2 x)^2), {x, 0, 50}], x] (* Vincenzo Librandi, Nov 21 2014 *)
  • PARI
    a(n)=if(n<0,0,(n-1)*2^n+1)
    
  • Python
    a=lambda n:((n-1)<<(n))+1 # Indranil Ghosh, Jan 05 2017
    

Formula

Binomial transform of A004273. Binomial transform of A008574 if the leading zero is dropped.
G.f.: x/((1-x)*(1-2*x)^2). - Simon Plouffe in his 1992 dissertation
E.g.f.: exp(x) - exp(2*x)*(1-2*x). a(n) = 4*a(n-1) - 4*a(n-2)+1, n>0. Series reversion of g.f. A(x) is x*A034015(-x). - Michael Somos
Binomial transform of n/(n+1) is a(n)/(n+1). - Paul Barry, Aug 19 2005
a(n) = A119258(n+1,n-1) for n>0. - Reinhard Zumkeller, May 11 2006
Convolution of "Number of fixed points in all 231-avoiding involutions in S_n" (A059570) with "The odd numbers" (A005408), treating the result as if offset=0. - Graeme McRae, Jul 12 2006
a(n) = Sum_{k=1..n} k*2^(k-1), partial sums of A001787. - Zerinvary Lajos, Oct 19 2006
a(n) = 5*a(n-1) - 8*a(n-2) + 4*a(n-3), n > 2. - Harvey P. Dale, Jun 21 2011
a(n) = Sum_{k=1..n} Sum_{i=1..n} i * C(k,i). - Wesley Ivan Hurt, Sep 19 2017
a(n) = A000295(n+1)^2 - A000295(n)*A000295(n+2). - Gregory Gerard Wojnar, Oct 23 2018

A016031 De Bruijn's sequence: 2^(2^(n-1) - n): number of ways of arranging 2^n bits in circle so all 2^n consecutive strings of length n are distinct.

Original entry on oeis.org

1, 1, 2, 16, 2048, 67108864, 144115188075855872, 1329227995784915872903807060280344576, 226156424291633194186662080095093570025917938800079226639565593765455331328
Offset: 1

Views

Author

Keywords

Comments

Sequence corresponds also to the largest number that may be determined by asking no more than 2^(n-1) - 1 questions("Yes"-or-"No" answerable) with lying allowed at most once. - Lekraj Beedassy, Jul 15 2002
Number of Ouroborean rings for binary n-tuplets. - Lekraj Beedassy, May 06 2004
Also the number of games of Nim that are wins for the second player when the starting position is either the empty heap or heaps of sizes 1 <= t_1 < t_2 < ... < t_k < 2^(n-1) (if n is 1, the only starting position is the empty heap). E.g.: a(4) = 16: the winning positions for the second player when all the heap sizes are different and less than 2^3: (4,5,6,7), (3,5,6), (3,4,7), (2,5,7), (2,4,6), (2,3,6,7), (2,3,4,5), (1,6,7), (1,4,5), (1,3,5,7), (1,3,4,6), (1,2,5,6), (1,2,4,7), (1,2,3), (1,2,3,4,5,6,7) and the empty heap. - Kennan Shelton (kennan.shelton(AT)gmail.com), Apr 14 2006
a(n + 1) = 2^(2^n-n-1) = 2^A000295(n) is also the number of set-systems on n vertices with no singletons. The case with singletons is A058891. The unlabeled case is A317794. The spanning/covering case is A323816. The antichain case is A006126. The connected case is A323817. The uniform case is A306021(n) - 1. The graphical case is A006125. The chain case is A005840. - Gus Wiseman, Feb 01 2019
Named after the Dutch mathematician Nicolaas Govert de Bruijn (1918-2012). - Amiram Eldar, Jun 20 2021

References

  • Jonathan L. Gross and Jay Yellen, eds., Handbook of Graph Theory, CRC Press, 2004, p. 255.
  • Frank Harary and Edgar M. Palmer, Graphical Enumeration, 1973, p. 31.
  • D. J. Newman, "A variation of the Game of Twenty Question", in: Murray S. Klamkin (ed.), Problems in Applied Mathematics, Philadelphia, PA: SIAM, 1990, Problem 66-20, pp. 121-122.
  • Richard P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Cor. 5.6.15.
  • Ian Stewart, Game, Set and Math, pp. 50, Penguin 1991.

Crossrefs

Cf. A000295, A003465, A006125, A058891 (set systems), A317794 (unlabeled case), A323816 (spanning case), A323817 (connected case), A331691 (alternating signs).

Programs

Formula

a(n) = 2^A000295(n-1). - Lekraj Beedassy, Jan 17 2007
Shifting once to the left gives the binomial transform of A323816. - Gus Wiseman, Feb 01 2019

A008949 Triangle read by rows of partial sums of binomial coefficients: T(n,k) = Sum_{i=0..k} binomial(n,i) (0 <= k <= n); also dimensions of Reed-Muller codes.

Original entry on oeis.org

1, 1, 2, 1, 3, 4, 1, 4, 7, 8, 1, 5, 11, 15, 16, 1, 6, 16, 26, 31, 32, 1, 7, 22, 42, 57, 63, 64, 1, 8, 29, 64, 99, 120, 127, 128, 1, 9, 37, 93, 163, 219, 247, 255, 256, 1, 10, 46, 130, 256, 382, 466, 502, 511, 512, 1, 11, 56, 176, 386, 638, 848, 968, 1013, 1023, 1024, 1, 12, 67, 232, 562, 1024, 1486, 1816, 1981, 2036, 2047, 2048
Offset: 0

Views

Author

Keywords

Comments

The second-left-from-middle column is A000346: T(2n+2, n) = A000346(n). - Ed Catmur (ed(AT)catmur.co.uk), Dec 09 2006
T(n,k) is the maximal number of regions into which n hyperplanes of co-dimension 1 divide R^k (the Cake-Without-Icing numbers). - Rob Johnson, Jul 27 2008
T(n,k) gives the number of vertices within distance k (measured along the edges) of an n-dimensional unit cube, (i.e., the number of vertices on the hypercube graph Q_n whose distance from a reference vertex is <= k). - Robert Munafo, Oct 26 2010
A triangle formed like Pascal's triangle, but with 2^n for n >= 0 on the right border instead of 1. - Boris Putievskiy, Aug 18 2013
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 04 2013
Consider each "1" as an apex of two sequences: the first is the set of terms in the same row as the "1", but the rightmost term in the row repeats infinitely. Example: the row (1, 4, 7, 8) becomes (1, 4, 7, 8, 8, 8, ...). The second sequence begins with the same "1" but is the diagonal going down and to the right, thus: (1, 5, 16, 42, 99, 219, 466, ...). It appears that for all such sequence pairs, the binomial transform of the first, (1, 4, 7, 8, 8, 8, ...) in this case; is equal to the second: (1, 5, 16, 42, 99, ...). - Gary W. Adamson, Aug 19 2015
Let T* be the infinite tree with root 0 generated by these rules: if p is in T*, then p+1 is in T* and x*p is in T*. Let q(n) be the sum of polynomials in the n-th generation of T*. For n >= 0, row n of A008949 gives the coefficients of q(n+1); e.g., (row 3) = (1, 4, 7, 8) matches x^3 + 4*x^2 + 7*x + 9, which is the sum of the 8 polynomials in the 4th generation of T*. - Clark Kimberling, Jun 16 2016
T(n,k) is the number of subsets of [n]={1,...,n} of at most size k. Equivalently, T(n,k) is the number of subsets of [n] of at least size n-k. Counting the subsets of at least size (n-k) by conditioning on the largest element m of the smallest (n-k) elements of such a subset provides the formula T(n,k) = Sum_{m=n-k..n} C(m-1,n-k-1)*2^(n-m), and, by letting j=m-n+k, we obtain T(n,k) = Sum_{j=0..k} C(n+j-k-1,j)*2^(k-j). - Dennis P. Walsh, Sep 25 2017
If the interval of integers 1..n is shifted up or down by k, making the new interval 1+k..n+k or 1-k..n-k, then T(n-1,n-1-k) (= 2^(n-1)-T(n-1,k-1)) is the number of subsets of the new interval that contain their own cardinal number as an element. - David Pasino, Nov 01 2018

Examples

			Triangle begins:
  1;
  1,  2;
  1,  3,  4;
  1,  4,  7,   8;
  1,  5, 11,  15,  16;
  1,  6, 16,  26,  31,  32;
  1,  7, 22,  42,  57,  63,  64;
  1,  8, 29,  64,  99, 120, 127, 128;
  1,  9, 37,  93, 163, 219, 247, 255,  256;
  1, 10, 46, 130, 256, 382, 466, 502,  511,  512;
  1, 11, 56, 176, 386, 638, 848, 968, 1013, 1023, 1024;
  ...
		

References

  • F. J. MacWilliams and N. J. A. Sloane, The Theory of Error-Correcting Codes, Elsevier-North Holland, 1978, p. 376.

Crossrefs

Row sums sequence is A001792.
T(n, m)= A055248(n, n-m).

Programs

  • GAP
    T:=Flat(List([0..11],n->List([0..n],k->Sum([0..k],j->Binomial(n+j-k-1,j)*2^(k-j))))); # Muniru A Asiru, Nov 25 2018
    
  • Haskell
    a008949 n k = a008949_tabl !! n !! k
    a008949_row n = a008949_tabl !! n
    a008949_tabl = map (scanl1 (+)) a007318_tabl
    -- Reinhard Zumkeller, Nov 23 2012
    
  • Magma
    [[(&+[Binomial(n,j): j in [0..k]]): k in [0..n]]: n in [0..12]]; // G. C. Greubel, Nov 25 2018
    
  • Maple
    A008949 := proc(n,k) local i; add(binomial(n,i),i=0..k) end; # Typo corrected by R. J. Mathar, Oct 26 2010
  • Mathematica
    Table[Length[Select[Subsets[n], (Length[ # ] <= k) &]], {n, 0, 12}, {k, 0, n}] // Grid (* Geoffrey Critzer, May 13 2009 *)
    Flatten[Accumulate/@Table[Binomial[n,i],{n,0,20},{i,0,n}]] (* Harvey P. Dale, Aug 08 2015 *)
    T[ n_, k_] := If[ n < 0 || k > n, 0, Binomial[n, k] Hypergeometric2F1[1, -k, n + 1 - k, -1]]; (* Michael Somos, Aug 05 2017 *)
  • PARI
    A008949(n)=T8949(t=sqrtint(2*n-sqrtint(2*n)),n-t*(t+1)/2)
    T8949(r,c)={ 2*c > r || return(sum(i=0,c,binomial(r,i))); 1<M. F. Hasler, May 30 2010
    
  • PARI
    {T(n, k) = if(k>n, 0, sum(i=0, k, binomial(n, i)))}; /* Michael Somos, Aug 05 2017 */
    
  • PARI
    row(n) = my(v=vector(n+1, k, binomial(n,k-1))); vector(#v, k, sum(i=1, k, v[i])); \\ Michel Marcus, Apr 13 2025
    
  • Sage
    [[sum(binomial(n,j) for j in range(k+1)) for k in range(n+1)] for n in range(12)] # G. C. Greubel, Nov 25 2018

Formula

From partial sums across rows of Pascal triangle A007318.
T(n, 0) = 1, T(n, n) = 2^n, T(n, k) = T(n-1, k-1) + T(n-1, k), 0 < k < n.
G.f.: (1 - x*y)/((1 - y - x*y)*(1 - 2*x*y)). - Antonio Gonzalez (gonfer00(AT)gmail.com), Sep 08 2009
T(2n,n) = A032443(n). - Philippe Deléham, Sep 16 2009
T(n,k) = 2 T(n-1,k-1) + binomial(n-1,k) = 2 T(n-1,k) - binomial(n-1,k). - M. F. Hasler, May 30 2010
T(n,k) = binomial(n,n-k)* 2F1(1, -k; n+1-k; -1). - Olivier Gérard, Aug 02 2012
For a closed-form formula for arbitrary left and right borders of Pascal like triangle see A228196. - Boris Putievskiy, Aug 18 2013
T(n,floor(n/2)) = A027306(n). - Reinhard Zumkeller, Nov 14 2014
T(n,n) = 2^n, otherwise for 0 <= k <= n-1, T(n,k) = 2^n - T(n,n-k-1). - Bob Selcoe, Mar 30 2017
For fixed j >= 0, lim_{n -> oo} T(n+1,n-j+1)/T(n,n-j) = 2. - Bob Selcoe, Apr 03 2017
T(n,k) = Sum_{j=0..k} C(n+j-k-1,j)*2^(k-j). - Dennis P. Walsh, Sep 25 2017

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Mar 23 2000

A014824 a(0) = 0; for n>0, a(n) = 10*a(n-1) + n.

Original entry on oeis.org

0, 1, 12, 123, 1234, 12345, 123456, 1234567, 12345678, 123456789, 1234567900, 12345679011, 123456790122, 1234567901233, 12345679012344, 123456790123455, 1234567901234566, 12345679012345677, 123456790123456788, 1234567901234567899, 12345679012345679010, 123456790123456790121
Offset: 0

Views

Author

Keywords

Comments

The square roots of these numbers have some remarkable properties - see the link to Schizophrenic numbers.
Partial sums of A002275. - Jonathan Vos Post, Apr 25 2010
This sequence is the particular case of a(0) = 0, a(n) = r*a(n-1) + n, when r = 10. If now the first N terms are computed for (r > N) then the resulting set of numbers is readable as the smallest k-digits permutations (1 <= k <= N): Those built from the concatenation of the first k digits in base-r (see links). - R. J. Cano, Jan 09 2013
There is also an interesting structure to the decimal expansion of 1/sqrt(a(2*n+1)), which has long strings of 0's (that gradually shorten in length until they disappear) interspersed with strings of what, at first sight, appear to be random digits. However, if we factorize these blocks of 'random' digits we find they are related to each other. An example illustrating this is given below. - Peter Bala, Sep 13 2015
From Peter Bala, Sep 15 2015: (Start)
Extending the previous empirical observation, it appears that the decimal expansions of the numbers 1/ (a(4*n+1))^(1/2), 1/a((4*n+3))^(1/4), 1/(10*a(4*n))^(1/2), 1/(10*a(4*n+2))^(1/4), and their powers, have the same pattern of beginning with long strings of 0's (that gradually shorten in length until they disappear) interlaced with strings of digits which, when read as integers and factorized, turn out to be related to each other.
The following result, which is a consequence of Bottomley's explicit formula for a(n), should be helpful in explaining these observations: the decimal expansion of 1/a(2*n-1) for n >= 5 begins with long strings of 0's interlaced successively with the digits of the numbers 81*(18*n + 1)^k for k going from 0 up to approximately n/log_10(18*n). For example, for n = 7 we have 1/a(13) = 0.00000000000081000000000102870000000130644900000 165919023..., with 10287 = 81*127, 1306449 = 81*127^2 and 165919023 = 81*127^3. A similar result holds for the decimal expansion of the number 1/a(2*n).
It appears that a(4*n+3)^(1/4), a(4*n+3)^(3/4), (10*a(4*n))^(1/2), (10*a(4*n+2))^(1/4) and (10*a(4*n+2))^(3/4) are further examples of Brown's schizophrenic numbers.
A theorem of Kuzmin in the measure theory of continued fractions says that for a random real number alpha, the probability that some given partial quotient of alpha is equal to a positive integer k is given by 1/log(2)*( log(1 + 1/k) - log(1 + 1/(k+1)) ). Thus large partial quotients are the exception in continued fraction expansions. Empirically, we observe the presence of unexpectedly large partial quotients early in the continued fraction expansions of the numbers (a(4*n+1))^(1/2), (a(4*n+3))^(1/4), (10*a(4*n))^(1/2), (10*a(4*n+2))^(1/4) and their powers. An example is given below. (End)
From Ya-Ping Lu, Dec 21 2024: (Start)
To get a(n), concatenate the first n digits in the cyclic string '123456790' and subtract the number of occurrences of '9' from the concatenated number. For example, a(8) = 12345679 - 1 = 12345678.
There are 2 prime terms for n <= 20000: a(2497) and a(3301). (End)

Examples

			From _Peter Bala_, Sep 13 2015: (Start)
The decimal expansion of 1/sqrt(a(51)) begins 9.0...0211050...07423683750...02901423065625000... x 10^(-26). The long strings of 0's gradually shorten in length and are interspersed with eleven blocks of digits  [9, 21105, 742368375, 2901423065625, 1190671490555859375, 5025824361636282421875, 216068565680679841787109375, 940978603539360710982861328125, 4137365297437126626102768402099609375, 18326229731370116994398540261077880859375, 816525165681195562685426961332324981689453125]. Read as ordinary integers these numbers factorize as [3^2, (3^2)*5*7*67, (3^3)*(5^3)*(7^2)*(67^2), (3^2)*(5^5)*(7^3)*(67^3), (3^2)*(5^8)*(7^5)*(67^4), (3^4)*(5^8)*(7^6)*(67^5), (3^3)*(5^10)*(7^7)*11*(67^6), (3^3)*(5^11)*(7^7)*11*13*(67^7), (3^4)*(5^16)*(7^8)*11*13*(67^8), (3^2)*(5^17)*(7^9)*11*13*17*(67^9), (3^2)*(5^18)*(7^10)*11*13*17*19*(67^10)]. (End)
From _Peter Bala_, Sep 15 2015: (Start)
The continued fraction expansion of 1/sqrt(a(51)) begins [0; 11111111111111111111111111, 9, 47382136934375740345889, 2, 21, 3, 1, 7, 2, 1, 101028010521057015662, 5, 14, 9, 1, 1, 2, 2, 8, 5, 1, 1, 1, 1, 215411536292232442, 5, 1, 5, 1, 1, 2, 1, 1, 8, 1, 4, 3, 1, 4, 2, 1, 8, 1, 1, 3, 10, 459299650942926, 4, 1, 1, 4, 1, 20, 64, 5, 9, 2, 2, 1, 2, 1, 1, 1, 1, 30, 1, 11, 3, 979316952969, 1, 2, 93, 1, 5, 1, 1, 11, 1, 1, 1, 1, 5, 1, 29, 1, 29, 1, 1, 1, 2, 4, 1, 37, 1, 1, 2, 8, 2, 2088095848, 12, 1, 3, 1, 3, 2, 2, 3, 1, 5, 6, 1, 3, 1, 4, 2, 2, 1, 2, 2, 14, 4, 1, 2, 1, 50, 2, 6, 1, 11, 135, 4452229, 1, ...] and has several unexpectedly large partial quotients early on. (End)
For n=5, a(5) = 1*15 + 9*20 + 9^2*15 + 9^3*6 + 9^4*1 + 9^5*0 = 12345. - _Bruno Berselli_, Nov 13 2015
		

Crossrefs

Cf. A060011.
Cf. A002275. - Jonathan Vos Post, Apr 25 2010
Similar sequences in other bases are: (base-2) A000295, (base-3) A000340, (base-4) A014825, (base-5) A014827, (base-6) A014829. - R. J. Cano, Jan 11 2013
Differs from A007908, A035239, A057137, A060555, A138957 from n=10 on. - M. F. Hasler, Jan 17 2013
Cf. A030512.

Programs

  • Magma
    [(10^n-1)*(10/81)-n/9: n in [0..20]]; // Vincenzo Librandi, Aug 23 2011
    
  • Maple
    a:=n->sum((10^(n-j)-1^(n-j))/9,j=0..n): seq(a(n), n=0..17); # Zerinvary Lajos, Jan 15 2007
    a:=n->sum(10^(n-j)*j,j=0..n): seq(a(n), n=0..16); # Zerinvary Lajos, Jun 05 2008
  • Mathematica
    Table[Sum[10^i - 1, {i, n}]/9, {n, 18}] (* Robert G. Wilson v, Nov 20 2004 *)
    CoefficientList[Series[x/(1 - 12*x + 21*x^2 - 10*x^3), {x, 0, 20}], x] (* Wesley Ivan Hurt, Sep 15 2015 *)
  • PARI
    linrec01(p,u,base)={my(r=!p,A=1);for(j=2,u,A=A*base+r+p*j); A};
    a(n)=(n!=0)*linrec01(1, n, 10); \\ R. J. Cano, Jan 09 2011; With (0, n, 10) it generates repunit numbers.
    
  • PARI
    A014824(n)=(10^(n+1)\9-n)\9  \\ M. F. Hasler, Jan 17 2013
    
  • Python
    def A014824(n): s = ''.join('123456790'[i%9] for i in range(n)); q, r = divmod(n, 9); return int(s) - q - r//8 # Ya-Ping Lu, Dec 21 2024

Formula

a(n) = (10^n-1)*(10/81) - n/9. - Henry Bottomley, Jul 04 2000
a(n)/10^n converges to 10/81 = 0.123456790123456790...
Let b(n) = if(n = 0, 1, if(n = 1, 10, 10*9^(n-2))). Then a(n) = Sum_{k=0..n} C(n, k)*b(k) (Binomial transform). - Paul Barry, Jan 29 2004
G.f.: x/(1-12*x+21*x^2-10*x^3). - Colin Barker, Jan 08 2012
a(n) = 12*a(n-1) - 21*a(n-2) + 10*a(n-3), n>2. - Wesley Ivan Hurt, Sep 15 2015
a(n) = Sum_{i=0..n} 9^i*binomial(n+1,n-1-i). - Bruno Berselli, Nov 13 2015
a(n) = Sum_{i=0..n} 10^(n-i)*i. - Ya-Ping Lu, Dec 21 2024
E.g.f.: exp(x)*(10*exp(9*x) - 9*x - 10)/81. - Elmo R. Oliveira, Mar 29 2025

A145271 Coefficients for expansion of (g(x)d/dx)^n g(x); refined Eulerian numbers for calculating compositional inverse of h(x) = (d/dx)^(-1) 1/g(x); iterated derivatives as infinitesimal generators of flows.

Original entry on oeis.org

1, 1, 1, 1, 1, 4, 1, 1, 11, 4, 7, 1, 1, 26, 34, 32, 15, 11, 1, 1, 57, 180, 122, 34, 192, 76, 15, 26, 16, 1, 1, 120, 768, 423, 496, 1494, 426, 294, 267, 474, 156, 56, 42, 22, 1, 1, 247, 2904, 1389, 4288, 9204, 2127, 496, 5946, 2829, 5142, 1206, 855, 768, 1344, 1038, 288, 56, 98, 64, 29, 1
Offset: 0

Views

Author

Tom Copeland, Oct 06 2008

Keywords

Comments

For more detail, including connections to Legendre transformations, rooted trees, A139605, A139002 and A074060, see Mathemagical Forests p. 9.
For connections to the h-polynomials associated to the refined f-polynomials of permutohedra see my comments in A008292 and A049019.
From Tom Copeland, Oct 14 2011: (Start)
Given analytic functions F(x) and FI(x) such that F(FI(x))=FI(F(x))=x about 0, i.e., they are compositional inverses of each other, then, with g(x) = 1/dFI(x)/dx, a flow function W(s,x) can be defined with the following relations:
W(s,x) = exp(s g(x)d/dx)x = F(s+FI(x)) ,
W(s,0) = F(s) ,
W(0,x) = x ,
dW(0,x)/ds = g(x) = F'[FI(x)] , implying
dW(0,F(x))/ds = g(F(x)) = F'(x) , and
W(s,W(r,x)) = F(s+FI(F(r+FI(x)))) = F(s+r+FI(x)) = W(s+r,x) . (See MF link below.) (End)
dW(s,x)/ds - g(x)dW(s,x)/dx = 0, so (1,-g(x)) are the components of a vector orthogonal to the gradient of W and, therefore, tangent to the contour of W, at (s,x) . - Tom Copeland, Oct 26 2011
Though A139605 contains A145271, the op. of A145271 contains that of A139605 in the sense that exp(s g(x)d/dx) w(x) = w(F(s+FI(x))) = exp((exp(s g(x)d/dx)x)d/du)w(u) evaluated at u=0. This is reflected in the fact that the forest of rooted trees assoc. to (g(x)d/dx)^n, FOR_n, can be generated by removing the single trunk of the planted rooted trees of FOR_(n+1). - Tom Copeland, Nov 29 2011
Related to formal group laws for elliptic curves (see Hoffman). - Tom Copeland, Feb 24 2012
The functional equation W(s,x) = F(s+FI(x)), or a restriction of it, is sometimes called the Abel equation or Abel's functional equation (see Houzel and Wikipedia) and is related to Schröder's functional equation and Koenigs functions for compositional iterates (Alexander, Goryainov and Kudryavtseva). - Tom Copeland, Apr 04 2012
g(W(s,x)) = F'(s + FI(x)) = dW(s,x)/ds = g(x) dW(s,x)/dx, connecting the operators here to presentations of the Koenigs / Königs function and Loewner / Löwner evolution equations of the Contreras et al. papers. - Tom Copeland, Jun 03 2018
The autonomous differential equation above also appears with a change in variable of the form x = log(u) in the renormalization group equation, or Beta function. See Wikipedia, Zinn-Justin equations 2.10 and 3.11, and Krajewski and Martinetti equation 21. - Tom Copeland, Jul 23 2020
A variant of these partition polynomials appears on p. 83 of Petreolle et al. with the indeterminates e_n there related to those given in the examples below by e_n = n!*(n'). The coefficients are interpreted as enumerating certain types of trees. See also A190015. - Tom Copeland, Oct 03 2022

Examples

			From _Tom Copeland_, Sep 19 2014: (Start)
Let h(x) = log((1+a*x)/(1+b*x))/(a-b); then, g(x) = 1/(dh(x)/dx) = (1+ax)(1+bx), so (0')=1, (1')=a+b, (2')=2ab, evaluated at x=0, and higher order derivatives of g(x) vanish. Therefore, evaluated at x=0,
R^0 g(x) =  1
R^1 g(x) =  a+b
R^2 g(x) = (a+b)^2 + 2ab = a^2 + 4 ab + b^2
R^3 g(x) = (a+b)^3 + 4*(a+b)*2ab = a^3 + 11 a^2*b + 11 ab^2 + b^3
R^4 g(x) = (a+b)^4 + 11*(a+b)^2*2ab + 4*(2ab)^2
         =  a^4 + 26 a^3*b + 66 a^2*b^2 + 26 ab^3 + b^4,
etc., and these bivariate Eulerian polynomials (A008292) are the first few coefficients of h^(-1)(x) = (e^(ax) - e^(bx))/(a*e^(bx) - b*e^(ax)), the inverse of h(x). (End)
Triangle starts:
  1;
  1;
  1,   1;
  1,   4,    1;
  1,  11,    4,    7,    1;
  1,  26,   34,   32,   15,   11,    1;
  1,  57,  180,  122,   34,  192,   76,  15,   26,   16,    1;
  1, 120,  768,  423,  496, 1494,  426, 294,  267,  474,  156,   56,  42,  22,    1;
  1, 247, 2904, 1389, 4288, 9204, 2127, 496, 5946, 2829, 5142, 1206, 855, 768, 1344, 1038, 288, 56, 98, 64, 29, 1;
		

References

  • D. S. Alexander, A History of Complex Dynamics: From Schröder to Fatou to Julia, Friedrich Vieweg & Sohn, 1994.
  • T. Mansour and M. Schork, Commutation Relations, Normal Ordering, and Stirling Numbers, Chapman and Hall/CRC, 2015.

Crossrefs

Cf. (A133437, A086810, A181289) = (LIF, reduced LIF, associated g(x)), where LIF is a Lagrange inversion formula. Similarly for (A134264, A001263, A119900), (A134685, A134991, A019538), (A133932, A111999, A007318).
Second column is A000295, subdiagonal is A000124, row sums are A000142, row lengths are A000041. - Peter Luschny, Jul 21 2016

Programs

  • Maple
    with(LinearAlgebra): with(ListTools):
    A145271_row := proc(n) local b, M, V, U, G, R, T;
    if n < 2 then return 1 fi;
    b := (n,k) -> `if`(k=1 or k>n+1,0,binomial(n-1,k-2)*g[n-k+1]);
    M := n -> Matrix(n, b):
    V := n -> Vector[row]([1, seq(0,i=2..n)]):
    U := n -> VectorMatrixMultiply(V(n), M(n)^(n-1)):
    G := n -> Vector([seq(g[i], i=0..n-1)]);
    R := n -> VectorMatrixMultiply(U(n), G(n)):
    T := Reverse([op(sort(expand(R(n+1))))]);
    seq(subs({seq(g[i]=1, i=0..n)},T[j]),j=1..nops(T)) end:
    for n from 0 to 9 do A145271_row(n) od; # Peter Luschny, Jul 20 2016

Formula

Let R = g(x)d/dx; then
R^0 g(x) = 1 (0')^1
R^1 g(x) = 1 (0')^1 (1')^1
R^2 g(x) = 1 (0')^1 (1')^2 + 1 (0')^2 (2')^1
R^3 g(x) = 1 (0')^1 (1')^3 + 4 (0')^2 (1')^1 (2')^1 + 1 (0')^3 (3')^1
R^4 g(x) = 1 (0')^1 (1')^4 + 11 (0')^2 (1')^2 (2')^1 + 4 (0')^3 (2')^2 + 7 (0')^3 (1')^1 (3')^1 + 1 (0')^4 (4')^1
R^5 g(x) = 1 (0') (1')^5 + 26 (0')^2 (1')^3 (2') + (0')^3 [34 (1') (2')^2 + 32 (1')^2 (3')] + (0')^4 [ 15 (2') (3') + 11 (1') (4')] + (0')^5 (5')
R^6 g(x) = 1 (0') (1')^6 + 57 (0')^2 (1')^4 (2') + (0')^3 [180 (1')^2 (2')^2 + 122 (1')^3 (3')] + (0')^4 [ 34 (2')^3 + 192 (1') (2') (3') + 76 (1')^2 (4')] + (0')^5 [15 (3')^2 + 26 (2') (4') + 16 (1') (5')] + (0')^6 (6')
where (j')^k = ((d/dx)^j g(x))^k. And R^(n-1) g(x) evaluated at x=0 is the n-th Taylor series coefficient of the compositional inverse of h(x) = (d/dx)^(-1) 1/g(x), with the integral from 0 to x.
The partitions are in reverse order to those in Abramowitz and Stegun p. 831. Summing over coefficients with like powers of (0') gives A008292.
Confer A190015 for another way to compute numbers for the array for each partition. - Tom Copeland, Oct 17 2014
Equivalent matrix computation: Multiply the n-th diagonal (with n=0 the main diagonal) of the lower triangular Pascal matrix by g_n = (d/dx)^n g(x) to obtain the matrix VP with VP(n,k) = binomial(n,k) g_(n-k). Then R^n g(x) = (1, 0, 0, 0, ...) [VP * S]^n (g_0, g_1, g_2, ...)^T, where S is the shift matrix A129185, representing differentiation in the divided powers basis x^n/n!. - Tom Copeland, Feb 10 2016 (An evaluation removed by author on Jul 19 2016. Cf. A139605 and A134685.)
Also, R^n g(x) = (1, 0, 0, 0, ...) [VP * S]^(n+1) (0, 1, 0, ...)^T in agreement with A139605. - Tom Copeland, Jul 21 2016
A recursion relation for computing each partition polynomial of this entry from the lower order polynomials and the coefficients of the cycle index polynomials of A036039 is presented in the blog entry "Formal group laws and binomial Sheffer sequences". - Tom Copeland, Feb 06 2018
A formula for computing the polynomials of each row of this matrix is presented as T_{n,1} on p. 196 of the Ihara reference in A139605. - Tom Copeland, Mar 25 2020
Indeterminate substitutions as illustrated in A356145 lead to [E] = [L][P] = [P][E]^(-1)[P] = [P][RT] and [E]^(-1) = [P][L] = [P][E][P] = [RT][P], where [E] contains the refined Eulerian partition polynomials of this entry; [E]^(-1), A356145, the inverse set to [E]; [P], the permutahedra polynomials of A133314; [L], the classic Lagrange inversion polynomials of A134685; and [RT], the reciprocal tangent polynomials of A356144. Since [L]^2 = [P]^2 = [RT]^2 = [I], the substitutional identity, [L] = [E][P] = [P][E]^(-1) = [RT][P], [RT] = [E]^(-1)[P] = [P][L][P] = [P][E], and [P] = [L][E] = [E][RT] = [E]^(-1)[L] = [RT][E]^(-1). - Tom Copeland, Oct 05 2022

Extensions

Title amplified by Tom Copeland, Mar 17 2014
R^5 and R^6 formulas and terms a(19)-a(29) added by Tom Copeland, Jul 11 2016
More terms from Peter Luschny, Jul 20 2016

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

Original entry on oeis.org

0, 0, 0, 1, 5, 16, 42, 99, 219, 466, 968, 1981, 4017, 8100, 16278, 32647, 65399, 130918, 261972, 524097, 1048365, 2096920, 4194050, 8388331, 16776915, 33554106, 67108512, 134217349, 268435049, 536870476, 1073741358, 2147483151, 4294966767, 8589934030
Offset: 0

Views

Author

Keywords

Comments

Number of subsets with at least 3 elements of an n-element set.
For n>4, number of simple rank-(n-1) matroids over S_n.
Number of non-interval subsets of {1,2,3,...,n} (cf. A000124). - Jose Luis Arregui (arregui(AT)unizar.es), Jun 27 2006
The partial sums of the second diagonal of A008292 or third column of A123125. - Tom Copeland, Sep 09 2008
a(n) is the number of binary sequences of length n having at least three 0's. - Geoffrey Critzer, Feb 11 2009
Starting with "1" = eigensequence of a triangle with the tetrahedral numbers (1, 4, 10, 20, ...) as the left border and the rest 1's. - Gary W. Adamson, Jul 24 2010
a(n) is also the number of crossing set partitions of [n+1] with two blocks. - Peter Luschny, Apr 29 2011
The Kn24 sums, see A180662, of triangle A065941 equal the terms (doubled) of this sequence minus the three leading zeros. - Johannes W. Meijer, Aug 14 2011
From L. Edson Jeffery, Dec 28 2011: (Start)
Nonzero terms of this sequence can be found from the row sums of the fourth sub-triangle extracted from Pascal's triangle as indicated below by braces:
1;
1, 1;
1, 2, 1;
{1}, 3, 3, 1;
{1, 4}, 6, 4, 1;
{1, 5, 10}, 10, 5, 1;
{1, 6, 15, 20}, 15, 6, 1;
... (End)
Partial sums of A000295 (Eulerian Numbers, Column 2).
Second differences equal 2^(n-2) - 1, for n >= 4. - Richard R. Forberg, Jul 11 2013
Starting (0, 0, 1, 5, 16, ...) is the binomial transform of (0, 0, 1, 2, 2, 2, ...). - Gary W. Adamson, Jul 27 2015
a(n - 1) is the rank of the divisor class group of the moduli space of stable rational curves with n marked points, see Keel p. 550. - Harry Richman, Aug 10 2024

Examples

			a(4) = 5 is the number of crossing set partitions of {1,2,..,5}, card{13|245, 14|235, 24|135, 25|134, 35|124}. - _Peter Luschny_, Apr 29 2011
		

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n) = A055248(n,3).
First differences are A000295.
Cf. also A000290, A001045.

Programs

Formula

G.f.: x^3/((1-2*x)*(1-x)^3).
a(n) = Sum_{k=0..n} binomial(n,k+3) = Sum_{k=3..n} binomial(n,k). - Paul Barry, Jul 30 2004
a(n+1) = 2*a(n) + binomial(n,2). - Paul Barry, Aug 23 2004
(1, 5, 16, 42, 99, ...) = binomial transform of (1, 4, 7, 8, 8, 8, ...). - Gary W. Adamson, Sep 30 2007
E.g.f.: exp(x)*(exp(x)-x^2/2-x-1). - Geoffrey Critzer, Feb 11 2009
a(n) = n - 2 + 3*a(n-1) - 2*a(n-2), for n >= 2. - Richard R. Forberg, Jul 11 2013
For n>1, a(n) = (1/4)*Sum_{k=1..n-2} 2^k*(n-k-1)*(n-k). For example, (1/4)*(2^1*(4*5) + 2^2*(3*4) + 2^3*(2*3) + 2^4*(1*2)) = 168/4 = 42. - J. M. Bergot, May 27 2014 [edited by Danny Rorabaugh, Apr 19 2015]
Convolution of A001045 and (A000290 shifted by one place). - Oboifeng Dira, Aug 16 2016
a(n) = Sum_{k=1..n-2} Sum_{i=1..n} (n-k-1) * C(k,i). - Wesley Ivan Hurt, Sep 19 2017
a(n) = 5*a(n-1) - 9*a(n-2) + 7*a(n-3) - 2*a(n-4) for n > 3. - Chai Wah Wu, Apr 03 2021
a(n) = a(n-1) + 1 + A000247(n-1). - Harry Richman, Aug 13 2024

A001891 Hit polynomials; convolution of natural numbers with Fibonacci numbers F(2), F(3), F(4), ....

Original entry on oeis.org

0, 1, 4, 10, 21, 40, 72, 125, 212, 354, 585, 960, 1568, 2553, 4148, 6730, 10909, 17672, 28616, 46325, 74980, 121346, 196369, 317760, 514176, 831985, 1346212, 2178250, 3524517, 5702824, 9227400, 14930285, 24157748, 39088098, 63245913, 102334080, 165580064
Offset: 0

Views

Author

Keywords

Comments

a(n) is the sum of the n-th row of the triangle in A119457 for n > 0. - Reinhard Zumkeller, May 20 2006
Convolution of odds (A005408) with Fibonacci numbers (A000045). - Graeme McRae, Jun 06 2006
Equals row sums of triangle A152203. - Gary W. Adamson, Nov 29 2008
Define a triangle by T(n,0) = n*(n+1)+1, T(n,n) = 1, and T(r,c) = T(r-1,c) + T(r-2,c-1). This triangle starts: 1; 3,1; 7,2,1; 13,5,2,1; 21,12,4,2,1; the sum of terms in row n is a(n+1). - J. M. Bergot, Apr 23 2013
a(n) = number of k-tuples (u(1), u(2), ..., u(k)) with 1 <= u(1) < u(2) < ... < u(k) <= n such that u(i) - u(i-1) <= 2 for i = 2,...,k. Changing the bound from 2 to 3, then 4, then 5, yields A356619, A356620, A356621. The patterns suggest that the limiting sequence as the bound increases is A000295. - Clark Kimberling, Aug 24 2022

References

  • J. Riordan, The enumeration of permutations with three-ply staircase restrictions, unpublished memorandum, Bell Telephone Laboratories, Murray Hill, NJ, Oct 1963. (See A001883)
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Partial sums of A001911.
A diagonal of triangle in A080061.
Right-hand column 5 of triangle A011794.

Programs

  • GAP
    List([0..40], n-> Fibonacci(n+5) -2*n-5); # G. C. Greubel, Jul 06 2019
  • Magma
    [Fibonacci(n+5)-(5+2*n): n in [0..40]]; // Vincenzo Librandi, Jun 07 2013
    
  • Mathematica
    LinearRecurrence[{3,-2,-1,1}, {0,1,4,10}, 40] (* Vladimir Joseph Stephan Orlovsky, Feb 16 2012 *)
    Table[Fibonacci[n+5] -(2*n+5), {n,0,40}] (* G. C. Greubel, Jul 06 2019 *)
    maxDiff = 2;
    Map[Length[Select[Map[{#, Max[Differences[#]]} &,
      Drop[Subsets[Range[#]], # + 1]], #[[2]] <= maxDiff &]] &,
      Range[16]] (* Peter J. C. Moses, Aug 14 2022 *)
  • PARI
    a(n)=([0,1,0,0; 0,0,1,0; 0,0,0,1; 1,-1,-2,3]^n*[0;1;4;10])[1,1] \\ Charles R Greathouse IV, Apr 08 2016
    
  • Sage
    [fibonacci(n+5) -2*n-5 for n in (0..40)] # G. C. Greubel, Jul 06 2019
    

Formula

G.f.: x*(1+x)/((1-x-x^2)*(1-x)^2). - Simon Plouffe in his 1992 dissertation
a(n) = Fibonacci(n+5) - (5+2*n). - Wolfdieter Lang
a(n) = a(n-1) + a(n-2) + (2n+1); a(-x)=0. - Barry E. Williams, Mar 27 2000
a(n) = 3*a(n-1) - 2*a(n-2) - a(n-3) + a(n-4). - Sam Lachterman (slachterman(AT)fuse.net), Sep 22 2003
a(n) - a(n-1) = A101220(2,1,n). - Ross La Haye, May 31 2006
a(n) = (-3 + (2^(-1-n)*((1-sqrt(5))^n*(-11+5*sqrt(5)) + (1+sqrt(5))^n*(11+5*sqrt(5)))) / sqrt(5) - 2*(1+n)). - Colin Barker, Mar 11 2017

A055248 Triangle of partial row sums of triangle A007318(n,m) (Pascal's triangle). Triangle A008949 read backwards. Riordan (1/(1-2x), x/(1-x)).

Original entry on oeis.org

1, 2, 1, 4, 3, 1, 8, 7, 4, 1, 16, 15, 11, 5, 1, 32, 31, 26, 16, 6, 1, 64, 63, 57, 42, 22, 7, 1, 128, 127, 120, 99, 64, 29, 8, 1, 256, 255, 247, 219, 163, 93, 37, 9, 1, 512, 511, 502, 466, 382, 256, 130, 46, 10, 1, 1024, 1023, 1013, 968, 848, 638, 386, 176, 56, 11, 1
Offset: 0

Views

Author

Wolfdieter Lang, May 26 2000

Keywords

Comments

In the language of the Shapiro et al. reference (also given in A053121) such a lower triangular (ordinary) convolution array, considered as matrix, belongs to the Riordan-group. The g.f. for the row polynomials p(n,x) (increasing powers of x) is 1/((1-2*z)*(1-x*z/(1-z))).
Binomial transform of the all 1's triangle: as a Riordan array, it factors to give (1/(1-x),x/(1-x))(1/(1-x),x). Viewed as a number square read by antidiagonals, it has T(n,k) = Sum_{j=0..n} binomial(n+k,n-j) and is then the binomial transform of the Whitney square A004070. - Paul Barry, Feb 03 2005
Riordan array (1/(1-2x), x/(1-x)). Antidiagonal sums are A027934(n+1), n >= 0. - Paul Barry, Jan 30 2005; edited by Wolfdieter Lang, Jan 09 2015
Eigensequence of the triangle = A005493: (1, 3, 10, 37, 151, 674, ...); row sums of triangles A011971 and A159573. - Gary W. Adamson, Apr 16 2009
Read as a square array, this is the generalized Riordan array ( 1/(1 - 2*x), 1/(1 - x) ) as defined in the Bala link (p. 5), which factorizes as ( 1/(1 - x), x/(1 - x) )*( 1/(1 - x), x )*( 1, 1 + x ) = P*U*transpose(P), where P denotes Pascal's triangle, A007318, and U is the lower unit triangular array with 1's on or below the main diagonal. - Peter Bala, Jan 13 2016

Examples

			The triangle a(n,m) begins:
n\m    0    1    2   3   4   5   6   7  8  9 10 ...
0:     1
1:     2    1
2:     4    3    1
3:     8    7    4   1
4:    16   15   11   5   1
5:    32   31   26  16   6   1
6:    64   63   57  42  22   7   1
7:   128  127  120  99  64  29   8   1
8:   256  255  247 219 163  93  37   9  1
9:   512  511  502 466 382 256 130  46 10  1
10: 1024 1023 1013 968 848 638 386 176 56 11  1
... Reformatted. - _Wolfdieter Lang_, Jan 09 2015
Fourth row polynomial (n=3): p(3,x)= 8 + 7*x + 4*x^2 + x^3.
The matrix inverse starts
   1;
  -2,   1;
   2,  -3,   1;
  -2,   5,  -4,    1;
   2,  -7,   9,   -5,    1;
  -2,   9, -16,   14,   -6,    1;
   2, -11,  25,-  30,   20,   -7,    1;
  -2,  13, -36,   55,  -50,   27,   -8,    1;
   2, -15,  49,  -91,  105,  -77,   35,   -9,  1;
  -2,  17, -64,  140, -196,  182, -112,   44, -10,   1;
   2, -19,  81, -204,  336, -378,  294, -156,  54, -11, 1;
   ...
which may be related to A029653. - _R. J. Mathar_, Mar 29 2013
From _Peter Bala_, Dec 23 2014: (Start)
With the array M(k) as defined in the Formula section, the infinite product M(0)*M(1)*M(2)*... begins
/1      \ /1        \ /1       \       /1       \
|2 1     ||0 1       ||0 1      |      |2  1     |
|4 3 1   ||0 2 1     ||0 0 1    |... = |4  5 1   |
|8 7 4 1 ||0 4 3 1   ||0 0 2 1  |      |8 19 9 1 |
|...     ||0 8 7 4 1 ||0 0 4 3 1|      |...      |
|...     ||...       ||...      |      |         |
= A143494. (End)
Matrix factorization of square array as P*U*transpose(P):
/1      \ /1        \ /1 1 1 1 ...\    /1  1  1  1 ...\
|1 1     ||1 1       ||0 1 2 3 ... |   |2  3  4  5 ... |
|1 2 1   ||1 1 1     ||0 0 1 3 ... | = |4  7 11 16 ... |
|1 3 3 1 ||1 1 1 1   ||0 0 0 1 ... |   |8 15 26 42 ... |
|...     ||...       ||...         |   |...            |
- _Peter Bala_, Jan 13 2016
		

Crossrefs

Column sequences: A000079 (powers of 2, m=0), A000225 (m=1), A000295 (m=2), A002662 (m=3), A002663 (m=4), A002664 (m=5), A035038 (m=6), A035039 (m=7), A035040 (m=8), A035041 (m=9), A035042 (m=10).
Row sums: A001792(n) = A055249(n, 0).
Alternating row sums: A011782.
Cf. A011971, A159573. - Gary W. Adamson, Apr 16 2009

Programs

  • Haskell
    a055248 n k = a055248_tabl !! n !! k
    a055248_row n = a055248_tabl !! n
    a055248_tabl = map reverse a008949_tabl
    -- Reinhard Zumkeller, Jun 20 2015
  • Maple
    T := (n,k) -> 2^n - (1/2)*binomial(n, k-1)*hypergeom([1, n + 1], [n-k + 2], 1/2).
    seq(seq(simplify(T(n,k)), k=0..n),n=0..10); # Peter Luschny, Oct 10 2019
  • Mathematica
    a[n_, m_] := Sum[ Binomial[n, m + j], {j, 0, n}]; Table[a[n, m], {n, 0, 10}, {m, 0, n}] // Flatten (* Jean-François Alcover, Jul 05 2013, after Paul Barry *)
    T[n_, k_] := Binomial[n, k] * Hypergeometric2F1[1, k - n, k + 1, -1];
    Flatten[Table[T[n, k], {n, 0, 7}, {k, 0, n}]]  (* Peter Luschny, Oct 06 2023 *)

Formula

a(n, m) = A008949(n, n-m), if n > m >= 0.
a(n, m) = Sum_{k=m..n} A007318(n, k) (partial row sums in columns m).
Column m recursion: a(n, m) = Sum_{j=m..n-1} a(j, m) + A007318(n, m) if n >= m >= 0, a(n, m) := 0 if n
G.f. for column m: (1/(1-2*x))*(x/(1-x))^m, m >= 0.
a(n, m) = Sum_{j=0..n} binomial(n, m+j). - Paul Barry, Feb 03 2005
Inverse binomial transform (by columns) of A112626. - Ross La Haye, Dec 31 2006
T(2n,n) = A032443(n). - Philippe Deléham, Sep 16 2009
From Peter Bala, Dec 23 2014: (Start)
Exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(8 + 7*x + 4*x^2/2! + x^3/3!) = 8 + 15*x + 26*x^2/2! + 42*x^3/3! + 64*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ).
Let M denote the present triangle. For k = 0,1,2,... define M(k) to be the lower unit triangular block array
/I_k 0\
\ 0 M/ having the k X k identity matrix I_k as the upper left block; in particular, M(0) = M. The infinite product M(0)*M(1)*M(2)*..., which is clearly well-defined, is equal to A143494 (but with a different offset). See the Example section. Cf. A106516. (End)
a(n,m) = Sum_{p=m..n} 2^(n-p)*binomial(p-1,m-1), n >= m >= 0, else 0. - Wolfdieter Lang, Jan 09 2015
T(n, k) = 2^n - (1/2)*binomial(n, k-1)*hypergeom([1, n+1], [n-k+2], 1/2). - Peter Luschny, Oct 10 2019
T(n, k) = binomial(n, k)*hypergeom([1, k - n], [k + 1], -1). - Peter Luschny, Oct 06 2023
n-th row polynomial R(n, x) = (2^n - x*(1 + x)^n)/(1 - x). These polynomials can be used to find series acceleration formulas for the constants log(2) and Pi. - Peter Bala, Mar 03 2025
Previous Showing 11-20 of 205 results. Next