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 19 results. Next

A133620 Binomial(n+p,n) mod n where p=10.

Original entry on oeis.org

0, 0, 1, 1, 3, 4, 2, 6, 2, 6, 1, 2, 1, 10, 5, 7, 1, 12, 1, 15, 18, 12, 1, 12, 21, 14, 4, 12, 1, 28, 1, 29, 1, 18, 6, 5, 1, 20, 14, 10, 1, 14, 1, 34, 15, 24, 1, 3, 8, 16, 18, 27, 1, 34, 23, 16, 1, 30, 1, 16, 1, 32, 17, 57, 40, 56, 1, 1, 47, 60, 1, 54, 1, 38, 36, 58, 12, 66, 1, 63, 10, 42, 1
Offset: 1

Views

Author

Hieronymus Fischer, Sep 30 2007

Keywords

Comments

Let d(m)...d(2)d(1)d(0) be the base-n representation of n+p. The relation a(n)=d(1) holds, if n is a prime index. For this reason there are infinitely many terms which are equal to 1.

Crossrefs

Programs

  • Mathematica
    Table[Mod[Binomial[n + 10, n], n], {n, 90}] (* Harvey P. Dale, Apr 04 2015 *)
  • PARI
    a(n) = binomial(n+10, n) % n \\ Michel Marcus, Jul 15 2013

Formula

a(n) = binomial(n+p,p) mod n.
a(n) = 1 if n is a prime > p, since binomial(n+p,n)==(1+floor(p/n))(mod n), provided n is a prime.
a(n) = A001287(n+10) mod n. - Michel Marcus, Jul 15 2013; corrected by Michel Marcus, Jan 27 2020
For n > 58060802, a(n) = 2*a(n-29030400) - a(n-58060800). - Ray Chandler, Apr 29 2023

A135278 Triangle read by rows, giving the numbers T(n,m) = binomial(n+1, m+1); or, Pascal's triangle A007318 with its left-hand edge removed.

Original entry on oeis.org

1, 2, 1, 3, 3, 1, 4, 6, 4, 1, 5, 10, 10, 5, 1, 6, 15, 20, 15, 6, 1, 7, 21, 35, 35, 21, 7, 1, 8, 28, 56, 70, 56, 28, 8, 1, 9, 36, 84, 126, 126, 84, 36, 9, 1, 10, 45, 120, 210, 252, 210, 120, 45, 10, 1, 11, 55, 165, 330, 462, 462, 330, 165, 55, 11, 1, 12, 66, 220, 495, 792, 924, 792
Offset: 0

Views

Author

Zerinvary Lajos, Dec 02 2007

Keywords

Comments

