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

A077445 Numbers k such that (k^2 - 8)/2 is a square.

Original entry on oeis.org

4, 20, 116, 676, 3940, 22964, 133844, 780100, 4546756, 26500436, 154455860, 900234724, 5246952484, 30581480180, 178241928596, 1038870091396, 6054978619780, 35291001627284, 205691031143924, 1198855185236260
Offset: 1

Views

Author

Gregory V. Richardson, Nov 09 2002

Keywords

Comments

The equation "(k^2 - 8)/2 is a square" is a version of the generalized Pell Equation "x^2 - D*y^2 = C".
a(n)^2 - 2*A077444(n) = 8.
From Wolfdieter Lang, Jan 18 2013: (Start)
4*(1-z)/(1-6*z+z^2) = Sum_{n>=0} a(n+1)*z^n is the formal power series for tan(4*x)/tan(x) if one lets
z = (tan(x))^2. For the numerator and denominator of this o.g.f. see A034867 and A034839, respectively. Convergence holds for 0 <= z < 3 - 2*sqrt(2), approximately 0.1715728753. This means for |x| < Pi/8, approximately 0.3926990818.
See also the o.g.f. given by Johannes W. Meijer, Aug 01 2010, in the formula section of A001653 = (this sequence)/4.
(End)
Positive values of x (or y) satisfying x^2 - 6*x*y + y^2 + 64 = 0. - Colin Barker, Feb 13 2014

References

  • A. H. Beiler, "The Pellian." Ch. 22 in Recreations in the Theory of Numbers: The Queen of Mathematics Entertains. Dover, New York, New York, pp. 248-268, 1966.
  • L. E. Dickson, History of the Theory of Numbers, Vol. II, Diophantine Analysis. AMS Chelsea Publishing, Providence, Rhode Island, 1999, pp. 341-400.
  • Peter G. L. Dirichlet, Lectures on Number Theory (History of Mathematics Source Series, V. 16); American Mathematical Society, Providence, Rhode Island, 1999, pp. 139-147.

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[4 (1 - x)/(1 - 6 x + x^2), {x, 0, 40}], x] (* Vincenzo Librandi, Feb 14 2014 *)
  • PARI
    a(n)=if(n<1,0,subst(poltchebi(n)+poltchebi(n-1),x,3))

Formula