T(n,m) is the number of m-faces of a regular n-simplex.
An n-simplex is the n-dimensional analog of a triangle. Specifically, a simplex is the convex hull of a set of (n + 1) affinely independent points in some Euclidean space of dimension n or higher, i.e., a set of points such that no m-plane contains more than (m + 1) of them. Such points are said to be in general position.
Reversing the rows gives A074909, which as a linear sequence is essentially the same as this.
From Tom Copeland, Dec 07 2007: (Start)
T(n,k) * (k+1)! = A068424. The comment on permuted words in A068424 shows that T is related to combinations of letters defined by connectivity of regular polytope simplexes.
If T is the diagonally-shifted Pascal matrix, binomial(n+m, k+m), for m=1, then T is a fundamental type of matrix that is discussed in A133314 and the following hold.
The infinitesimal matrix generator is given by A132681, so T = LM(1) of A132681 with inverse LM(-1).
With a(k) = (-x)^k / k!, T * a = [ Laguerre(n,x,1) ], a vector array with index n for the Laguerre polynomials of order 1. Other formulas for the action of T are given in A132681.
T(n,k) = (1/n!) (D_x)^n (D_t)^k Gf(x,t) evaluated at x=t=0 with Gf(x,t) = exp[ t * x/(1-x) ] / (1-x)^2.
[O.g.f. for T ] = 1 / { [ 1 - t * x/(1-x) ] * (1-x)^2 }. [ O.g.f. for row sums ] = 1 / { (1-x) * (1-2x) }, giving A000225 (without a leading zero) for the row sums. Alternating sign row sums are all 1. [Sign correction noted by Vincent J. Matsko, Jul 19 2015]
O.g.f. for row polynomials = [ (1+q)**(n+1) - 1 ] / [ (1+q) -1 ] = A(1,n+1,q) on page 15 of reference on Grassmann cells in A008292. (End)
Given matrices A and B with A(n,k) = T(n,k)*a(n-k) and B(n,k) = T(n,k)*b(n-k), then A*B = C where C(n,k) = T(n,k)*[a(.)+b(.)]^(n-k), umbrally. The e.g.f. for the row polynomials of A is {(a+t) exp[(a+t)x] - a exp(a x)}/t, umbrally. - Tom Copeland, Aug 21 2008
A007318*A097806 as infinite lower triangular matrices. - Philippe Deléham, Feb 08 2009
Riordan array (1/(1-x)^2, x/(1-x)). - Philippe Deléham, Feb 22 2012
The elements of the matrix inverse are T^(-1)(n,k)=(-1)^(n+k)*T(n,k). - R. J. Mathar, Mar 12 2013
Relation to K-theory: T acting on the column vector (-0,d,-d^2,d^3,...) generates the Euler classes for a hypersurface of degree d in CP^n. Cf. Dugger p. 168 and also A104712, A111492, and A238363. - Tom Copeland, Apr 11 2014
Number of walks of length p>0 between any two distinct vertices of the complete graph K_(n+2) is W(n+2,p)=(-1)^(p-1)*Sum_{k=0..p-1} T(p-1,k)*(-n-2)^k = ((n+1)^p - (-1)^p)/(n+2) = (-1)^(p-1)*Sum_{k=0..p-1} (-n-1)^k. This is equal to (-1)^(p-1)*Phi(p,-n-1), where Phi is the cyclotomic polynomial when p is an odd prime. For K_3, see A001045; for K_4, A015518; for K_5, A015521; for K_6, A015531; for K_7, A015540. - Tom Copeland, Apr 14 2014
Consider the transformation 1 + x + x^2 + x^3 + ... + x^n = A_0*(x-1)^0 + A_1*(x-1)^1 + A_2*(x-1)^2 + ... + A_n*(x-1)^n. This sequence gives A_0, ..., A_n as the entries in the n-th row of this triangle, starting at n = 0. - Derek Orr, Oct 14 2014
See A074909 for associations among this array, the Bernoulli polynomials and their umbral compositional inverses, and the face polynomials of permutahedra and their duals (cf. A019538). - Tom Copeland, Nov 14 2014
From Wolfdieter Lang, Dec 10 2015: (Start)
A(r, n) = T(n+r-2, r-1) = risefac(n,r)/r! = binomial(n+r-1, r), for n >= 1 and r >= 1, gives the array with the number of independent components of a symmetric tensors of rank r (number of indices) and dimension n (indices run from 1 to n). Here risefac(n, k) is the rising factorial.
As(r, n) = T(n+1, r+1) = fallfac(n, r)/r! = binomial(n, r), r >= 1 and n >= 1 (with the triangle entries T(n, k) = 0 for n < k) gives the array with the number of independent components of an antisymmetric tensor of rank r and dimension n. Here fallfac is the falling factorial. (End)
The h-vectors associated to these f-vectors are given by A000012 regarded as a lower triangular matrix. Read as bivariate polynomials, the h-polynomials are the complete homogeneous symmetric polynomials in two variables, found in the compositional inverse of an e.g.f. for A008292, the h-vectors of the permutahedra. - Tom Copeland, Jan 10 2017
For a correlation between the states of a quantum system and the combinatorics of the n-simplex, see Boya and Dixit. - Tom Copeland, Jul 24 2017

Examples

			The triangle T(n, k) begins:
   n\k  0  1   2   3   4   5   6   7   8  9 10 11 ...
   0:   1
   1:   2  1
   2:   3  3   1
   3:   4  6   4   1
   4:   5 10  10   5   1
   5:   6 15  20  15   6   1
   6:   7 21  35  35  21   7   1
   7:   8 28  56  70  56  28   8   1
   8:   9 36  84 126 126  84  36   9   1
   9:  10 45 120 210 252 210 120  45  10  1
  10:  11 55 165 330 462 462 330 165  55 11  1
  11:  12 66 220 495 792 924 792 495 220 66 12  1
  ... reformatted by _Wolfdieter Lang_, Mar 23 2015
Production matrix begins
   2   1
  -1   1   1
   1   0   1   1
  -1   0   0   1   1
   1   0   0   0   1   1
  -1   0   0   0   0   1   1
   1   0   0   0   0   0   1   1
  -1   0   0   0   0   0   0   1   1
   1   0   0   0   0   0   0   0   1   1
- _Philippe Deléham_, Jan 29 2014
From _Wolfdieter Lang_, Nov 08 2018: (Start)
Recurrence [_Philippe Deléham_]: T(7, 3) = 2*35 + 35 - 15 - 20 = 70.
Recurrence from Riordan A- and Z-sequences: [1,1,repeat(0)] and [2, repeat(-1, +1)]: From Z: T(5, 0) = 2*5 - 10 + 10 - 5 + 1 = 6. From A: T(7, 3) = 35 + 35 = 70.
Boas-Buck column k=3 recurrence: T(7, 3) = (5/4)*(1 + 5 + 15 + 35) = 70. (End)
		

Crossrefs

Programs

  • Maple
    for i from 0 to 12 do seq(binomial(i, j)*1^(i-j), j = 1 .. i) od;
  • Mathematica
    Flatten[Table[CoefficientList[D[1/x ((x + 1) Exp[(x + 1) z] - Exp[z]), {z, k}] /. z -> 0, x], {k, 0, 11}]]
    CoefficientList[CoefficientList[Series[1/((1 - x)*(1 - x - x*y)), {x, 0, 10}, {y, 0, 10}], x], y] // Flatten (* G. C. Greubel, Nov 22 2017 *)
  • PARI
    for(n=0, 20, for(k=0, n, print1(1/k!*sum(i=0, n, (prod(j=0, k-1, i-j))), ", "))) \\ Derek Orr, Oct 14 2014
    
  • Sage
    Trow = lambda n: sum((x+1)^j for j in (0..n)).list()
    for n in (0..10): print(Trow(n)) # Peter Luschny, Jul 09 2019

Formula

T(n, k) = Sum_{j=k..n} binomial(j,k) = binomial(n+1, k+1), n >= k >= 0, else 0. (Partial sum of column k of A007318 (Pascal), or summation on the upper binomial index (Graham et al. (GKP), eq. (5.10). For the GKP reference see A007318.) - Wolfdieter Lang, Aug 22 2012
E.g.f.: 1/x*((1 + x)*exp(t*(1 + x)) - exp(t)) = 1 + (2 + x)*t + (3 + 3*x + x^2)*t^2/2! + .... The infinitesimal generator for this triangle has the sequence [2,3,4,...] on the main subdiagonal and 0's elsewhere. - Peter Bala, Jul 16 2013
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k) - T(n-2,k-1), T(0,0)=1, T(1,0)=2, T(1,1)=1, T(n,k)=0 if k<0 or if k>n. - Philippe Deléham, Dec 27 2013
T(n,k) = A193862(n,k)/2^k. - Philippe Deléham, Jan 29 2014
G.f.: 1/((1-x)*(1-x-x*y)). - Philippe Deléham, Mar 13 2014
From Tom Copeland, Mar 26 2014: (Start)
[From Copeland's 2007 and 2008 comments]
A) O.g.f.: 1 / { [ 1 - t * x/(1-x) ] * (1-x)^2 } (same as Deleham's).
B) The infinitesimal generator for T is given in A132681 with m=1 (same as Bala's), which makes connections to the ubiquitous associated Laguerre polynomials of integer orders, for this case the Laguerre polynomials of order one L(n,-t,1).
C) O.g.f. of row e.g.f.s: Sum_{n>=0} L(n,-t,1) x^n = exp[t*x/(1-x)]/(1-x)^2 = 1 + (2+t)x + (3+3*t+t^2/2!)x^2 + (4+6*t+4*t^2/2!+t^3/3!)x^3+ ... .
D) E.g.f. of row o.g.f.s: ((1+t)*exp((1+t)*x)-exp(x))/t (same as Bala's).
E) E.g.f. for T(n,k)*a(n-k): {(a+t) exp[(a+t)x] - a exp(a x)}/t, umbrally. For example, for a(k)=2^k, the e.g.f. for the row o.g.f.s is {(2+t) exp[(2+t)x] - 2 exp(2x)}/t.
(End)
From Tom Copeland, Apr 28 2014: (Start)
With different indexing
A) O.g.f. by row: [(1+t)^n-1]/t.
B) O.g.f. of row o.g.f.s: {1/[1-(1+t)*x] - 1/(1-x)}/t.
C) E.g.f. of row o.g.f.s: {exp[(1+t)*x]-exp(x)}/t.
These generating functions are related to row e.g.f.s of A111492. (End)
From Tom Copeland, Sep 17 2014: (Start)
A) U(x,s,t)= x^2/[(1-t*x)(1-(s+t)x)] = Sum_{n >= 0} F(n,s,t)x^(n+2) is a generating function for bivariate row polynomials of T, e.g., F(2,s,t)= s^2 + 3s*t + 3t^2 (Buchstaber, 2008).
B) dU/dt=x^2 dU/dx with U(x,s,0)= x^2/(1-s*x) (Buchstaber, 2008).
C) U(x,s,t) = exp(t*x^2*d/dx)U(x,s,0) = U(x/(1-t*x),s,0).
D) U(x,s,t) = Sum[n >= 0, (t*x)^n L(n,-:xD:,-1)] U(x,s,0), where (:xD:)^k=x^k*(d/dx)^k and L(n,x,-1) are the Laguerre polynomials of order -1, related to normalized Lah numbers. (End)
E.g.f. satisfies the differential equation d/dt(e.g.f.(x,t)) = (x+1)*e.g.f.(x,t) + exp(t). - Vincent J. Matsko, Jul 18 2015
The e.g.f. of the Norlund generalized Bernoulli (Appell) polynomials of order m, NB(n,x;m), is given by exponentiation of the e.g.f. of the Bernoulli numbers, i.e., multiple binomial self-convolutions of the Bernoulli numbers, through the e.g.f. exp[NB(.,x;m)t] = (t/(e^t - 1))^(m+1) * e^(xt). Norlund gave the relation to the factorials (x-1)!/(x-1-n)! = (x-1) ... (x-n) = NB(n,x;n), so T(n,m) = NB(m+1,n+2;m+1)/(m+1)!. - Tom Copeland, Oct 01 2015
From Wolfdieter Lang, Nov 08 2018: (Start)
Recurrences from the A- and Z- sequences for the Riordan triangle (see the W. Lang link under A006232 with references), which are A(n) = A019590(n+1), [1, 1, repeat (0)] and Z(n) = (-1)^(n+1)*A054977(n), [2, repeat(-1, 1)]:
T(0, 0) = 1, T(n, k) = 0 for n < k, and T(n, 0) = Sum_{j=0..n-1} Z(j)*T(n-1, j), for n >= 1, and T(n, k) = T(n-1, k-1) + T(n-1, k), for n >= m >= 1.
Boas-Buck recurrence for columns (see the Aug 10 2017 remark in A036521 also for references):
T(n, k) = ((2 + k)/(n - k))*Sum_{j=k..n-1} T(j, k), for n >= 1, k = 0, 1, ..., n-1, and input T(n, n) = 1, for n >= 0, (the BB-sequences are alpha(n) = 2 and beta(n) = 1). (End)
T(n, k) = [x^k] Sum_{j=0..n} (x+1)^j. - Peter Luschny, Jul 09 2019