a(n) = (((3+2*sqrt(2))^n + (3-2*sqrt(2))^n) + ((3+2*sqrt(2))^(n-1) + (3-2*sqrt(2))^(n-1))) / 2.
a(n) = 6*a(n-1) - a(n-2) = 4*A001653(n).
G.f.: 4*(x-x^2)/(1-6*x+x^2).
With a=3+2*sqrt(2), b=3-2*sqrt(2): a(n) = sqrt(2)*(a^((2*n-1)/2) + b^((2*n-1)/2)). a(n) = sqrt(2*A003499(2*n-1)+4). - Mario Catalani (mario.catalani(AT)unito.it), Mar 24 2003
a(n) = (A003499(n+1) + A003499(n))/2. - Mario Catalani (mario.catalani(AT)unito.it), Mar 31 2003
a(n) = (2 + sqrt(2))*(3 + 2*sqrt(2))^n + (2 - sqrt(2))*(3- 2*sqrt(2))^n. - Antonio Alberto Olivares, Feb 23 2006
a(n) = 2*A075870(n). - Bruno Berselli, Nov 27 2013
G.f.: 2*Q(0)*x*(1-x)/(1-3*x), where Q(k) = 1 + 1/( 1 - x*(8*k-9)/( x*(8*k-1) - 3/Q(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Dec 10 2013

A202064 Triangle T(n,k), read by rows, given by (2, -1/2, 1/2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, 1/2, -1/2, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 2, 0, 3, 1, 0, 4, 4, 0, 0, 5, 10, 1, 0, 0, 6, 20, 6, 0, 0, 0, 7, 35, 21, 1, 0, 0, 0, 8, 56, 56, 8, 0, 0, 0, 0, 9, 84, 126, 36, 1, 0, 0, 0, 0, 10, 120, 252, 120, 10, 0, 0, 0, 0, 0, 11, 165, 462, 330, 55, 1, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Philippe Deléham, Dec 10 2011

Keywords

Comments

Riordan array (x/(1-x)^2, x^2/(1-x)^2).
Mirror image of triangle in A119900.
A203322*A130595 as infinite lower triangular matrices. - Philippe Deléham, Jan 05 2011
From Gus Wiseman, Jul 07 2025: (Start)
Also the number of subsets of {1..n} containing n with k maximal runs (sequences of consecutive elements increasing by 1). For example, row n = 5 counts the following subsets:
{5} {1,5} {1,3,5}
{4,5} {2,5}
{3,4,5} {3,5}
{2,3,4,5} {1,2,5}
{1,2,3,4,5} {1,4,5}
{2,3,5}
{2,4,5}
{1,2,3,5}
{1,2,4,5}
{1,3,4,5}
For anti-runs instead of runs we have A053538.
Without requiring n see A210039, A202023, reverse A098158, A109446.
(End)

Examples

			Triangle begins :
1
2, 0
3, 1, 0
4, 4, 0, 0
5, 10, 1, 0, 0
6, 20, 6, 0, 0, 0
7, 35, 21, 1, 0, 0, 0
8, 56, 56, 8, 0, 0, 0, 0
		

Crossrefs

Cf. A007318, A005314 (antidiagonal sums), A119900, A084938, A130595, A203322.
Column k = 1 is A000027.
Row sums are A000079.
Column k = 2 is A000292.
Without zeros we have A034867.
Last nonzero term in each row appears to be A124625.
A034839 counts subsets by number of maximal runs, for anti-runs A384893.
A116674 counts strict partitions by number of maximal runs, for anti-runs A384905.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],MemberQ[#,n]&&Length[Split[#,#2==#1+1&]]==k&]],{n,12},{k,n}] (* Gus Wiseman, Jul 07 2025 *)

Formula

G.f.: 1/((1-x)^2-y*x^2).
Sum_{k, 0<=k<=n} T(n,k)*x^k = A000027(n+1), A000079(n), A000129(n+1), A002605(n+1), A015518(n+1), A063727(n), A002532(n+1), A083099(n+1), A015519(n+1), A003683(n+1), A002534(n+1), A083102(n), A015520(n+1), A091914(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 10, 11, 12, 13 respectively.
T(n,k) = binomial(n+1,2k+1).
T(n,k) = 2*T(n-1,k) + T(n-2,k-1) - T(n-2,k), T(0,0) = 1, T(1,0) = 2, T(1,1) = 0 and T(n,k) = 0 if k<0 or if k>n. - Philippe Deléham, Mar 15 2012

A105070 T(n,k) = 2^k*binomial(n,2k+1), where 0 <= k <= floor((n-1)/2), n >= 1.

Original entry on oeis.org

1, 2, 3, 2, 4, 8, 5, 20, 4, 6, 40, 24, 7, 70, 84, 8, 8, 112, 224, 64, 9, 168, 504, 288, 16, 10, 240, 1008, 960, 160, 11, 330, 1848, 2640, 880, 32, 12, 440, 3168, 6336, 3520, 384, 13, 572, 5148, 13728, 11440, 2496, 64, 14, 728, 8008, 27456, 32032, 11648, 896, 15, 910, 12012, 51480, 80080, 43680, 6720, 128
Offset: 1

Views

Author

Emeric Deutsch, Apr 05 2005

Keywords

Comments

Row n contains ceiling(n/2) terms. Row sums yield the Pell numbers (A000129). Column 1 yields A007290.
Eigenvector equals A118397, so that A118397(n) = Sum_{k=0..[n/2]} T(n+1,k)*A118397(k) for n >= 0. - Paul D. Hanna, May 08 2006
Essentially a triangle, read by rows, given by (2, -1/2, 1/2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 12 2011
Subtriangle of the triangle given by (1, 1, -1, 1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, 0, 2, -2, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Apr 07 2012

Examples

			Triangle begins:
  1;
  2;
  3,  2;
  4,  8;
  5, 20,  4;
  6, 40, 24;
(2, -1/2, 1/2, 0, 0, ...) DELTA (0, 1, -1, 0, 0, ...) begins:
  1;
  2,  0;
  3,  2,  0;
  4,  8,  0,  0;
  5, 20,  4,  0,  0;
  6, 40, 24,  0,  0,  0.
(1, 1, -1, 1, 0, 0, ...) DELTA (0, 0, 2, -2, 0, 0, ...) begins:
  1;
  1,  0;
  2,  0,  0;
  3,  2,  0,  0;
  4,  8,  0,  0,  0;
  5, 20,  4,  0,  0,  0;
  6, 40, 24,  0,  0,  0,  0. - _Philippe Deléham_, Apr 07 2012
		

Crossrefs

Cf. A118397 (eigenvector).

Programs

  • Magma
    [2^k*Binomial(n,2*k+1): k in [0..Floor((n-1)/2)], n in [1..15]]; // G. C. Greubel, Mar 15 2020
    
  • Maple
    T:=(n,k)->binomial(n,2*k+1)*2^k:for n from 1 to 15 do seq(T(n,k),k=0..floor((n-1)/2)) od; # yields sequence in triangular form
  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := u[n - 1, x] + 2 x*v[n - 1, x]
    v[n_, x_] := u[n - 1, x] + v[n - 1, x]
    Table[Factor[u[n, x]], {n, 1, z}]
    Table[Factor[v[n, x]], {n, 1, z}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%]  (* A207536 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]  (* A105070 *)
    (* Clark Kimberling, Feb 18 2010 *)
    Table[2^k*Binomial[n, 2*k+1], {n, 15}, {k,0,Floor[(n-1)/2]}]//Flatten (* G. C. Greubel, Mar 15 2020 *)
  • Sage
    [[2^k*binomial(n,2*k+1) for k in (0..floor((n-1)/2))] for n in (1..15)] # G. C. Greubel, Mar 15 2020

Formula

E.g.f.: exp(x)*sinh(x*sqrt(2*y))/sqrt(2*y), cf. A034867. - Vladeta Jovovic, Apr 06 2005
From Philippe Deléham, Apr 07 2012: (Start)
As DELTA-triangle T(n,k) with 0 <= k <= n:
G.f.: (1-x+x^2-y*x^2)/(1-2*x+x^2-2*y*x^2).
T(n,k) = 2*T(n-1,k) - T(n-2,k) + 2*T(n-2,k-1), T(0,0) = T(1,0) = 1, T(1,1) = T(2,1) = T(2,2) = 0, T(2,0) = 2 and T(n,k) = 0 if k<0 or if k>n. (End)
Sum_{k=0..floor((n-1)/2)} T(n,k) = { P(n) (A000129(n)), A215928(n), (-1)^(n-1) *A077985(n-1), -A176981(n+1), (-1)^(n-1)*A215936(n+2) }, for n >= 1. - G. C. Greubel, Mar 15 2020

A162590 Polynomials with e.g.f. exp(x*t)/csch(t), triangle of coefficients read by rows.

Original entry on oeis.org

0, 1, 0, 0, 2, 0, 1, 0, 3, 0, 0, 4, 0, 4, 0, 1, 0, 10, 0, 5, 0, 0, 6, 0, 20, 0, 6, 0, 1, 0, 21, 0, 35, 0, 7, 0, 0, 8, 0, 56, 0, 56, 0, 8, 0, 1, 0, 36, 0, 126, 0, 84, 0, 9, 0, 0, 10, 0, 120, 0, 252, 0, 120, 0, 10, 0, 1, 0, 55, 0, 330, 0, 462, 0, 165, 0, 11, 0, 0, 12, 0, 220, 0, 792, 0, 792, 0
Offset: 0

Views

Author

Peter Luschny, Jul 07 2009

Keywords

Comments

Comment from Peter Bala (Dec 06 2011): "Let P denote Pascal's triangle A070318 and put M = 1/2*(P-P^-1). M is A162590 (see also A131047). Then the first column of (I-t*M)^-1 (apart from the initial 1) lists the row polynomials for" A196776(n,k), which gives the number of ordered partitions of an n set into k odd-sized blocks. - Peter Luschny, Dec 06 2011
The n-th row of the triangle is formed by multiplying by 2^(n-1) the elements of the first row of the limit as k approaches infinity of the stochastic matrix P^(2k-1) where P is the stochastic matrix associated with the Ehrenfest model with n balls. The elements of a stochastic matrix P give the probability of arriving in a state j given the previous state i. In particular the sum of every row of the matrix must be 1, and so the sum of the terms in the n-th row of this triangle is 2^(n-1). Furthermore, by the properties of Markov chains, we can interpret P^(2k) as the (2k)-step transition matrix of the Ehrenfest model and its limit exists and it is again a stochastic matrix. The rows of the triangle divided by 2^(n-1) are the even rows (second, fourth, ...) and the odd rows (first, third, ...) of the limit matrix P^(2k). - Luca Onnis, Oct 29 2023

Examples

			Triangle begins:
  0
  1,  0
  0,  2,  0
  1,  0,  3,  0
  0,  4,  0,  4,  0
  1,  0, 10,  0,  5,  0
  0,  6,  0, 20,  0,  6,  0
  1,  0, 21,  0, 35,  0,  7,  0
  ...
  p[0](x) = 0;
  p[1](x) = 1
  p[2](x) = 2*x
  p[3](x) = 3*x^2 +  1
  p[4](x) = 4*x^3 +  4*x
  p[5](x) = 5*x^4 + 10*x^2 +  1
  p[6](x) = 6*x^5 + 20*x^3 +  6*x
  p[7](x) = 7*x^6 + 35*x^4 + 21*x^2 + 1
  p[8](x) = 8*x^7 + 56*x^5 + 56*x^3 + 8*x
.
Cf. the triangle of odd-numbered terms in rows of Pascal's triangle (A034867).
p[n] (k), n=0,1,...
k=0:  0, 1,  0,   1,    0,     1, ... A000035, (A059841)
k=1:  0, 1,  2,   4,    8,    16, ... A131577, (A000079)
k=2:  0, 1,  4,  13,   40,   121, ... A003462
k=3:  0, 1,  6,  28,  120,   496, ... A006516
k=4:  0, 1,  8,  49,  272,  1441, ... A005059
k=5:  0, 1, 10,  76,  520,  3376, ... A081199, (A016149)
k=6:  0, 1, 12, 109,  888,  6841, ... A081200, (A016161)
k=7:  0, 1, 14, 148, 1400, 12496, ... A081201, (A016170)
k=8:  0, 1, 16, 193, 2080, 21121, ... A081202, (A016178)
k=9:  0, 1, 18, 244, 2952, 33616, ... A081203, (A016186)
k=10: 0, 1, 20, 301, 4040, 51001, ... ......., (A016190)
.
p[n] (k), k=0,1,...
p[0]: 0,  0,   0,    0,    0,     0, ... A000004
p[1]: 1,  1,   1,    1,    1,     1, ... A000012
p[2]: 0,  2,   4,    6,    8,    10, ... A005843
p[3]: 1,  4,  13,   28,   49,    76, ... A056107
p[4]: 0,  8,  40,  120,  272,   520, ... A105374
p[5]: 1, 16, 121,  496, 1441,  3376, ...
p[6]: 0, 32, 364, 2016, 7448, 21280, ...
		

Crossrefs

Cf. A119467.

Programs

  • Maple
    # Polynomials: p_n(x)
    p := proc(n,x) local k;
    pow := (n,k) -> `if`(n=0 and k=0,1,n^k);
    add((k mod 2)*binomial(n,k)*pow(x,n-k),k=0..n) end;
    # Coefficients: a(n)
    seq(print(seq(coeff(i!*coeff(series(exp(x*t)/csch(t), t,16),t,i),x,n), n=0..i)), i=0..8);
  • Mathematica
    p[n_, x_] := Sum[Binomial[n, 2*k-1]*x^(n-2*k+1), {k, 0, n+2}]; row[n_] := CoefficientList[p[n, x], x] // Append[#, 0]&; Table[row[n], {n, 0, 12}] // Flatten (* Jean-François Alcover, Jun 28 2013 *)
    n = 15; "n-th row"
    mat = Table[Table[0, {j, 1, n + 1}], {i, 1, n + 1}];
    mat[[1, 2]] = 1;
    mat[[n + 1, n]] = 1;
    For[i = 2, i <= n, i++, mat[[i, i - 1]] = (i - 1)/n ];
    For[i = 2, i <= n, i++, mat[[i, i + 1]] = (n - i + 1)/n];
    mat // MatrixForm;
    P2 = Dot[mat, mat];
    R1 = Simplify[
      Eigenvectors[Transpose[P2]][[1]]/
       Total[Eigenvectors[Transpose[P2]][[1]]]]
    R2 = Table[Dot[R1, Transpose[mat][[k]]], {k, 1, n + 1}]
    even = R1*2^(n - 1) (* Luca Onnis, Oct 29 2023 *)

Formula

p_n(x) = Sum_{k=0..n} (k mod 2)*binomial(n,k)*x^(n-k).
E.g.f.: exp(x*t)/csch(t) = 0*(t^0/0!) + 1*(t^1/1!) + (2*x)*(t^2/2!) + (3*x^2+1)*(t^3/3!) + ...
The 'co'-polynomials with generating function exp(x*t)*sech(t) are the Swiss-Knife polynomials (A153641).

A080923 First differences of A003946.

Original entry on oeis.org

1, 3, 8, 24, 72, 216, 648, 1944, 5832, 17496, 52488, 157464, 472392, 1417176, 4251528, 12754584, 38263752, 114791256, 344373768, 1033121304, 3099363912, 9298091736, 27894275208, 83682825624, 251048476872, 753145430616
Offset: 0

Views

Author

Paul Barry, Feb 26 2003

Keywords

Comments

Sum of consecutive pairs of elements of A025192.
The alternating sign sequence with g.f. (1-x^2)/(1+3x) gives the diagonal sums of A110168. - Paul Barry, Jul 14 2005
Let M = an infinite lower triangular matrix with the odd integers (1,3,5,...) in every column, with the leftmost column shifted up one row. Then A080923 = lim_{n->inf} M^n. - Gary W. Adamson, Feb 18 2010
a(n+1), n >= 0, with o.g.f. ((1-x^2)/(1-3*x)-1)/x = (3-x)/(1-3*x) provides the coefficients in the formal power series for tan(3*x)/tan(x) = (3-z)/(1-3*z) = Sum_{n>=0} a(n+1)*z^n, with z = tan(x)^2. Convergence holds for 0 <= z < 1/3, i.e., |x| < Pi/6, approximately 0.5235987758. For the numerator and denominator of this o.g.f. see A034867 and A034839, respectively. - Wolfdieter Lang, Jan 18 2013

Crossrefs

Essentially the same as A005051, A026097 and A083583.

Programs

  • Mathematica
    CoefficientList[Series[(1 - x^2) / (1 - 3 x), {x, 0, 20}], x] (* Vincenzo Librandi, Aug 05 2013 *)

Formula

G.f.: (1-x^2)/(1-3*x).
G.f.: 1/(1 - 3*x + x^2 - 3*x^3 + x^4 - 3*x^5 + ...). - Gary W. Adamson, Jan 06 2011
a(n) = 2^3*3^(n-2), n >= 2, a(0) = 1, a(1) = 3. - Wolfdieter Lang, Jan 18 2013

A202023 Triangle T(n,k), read by rows, given by (1, 0, 1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 3, 0, 0, 1, 6, 1, 0, 0, 1, 10, 5, 0, 0, 0, 1, 15, 15, 1, 0, 0, 0, 1, 21, 35, 7, 0, 0, 0, 0, 1, 28, 70, 28, 1, 0, 0, 0, 0, 1, 36, 126, 84, 9, 0, 0, 0, 0, 0, 1, 45, 210, 210, 45, 1, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Philippe Deléham, Dec 10 2011

Keywords

Comments

Riordan array (1/(1-x), x^2/(1-x)^2).
A skewed version of triangular array A085478.
Mirror image of triangle in A098158.
Sum_{k, 0<=k<=n} T(n,k)*x^k = A138229(n), A006495(n), A138230(n),A087455(n), A146559(n), A000012(n), A011782(n), A001333(n),A026150(n), A046717(n), A084057(n), A002533(n), A083098(n),A084058(n), A003665(n), A002535(n), A133294(n), A090042(n),A125816(n), A133343(n), A133345(n), A120612(n), A133356(n), A125818(n) for x = -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18 respectively.
Sum_{k, 0<=k<=n} T(n,k)*x^(n-k) = A009116(n), A000007(n), A011782(n), A006012(n), A083881(n), A081335(n), A090139(n), A145301(n), A145302(n), A145303(n), A143079(n) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively.
From Gus Wiseman, Jul 08 2025: (Start)
After the first row this is also the number of subsets of {1..n-1} with k maximal runs (sequences of consecutive elements increasing by 1) for k = 0..n. For example, row n = 5 counts the following subsets:
{} {1} {1,3} . . .
{2} {1,4}
{3} {2,4}
{4} {1,2,4}
{1,2} {1,3,4}
{2,3}
{3,4}
{1,2,3}
{2,3,4}
{1,2,3,4}
Requiring n-1 gives A202064.
For anti-runs instead of runs we have A384893.
(End)

Examples

			Triangle begins :
1
1, 0
1, 1, 0
1, 3, 0, 0
1, 6, 1, 0, 0
1, 10, 5, 0, 0, 0
1, 15, 15, 1, 0, 0, 0
1, 21, 35, 7, 0, 0, 0, 0
1, 28, 70, 28, 1, 0, 0, 0, 0
		

Crossrefs

Column k = 1 is A000217.
Column k = 2 is A000332.
Row sums are A011782 (or A000079 shifted right).
Removing all zeros gives A034839 (requiring n-1 A034867).
Last nonzero term in each row appears to be A093178, requiring n-1 A124625.
Reversing rows gives A098158, without zeros A109446.
Without the k = 0 column we get A210039.
Row maxima appear to be A214282.
A116674 counts strict partitions by number of maximal runs, for anti-runs A384905.
A268193 counts integer partitions by number of maximal runs, for anti-runs A384881.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n-1]],Length[Split[#,#2==#1+1&]]==k&]],{n,0,10},{k,0,n}] (* Gus Wiseman, Jul 08 2025 *)

Formula

T(n,k) = binomial(n,2k).
G.f.: (1-x)/((1-x)^2-y*x^2).
T(n,k)= Sum_{j, j>=0} T(n-1-j,k-1)*j with T(n,0)=1 and T(n,k)= 0 if k<0 or if n
T(n,k) = 2*T(n-1,k) + T(n-2,k-1) - T(n-2,k) for n>1, T(0,0) = T(1,0) = 1, T(1,1) = 0, T(n,k) = 0 if k>n or if k<0. - Philippe Deléham, Nov 10 2013

A095704 Triangle read by rows giving coefficients of the trigonometric expansion of sin(n*x).

Original entry on oeis.org

1, 2, 0, 3, 0, -1, 4, 0, -4, 0, 5, 0, -10, 0, 1, 6, 0, -20, 0, 6, 0, 7, 0, -35, 0, 21, 0, -1, 8, 0, -56, 0, 56, 0, -8, 0, 9, 0, -84, 0, 126, 0, -36, 0, 1, 10, 0, -120, 0, 252, 0, -120, 0, 10, 0, 11, 0, -165, 0, 462, 0, -330, 0, 55, 0, -1, 12, 0, -220, 0, 792, 0, -792, 0, 220, 0, -12, 0, 13, 0, -286, 0, 1287, 0
Offset: 1

Author

Robert G. Wilson v, Jul 06 2004

Keywords

Examples

			The trigonometric expansion of sin(4x) is 4*cos(x)^3*sin(x) - 4*cos(x)*sin(x)^3, so the fourth row is 4, 0, -4, 0.
Triangle begins:
1
2 0
3 0 -1
4 0 -4 0
5 0 -10 0 1
6 0 -20 0 6 0
7 0 -35 0 21 0 -1
8 0 -56 0 56 0 -8 0
		

Crossrefs

First column is A000027 = C(n, 1), third column is A000292 = C(n, 3), fifth column is A000389 = C(n, 5), seventh column is A000580 = C(n, 7), ninth column is A000582 = C(n, 9).
A001288 = C(n, 11), A010966 = C(n, 13), A010968 = C(n, 15), A010970 = C(n, 17), A010972 = C(n, 19),
A010974 = C(n, 21), A010976 = C(n, 23), A010978 = C(n, 25), A010980 = C(n, 27), A010982 = C(n, 29),
A010984 = C(n, 31), A010986 = C(n, 33), A010988 = C(n, 35), A010990 = C(n, 37), A010992 = C(n, 39),
A010994 = C(n, 41), A010996 = C(n, 43), A010998 = C(n, 45), A011000 = C(n, 47), A017713 = C(n, 49)
Another version of the triangle in A034867. Cf. A096754.
A017715 = C(n, 51), A017717 = C(n, 53), A017719 = C(n, 55), A017721 = C(n, 57), etc.

Programs

  • Mathematica
    Flatten[ Table[ Plus @@ CoefficientList[ TrigExpand[ Sin[n*x]], {Sin[x], Cos[x]}], {n, 13}]]

Formula

T(n,k) = C(n+1,k+1)*sin(Pi*(k+1)/2). - Paul Barry, May 21 2006

A123162 Triangle read by rows: T(n,k) = binomial(2*n - 1, 2*k - 1) for 0 < k <= n and T(n,0) = 1.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 10, 1, 1, 7, 35, 21, 1, 1, 9, 84, 126, 36, 1, 1, 11, 165, 462, 330, 55, 1, 1, 13, 286, 1287, 1716, 715, 78, 1, 1, 15, 455, 3003, 6435, 5005, 1365, 105, 1, 1, 17, 680, 6188, 19448, 24310, 12376, 2380, 136, 1, 1, 19, 969, 11628, 50388, 92378, 75582, 27132, 3876, 171, 1
Offset: 0

Author

Roger L. Bagula, Oct 02 2006

Keywords

Examples

			Triangle begins:
     1;
     1,  1;
     1,  3,   1;
     1,  5,  10,    1;
     1,  7,  35,   21,    1;
     1,  9,  84,  126,   36,    1;
     1, 11, 165,  462,  330,   55,    1;
     1, 13, 286, 1287, 1716,  715,   78,  1;
     1, 15, 455, 3003, 6435, 5005, 1365, 105, 1;
     ...
		

Programs

  • GAP
    Flat(Concatenation([1],List([1..10],n->Concatenation([1],List([1..n],m->Binomial(2*n-1,2*m-1)))))); # Muniru A Asiru, Oct 11 2018
    
  • Magma
    A123162:= func< n,k | k eq 0 select 1 else Binomial(2*n-1, 2*k-1) >;
    [A123162(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Jul 18 2023
  • Mathematica
    T[n_, k_]= If [k==0, 1, Binomial[2*n-1, 2*k-1]];
    Table[T[n, k], {n,0,12}, {k,0,n}]//Flatten
  • Maxima
    T(n, k) := if k = 0 then 1 else binomial(2*n - 1, 2*k - 1)$
    tabl(nn) := for n:0 thru nn do print(makelist(T(n, k), k, 0, n))$ /* Franck Maminirina Ramaharo, Oct 10 2018 */
    
  • SageMath
    def A123162(n,k): return binomial(2*n-1, 2*k-1) + int(k==0)
    flatten([[A123162(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Jun 01 2022
    

Formula

From Paul Barry, May 26 2008: (Start)
T(n,k) = binomial(2*n - 1, 2*k - 1) + 0^k.
Column k has g.f. (x^k/(1 - x)^(2*k + 0^k))*Sum_{j=0..k} binomial(2*k, 2*j)*x^j. (End)
From Franck Maminirina Ramaharo, Oct 10 2018: (Start)
Row n = coefficients in the expansion of ((x + sqrt(x))*(sqrt(x) - 1)^(2*n) + (x - sqrt(x))*(sqrt(x) + 1)^(2*n) + 2*x - 2)/(2*x - 2).
G.f.: (1 - (2 + x)*y + (1 - 2*x)*y^2 - (x - x^2)*y^3)/(1 - (3 + 2*x)*y + (3 + x^2)*y^2 - (1 - 2*x + x^2)*y^3).
E.g.f.: ((x + sqrt(x))*exp(y*(sqrt(x) - 1)^2) + (x - sqrt(x))*exp(y*(sqrt(x) + 1)^2) + (2*x - 2)*exp(y) - 2*x)/(2*x - 2). (End)
From G. C. Greubel, Jul 18 2023: (Start)
Sum_{k=0..n} T(n,k) = A123166(n).
T(n, n-1) = (n-1)*T(n, 1), n >= 2.
T(2*n, n) = A259557(n).
T(2*n+1, n+1) = A002458(n). (End)

Extensions

Edited by N. J. A. Sloane, Oct 04 2006
Partially edited and offset corrected by Franck Maminirina Ramaharo, Oct 10 2018

A131980 A coefficient tree from the list partition transform relating A000129, A000142, A000165, A110327, and A110330.

Original entry on oeis.org

1, 2, 6, 2, 24, 24, 120, 240, 24, 720, 2400, 720, 5040, 25200, 15120, 720, 40320, 282240, 282240, 40320, 362880, 3386880, 5080320, 1451520, 40320, 3628800, 43545600, 91445760, 43545600, 3628800, 39916800, 598752000, 1676505600, 1197504000, 199584000, 3628800
Offset: 0

Author

Tom Copeland, Oct 30 2007, Nov 29 2007, Nov 30 2007

Keywords

Comments

Construct the infinite array of polynomials
a(0,t) = 1
a(1,t) = 2
a(2,t) = 6 + 2 t
a(3,t) = 24 + 24 t
a(4,t) = 120 + 240 t + 24 t^2
a(5,t) = 720 + 2400 t + 720 t^2
a(6,t) = 5040 + 25200 t + 15120 t^2 + 720 t^3
This array is the reciprocal array of the following array b(n,t) under the list partition transform and its associated operations described in A133314.
b(0,t) = 1, b(1,t) = -2, b(2,t) = -2*(t-1), b(n,t) = 0 for n>2.
Then A000165(n) = a(n,1).
Lower triangular matrix A110327 = binomial(n,k)*a(n-k,2).
n! * A000129(n+1) = a(n,2) = A110327(n,0).
A110330 = matrix inverse of binomial(n,k)*a(n-k,2) = binomial(n,k)*b(n-k,2).
A000142(n+1) = a(n,0).
From Peter Bala, Sep 09 2013: (Start)
Let {P(n,x)}n>=0 be a polynomial sequence. Koutras has defined generalized Eulerian numbers associated with the sequence P(n,x) as the coefficients A(n,k) in the expansion of P(n,x) in a series of factorials of degree n, namely P(n,x) = Sum_{k=0..n} A(n,k)* binomial(x+n-k,n). The choice P(n,x) = x^n produces the classical Eulerian numbers of A008292. Let now P(n,x) = x*(x + 1)*...*(x + n - 1) denote the n-th rising factorial polynomial. Then the present table is the generalized Eulerian numbers associated with the polynomial sequence P(n,2*x). See A228955 for the generalized Eulerian numbers associated with the polynomial sequence P(n,2*x + 1). (End)

Examples

			Triangle begins as:
        1;
        2;
        6,        2;
       24,       24;
      120,      240,       24;
      720,     2400,      720;
     5040,    25200,    15120,      720;
    40320,   282240,   282240,    40320;
   362880,  3386880,  5080320,  1451520,   40320;
  3628800, 43545600, 91445760, 43545600, 3628800;
		

Crossrefs

Cf. A228955.

Programs

  • GAP
    Flat(List([0..10], n-> List([0..Int(n/2)], k-> Factorial(n)*Binomial(n+1, 2*k+1) ))); # G. C. Greubel, Dec 30 2019
  • Magma
    [Factorial(n)*Binomial(n+1, 2*k+1): k in [0..Floor(n/2)], n in [0..10]]; // G. C. Greubel, Dec 30 2019
    
  • Maple
    for n from 0 to 10 do
    seq( n!*binomial(n+1,2*k+1), k = 0..floor(n/2) )
    end do; # Peter Bala, Sep 09 2013
  • Mathematica
    Table[n!*Binomial[n+1, 2*k+1], {n,0,10}, {k,0,Floor[n/2]}]//Flatten (* G. C. Greubel, Dec 30 2019 *)
  • PARI
    T(n,k) = n!*binomial(n+1, 2*k+1);
    for(n=0,10, for(k=0, n\2, print1(T(n,k), ", "))) \\ G. C. Greubel, Dec 30 2019
    
  • Sage
    [[factorial(n)*binomial(n+1, 2*k+1) for k in (0..floor(n/2))] for n in (0..10)] # G. C. Greubel, Dec 30 2019
    

Formula

E.g.f. for the polynomials b(.,t), introduced above, is 1 - 2x - (t-1) * x^2; therefore e.g.f. for the polynomials a(.,t), which are the row polynomials of this array, is 1 / ( 1 - 2x - (t-1) * x^2 ) = (t-1) / ( t - ( 1 + x*(t-1) )^2 ).
Also, a(n,t) = (1 - t*u^2)^(n+1) (D_u)^n [ 1 / (1 - t*u^2) ] with eval. at u = 1/t. Compare A076743.
a(n,t) = n!*Sum_{k>=0} binomial(n+1,2k+1) * t^k = n!*Sum_{k>=0} A034867(n,k) * t^k.
Additional relations are given by formulas in A133314.
From Peter Bala, Sep 09 2013: (Start)
Recurrence equation: T(n+1,k) = (n+2 +2*k)T(n,k) + (n +2 -2*k)T(n,k-1).
Let P(n,x) = x*(x + 1)*...*(x + n - 1) denote the n-th rising factorial.
T(n,k) = Sum_{j=0..k+1} (-1)^(k+1-j)*binomial(n+1,k+1-j)*P(n,2*j) for n >= 1.
The row polynomial a(n,t) satisfies t*a(n,t)/(1 - t)^(n+1) = Sum_{j>=1} P(n,2*j)*t^j. For example, for n = 3 we have t*(24 + 24*t)/(1 - t)^4 = 2*3*4*t + (4*5*6)*t^2 + (6*7*8)*t^3 + ..., while for n = 4 we have t*(120 + 240*t + 24*t^2)/(1 - t)^5 = (2*3*4*5)*t + (4*5*6*7)*t^2 + (6*7*8*9)*t^3 + .... (End)

Extensions

Removed erroneous and duplicate statements. - Tom Copeland, Dec 03 2013

A201509 Irregular triangle read by rows: T(n,k) = 2*T(n-1,k) + T(n-2,k-1) with T(0,0) = 0, T(n,0) = T(1,1) = 1 and T(n,k) = 0 if k < 0 or if n < k.

Original entry on oeis.org

1, 1, 2, 2, 4, 5, 1, 8, 12, 4, 16, 28, 13, 1, 32, 64, 38, 6, 64, 144, 104, 25, 1, 128, 320, 272, 88, 8, 256, 704, 688, 280, 41, 1, 512, 1536, 1696, 832, 170, 10, 1024, 3328, 4096, 2352, 620, 61, 1, 2048, 7168
Offset: 0

Author

Paul Curtz, Dec 02 2011

Keywords

Comments

This is the pseudo-triangle whose successive lines are of the type T(n,0), T(n,1)+T(n-1,0), T(n,2)+T(n-1,1), ... T(n,k)+T(n-1,k-1), without 0's, with T=A201701. [e-mail, Philippe Deléham, Dec 04 2011]

Examples

			Triangle starts:
    1   1
    2   2
    4   5   1
    8  12   4
   16  28  13  1
   32  64  38  6
   64 144 104 25 1
  128 320 272 88 8
  ...
Triangle begins (full version):
    0
    1,   1
    2,   2,   0
    4,   5,   1,  0
    8,  12,   4,  0, 0
   16,  28,  13,  1, 0, 0
   32,  64,  38,  6, 0, 0, 0
   64, 144, 104, 25, 1, 0, 0, 0
  128, 320, 272, 88, 8, 0, 0, 0, 0
		

Crossrefs

Cf. A052542 (row sums).

Formula

T(n,k) = 2*T(n-1,k) + T(n-2,k-1) with T(0,0) = 0, T(n,0) = T(1,1) = 1 and T(n,k) = 0 if k < 0 or if n < k. - Philippe Deléham, Dec 05 2011
The n-th row polynomial appears to equal Sum_{k = 1..floor((n+1)/2)} binomial(n,2*k-1)*(1+t)^k. Cf. A034867. - Peter Bala, Sep 10 2012
Aside from the first two rows below, the signed coefficients appear in the expansion (b*x - 1)^2 / (a*b*x^2 - 2a*x + 1) = 1 + (2 a - 2 b)x + (4 a^2 - 5 a b + b^2)x^2 + (8 a^3 - 12 a^2b + 4 ab^2)x^3 + ..., the reciprocal of the derivative of x*(1-a*x) / (1-b*x). This is related to A263633 via the expansion (a*b*x^2 - 2a*x + 1) / (b*x - 1)^2 = 1 + (b - a) (2x + 3b x^2 + 4b^2 x^3 + ...). See also A201780. - Tom Copeland, Oct 30 2023

Extensions

Edited and new name using Philippe Deléham's formula, Joerg Arndt, Dec 13 2023
Previous Showing 11-20 of 24 results. Next