Extensions

Edited by Tom Copeland and N. J. A. Sloane, Dec 11 2007

A110555 Triangle of partial sums of alternating binomial coefficients: T(n, k) = Sum_{j=0..k} binomial(n, j)*(-1)^j, for n >= 0, 0 <= k <= n.

Original entry on oeis.org

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

Views

Author

Reinhard Zumkeller, Jul 27 2005

Keywords

Examples

			Triangle T(n, k) starts:
  [0] 1;
  [1] 1,  0;
  [2] 1, -1,  0;
  [3] 1, -2,  1,   0;
  [4] 1, -3,  3,  -1,  0;
  [5] 1, -4,  6,  -4,  1,   0;
  [6] 1, -5, 10, -10,  5,  -1,  0;
  [7] 1, -6, 15, -20, 15,  -6,  1,  0;
  [8] 1, -7, 21, -35, 35, -21,  7, -1,  0.
		

Crossrefs

T(n,1) = -n + 1 for n>0;
T(n,2) = A000217(n-2) for n > 1;
T(n,3) = -A000292(n-4) for n > 2;
T(n,4) = A000332(n-1) for n > 3;
T(n,5) = -A000389(n-1) for n > 5;
T(n,6) = A000579(n-1) for n > 6;
T(n,7) = -A000580(n-1) for n > 7;
T(n,8) = A000581(n-1) for n > 8;
T(n,9) = -A000582(n-1) for n > 9;
T(n,10) = A001287(n-1) for n > 10;
T(n,11) = -A001288(n-1) for n > 11;
T(n,12) = A010965(n-1) for n > 12;
T(n,13) = -A010966(n-1) for n > 13;
T(n,14) = A010967(n-1) for n > 14;
T(n,15) = -A010968(n-1) for n > 15;
T(n,16) = A010969(n-1) for n > 16.
Cf. A071919 (variant), A000007 (row sums), A110556 (central terms).

Programs

  • Maple
    T := (n, k) -> (-1)^k * binomial(n-1, k):
    seq(print(seq(T(n, k), k = 0..n)), n = 0..7); # Peter Luschny, Apr 13 2023
  • Mathematica
    T[0, 0] := 1;  T[n_, n_] := 0; T[n_, k_] := (-1)^k*Binomial[n - 1, k]; Table[T[n, k], {n, 0, 20}, {k, 0, n}] // Flatten (* G. C. Greubel, Aug 31 2017 *)
  • PARI
    concat(1, for(n=1,10, for(k=0,n, print1(if(k != n, (-1)^k*binomial(n-1,k), 0), ", ")))) \\ G. C. Greubel, Aug 31 2017

Formula

T(n, 0) = 1, T(n, n) = 0^n, T(n, k) = -T(n-1, k-1) + T(n-1, k), for 0 < k < n.
T(n, k) = binomial(n-1, k)*(-1)^k, 0 <= k < n, T(n, n) = 0^n.
T(n, n-k-1) = -T(n, k), for 0 < k < n.
T(n, k) = A071919(n, k)*(-1)^k and A071919(n, k) = abs(T(n, k)).
Triangle T(n,k), 0 <= k <= n, read by rows, given by [1, 0, 0, 0, 0, 0, 0, 0, ...] DELTA [0, -1, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Sep 05 2005
G.f.: (1 + x*y) / (1 + x*y - x). - R. J. Mathar, Aug 11 2015

Extensions

Typo in name corrected by Andrey Zabolotskiy, Feb 22 2022
Offset corrected by Peter Luschny, Apr 13 2023

A221857 Number A(n,k) of shapes of balanced k-ary trees with n nodes, where a tree is balanced if the total number of nodes in subtrees corresponding to the branches of any node differ by at most one; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 2, 1, 0, 1, 1, 3, 1, 1, 0, 1, 1, 4, 3, 4, 1, 0, 1, 1, 5, 6, 1, 4, 1, 0, 1, 1, 6, 10, 4, 9, 4, 1, 0, 1, 1, 7, 15, 10, 1, 27, 1, 1, 0, 1, 1, 8, 21, 20, 5, 16, 27, 8, 1, 0, 1, 1, 9, 28, 35, 15, 1, 96, 81, 16, 1, 0, 1, 1, 10, 36, 56, 35, 6, 25, 256, 81, 32, 1, 0
Offset: 0

Views

Author

Alois P. Heinz, Apr 10 2013

Keywords

Examples

			: A(2,2) = 2  : A(2,3) = 3      : A(3,3) = 3          :
:   o     o   :   o    o    o   :   o      o      o   :
:  / \   / \  :  /|\  /|\  /|\  :  /|\    /|\    /|\  :
: o         o : o      o      o : o o    o   o    o o :
:.............:.................:.....................:
: A(3,4) = 6                                          :
:    o        o        o        o       o        o    :
:  /( )\    /( )\    /( )\    /( )\   /( )\    /( )\  :
: o o      o   o    o     o    o o     o   o      o o :
Square array A(n,k) begins:
  1, 1, 1,  1,   1,   1,  1,  1,  1,   1,   1, ...
  1, 1, 1,  1,   1,   1,  1,  1,  1,   1,   1, ...
  0, 1, 2,  3,   4,   5,  6,  7,  8,   9,  10, ...
  0, 1, 1,  3,   6,  10, 15, 21, 28,  36,  45, ...
  0, 1, 4,  1,   4,  10, 20, 35, 56,  84, 120, ...
  0, 1, 4,  9,   1,   5, 15, 35, 70, 126, 210, ...
  0, 1, 4, 27,  16,   1,  6, 21, 56, 126, 252, ...
  0, 1, 1, 27,  96,  25,  1,  7, 28,  84, 210, ...
  0, 1, 8, 81, 256, 250, 36,  1,  8,  36, 120, ...
		

Crossrefs

Rows n=0+1, 2-3, give: A000012, A001477, A179865.
Diagonal and upper diagonals give: A028310, A000217, A000292, A000332, A000389, A000579, A000580, A000581, A000582, A001287, A001288.
Lower diagonals give: A000012, A000290, A092364(n) for n>1.

Programs

  • Maple
    A:= proc(n, k) option remember; local m, r; if n<2 or k=1 then 1
          elif k=0 then 0 else r:= iquo(n-1, k, 'm');
          binomial(k, m)*A(r+1, k)^m*A(r, k)^(k-m) fi
        end:
    seq(seq(A(n, d-n), n=0..d), d=0..12);
  • Mathematica
    a[n_, k_] := a[n, k] = Module[{m, r}, If[n < 2 || k == 1, 1, If[k == 0, 0, {r, m} = QuotientRemainder[n-1, k]; Binomial[k, m]*a[r+1, k]^m*a[r, k]^(k-m)]]]; Table[a[n, d-n], {d, 0, 12}, {n, 0, d}] // Flatten (* Jean-François Alcover, Apr 17 2013, translated from Maple *)

A254142 a(n) = (9*n+10)*binomial(n+9,9)/10.

Original entry on oeis.org

1, 19, 154, 814, 3289, 11011, 32032, 83512, 199342, 442442, 923780, 1830764, 3468374, 6317234, 11113784, 18958808, 31461815, 50930165, 80613390, 125014890, 190285095, 284712285, 419329560, 608658960, 871616460, 1232604516, 1722822024, 2381824984
Offset: 0

Views

Author

Bruno Berselli, Jan 26 2015

Keywords

Comments

Partial sums of A056003.
If n is of the form 8*k+2*(-1)^k-1 or 8*k+2*(-1)^k-2 then a(n) is odd.

Crossrefs

Cf. sequences of the type (k*n+k+1)*binomial(n+k,k)/(k+1): A000217 (k=1), A000330 (k=2), A001296 (k=3), A034263 (k=4), A051946 (k=5), A034265 (k=6), A034266 (k=7), A056122 (k=8), this sequence (k=9).

Programs

  • GAP
    List([0..30], n-> (9*n+10)*Binomial(n+9,9)/10); # G. C. Greubel, Aug 28 2019
  • Magma
    [(9*n+10)*Binomial(n+9,9)/10: n in [0..30]];
    
  • Maple
    seq((9*n+10)*binomial(n+9,9)/10, n=0..30); # G. C. Greubel, Aug 28 2019
  • Mathematica
    Table[(9n+10)Binomial[n+9, 9]/10, {n, 0, 30}]
  • PARI
    vector(30, n, n--; (9*n+10)*binomial(n+9, 9)/10)
    
  • Sage
    [(9*n+10)*binomial(n+9,9)/10 for n in (0..30)]
    

Formula

G.f.: (1 + 8*x)/(1-x)^11.
a(n) = Sum_{i=0..n} (i+1)*A000581(i+8).
a(n+1) = 8*A001287(n+10) + A001287(n+11).

A258993 Triangle read by rows: T(n,k) = binomial(n+k,n-k), k = 0..n-1.

Original entry on oeis.org

1, 1, 3, 1, 6, 5, 1, 10, 15, 7, 1, 15, 35, 28, 9, 1, 21, 70, 84, 45, 11, 1, 28, 126, 210, 165, 66, 13, 1, 36, 210, 462, 495, 286, 91, 15, 1, 45, 330, 924, 1287, 1001, 455, 120, 17, 1, 55, 495, 1716, 3003, 3003, 1820, 680, 153, 19, 1, 66, 715, 3003, 6435, 8008, 6188, 3060, 969, 190, 21
Offset: 1

Views

Author

Reinhard Zumkeller, Jun 22 2015

Keywords

Comments

T(n,k) = A085478(n,k) = A007318(A094727(n),A004736(k)), k = 0..n-1;
rounded(T(n,k)/(2*k+1)) = A258708(n,k);
rounded(sum(T(n,k)/(2*k+1)): k = 0..n-1) = A000967(n).

Examples

			.  n\k |  0  1    2    3     4     5     6     7    8    9  10 11
. -----+-----------------------------------------------------------
.   1  |  1
.   2  |  1  3
.   3  |  1  6    5
.   4  |  1 10   15    7
.   5  |  1 15   35   28     9
.   6  |  1 21   70   84    45    11
.   7  |  1 28  126  210   165    66    13
.   8  |  1 36  210  462   495   286    91    15
.   9  |  1 45  330  924  1287  1001   455   120   17
.  10  |  1 55  495 1716  3003  3003  1820   680  153   19
.  11  |  1 66  715 3003  6435  8008  6188  3060  969  190  21
.  12  |  1 78 1001 5005 12870 19448 18564 11628 4845 1330 231 23  .
		

Crossrefs

If a diagonal of 1's is added on the right, this becomes A085478.
Essentially the same as A143858.
Cf. A027941 (row sums), A117671 (central terms), A143858, A000967, A258708.
T(n,k): A000217 (k=1), A000332 (k=2), A000579 (k=3), A000581 (k=4), A001287 (k=5), A010965 (k=6), A010967 (k=7), A010969 (k=8), A010971 (k=9), A010973 (k=10), A010975 (k=11), A010977 (k=12), A010979 (k=13), A010981 (k=14), A010983 (k=15), A010985 (k=16), A010987 (k=17), A010989 (k=18), A010991 (k=19), A010993 (k=20), A010995 (k=21), A010997 (k=22), A010999 (k=23), A011001 (k=24), A017714 (k=25), A017716 (k=26), A017718 (k=27), A017720 (k=28), A017722 (k=29), A017724 (k=30), A017726 (k=31), A017728 (k=32), A017730 (k=33), A017732 (k=34), A017734 (k=35), A017736 (k=36), A017738 (k=37), A017740 (k=38), A017742 (k=39), A017744 (k=40), A017746 (k=41), A017748 (k=42), A017750 (k=43), A017752 (k=44), A017754 (k=45), A017756 (k=46), A017758 (k=47), A017760 (k=48), A017762 (k=49), A017764 (k=50).
T(n+k,n): A005408 (k=1), A000384 (k=2), A000447 (k=3), A053134 (k=4), A002299 (k=5), A053135 (k=6), A053136 (k=7), A053137 (k=8), A053138 (k=9), A196789 (k=10).
Cf. A165253.

Programs

  • GAP
    Flat(List([1..12], n-> List([0..n-1], k-> Binomial(n+k,n-k) ))); # G. C. Greubel, Aug 01 2019
  • Haskell
    a258993 n k = a258993_tabl !! (n-1) !! k
    a258993_row n = a258993_tabl !! (n-1)
    a258993_tabl = zipWith (zipWith a007318) a094727_tabl a004736_tabl
    
  • Magma
    [Binomial(n+k,n-k): k in [0..n-1], n in [1..12]]; // G. C. Greubel, Aug 01 2019
    
  • Mathematica
    Table[Binomial[n+k,n-k], {n,1,12}, {k,0,n-1}]//Flatten (* G. C. Greubel, Aug 01 2019 *)
  • PARI
    T(n,k) = binomial(n+k,n-k);
    for(n=1, 12, for(k=0,n-1, print1(T(n,k), ", "))) \\ G. C. Greubel, Aug 01 2019
    
  • Sage
    [[binomial(n+k,n-k) for k in (0..n-1)] for n in (1..12)] # G. C. Greubel, Aug 01 2019
    

Formula

T(n,k) = A085478(n,k) = A007318(A094727(n),A004736(k)), k = 0..n-1;
rounded(T(n,k)/(2*k+1)) = A258708(n,k);
rounded(sum(T(n,k)/(2*k+1)): k = 0..n-1) = A000967(n).

A039948 A triangle related to A000045 (Fibonacci numbers).

Original entry on oeis.org

1, 1, 1, 4, 2, 1, 18, 12, 3, 1, 120, 72, 24, 4, 1, 960, 600, 180, 40, 5, 1, 9360, 5760, 1800, 360, 60, 6, 1, 105840, 65520, 20160, 4200, 630, 84, 7, 1, 1370880, 846720, 262080, 53760, 8400, 1008, 112, 8, 1, 19958400, 12337920, 3810240, 786240, 120960, 15120, 1512, 144, 9, 1
Offset: 0

Views

Author

Keywords

Examples

			Triangle begins :
    1;
    1,   1;
    4,   2,   1;
   18,  12,   3,  1;
  120,  72,  24,  4, 1;
  960, 600, 180, 40, 5, 1;
... - _Philippe Deléham_, Nov 08 2011
		

Crossrefs

Programs

  • Magma
    [(Factorial(n)/Factorial(k))*Fibonacci(n-k+1): k in [0..n], n in [0..12]]; // G. C. Greubel, Nov 20 2022
    
  • Mathematica
    T[n_,k_]:= (n!/k!)*Fibonacci[n-k+1];
    Table[T[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Nov 20 2022 *)
  • SageMath
    def A039948(n, k): return factorial(n-k)*binomial(n,k)*fibonacci(n-k+1)
    flatten([[A039948(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Nov 20 2022

Formula

T(n, m) = n!*Fibonacci(n-m+1)/m!, n >= m >= 0.
T(n, 0) = A005442(n).
T(n, 1) = A005443(n).
E.g.f. for column m: x^m/(m!*(1-x-x^2)), m >= 0.
From G. C. Greubel, Nov 20 2022: (Start)
T(n, n-1) = A000027(n).
T(n, n-2) = 4*A000217(n-1), n >= 2.
T(n, n-3) = 18*A000292(n-2), n >= 3.
T(n, n-4) = 5! * A000332(n), n >= 4.
T(n, n-5) = 8 * 5! * A000389(n), n >= 5.
T(n, n-6) = 13 * 6! * A000579(n), n >= 6.
T(n, n-7) = 21 * 7! * A000580(n), n >= 7.
T(n, n-8) = 34 * 8! * A000581(n), n >= 8.
T(n, n-9) = 55 * 9! * A000582(n), n >= 9.
T(n, n-10) = 89 * 10! * A001287(n), n >= 10.
T(n, n-11) = 12 * 12! * A001288(n), n >= 11.
T(n, n-12) = 233 * 12! * A010965(n), n >= 12.
T(n, n-13) = 89 * 13! * A010966(n), n >= 13.
Sum_{k=0..n} T(n, k) = A110313(n). (End)

A101095 Fourth difference of fifth powers (A000584).

Original entry on oeis.org

1, 28, 121, 240, 360, 480, 600, 720, 840, 960, 1080, 1200, 1320, 1440, 1560, 1680, 1800, 1920, 2040, 2160, 2280, 2400, 2520, 2640, 2760, 2880, 3000, 3120, 3240, 3360, 3480, 3600, 3720, 3840, 3960, 4080, 4200, 4320, 4440, 4560, 4680, 4800, 4920, 5040, 5160, 5280
Offset: 1

Views

Author

Cecilia Rossiter, Dec 15 2004

Keywords

Comments

Original Name: Shells (nexus numbers) of shells of shells of shells of the power of 5.
The (Worpitzky/Euler/Pascal Cube) "MagicNKZ" algorithm is: MagicNKZ(n,k,z) = Sum_{j=0..k+1} (-1)^j*binomial(n + 1 - z, j)*(k - j + 1)^n, with k>=0, n>=1, z>=0. MagicNKZ is used to generate the n-th accumulation sequence of the z-th row of the Euler Triangle (A008292). For example, MagicNKZ(3,k,0) is the 3rd row of the Euler Triangle (followed by zeros) and MagicNKZ(10,k,1) is the partial sums of the 10th row of the Euler Triangle. This sequence is MagicNKZ(5,k-1,2).

Crossrefs

Fourth differences of A000584, third differences of A022521, second differences of A101098, and first differences of A101096.
For other sequences based upon MagicNKZ(n,k,z):
...... | n = 1 | n = 2 | n = 3 | n = 4 | n = 5 | n = 6 | n = 7 | n = 8
--------------------------------------------------------------------------------------
z = 0 | A000007 | A019590 | ....... MagicNKZ(n,k,0) = T(n,k+1) from A008292 .......
z = 1 | A000012 | A040000 | A101101 | A101104 | A101100 | ....... | ....... | .......
z = 2 | A000027 | A005408 | A008458 | A101103 | thisSeq | ....... | ....... | .......
z = 3 | A000217 | A000290 | A003215 | A005914 | A101096 | ....... | ....... | .......
z = 4 | A000292 | A000330 | A000578 | A005917 | A101098 | ....... | ....... | .......
z = 5 | A000332 | A002415 | A000537 | A000583 | A022521 | ....... | A255181 | .......
z = 12 | A001288 | A057788 | ....... | A254870 | A254471 | A254683 | A254646 | A254642
z = 13 | A010965 | ....... | ....... | ....... | A254871 | A254472 | A254684 | A254647
z = 14 | A010966 | ....... | ....... | ....... | ....... | A254872 | ....... | .......
--------------------------------------------------------------------------------------
Cf. A047969.

Programs

  • Magma
    I:=[1,28,121,240,360]; [n le 5 select I[n] else 2*Self(n-1)-Self(n-2): n in [1..50]]; // Vincenzo Librandi, May 07 2015
    
  • Mathematica
    MagicNKZ=Sum[(-1)^j*Binomial[n+1-z, j]*(k-j+1)^n, {j, 0, k+1}];Table[MagicNKZ, {n, 5, 5}, {z, 2, 2}, {k, 0, 34}]
    CoefficientList[Series[(1 + 26 x + 66 x^2 + 26 x^3 + x^4)/(1 - x)^2, {x, 0, 50}], x] (* Vincenzo Librandi, May 07 2015 *)
    Join[{1,28,121,240},Differences[Range[50]^5,4]] (* or *) LinearRecurrence[{2,-1},{1,28,121,240,360},50] (* Harvey P. Dale, Jun 11 2016 *)
  • PARI
    a(n)=if(n>3, 120*n-240, 33*n^2-72*n+40) \\ Charles R Greathouse IV, Oct 11 2015
  • Sage
    [1,28,121]+[120*(k-2) for k in range(4,36)] # Danny Rorabaugh, Apr 23 2015
    

Formula

a(k+1) = Sum_{j=0..k+1} (-1)^j*binomial(n + 1 - z, j)*(k - j + 1)^n; n = 5, z = 2.
For k>3, a(k) = Sum_{j=0..4} (-1)^j*binomial(4, j)*(k - j)^5 = 120*(k - 2).
a(n) = 2*a(n-1) - a(n-2), n>5. G.f.: x*(1+26*x+66*x^2+26*x^3+x^4) / (1-x)^2. - Colin Barker, Mar 01 2012

Extensions

MagicNKZ material edited, Crossrefs table added, SeriesAtLevelR material removed by Danny Rorabaugh, Apr 23 2015
Name changed and keyword 'uned' removed by Danny Rorabaugh, May 06 2015

A155856 Triangle T(n,k) = binomial(2*n-k, k)*(n-k)!, read by rows.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 6, 10, 6, 1, 24, 42, 30, 10, 1, 120, 216, 168, 70, 15, 1, 720, 1320, 1080, 504, 140, 21, 1, 5040, 9360, 7920, 3960, 1260, 252, 28, 1, 40320, 75600, 65520, 34320, 11880, 2772, 420, 36, 1, 362880, 685440, 604800, 327600, 120120, 30888, 5544, 660, 45, 1
Offset: 0

Views

Author

Paul Barry, Jan 29 2009

Keywords

Comments

Row sums of B^{-1}*A155856*B^{-1} are A000166 with B=A007318.
Downward diagonals T(n+j, n) = j!*binomial(n+j, n) = j!*seq(j), where seq(j) are sequences A010965, A010967, ..., A011101, A017714, A017716, ..., A017764, for 6 <= j <= 50, respectively. - G. C. Greubel, Jun 04 2021

Examples

			Triangle begins:
     1;
     1,    1;
     2,    3,    1;
     6,   10,    6,    1;
    24,   42,   30,   10,    1;
   120,  216,  168,   70,   15,   1;
   720, 1320, 1080,  504,  140,  21,  1;
  5040, 9360, 7920, 3960, 1260, 252, 28, 1;
		

Crossrefs

Cf. A155857 (row sums), A155858 (diagonal sums).

Programs

  • Mathematica
    Table[Binomial[2n-k,k](n-k)!,{n,0,10},{k,0,n}]//Flatten (* Harvey P. Dale, Mar 24 2017 *)
  • Sage
    flatten([[factorial(n-k)*binomial(2*n-k, k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Jun 04 2021

Formula

T(n,k) = binomial(2*n-k, k)*(n-k)!.
Sum_{k=0..n} T(n, k) = A155857(n)
Sum_{k=0..floor(n/2)} T(n-k, k) = A155858(n) (diagonal sums).
G.f.: 1/(1-xy-x/(1-xy-x/(1-xy-2x/(1-xy-2x/(1-xy-3x/(1-.... (continued fraction).
From G. C. Greubel, Jun 04 2021: (Start)
T(n, 0) = A000142(n). T(n+1, n) = A000217(n+1).
T(n+1, 1) = A007680(n). T(n+2, n) = A034827(n+4).
T(n+2, 2) = A175925(n). T(n+3, n) = A253946(n).
T(2*n, n) = A064352(n) T(n+4, n) = 4!*A000581(n).
T(n+1, n) = A000217(n+1). T(n+5, n) = 5!*A001287(n). (End)

A017764 a(n) = binomial coefficient C(n,100).

Original entry on oeis.org

1, 101, 5151, 176851, 4598126, 96560646, 1705904746, 26075972546, 352025629371, 4263421511271, 46897636623981, 473239787751081, 4416904685676756, 38393094575497956, 312629484400483356, 2396826047070372396, 17376988841260199871, 119594570260437846171
Offset: 100

Views

Author

Keywords

Comments

More generally, the ordinary generating function for the binomial coefficients C(n,k) is x^k/(1 - x)^(k+1). - Ilya Gutkovskiy, Mar 21 2016

Crossrefs

Cf. similar sequences of the binomial coefficients C(n,k): A000012 (k = 0), A001477 (k = 1), A000217 (k = 2), A000292 (k = 3), A000332 (k = 4), A000389 (k = 5), A000579-A000582 (k = 6..9) A001287 (k = 10), A001288 (k = 11), A010965-A011001 (k = 12..48), A017713-A017763 (k = 49..99), this sequence (k = 100).

Programs

Formula

G.f.: x^100/(1 - x)^101. - Ilya Gutkovskiy, Mar 21 2016
E.g.f.: x^100 * exp(x)/(100)!. - G. C. Greubel, Nov 24 2017
From Amiram Eldar, Dec 20 2020: (Start)
Sum_{n>=100} 1/a(n) = 100/99.
Sum_{n>=100} (-1)^n/a(n) = A001787(100)*log(2) - A242091(100)/99! = 63382530011411470074835160268800*log(2) - 1914409165727592211172313915606932788039791776845041612575266508424929 / 43575234518570298227833630584570189723 = 0.9902877001... (End)
Showing 1-10 of 19 results. Next