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

A010892 Inverse of 6th cyclotomic polynomial. A period 6 sequence.

Original entry on oeis.org

1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0
Offset: 0

Views

Author

Keywords

Comments

Any sequence b(n) satisfying the recurrence b(n) = b(n-1) - b(n-2) can be written as b(n) = b(0)*a(n) + (b(1)-b(0))*a(n-1).
a(n) is the determinant of the n X n matrix M with m(i,j)=1 if |i-j| <= 1 and 0 otherwise. - Mario Catalani (mario.catalani(AT)unito.it), Jan 25 2003
Also row sums of triangle in A108299; a(n)=L(n-1,1), where L is also defined as in A108299; see A061347 for L(n,-1). - Reinhard Zumkeller, Jun 01 2005
Pisano period lengths: 1, 3, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, ... - R. J. Mathar, Aug 10 2012
Periodic sequences of this type can also be calculated as a(n) = c + floor(q/(p^m-1)*p^n) mod p, where c is a constant, q is the number representing the periodic digit pattern and m is the period. c, p and q can be calculated as follows: Let D be the array representing the number pattern to be repeated, m = size of D, max = maximum value of elements in D, min = minimum value of elements in D. Then c := min, p := max - min + 1 and q := p^m*Sum_{i=1..m} (D(i)-min)/p^i. Example: D = (1, 1, 0, -1, -1, 0), c = -1, m = 6, p = 3 and q = 676 for this sequence. - Hieronymus Fischer, Jan 04 2013
B(n) = a(n+5) = S(n-1, 1) appears, together with a(n) = A057079(n+1), in the formula 2*exp(Pi*n*i/3) = A(n) + B(n)*sqrt(3)*i with i = sqrt(-1). For S(n, x) see A049310. See also a Feb 27 2014 comment on A099837. - Wolfdieter Lang, Feb 27 2014
a(n) (for n>=1) is the difference between numbers of even and odd permutations p of 1,2,...,n such that |p(i)-i|<=1 for i=1,2,...,n. - Dmitry Efimov, Jan 08 2016
From Tom Copeland, Jan 31 2016: (Start)
Specialization of the o.g.f. 1 / ((x - w1)(x-w2)) = (1/(w1-w2)) ((w1-w2) + (w1^2 - w2^2) x + (w1^3-w2^3) x^2 + ...) with w1*w2 = (1/w1) + (1/w2) = 1. Then w1 = q = e^(i*Pi/3) and w2 = 1/q = e^(-i*Pi/3), giving the o.g.f. 1 /(1-x+x^2) for this entry with a(n) = (2/sqrt(3)) sin((n+1)Pi/3). See the Copeland link for more relations.
a(n) = (q^(n+1) - q^(-(n+1))) / (q - q^(-1)), so this entry gives the o.g.f. for an instance of the quantum integers denoted by [m]_q in Morrison et al. and Tingley. (End)

Examples

			G.f. = 1 + x - x^3 - x^4 + x^6 + x^7 - x^9 - x^10 + x^12 + x^13 - x^15 + ...
		

Crossrefs

a(n) = row sums of signed triangle A049310.
Differs only by a shift from A128834.
a(n+1) = row sums of triangle A130777: repeat(1,0,-1,-1,0,1).

Programs

  • Magma
    &cat[[1,1,0,-1,-1,0]: n in [0..20]]; // Vincenzo Librandi, Apr 03 2014
  • Maple
    a:=n->coeftayl(1/(x^2-x+1), x=0, n);
    a:=n->2*sin(Pi*(n+1)/3)/sqrt(3);
    A010892:=n->[1,1,0,-1,-1,0][irem(n,6)+1];
    A010892:=n->Array(0..5,[1,1,0,-1,-1,0])[irem(n,6)];
    A010892:=n->table([0=1,1=1,2=0,3=-1,4=-1,5=0])[irem(n,6)];
    with(numtheory,cyclotomic); c := series(1/cyclotomic(6,x),x,102): seq(coeff(c,x,n),n=0..101); # Rainer Rosenthal, Jan 01 2007
  • Mathematica
    a[n_] := {1, 1, 0, -1, -1, 0}[[Mod[n, 6] + 1]]; Table[a[n], {n, 0, 101}] (* Jean-François Alcover, Jul 19 2013 *)
    CoefficientList[Series[1/Cyclotomic[6, x], {x, 0, 100}], x] (* Vincenzo Librandi, Apr 03 2014 *)
    PadRight[{},120,{1,1,0,-1,-1,0}] (* Harvey P. Dale, Jul 07 2020 *)
  • PARI
    {a(n) = (-1)^(n\3) * sign((n + 1)%3)}; /* Michael Somos, Sep 23 2005 */
    
  • PARI
    {a(n) = subst( poltchebi(n) + poltchebi(n-1), 'x, 1/2) * 2/3}; /* Michael Somos, Sep 23 2005 */
    
  • PARI
    {a(n) = [1, 1, 0, -1, -1, 0][n%6 + 1]}; /* Michael Somos, Feb 14 2006 */
    
  • PARI
    {a(n) = my(A, p, e); if( n<0, 0, n++; A = factor(n); prod( k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==2, -(-1)^e, p==3, 0, p%6 == 1, 1, (-1)^e)))}; /* Michael Somos, Oct 29 2006 */
    
  • Python
    def A010892(n): return [1,1,0,-1,-1,0][n%6] # Alec Mihailovs, Jan 01 2007
    
  • Sage
    [lucas_number1(n,1,+1) for n in range(-5, 97)] # Zerinvary Lajos, Apr 22 2009
    
  • Sage
    def A010892():
        x, y = -1, -1
        while True:
            yield -x
            x, y = y, -x + y
    a = A010892()
    [next(a) for i in range(40)]  # Peter Luschny, Jul 11 2013
    

Formula

G.f.: 1 / (1 - x + x^2).
a(n) = a(n-1) - a(n-2), a(0)=1, a(1)=1.
a(n) = ((-1)^floor(n/3) + (-1)^floor((n+1)/3))/2.
a(n) = 0 if n mod 6 = 2 or 5, a(n) = +1 if n mod 6 = 0 or 1, a(n) = -1 otherwise. a(n) = S(n, 1) = U(n, 1/2) (Chebyshev U(n, x) polynomials).
a(n) = sqrt(4/3)*Im((1/2 + i*sqrt(3/4))^(n+1)). - Henry Bottomley, Apr 12 2000
Binomial transform of A057078. a(n) = Sum_{k=0..n} C(k, n-k)*(-1)^(n-k). - Paul Barry, Sep 13 2003
a(n) = 2*sin(Pi*n/3 + Pi/3)/sqrt(3). - Paul Barry, Jan 28 2004
a(n) = Sum_{k=0..floor(n/2)} C(n-k, k)*(-1)^k. - Paul Barry, Jul 28 2004
Euler transform of length 6 sequence [1, -1, -1, 0, 0, 1]. - Michael Somos, Sep 23 2005
a(n) = a(1 - n) = -a(-2 - n) for all n in Z. - Michael Somos, Feb 14 2006
a(n) = Sum_{k=0..n} (-2)^(n-k) * A085838(n,k). - Philippe Deléham, Oct 26 2006
a(n) = b(n+1) where b(n) is multiplicative with b(2^e) = -(-1)^e if e>0, b(3^e) = 0^e, b(p^e) = 1 if p == 1 (mod 6), b(p^e) = (-1)^e if p == 5 (mod 6). - Michael Somos, Oct 29 2006
Given g.f. A(x), then, B(x) = x * A(x) satisfies 0 = f(B(x), B(x^2)) where f(u, v) = u^2 - v - 2*u*v * (1 - u). - Michael Somos, Oct 29 2006
a(2*n) = A057078(n), a(2*n+1) = A049347(n).
a(n) = Sum_{k=0..n} A109466(n,k). - Philippe Deléham, Nov 14 2006
a(n) = Sum_{k=0..n} A133607(n,k). - Philippe Deléham, Dec 30 2007
a(n) = A128834(n+1). - Jaume Oliver Lafont, Dec 05 2008
a(n) = Sum_{k=0..n} C(n+k+1,2k+1) * (-1)^k. - Paul Barry, Jun 03 2009
a(n) = A101950(n,0) = (-1)^n * A049347(n). - Philippe Deléham, Feb 10 2012
a(n) = Product_{k=1..floor(n/2)} 1 - 4*(cos(k*Pi/(n+1)))^2. - Mircea Merca, Apr 01 2012
G.f.: 1 / (1 - x / (1 + x / (1 - x))). - Michael Somos, Apr 02 2012
a(n) = -1 + floor(181/819*10^(n+1)) mod 10. - Hieronymus Fischer, Jan 03 2013
a(n) = -1 + floor(13/14*3^(n+1)) mod 3. - Hieronymus Fischer, Jan 04 2013
a(n) = 1/(1+r2)*(1/r1)^n + 1/(1+r1)*(1/r2)^n, with r1=(1-i*sqrt(3))/2 and r2=(1+i*sqrt(3))/2. - Ralf Stephan, Jul 19 2013
a(n) = ((n+1)^2 mod 3) * (-1)^floor((n+1)/3). - Wesley Ivan Hurt, Mar 15 2015
a(n-1) = n - Sum_{i=1..n-1} i*a(n-i). - Derek Orr, Apr 28 2015
a(n) = S(2*n+1, sqrt(3))/sqrt(3) = S(n, 1) with S(n, x) coefficients given in A049310. The S(n, 1) formula appeared already above. S(2*n, sqrt(3)) = A057079(n). See also a Feb 27 2014 comment above. - Wolfdieter Lang, Jan 16 2018
E.g.f.: sqrt(exp(x)*4/3) * cos(x*sqrt(3/4) - Pi/6). - Michael Somos, Jul 05 2018
a(n) = Determinant(Tri(n)), for n >= 1, with Tri(n) the n X n tridiagonal matrix with entries 1 (a special Toeplitz matrix). - Wolfdieter Lang, Sep 20 2019
a(n) = Product_{k=1..n}(1 + 2*cos(k*Pi/(n+1))). - Peter Luschny, Nov 28 2019

Extensions

Entry revised by N. J. A. Sloane, Jul 16 2004

A061347 Period 3: repeat [1, 1, -2].

Original entry on oeis.org

1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2, 1, 1, -2
Offset: 1

Views

Author

Jason Earls, Jun 07 2001

Keywords

Comments

WARNING: It is unclear whether this sequence should start at offset 1 (as written) or offset 0 (in analogy to many similar sequences, which seems to be assumed in many of the given formulas).
Inverse binomial transform of A057079. - Paul Barry, May 15 2003
The unsigned version, with g.f. (1 + x + 2*x^2)/(1 - x^3), has a(n) = 4/3 -cos(2*Pi*n/3)/3 - sqrt(3)*sin(2*Pi*n/3)/3 = gcd(Fib(n+4), Fib(n+1)). - Paul Barry, Apr 02 2004
a(n) = L(n-2,-1), where L is defined as in A108299; see also A010892 for L(n,+1). - Reinhard Zumkeller, Jun 01 2005
From the Taylor expansion of log(1 + x + x^2) at x = 1, Sum_{k > 0} a(k)/k = log(3) = A002391. This is case n = 3 of the general expression Sum_{k > 0} (1-n*!(k mod n))/k = log(n). - Jaume Oliver Lafont, Oct 16 2009
If used with offset zero, a non-simple continued fraction representation of 2+sqrt(2). - R. J. Mathar, Mar 08 2012
Periodic sequences of this type can be also calculated by a(n) = c + floor(q/(p^m-1)*p^n) mod p, where c is a constant, q is the number representing the periodic digit pattern and m is the period length. c, p and q can be calculated as follows: Let D be the array representing the number pattern to be repeated, m = size of D, max = maximum value of elements in D, min = minimum value of elements in D. Than c := min, p := max - min + 1 and q := p^m*Sum_{i=1..m} (D(i)-min)/p^i. Example: D = (1, 1, -2), c = -2, p = 4 and q = 60 for this sequence. - Hieronymus Fischer, Jan 04 2013
This is the Dirichlet inverse of A117997. - Petros Hadjicostas, Jul 25 2020

Examples

			G.f.: x + x^2 - 2*x^3 + x^4 + x^5 - 2*x^6 + x^7 + x^8 - 2*x^9 + ... - _Michael Somos_, Nov 27 2019
		

Crossrefs

Apart from signs, same as A057079 (also bin. Transf), A100063. Cf. A000045, A010892 for the rules a(n) = a(n - 1) + a(n - 2), a(n) = a(n - 1) - a(n - 2). a(n) = - a(n - 1) + a(n - 2) gives a signed version of Fibonacci numbers.
Alternating row sums of A130777: repeat(1,-2,1).

Programs

  • GAP
    Flat(List([1..50],n->[1,1,-2])); # Muniru A Asiru, Aug 02 2018
  • Magma
    &cat [[1, 1, -2]^^30]; // Wesley Ivan Hurt, Jul 01 2016
    
  • Maple
    seq(op([1, 1, -2]), n=1..50); # Wesley Ivan Hurt, Jul 01 2016
  • Mathematica
    a[n_] := {1, 1, -2}[[Mod[n - 1, 3] + 1]]; Table[a[n], {n, 108}] (* Jean-François Alcover, Jul 19 2013 *)
    PadRight[{}, 90, {1, 1, -2}] (* After Harvey P. Dale, or *)
    CoefficientList[ Series[(2x + 1)/(x^2 + x + 1), {x, 0, 89}], x]  (* or *)
    LinearRecurrence[{-1, -1}, {1, 1}, 90] (* Robert G. Wilson v, Jul 30 2018 *)
  • PARI
    a(n)=1-3*!(n%3) \\ Jaume Oliver Lafont, Oct 16 2009
    
  • Sage
    def A061347():
        x, y = -1, -1
        while True:
            yield -x
            x, y = y, -x -y
    a = A061347(); [next(a) for i in range(40)] # Peter Luschny, Jul 11 2013
    

Formula

With offset zero, a(n) = A057079(2n). a(n) = -a(n-1) - a(n-2) with a(0) = a(1) = 1.
From Mario Catalani (mario.catalani(AT)unito.it), Jan 07 2003: (Start)
G.f.: x*(1 + 2*x)/(1 + x + x^2).
a(n) = (-1)^floor(2n/3) + ((-1)^floor((2n-1)/3) + (-1)^floor((2n+1)/3))/2. (End)
a(n) = -2*cos(2*Pi*n/3). - Jaume Oliver Lafont, May 06 2008
Dirichlet g.f.: zeta(s)*(1-1/3^(s-1)). - R. J. Mathar, Feb 09 2011
a(n) = n * Sum_{k=1..n} binomial(k,n-k)/k*(-1)^(k+1). - Dmitry Kruchinin, Jun 03 2011
a(n) = -2 + floor(110/333*10^(n+1)) mod 10. - Hieronymus Fischer, Jan 04 2013
a(n) = -2 + floor(20/21*4^(n+1)) mod 4. - Hieronymus Fischer, Jan 04 2013
a(n) = a(n-3) for n > 3. - Wesley Ivan Hurt, Jul 01 2016
E.g.f.: 2 - 2*cos(sqrt(3)*x/2)*exp(-x/2). - Ilya Gutkovskiy, Jul 01 2016
a(n) = (-1)^n*hypergeom([-n/2-1, -n/2-3/2], [-n-2], 4). - Peter Luschny, Dec 17 2016
a(n) = A000032(n) - A007040(n), for n > 1. - Wojciech Florek, Feb 20 2018

Extensions

Better definition from M. F. Hasler, Jan 13 2013

A046854 Triangle read by rows: T(n, k) = binomial(floor((n+k)/2), k), n >= k >= 0.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Row sums are Fibonacci(n+2). Diagonal sums are A016116. - Paul Barry, Jul 07 2004
Riordan array (1/(1-x), x/(1-x^2)). Matrix inverse is A106180. - Paul Barry, Apr 24 2005
As an infinite lower triangular matrix * [1,2,3,...] = A055244. - Gary W. Adamson, Dec 23 2008
From Emeric Deutsch, Jun 18 2010: (Start)
T(n,k) is the number of alternating parity increasing subsequences of {1,2,...,n} of size k, starting with an odd number (Terquem's problem, see the Riordan reference, p. 17). Example: T(8,5)=6 because we have 12345, 12347, 12367, 12567, 14567, and 34567.
T(n,k) is the number of alternating parity increasing subsequences of {1,2,...,n,n+1} of size k, starting with an even number. Example: T(7,4)=5 because we have 2345, 2347, 2367, 2567, and 4567. (End)
From L. Edson Jeffery, Mar 01 2011: (Start)
This triangle can be constructed as follows. Interlace two copies of the table of binomial coefficients to get the preliminary table
1
1
1 1
1 1
1 2 1
1 2 1
1 3 3 1
1 3 3 1
...,
then shift each entire r-th column up r rows, r=0,1,2,.... Also, a signed version of this sequence (A187660 in tabular form) begins with
1;
1, -1;
1, -1, -1;
1, -2, -1, 1;
1, -2, -3, 1, 1;
...
(compare with A066170, A130777). Let T(N,k) denote the k-th entry in row N of the signed table. Then, for N>1, row N gives the coefficients of the characteristic function p_N(x) = Sum_{k=0..N} T(N,k)*x^(N-k) = 0 of the N X N matrix U_N=[(0 ... 0 1);(0 ... 0 1 1);...;(0 1 ... 1);(1 ... 1)]. Now let Q_r(t) be a polynomial with recurrence relation Q_r(t)=t*Q_(r-1)(t)-Q_(r-2)(t) (r>1), with Q_0(t)=1 and Q_1(t)=t. Then p_N(x)=0 has solutions Q_(N-1)(phi_j), where phi_j=2*(-1)^(j-1)*cos(j*Pi/(2*N+1)), j=1,2,...,N.
For example, row N=3 is {1,-2,-1,1}, giving the coefficients of the characteristic function p_3(x) = x^3-2*x^2-x+1 = 0 for the 3 X 3 matrix U_3=[(0 0 1);(0 1 1);(1 1 1)], with eigenvalues Q_2(phi_j)=[2*(-1)^(j-1)*cos(j*Pi/7)]^2-1, j=1,2,3. (End)
Given the signed polynomials (+--++--,...) of the triangle, the largest root of the n-th row polynomial is the longest (2n+1) regular polygon diagonal length, with edge = 1. Example: the largest root to x^3 - 2x^2 - x + 1 = 0 is 2.24697...; the longest heptagon diagonal, sin(3*Pi/7)/sin(Pi/7). - Gary W. Adamson, Sep 06 2011
Given the signed polynomials from Gary W. Adamson's comment, the largest root of the n-th polynomial also equals the length from the center to a corner (vertex) of a regular 2*(2*n+1)-sided polygon with side (edge) length = 1. - L. Edson Jeffery, Jan 01 2012
Put f(x,0) = 1 and f(x,n) = x + 1/f(x,n-1). Then f(x,n) = u(x,n)/v(x,n), where u(x,n) and v(x,n) are polynomials. The flattened triangles of coefficients of u and v are both essentially A046854, as indicated by the Mathematica program headed "Polynomials". - Clark Kimberling, Oct 12 2014
From Jeremy Dover, Jun 07 2016: (Start)
T(n,k) is the number of binary strings of length n+1 starting with 0 that have exactly k pairs of consecutive 0's and no pairs of consecutive 1's.
T(n,k) is the number of binary strings of length n+2 starting with 1 that have exactly k pairs of consecutive 0's and no pairs of consecutive 1's. (End)

Examples

			Triangle begins:
  1;
  1 1;
  1 1 1;
  1 2 1 1;
  1 2 3 1 1;
  1 3 3 4 1 1;
  ...
		

References

  • J. Riordan, An Introduction to Combinatorial Analysis, Princeton University Press, 1978. [Emeric Deutsch, Jun 18 2010]

Crossrefs

Reflected version of A065941, which is considered the main entry. A deficient version is in A030111.
Cf. A055244. - Gary W. Adamson, Dec 23 2008

Programs

  • GAP
    Flat(List([0..16], n-> List([0..n], k-> Binomial(Int((n+k)/2), k) ))); # G. C. Greubel, Jul 13 2019
  • Haskell
    a046854 n k = a046854_tabl !! n !! k
    a046854_row n = a046854_tabl !! n
    a046854_tabl = [1] : f [1] [1,1] where
       f us vs = vs : f vs  (zipWith (+) (us ++ [0,0]) ([0] ++ vs))
    -- Reinhard Zumkeller, Apr 24 2013
    
  • Magma
    [Binomial(Floor((n+k)/2), k): k in [0..n], n in [0..16]]; // G. C. Greubel, Jul 13 2019
    
  • Maple
    A046854:= proc(n,k): binomial(floor(n/2+k/2), k) end: seq(seq(A046854(n,k),k=0..n),n=0..16); # Nathaniel Johnston, Jun 30 2011
  • Mathematica
    Table[Binomial[Floor[(n+k)/2], k], {n,0,16}, {k,0,n}]//Flatten
    (* next program: Polynomials *)
    z = 12; f[x_, n_] := x + 1/f[x, n - 1]; f[x_, 1] = 1;
    t = Table[Factor[f[x, n]], {n, 1, z}]
    u = Flatten[CoefficientList[Numerator[t], x]] (* this sequence *)
    v = Flatten[CoefficientList[Denominator[t], x]]
    (* Clark Kimberling, Oct 13 2014 *)
  • PARI
    T(n,k) = binomial((n+k)\2, k); \\ G. C. Greubel, Jul 13 2019
    
  • Sage
    [[binomial(floor((n+k)/2), k) for k in (0..n)] for n in (0..16)] # G. C. Greubel, Jul 13 2019
    

Formula

T(n,k) = binomial(floor((n+k)/2), k).
G.f.: (1+x)/(1-x*y-x^2). - Ralf Stephan, Feb 13 2005
Triangle = A097806 * A168561, as infinite lower triangular matrices. - Gary W. Adamson, Oct 28 2007
T(n,k) = A065941(n,n-k) = abs(A130777(n,k)) = abs(A066170(n,k)) = abs(A187660(n,k)). - Johannes W. Meijer, Aug 08 2011
For n > 1: T(n, k) = T(n-1, k-1) + T(n-2, k), 0 < k < n-1. - Reinhard Zumkeller, Apr 24 2013
T(n,k) = A168561(n,k) + A168561(n-1,k). - R. J. Mathar, Feb 10 2024

A108299 Triangle read by rows, 0 <= k <= n: T(n,k) = binomial(n-[(k+1)/2],[k/2])*(-1)^[(k+1)/2].

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
Offset: 0

Views

Author

Reinhard Zumkeller, Jun 01 2005

Keywords

Comments

Matrix inverse of A124645.
Let L(n,x) = Sum_{k=0..n} T(n,k)*x^(n-k) and Pi=3.14...:
L(n,x) = Product_{k=1..n} (x - 2*cos((2*k-1)*Pi/(2*n+1)));
Sum_{k=0..n} T(n,k) = L(n,1) = A010892(n+1);
Sum_{k=0..n} abs(T(n,k)) = A000045(n+2);
abs(T(n,k)) = A065941(n,k), T(n,k) = A065941(n,k)*A087960(k);
T(2*n,k) + T(2*n+1,k+1) = 0 for 0 <= k <= 2*n;
T(n,0) = A000012(n) = 1; T(n,1) = -1 for n > 0;
T(n,2) = -(n-1) for n > 1; T(n,3) = A000027(n)=n for n > 2;
T(n,4) = A000217(n-3) for n > 3; T(n,5) = -A000217(n-4) for n > 4;
T(n,6) = -A000292(n-5) for n > 5; T(n,7) = A000292(n-6) for n > 6;
T(n,n-3) = A058187(n-3)*(-1)^floor(n/2) for n > 2;
T(n,n-2) = A008805(n-2)*(-1)^floor((n+1)/2) for n > 1;
T(n,n-1) = A008619(n-1)*(-1)^floor(n/2) for n > 0;
T(n,n) = L(n,0) = (-1)^floor((n+1)/2);
L(n,1) = A010892(n+1); L(n,-1) = A061347(n+2);
L(n,2) = 1; L(n,-2) = A005408(n)*(-1)^n;
L(n,3) = A001519(n); L(n,-3) = A002878(n)*(-1)^n;
L(n,4) = A001835(n+1); L(n,-4) = A001834(n)*(-1)^n;
L(n,5) = A004253(n); L(n,-5) = A030221(n)*(-1)^n;
L(n,6) = A001653(n); L(n,-6) = A002315(n)*(-1)^n;
L(n,7) = A049685(n); L(n,-7) = A033890(n)*(-1)^n;
L(n,8) = A070997(n); L(n,-8) = A057080(n)*(-1)^n;
L(n,9) = A070998(n); L(n,-9) = A057081(n)*(-1)^n;
L(n,10) = A072256(n+1); L(n,-10) = A054320(n)*(-1)^n;
L(n,11) = A078922(n+1); L(n,-11) = A097783(n)*(-1)^n;
L(n,12) = A077417(n); L(n,-12) = A077416(n)*(-1)^n;
L(n,13) = A085260(n);
L(n,14) = A001570(n); L(n,-14) = A028230(n)*(-1)^n;
L(n,n) = A108366(n); L(n,-n) = A108367(n).
Row n of the matrix inverse (A124645) has g.f.: x^floor(n/2)*(1-x)^(n-floor(n/2)). - Paul D. Hanna, Jun 12 2005
From L. Edson Jeffery, Mar 12 2011: (Start)
Conjecture: Let N=2*n+1, with n > 2. Then T(n,k) (0 <= k <= n) gives the k-th coefficient in the characteristic function p_N(x)=0, of degree n in x, for the n X n tridiagonal unit-primitive matrix G_N (see [Jeffery]) of the form
G_N=A_{N,1}=
(0 1 0 ... 0)
(1 0 1 0 ... 0)
(0 1 0 1 0 ... 0)
...
(0 ... 0 1 0 1)
(0 ... 0 1 1),
with solutions phi_j = 2*cos((2*j-1)*Pi/N), j=1,2,...,n. For example, for n=3,
G_7=A_{7,1}=
(0 1 0)
(1 0 1)
(0 1 1).
We have {T(3,k)}=(1,-1,-2,1), while the characteristic function of G_7 is p(x) = x^3-x^2-2*x+1 = 0, with solutions phi_j = 2*cos((2*j-1)*Pi/7), j=1,2,3. (End)
The triangle sums, see A180662 for their definitions, link A108299 with several sequences, see the crossrefs. - Johannes W. Meijer, Aug 08 2011
The roots to the polynomials are chaotic using iterates of the operation (x^2 - 2), with cycle lengths L and initial seeds returning to the same term or (-1)* the seed. Periodic cycle lengths L are shown in A003558 such that for the polynomial represented by row r, the cycle length L is A003558(r-1). The matrices corresponding to the rows as characteristic polynomials are likewise chaotic [cf. Kappraff et al., 2005] with the same cycle lengths but substituting 2*I for the "2" in (x^2 - 2), where I = the Identity matrix. For example, the roots to x^3 - x^2 - 2x + 1 = 0 are 1.801937..., -1.246979..., and 0.445041... With 1.801937... as the initial seed and using (x^2 - 2), we obtain the 3-period trajectory of 8.801937... -> 1.246979... -> -0.445041... (returning to -1.801937...). We note that A003558(2) = 3. The corresponding matrix M is: [0,1,0; 1,0,1; 0,1,1,]. Using seed M with (x^2 - 2*I), we obtain the 3-period with the cycle completed at (-1)*M. - Gary W. Adamson, Feb 07 2012

Examples

			Triangle begins:
  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;
  ...
		

References

  • Friedrich L. Bauer, 'De Moivre und Lagrange: Cosinus eines rationalen Vielfachen von Pi', Informatik Spektrum 28 (Springer, 2005).
  • Jay Kappraff, S. Jablan, G. Adamson, & R. Sazdonovich: "Golden Fields, Generalized Fibonacci Sequences, & Chaotic Matrices"; FORMA, Vol 19, No 4, (2005).

Crossrefs

Cf. A049310, A039961, A124645 (matrix inverse).
Triangle sums (see the comments): A193884 (Kn11), A154955 (Kn21), A087960 (Kn22), A000007 (Kn3), A010892 (Fi1), A134668 (Fi2), A078031 (Ca2), A193669 (Gi1), A001519 (Gi3), A193885 (Ze1), A050935 (Ze3). - Johannes W. Meijer, Aug 08 2011
Cf. A003558.

Programs

  • Haskell
    a108299 n k = a108299_tabl !! n !! k
    a108299_row n = a108299_tabl !! n
    a108299_tabl = [1] : iterate (\row ->
       zipWith (+) (zipWith (*) ([0] ++ row) a033999_list)
                   (zipWith (*) (row ++ [0]) a059841_list)) [1,-1]
    -- Reinhard Zumkeller, May 06 2012
  • Maple
    A108299 := proc(n,k): binomial(n-floor((k+1)/2), floor(k/2))*(-1)^floor((k+1)/2) end: seq(seq(A108299 (n,k), k=0..n), n=0..11); # Johannes W. Meijer, Aug 08 2011
  • Mathematica
    t[n_, k_?EvenQ] := I^k*Binomial[n-k/2, k/2]; t[n_, k_?OddQ] := -I^(k-1)*Binomial[n+(1-k)/2-1, (k-1)/2]; Table[t[n, k], {n, 0, 12}, {k, 0, n}] // Flatten (* Jean-François Alcover, May 16 2013 *)
  • PARI
    {T(n,k)=polcoeff(polcoeff((1-x*y)/(1-x+x^2*y^2+x^2*O(x^n)),n,x)+y*O(y^k),k,y)} (Hanna)
    

Formula

T(n,k) = binomial(n-floor((k+1)/2),floor(k/2))*(-1)^floor((k+1)/2).
T(n+1, k) = if sign(T(n, k-1))=sign(T(n, k)) then T(n, k-1)+T(n, k) else -T(n, k-1) for 0 < k < n, T(n, 0) = 1, T(n, n) = (-1)^floor((n+1)/2).
G.f.: A(x, y) = (1 - x*y)/(1 - x + x^2*y^2). - Paul D. Hanna, Jun 12 2005
The generating polynomial (in z) of row n >= 0 is (u^(2*n+1) + v^(2*n+1))/(u + v), where u and v are defined by u^2 + v^2 = 1 and u*v = z. - Emeric Deutsch, Jun 16 2011
From Johannes W. Meijer, Aug 08 2011: (Start)
abs(T(n,k)) = A065941(n,k) = abs(A187660(n,n-k));
T(n,n-k) = A130777(n,k); abs(T(n,n-k)) = A046854(n,k) = abs(A066170(n,k)). (End)

Extensions

Corrected and edited by Philippe Deléham, Oct 20 2008

A111125 Triangle read by rows: T(k,s) = ((2*k+1)/(2*s+1))*binomial(k+s,2*s), 0 <= s <= k.

Original entry on oeis.org

1, 3, 1, 5, 5, 1, 7, 14, 7, 1, 9, 30, 27, 9, 1, 11, 55, 77, 44, 11, 1, 13, 91, 182, 156, 65, 13, 1, 15, 140, 378, 450, 275, 90, 15, 1, 17, 204, 714, 1122, 935, 442, 119, 17, 1, 19, 285, 1254, 2508, 2717, 1729, 665, 152, 19, 1, 21, 385, 2079, 5148, 7007, 5733, 2940, 952, 189, 21, 1
Offset: 0

Views

Author

N. J. A. Sloane, Oct 16 2005

Keywords

Comments

Riordan array ((1+x)/(1-x)^2, x/(1-x)^2). Row sums are A002878. Diagonal sums are A003945. Inverse is A113187. An interesting factorization is (1/(1-x), x/(1-x))(1+2*x, x*(1+x)). - Paul Barry, Oct 17 2005
Central coefficients of rows with odd numbers of term are A052227.
From Wolfdieter Lang, Jun 26 2011: (Start)
T(k,s) appears as T_s(k) in the Knuth reference, p. 285.
This triangle is related to triangle A156308(n,m), appearing in this reference as U_m(n) on p. 285, by T(k,s) - T(k-1,s) = A156308(k,s), k>=s>=1 (identity on p. 286). T(k,s) = A156308(k+1,s+1) - A156308(k,s+1), k>=s>=0 (identity on p. 286).
(End)
A111125 is jointly generated with A208513 as an array of coefficients of polynomials v(n,x): initially, u(1,x)= v(1,x)= 1; for n>1, u(n,x)= u(n-1,x) +x*(x+1)*v(n-1) and v(n,x)= u(n-1,x) +x*v(n-1,x) +1. See the Mathematica section. The columns of A111125 are identical to those of A208508. Here, however, the alternating row sums are periodic (with period 1,2,1,-1,-2,-1). - Clark Kimberling, Feb 28 2012
This triangle T(k,s) (with signs and columns scaled with powers of 5) appears in the expansion of Fibonacci numbers F=A000045 with multiples of odd numbers as indices in terms of odd powers of F-numbers. See the Jennings reference, p. 108, Theorem 1. Quoted as Lemma 3 in the Ozeki reference given in A111418. The formula is: F_{(2*k+1)*n} = Sum_{s=0..k} ( T(k,s)*(-1)^((k+s)*n)*5^s*F_{n}^(2*s+1) ), k >= 0, n >= 0. - Wolfdieter Lang, Aug 24 2012
From Wolfdieter Lang, Oct 18 2012: (Start)
This triangle T(k,s) appears in the formula x^(2*k+1) - x^(-(2*k+1)) = Sum_{s=0..k} ( T(k,s)*(x-x^(-1))^(2*s+1) ), k>=0. Prove the inverse formula (due to the Riordan property this will suffice) with the binomial theorem. Motivated to look into this by the quoted paper of Wang and Zhang, eq. (1.4).
Alternating row sums are A057079.
The Z-sequence of this Riordan array is A217477, and the A-sequence is (-1)^n*A115141(n). For the notion of A- and Z-sequences for Riordan triangles see a W. Lang link under A006232. (End)
The signed triangle ((-1)^(k-s))*T(k,s) gives the coefficients of (x^2)^s of the polynomials C(2*k+1,x)/x, with C the monic integer Chebyshev T-polynomials whose coefficients are given in A127672 (C is there called R). See the odd numbered rows there. This signed triangle is the Riordan array ((1-x)/(1+x)^2, x/(1+x)^2). Proof by comparing the o.g.f. of the row polynomials where x is replaced by x^2 with the odd part of the bisection of the o.g.f. for C(n,x)/x. - Wolfdieter Lang, Oct 23 2012
From Wolfdieter Lang, Oct 04 2013: (Start)
The signed triangle S(k,s) := ((-1)^(k-s))*T(k,s) (see the preceding comment) is used to express in a (4*(k+1))-gon the length ratio s(4*(k+1)) = 2*sin(Pi/4*(k+1)) = 2*cos((2*k+1)*Pi/(4*(k+1))) of a side/radius as a polynomial in rho(4*(k+1)) = 2*cos(Pi/4*(k+1)), the length ratio (smallest diagonal)/side:
s(4*(k+1)) = Sum_{s=0..k} ( S(k,s)*rho(4*(k+1))^(2*s+1) ).
This is to be computed modulo C(4*(k+1), rho(4*(k+1)) = 0, the minimal polynomial (see A187360) in order to obtain s(4*(k+1)) as an integer in the algebraic number field Q(rho(4*(k+1))) of degree delta(4*(k+1)) (see A055034). Thanks go to Seppo Mustonen for asking me to look into the problem of the square of the total length in a regular n-gon, where this formula is used in the even n case. See A127677 for the formula in the (4*k+2)-gon. (End)
From Wolfdieter Lang, Aug 14 2014: (Start)
The row polynomials for the signed triangle (see the Oct 23 2012 comment above), call them todd(k,x) = Sum_{s=0..k} ( (-1)^(k-s)*T(k,s)*x^s ) = S(k, x-2) - S(k-1, x-2), k >= 0, with the Chebyshev S-polynomials (see their coefficient triangle (A049310) and S(-1, x) = 0), satisfy the recurrence todd(k, x) = (-1)^(k-1)*((x-4)/2)*todd(k-1, 4-x) + ((x-2)/2)*todd(k-1, x), k >= 1, todd(0, x) = 1. From the Aug 03 2014 comment on A130777.
This leads to a recurrence for the signed triangle, call it S(k,s) as in the Oct 04 2013 comment: S(k,s) = (1/2)*(1 + (-1)^(k-s))*S(k-1,s-1) + (2*(s+1)*(-1)^(k-s) - 1)*S(k-1,s) + (1/2)*(-1)^(k-s)*Sum_{j=0..k-s-2} ( binomial(j+s+2,s)*4^(j+2)* S(k-1, s+1+j) ) for k >= s >= 1, and S(k,s) = 0 if k < s and S(k,0) = (-1)^k*(2*k+1). Note that the recurrence derived from the Riordan A-sequence A115141 is similar but has simpler coefficients: S(k,s) = sum(A115141(j)*S(k-1,s-1+j), j=0..k-s), k >= s >=1.
(End)
From Tom Copeland, Nov 07 2015: (Start)
Rephrasing notes here: Append an initial column of zeros, except for a 1 at the top, to A111125 here. Then the partial sums of the columns of this modified entry are contained in A208513. Append an initial row of zeros to A208513. Then the difference of consecutive pairs of rows of the modified A208513 generates the modified A111125. Cf. A034807 and A127677.
For relations among the characteristic polynomials of Cartan matrices of the Coxeter root groups, Chebyshev polynomials, cyclotomic polynomials, and the polynomials of this entry, see Damianou (p. 20 and 21) and Damianou and Evripidou (p. 7).
As suggested by the equations on p. 7 of Damianou and Evripidou, the signed row polynomials of this entry are given by (p(n,x))^2 = (A(2*n+1, x) + 2)/x = (F(2*n+1, (2-x), 1, 0, 0, ... ) + 2)/x = F(2*n+1, -x, 2*x, -3*x, ..., (-1)^n n*x)/x = -F(2*n+1, x, 2*x, 3*x, ..., n*x)/x, where A(n,x) are the polynomials of A127677 and F(n, ...) are the Faber polynomials of A263196. Cf. A127672 and A127677.
(End)
The row polynomials P(k, x) of the signed triangle S(k, s) = ((-1)^(k-s))*T(k, s) are given from the row polynomials R(2*k+1, x) of triangle A127672 by
P(k, x) = R(2*k+1, sqrt(x))/sqrt(x). - Wolfdieter Lang, May 02 2021

Examples

			Triangle T(k,s) begins:
k\s  0    1     2     3     4     5     6    7    8   9 10
0:   1
1:   3    1
2:   5    5     1
3:   7   14     7     1
4:   9   30    27     9     1
5:  11   55    77    44    11     1
6:  13   91   182   156    65    13     1
7:  15  140   378   450   275    90    15    1
8:  17  204   714  1122   935   442   119   17    1
9:  19  285  1254  2508  2717  1729   665  152   19   1
10: 21  385  2079  5148  7007  5733  2940  952  189  21  1
... Extended and reformatted by _Wolfdieter Lang_, Oct 18 2012
Application for Fibonacci numbers F_{(2*k+1)*n}, row k=3:
F_{7*n} = 7*(-1)^(3*n)*F_n + 14*(-1)^(4*n)*5*F_n^3 + 7*(-1)^(5*n)*5^2*F_n^5 + 1*(-1)^(6*n)*5^3*F_n^7, n>=0. - _Wolfdieter Lang_, Aug 24 2012
Example for the  Z- and A-sequence recurrences  of this Riordan triangle: Z = A217477 = [3,-4,12,-40,...]; T(4,0) = 3*7 -4*14 +12*7 -40*1 = 9. A =  [1, 2, -1, 2, -5, 14, ..]; T(5,2) = 1*30 + 2*27 - 1*9 + 2*1= 77. _Wolfdieter Lang_, Oct 18 2012
Example for the (4*(k+1))-gon length ratio s(4*(k+1))(side/radius) as polynomial in the ratio rho(4*(k+1)) ((smallest diagonal)/side): k=0, s(4) = 1*rho(4) = sqrt(2); k=1, s(8) = -3*rho(8) + rho(8)^3 = sqrt(2-sqrt(2)); k=2, s(12) = 5*rho(12) - 5*rho(12)^3 + rho(12)^5, and C(12,x) = x^4 - 4*x^2 + 1, hence rho(12)^5 = 4*rho(12)^3 - rho(12), and s(12) = 4*rho(12) - rho(12)^3 = sqrt(2 - sqrt(3)). - _Wolfdieter Lang_, Oct 04 2013
Example for the recurrence for the signed triangle S(k,s)= ((-1)^(k-s))*T(k,s) (see the Aug 14 2014 comment above):
S(4,1) = 0 + (-2*2 - 1)*S(3,1) - (1/2)*(3*4^2*S(3,2) + 4*4^3*S(3,3)) = - 5*14 - 3*8*(-7) - 128*1 = -30. The recurrence from the Riordan A-sequence A115141 is S(4,1) = -7 -2*14 -(-7) -2*1 = -30. - _Wolfdieter Lang_, Aug 14 2014
		

Crossrefs

Mirror image of A082985, which see for further references, etc.
Also closely related to triangles in A098599 and A100218.

Programs

  • Magma
    [((2*n+1)/(n+k+1))*Binomial(n+k+1, 2*k+1): k in [0..n], n in [0..12]];  // G. C. Greubel, Feb 01 2022
  • Mathematica
    (* First program *)
    u[1, x_]:=1; v[1, x_]:=1; z=16;
    u[n_, x_]:= u[n-1, x] + x*v[n-1, x];
    v[n_, x_]:= u[n-1, x] + (x+1)*v[n-1, x] + 1;
    Table[Expand[u[n, x]], {n, 1, z/2}]
    Table[Expand[v[n, x]], {n, 1, z/2}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%]  (* A208513 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]  (* A111125 *) (* Clark Kimberling, Feb 28 2012 *)
    (* Second program *)
    T[n_, k_]:= ((2*n+1)/(2*k+1))*Binomial[n+k, 2*k];
    Table[T[n, k], {n,0,15}, {k,0,n}]//Flatten (* G. C. Greubel, Feb 01 2022 *)
  • Sage
    @CachedFunction
    def T(n,k):
        if n< 0: return 0
        if n==0: return 1 if k == 0 else 0
        h = 3*T(n-1,k) if n==1 else 2*T(n-1,k)
        return T(n-1,k-1) - T(n-2,k) - h
    A111125 = lambda n,k: (-1)^(n-k)*T(n,k)
    for n in (0..9): [A111125(n,k) for k in (0..n)] # Peter Luschny, Nov 20 2012
    

Formula

T(k,s) = ((2*k+1)/(2*s+1))*binomial(k+s,2*s), 0 <= s <= k.
From Peter Bala, Apr 30 2012: (Start)
T(n,k) = binomial(n+k,2*k) + 2*binomial(n+k,2*k+1).
The row generating polynomials P(n,x) are a generalization of the Morgan-Voyce polynomials b(n,x) and B(n,x). They satisfy the recurrence equation P(n,x) = (x+2)*P(n-1,x) - P(n-2,x) for n >= 2, with initial conditions P(0,x) = 1, P(1,x) = x+r+1 and with r = 2. The cases r = 0 and r = 1 give the Morgan-Voyce polynomials A085478 and A078812 respectively. Andre-Jeannin has considered the case of general r.
P(n,x) = U(n+1,1+x/2) + U(n,1+x/2), where U(n,x) denotes the Chebyshev polynomial of the second kind - see A053117. P(n,x) = (2/x)*(T(2*n+2,u)-T(2*n,u)), where u = sqrt((x+4)/4) and T(n,x) denotes the Chebyshev polynomial of the first kind - see A053120. P(n,x) = Product_{k = 1..n} ( x + 4*(sin(k*Pi/(2*n+1)))^2 ). P(n,x) = 1/x*(b(n+1,x) - b(n-1,x)) and P(n,x) = 1/x*{(b(2*n+2,x)+1)/b(n+1,x) - (b(2*n,x)+1)/b(n,x)}, where b(n,x) := Sum_{k = 0..n} binomial(n+k,2*k)*x^k are the Morgan-Voyce polynomials of A085478. Cf. A211957.
(End)
From Wolfdieter Lang, Oct 18 2012 (Start)
O.g.f. column No. s: ((1+x)/(1-x)^2)*(x/(1-x)^2)^s, s >= 0. (from the Riordan data given in a comment above).
O.g.f. of the row polynomials R(k,x):= Sum_{s=0..k} ( T(k,s)*x^s ), k>=0: (1+z)/(1-(2+x)*z+z^2) (from the Riordan property).
(End)
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k), T(0,0) = 1, T(1,0) = 3, T(1,1) = 1, T(n,k) = 0 if k<0 or if k>n. - Philippe Deléham, Nov 12 2013

Extensions

More terms from Paul Barry, Oct 17 2005

A063886 Number of n-step walks on a line starting from the origin but not returning to it.

Original entry on oeis.org

1, 2, 2, 4, 6, 12, 20, 40, 70, 140, 252, 504, 924, 1848, 3432, 6864, 12870, 25740, 48620, 97240, 184756, 369512, 705432, 1410864, 2704156, 5408312, 10400600, 20801200, 40116600, 80233200, 155117520, 310235040, 601080390, 1202160780, 2333606220, 4667212440
Offset: 0

Views

Author

Henry Bottomley, Aug 28 2001

Keywords

Comments

A Chebyshev transform of A007877(n+1). The g.f. is transformed to (1+x)/((1-x)(1+x^2)) under the mapping G(x)->(1/(1+x^2))G(1/(1+x^2)). - Paul Barry, Oct 12 2004
a(n-1) = 2*C(n-2, floor((n-2)/2)) is also the number of bit strings of length n in which the number of 00 substrings is equal to the number of 11 substrings. For example, when n = 4 we have 4 such bit strings: 0011, 0101, 1010, and 1100. - Angel Plaza, Apr 23 2009
Hankel transform is A120617. - Paul Barry, Aug 10 2009
The Hankel transform of a(n) is (-2)^C(n+1,2). The Hankel transform of (-1)^C(n+1,2)*a(n) is (-1)^C(n+1,2)*A164584(n). - Paul Barry, Aug 17 2009
For n > 1, a(n) is also the number of n-step walks starting from the origin and returning to it exactly once. - Geoffrey Critzer, Jan 24 2010
-a(n) is the Z-sequence for the Riordan array A130777. (See the W. Lang link under A006232 for A- and Z-sequences for Riordan matrices). - Wolfdieter Lang, Jul 12 2011
Number of subsets of {1,...,n} in which the even elements appear as often at even positions as at odd positions. - Gus Wiseman, Mar 17 2018

Examples

			a(4) = 6 because there are six length four walks that do not return to the origin: {-1, -2, -3, -4}, {-1, -2, -3, -2}, {-1, -2, -1, -2}, {1, 2, 1, 2}, {1, 2, 3, 2}, {1, 2, 3, 4}. There are also six such walks that return exactly one time: {-1, -2, -1, 0}, {-1, 0, -1, -2}, {-1, 0, 1, 2}, {1, 0, -1, -2}, {1, 0, 1, 2}, {1, 2, 1, 0}. - _Geoffrey Critzer_, Jan 24 2010
The a(5) = 12 subsets in which the even elements appear as often at even positions as at odd positions: {}, {1}, {3}, {5}, {1,3}, {1,5}, {2,4}, {3,5}, {1,2,4}, {1,3,5}, {2,4,5}, {1,2,4,5}. - _Gus Wiseman_, Mar 17 2018
		

Crossrefs

Programs

  • Magma
    [1] cat [2*Binomial(n-1, Floor((n-1)/2)): n in [1..40]]; // G. C. Greubel, Jun 07 2023
    
  • Maple
    seq(seq(binomial(2*j,j)*i, i=1..2),j=0..16); # Zerinvary Lajos, Apr 28 2007
    # second Maple program:
    a:= proc(n) option remember; `if`(n<2, n+1,
           4*a(n-2) +2*(a(n-1) -4*a(n-2))/n)
        end:
    seq(a(n), n=0..40);  # Alois P. Heinz, Feb 10 2014
    # third program:
    A063886 := series(BesselI(0, 2*x)*(1 + x*2 + x*Pi*StruveL(1, 2*x)) - Pi*x*BesselI(1, 2*x)*StruveL(0, 2*x), x = 0, 34): seq(n!*coeff(A063886, x, n), n = 0 .. 33); # Mélika Tebni, Jun 17 2024
  • Mathematica
    Table[Length[Select[Map[Accumulate, Strings[{-1, 1}, n]], Count[ #, 0] == 0 &]], {n, 0, 20}] (* Geoffrey Critzer, Jan 24 2010 *)
    CoefficientList[Series[Sqrt[(1+2x)/(1-2x)],{x,0,40}],x] (* Harvey P. Dale, Apr 28 2016 *)
  • PARI
    a(n)=(n==0)+2*binomial(n-1,(n-1)\2)
    
  • PARI
    a(n) = 2^n*prod(k=0,n-1,(k/n+1/n)^((-1)^k)); \\ Michel Marcus, Dec 03 2013
    
  • Python
    from math import ceil
    from sympy import binomial
    def a(n):
        if n==0: return 1
        return 2*binomial(n-1,(n-1)//2)
    print([a(n) for n in range(18)])
    # David Nacin, Feb 29 2012
    
  • SageMath
    [2*binomial(n-1, (n-1)//2) + int(n==0) for n in range(41)] # G. C. Greubel, Jun 07 2023

Formula

G.f.: sqrt((1+2*x)/(1-2*x)).
a(n+1) = 2*C(n, floor(n/2)) = 2*A001405(n); a(2n) = C(2n, n) = A000984(n) = 4*a(2n-2)-|A002420(n)| = 4*a(2n-2)-2*A000108(n-1) = 2*A001700(n-1); a(2n+1) = 2*a(2n) = A028329(n).
2*a(n) = A047073(n+1).
a(n) = Sum_{k=0..n} abs(A106180(n,k)). - Philippe Deléham, Oct 06 2006
a(n) = Sum_{k=0..n} (k+1)binomial(n, (n-k)/2) ( 1-cos((k+1)*Pi/2) (1+(-1)^(n-k))/(n+k+2) ). - Paul Barry, Oct 12 2004
G.f.: 1/(1-2*x/(1+x/(1+x/(1-x/(1-x/(1+x/(1+x/(1-x/(1-x/(1+ ... (continued fraction). - Paul Barry, Aug 10 2009
G.f.: 1 + 2*x/(G(0)-x+x^2) where G(k)= 1 - 2*x^2 - x^4/G(k+1); (continued fraction, 1-step). - Sergei N. Gladkovskii, Aug 10 2012
D-finite with recurrence: n*a(n) = 2*a(n-1) + 4*(n-2)*a(n-2). - R. J. Mathar, Dec 03 2012
From Sergei N. Gladkovskii, Jul 26 2013: (Start)
G.f.: 1/G(0), where G(k) = 1 - 2*x/(1 + 2*x/(1 + 1/G(k+1) )); (continued fraction).
G.f.: G(0), where G(k) = 1 + 2*x/(1 - 2*x/(1 + 1/G(k+1) )); (continued fraction).
G.f.: W(0)/2*(1+2*x), where W(k) = 1 + 1/(1 - 2*x/(2*x + (k+1)/(x*(2*k+1))/W(k+1) )), abs(x) < 1/2; (continued fraction). (End)
a(n) = 2^n*Product_{k=0..n-1} (k/n + 1/n)^((-1)^k). - Peter Luschny, Dec 02 2013
G.f.: G(0), where G(k) = 1 + 2*x*(4*k+1)/((2*k+1)*(1+2*x) - (2*k+1)*(4*k+3)*x*(1+2*x)/((4*k+3)*x + (k+1)*(1+2*x)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 19 2014
From Peter Bala, Mar 29 2024: (Start)
a(n) = 2^n * Sum_{k = 0..n} (-1)^(n+k)*binomial(1/2, k)*binomial(- 1/2, n-k) = 2^n * A000246(n)/n!.
a(n) = (1/2^n) * binomial(2*n, n) * hypergeom([-1/2, -n], [1/2 - n], -1). (End)
E.g.f.: BesselI(0, 2*x)*(1 + x*(2 + Pi)*StruveL(1, 2*x)) - Pi*x*BesselI(1, 2*x)*StruveL(0, 2*x). - Stefano Spezia, May 11 2024
a(n) = A089849(n) + A138364(n). - Mélika Tebni, Jun 17 2024
From Amiram Eldar, Aug 15 2025: (Start)
Sum_{n>=0} 1/a(n) = Pi/(3*sqrt(3)) + 2.
Sum_{n>=0} (-1)^n/a(n) = 2/3 + Pi/(9*sqrt(3)). (End)

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

Original entry on oeis.org

1, 1, 5, 49, 729, 14641, 371293, 11390625, 410338673, 16983563041, 794280046581, 41426511213649, 2384185791015625, 150094635296999121, 10260628712958602189, 756943935220796320321, 59938945498865420543457
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Comments

a(n+1) is the number of labeled incomplete ternary trees on n vertices in which each left child has a larger label than its parent. - Brian Drake, Jul 28 2008
Put a(0) = 1. For n > 0, let x(n,k) = 2*cos((2*k-1)*Pi/(2*n+1)), k=1..n. Define the recurrences S(n;0,x(n,k)) = 1, S(n;1,x(n,k)) = x(n,k), S(n;r,x(n,k)) = x(n,k)*S(n;r-1,x(n,k)) - S(n;r-2,x(n,k)), r > 1 an integer, k=1..n. CONJECTURE: For n > 0, a(n) = Product_{k=1..n} (Sum_{m=0..n-1} S(n;2*m,x(n,k))^2). - L. Edson Jeffery, Sep 11 2013
From Wolfdieter Lang, Dec 16 2013: (Start)
Discriminants of the first difference of Chebyshev S-polynomials.
The coefficient table for the first difference polynomials P(n, x) = S(n, x) - S(n-1, x), n >= 0, S(-1, x) = 0, with the Chebyshev S polynomials (see A049310), is given in A130777.
For the discriminant of a polynomial in terms of the square of a determinant of a Vandermonde matrix build from the zeros of the polynomial see, e.g., A127670.
For the proof that D(n) := discriminant(P(n,x)) = (2*n + 1)^(n - 1), n >= 1, use the formula given e.g., in the Rivlin reference, p. 218, Theorem 5.13, eq. (5.3), namely D(n) = (-1)^(n*(n-1)/2)*Product_{j=1..n} P'(n, x(n,j)), with the zeros x(n,j) = -2*cos(2*Pi*j/(2*n+1)) of P(n, x) (see A130777). P'(n, x(n,j)) = (2*n+1)*P(n-1, x(n,j))/(2*sin(Pi*j/(2*n+1))*2*cos(Pi*j/(2*n+1)))^2. P(n-1, x(n,j)) = (-1)^(n+j)*2*cos(Pi*j/(2*n+1)). Product_{j=1..n} 2*sin(Pi*j/(2*n+1)) = 2*n+1 (see the Oct 10 2013 formula in A005408. Product_{j=1..n} 2*cos(Pi*j/(2*n+1)) = 1, because S(2*n, 0) = (-1)^n.
(End)
a(n) is the number of labeled 2-trees with n+2 vertices, rooted at a given edge. - Nikos Apostolakis, Nov 30 2018
a(n) is also the number of 2-trees with n labeled triangles and with a distinguished oriented edge. - Nikos Apostolakis, Dec 14 2018

Examples

			Discriminant: n=4: P(4, x) = 1 + 2*x - 3*x^2 - x^3 + x^4 with the zeros x[1] = -2*cos((2/9)*Pi), x[2] = -2*cos((4/9)*Pi), x[3] = 1, x[4] = 2*cos((1/9)*Pi). D(4) = (Det(Vandermonde(4,[x[1],x[2],x[3],x[4]])))^2 = 729 = a(4). - _Wolfdieter Lang_, Dec 16 2013
		

References

  • L. W. Beineke, and J. W. Moon, Several proofs of the number of labelled 2-dimensional trees, In "Proof Techniques in Graph Theory" (F. Harary editor). Academic Press, New York, 1969, pp. 11-20.
  • Theodore J. Rivlin, Chebyshev polynomials: from approximation theory to algebra and number theory, 2. ed., Wiley, New York, 1990.

Crossrefs

Programs

  • GAP
    List([0..20],n->(2*n+1)^(n-1)); # Muniru A Asiru, Dec 05 2018
  • Magma
    [(2*n+1)^(n-1) : n in [0..20]]; // Wesley Ivan Hurt, Jan 20 2017
    
  • Maple
    spec := [S,{B=Prod(Z,S,S),S=Set(B)},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
  • Mathematica
    max = 16; (Series[Exp[-1/2*ProductLog[-2*x]], {x, 0, max}] // CoefficientList[#, x] & ) * Range[0, max]! (* Jean-François Alcover, Jun 20 2013 *)
    Table[(2*n+1)^(n-1),{n,0,20}] (* Harvey P. Dale, Jul 14 2025 *)
  • PARI
    a(n)=(2*n+1)^(n-1) \\ Charles R Greathouse IV, Nov 20 2011
    
  • PARI
    {a(n)=local(A=1+x);for(i=1,21,A=sqrt(1+2*sum(n=1,21,x^(2*n-1)/(2*n-1)!*A^(4*n-1))+x*O(x^n)));n!*polcoeff(A,n)} \\ Paul D. Hanna, Sep 07 2012
    
  • Python
    for n in range(0, 20): print((2*n + 1)**(n - 1), end=', ') # Stefano Spezia, Dec 01 2018
    

Formula

E.g.f.: exp(-1/2*W(-2*x)), where W is Lambert's W function.
E.g.f. satisfies: A(x) = sqrt(1 + 2*Sum_{n>=1} x^(2*n-1)/(2*n-1)! * A(x)^(4*n-1)). - Paul D. Hanna, Sep 07 2012
E.g.f. satisfies: A(x) = 1/A(-x*A(x)^4). - Paul D. Hanna, Sep 07 2012
a(n) = discriminant of P(n,x) = S(n,x) - S(n-1,x), n >= 1, with the Chebyshev S polynomials from A049310. For the proof see the comment above. a(n) is also the discriminant of S(n,x) + S(n-1,x) = (-1)^n*(S(n,-x) - S(n-1,-x)). - Wolfdieter Lang, Dec 16 2013
From Peter Bala, Dec 19 2013: (Start)
The e.g.f. A(x) = 1 + x + 5*x^2/2! + 49*x^3/3! + 729*x^4/4! + ... satisfies:
1) A(x*exp(-2*x)) = exp(x) = 1/A(-x*exp(2*x));
2) A^2(x) = 1/x*series reversion(x*exp(-2*x));
3) A(x^2) = 1/x*series reversion(x*exp(-x^2));
4) A(x) = exp(x*A(x)^2). (End)
E.g.f.: sqrt(-LambertW(-2*x)/(2*x)). - Vaclav Kotesovec, Dec 07 2014
Related to A001705 by Sum_{n >= 1} a(n)*x^n/n! = series reversion( 1/(1 + x)^2*log(1 + x) ) = series reversion(x - 5*x^2/2! + 26*x^3/3! - 154*x^4/4! + ...). Cf. A000272, A052752, A052774, A052782. - Peter Bala, Jun 15 2016
From Peter Bala, Dec 13 2022: (Start)
The e.g.f. A(x) = 1/x * series reversion of x^2/T(x), where the tree function T(x) = Sum_{n >= 1} n^(n-1)*x^n/n!. See A000169.
For c in C, A(x)^c = 1 + Sum_{n >= 1} c*(2*n + c)^(n-1)*x^n/n!.
First derivative A'(x) = A(x)^3/(1 - 2*x*A(x)^2).
Series reversion of (1 - A(-z)) = -log(1 - z)/(1 - z)^2 is the e.g.f. of A001705.
1/z * series reversion of z/A(z) = 1 + z + 7*z^2/2! + (10^2)*z^3/3! + (13^3)*z^4/4! + ... is the e.g.f. of A052752.
1/z * series reversion of z/A(z^2) = 1 + z^2 + 9*z^4/2! + (13^2)*z^6/3! + (17^3)*z^8/4! + ... = Sum_{n >= 0} A052774(n)*z^(2*n)/n!.
1/z * series reversion of z/A(z^3) = 1 + z^3 + 11*z^6/2! + (16^2)*z^9/3! + (21^3)*z^12/4! + ... = Sum_{n >= 0} A052782(n)*z^(3*n)/n!.
1/z * series reversion of z/A(z)^2 = A(2*z) = 2*Sum_{n >= 0} (4*n + 2)^(n-1)*z^n/n!.
1/z * series reversion of z/A(z)^k = k*Sum_{n >= 0} ((k+2)*n + k)^(n-1)*z^n/n!. (End)
a(n) = Sum_{k=1..n} (-1)^(n-k)*(n+k)^(n-1)*binomial(n,k-1), a(0)=1. - Vladimir Kruchinin, Aug 14 2025

Extensions

Better description from Vladeta Jovovic, Sep 02 2003

A084930 Triangle of coefficients of Chebyshev polynomials T_{2n+1} (x).

Original entry on oeis.org

1, -3, 4, 5, -20, 16, -7, 56, -112, 64, 9, -120, 432, -576, 256, -11, 220, -1232, 2816, -2816, 1024, 13, -364, 2912, -9984, 16640, -13312, 4096, -15, 560, -6048, 28800, -70400, 92160, -61440, 16384, 17, -816, 11424, -71808, 239360, -452608, 487424, -278528, 65536, -19, 1140, -20064, 160512, -695552
Offset: 0

Views

Author

Gary W. Adamson, Jun 12 2003

Keywords

Comments

From Herb Conn, Jan 28 2005: (Start)
"Letting x = 2 Cos 2A, we have the following trigonometric identities:
"Sin 3A = 3*Sin A - 4*Sin^3 A
"Sin 5A = 5*Sin A - 20*Sin^3 A + 16*Sin^5 A
"Sin 7A = 7*Sin A - 56*Sin^3 A + 112*Sin^5 A - 64*Sin^7 A
"Sin 9A = 9*Sin A - 120*Sin^3 A + 432*Sin^5 A - 576*Sin^7 A + 256*Sin^9 A, etc." (End)
Cayley (1876) states "Write sin u = x, then we have sin u = x, [...] sin 3u = 3x - 4x^3, [...] sin 5u = 5x - 20x^3 + 16 x^5, [...]". Since T_n(cos(u)) = cos(nu) for all integer n, sin(u) = cos(u - Pi/2), and sin(u + k*Pi) = (-1)^k sin(u) it follows that T_n(sin(u)) = (-1)^((n-1)/2) sin(nu) for all odd integer n. - Michael Somos, Jun 19 2012
From Wolfdieter Lang, Aug 05 2014: (Start)
The coefficient triangle t(n,k) for the row polynomials Todd(n, x) := T_{2*n+1}(sqrt(x))/sqrt(x) = sum(t(n,k)*x^k, k=0..n) is the Riordan triangle ((1-z)/(1+z)^2, 4*z/(1+z)^2) (rewrite the g.f. for the present triangle a(n,k) given in the formula section). The triangle entries t(n,k) = a(n,k), but the interpretation of the row polynomials is different for both cases.
From the relation Todd(n, x) = S(n, 2*(2*x-1)) - S(n-1, 2*(2*x-1)) with the Chebyshev S-polynomials (see A049310 and the formula section of A130777) follows the recurrence: Todd(n, x) = 2*(-1)^n*(1-x)*Todd(n-1, 1-x) + (2x-1)*Todd(n-1, x), n >= 1, Todd(0, x) = 1.
This corresponds to the triangle recurrence t(n,k) = (2*(k+1)*(-1)^(n-k) - 1)*t(n-1,k) + 2*(1 +(-1)^(n-k))*t(n-1,k-1) + 2*(-1)^(n-k)*sum(binomial(l+1,k)*t(n-1,l), l=k+1..n-1), n >= k >= 1, t(n,k) = 0 if n < k, t(n,0) = (-1)^n*(2*n+1). Compare this with the shorter recurrence involving the rational A-sequence for this Riordan triangle which has g.f. x^2/(2-x-2*sqrt(1-x)). t(n,k) = sum(A(j)*t(n-1,k-1+j), j=0..n-k), n >= k >= 1. The Z-sequence has g.f. -(1 + 2/sqrt(1-x)). For the A- and Z-sequence see a link under A006232. (End)

Examples

			The triangle a(n,k):
n   2n+1\k 0     1      2       3       4        5        6         7        8        9      10 ...
0    1:    1
1    3:   -3     4
2    5:    5   -20     16
3    7:   -7    56   -112      64
4    9:    9  -120    432    -576     256
5   11:  -11   220  -1232    2816   -2816     1024
6   13:   13  -364   2912   -9984   16640   -13312     4096
7   15:  -15   560  -6048   28800  -70400    92160   -61440     16384
8   17:   17  -816  11424  -71808  239360  -452608   487424   -278528    65536
9   19:  -19  1140 -20064  160512 -695552  1770496 -2723840   2490368 -1245184   262144
10  21:   21 -1540  33264 -329472 1793792 -5870592 12042240 -15597568 12386304 -5505024 1048576
... formatted and extended by _Wolfdieter Lang_, Aug 02 2014.
---------------------------------------------------------------------------------------------------
First few polynomials T_{2n+1}(x) are
1*x - 3*x + 4*x^3
5*x - 20*x^3 + 16*x^5
- 7*x + 56*x^3 - 112*x^5 + 64*x^7
9*x - 120*x^3 + 432*x^5 - 576*x^7  + 256*x^9
		

References

  • A. Cayley, On an Expression for 1 +- sin(2p+1)u in Terms of sin u, Messenger of Mathematics, 5 (1876), pp. 7-8 = Mathematical Papers Vol. 10, n. 630, pp. 1-2.
  • Theodore J. Rivlin, Chebyshev polynomials: from approximation theory to algebra and number theory, 2nd ed., Wiley, New York, 1990. p. 37, eq. (1.96) and p. 4, eq. (1.10).

Crossrefs

Cf. A053120 (coefficient triangle of T-polynomials), A127674 (even-indexed T polynomials).
Cf. A127675 (row reversed triangle with different signs).

Programs

  • Mathematica
    row[n_] := (T = ChebyshevT[2*n+1, x]; Coefficient[T, x, #]& /@ Range[1, Exponent[T, x], 2]); Table[row[n], {n, 0, 9} ] // Flatten (* Jean-François Alcover, Aug 06 2014 *)

Formula

Alternate rows of A008310.
a(n,k)=((-1)^(n-k))*(2^(2*k))*binomial(n+1+k,2*k+1)*(2*n+1)/(n+1+k) if n>=k>=0 else 0.
From Wolfdieter Lang, Aug 02 2014: (Start)
a(n,k) = [x^(2*k+1)] T_{2*n+1}(x), n >= k >= 0.
G.f. for row polynomials: x*(1-z)/(1 + 2*(1- 2*x^2)*z + z^2). (End)
The first column sequences are: A157142, 4*(-1)^(n+1)*A000330(n), 16*(-1)^n*A005585(n-1), 64*(-1)^(n+1)*A050486(n-3), 256*(-1)^n*A054333(n-4), ... - Wolfdieter Lang, Aug 05 2014

Extensions

More terms from Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Jun 26 2003
Edited; example rewritten (to conform with the triangle) by Wolfdieter Lang, Aug 02 2014

A187660 Triangle read by rows: T(n,k) = (-1)^(floor(3*k/2))*binomial(floor((n+k)/2),k), 0 <= k <= n.

Original entry on oeis.org

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

Views

Author

L. Edson Jeffery, Mar 12 2011

Keywords

Comments

Conjecture: (i) Let n > 1 and N=2*n+1. Row n of T gives the coefficients of the characteristic polynomial p_N(x)=Sum_{k=0..n} T(n,k)*x^(n-k) of the n X n Danzer matrix D_{N,n-1} = {{0,...,0,1}, {0,...,0,1,1}, ..., {0,1,...,1}, {1,...,1}}. (ii) Let S_0(t)=1, S_1(t)=t and S_r(t)=t*S_(r-1)(t)-S_(r-2)(t), r > 1 (cf. A049310). Then p_N(x)=0 has solutions w_{N,j}=S_(n-1)(phi_{N,j}), where phi_{N,j}=2*(-1)^(j+1)*cos(j*Pi/N), j = 1..n. - L. Edson Jeffery, Dec 18 2011

Examples

			Triangle begins:
  1;
  1,  -1;
  1,  -1,  -1;
  1,  -2,  -1,   1;
  1,  -2,  -3,   1,   1;
  1,  -3,  -3,   4,   1,  -1;
  1,  -3,  -6,   4,   5,  -1,  -1;
  1,  -4,  -6,  10,   5,  -6,  -1,   1;
  1,  -4, -10,  10,  15,  -6,  -7,   1,   1;
  1,  -5, -10,  20,  15, -21,  -7,   8,   1,  -1;
  1,  -5, -15,  20,  35, -21, -28,   8,   9,  -1,  -1;
  1,  -6, -15,  35,  35, -56, -28,  36,   9, -10,  -1,   1;
		

Crossrefs

Signed version of A046854.
Absolute values of a(n) form a reflected version of A065941, which is considered the main entry.

Programs

  • Maple
    A187660 := proc(n,k): (-1)^(floor(3*k/2))*binomial(floor((n+k)/2),k) end: seq(seq(A187660(n,k), k=0..n), n=0..11); # Johannes W. Meijer, Aug 08 2011
  • Mathematica
    t[n_, k_] := (-1)^Floor[3 k/2] Binomial[Floor[(n + k)/2], k]; Table[t[n, k], {n, 0, 11}, {k, 0, n}] (* L. Edson Jeffery, Oct 20 2017 *)

Formula

T(n,k) = (-1)^n*A066170(n,k).
abs(T(n,k)) = A046854(n,k) = abs(A066170(n,k)) = abs(A130777(n,k)).
abs(T(n,k)) = A065941(n,n-k) = abs(A108299(n,n-k)).

Extensions

Edited and corrected by L. Edson Jeffery, Oct 20 2017

A267482 Triangle of coefficients of Gaussian polynomials [2n+1,1]_q represented as finite sum of terms (1+q^2)^k*q^(g-k), where k = 0,1,...,g with g=n.

Original entry on oeis.org

1, 1, 1, -1, 1, 1, -1, -2, 1, 1, 1, -2, -3, 1, 1, 1, 3, -3, -4, 1, 1, -1, 3, 6, -4, -5, 1, 1, -1, -4, 6, 10, -5, -6, 1, 1, 1, -4, -10, 10, 15, -6, -7, 1, 1, 1, 5, -10, -20, 15, 21, -7, -8, 1, 1
Offset: 0

Views

Author

Stephen O'Sullivan, Jan 15 2016

Keywords

Comments

The entry a(n,k), n >= 0, k = 0,1,...,g, where g=n, of this irregular triangle is the coefficient of (1+q^2)^k*q^(g-k) in the representation of the Gaussian polynomial [2n+1,1]q = Sum{k=0..g) a(n,k)*(1+q^2)^k*q^(g-k).
The sequence arises in the formal derivation of the stability polynomial B(x) = Sum_{i=0..N} d_i T(iM,x) of rank N, and degree L, where T(iM,x) denotes the Chebyshev polynomial of the first kind of degree iM. The coefficients d_i are determined by order conditions on the stability polynomial.
Conjecture: More generally, the Gaussian polynomial [2*n+m+1-(m mod 2),m]q = Sum{k=0..g(m;n)} a(m;n,k)*(1+q^2)^k*q^(g(m;n)-k), for m >= 0, n >= 0, where g(m;n) = m*n if m is odd and (2*n+1)*m/2 if m is even, and the tabf array entries a(m;n,k) are the coefficients of the g.f. for the row n polynomials G(m;n,x) = (d^m/dt^m)G(m;n,t,x)/m!|{t=0}, with G(m;n,t,x) = (1+t)*Product{k=1..n+(m - m (mod 2))/2}(1 + t^2 + 2*t*T(k,x/2) (Chebyshev's T-polynomials). Hence a(m;n,k) = [x^k]G(m;n,x), for k=0..g(m;n). The present entry is the instance m = 2. (Thanks to Wolfdieter Lang for clarifying the text on the general prescription of a(m;n,k).)
Signed version of A046854, A130777.
Conjecture: row n is U(n, x/2) + U(n-1, x/2) where U is the sequence of Chebyshev polynomials of the second kind. - Thomas Baruchel, Jun 03 2018 [For a proof see the following comment.]
From Wolfdieter Lang, Oct 19 2019: (Start)
The row polynomial R(n, x) = Sum_{k=0..n} a(n, k)*x^k = [2*n+1]_q / q^n with the q-number [2*n+1]_q := (1 - q^n)/(1 - q), which for q = 1 becomes 2*n+1, and x = x(q) = q + q^(-1). See the simplified Name and the first comment. In terms of Chebyshev S polynomials (A049310) this q-number is written as [2*n+1]_q = q^n*S(2*n, q^(1/2) + q^(-1/2)), hence R(n, x) = S(2*n, sqrt(2+x)) = S(n, x) + S(n-1, x) (which proves the conjecture of the previous comment).
For the o.g.f. of R(n, x) see the formula section.
My motivation for looking at this sequence came from the Brändli and Beyne paper's recurrence for the polynomial P_m(s) which coincides with R(n, x), with m -> n and s -> x. (End)
A294099(n, k) = Sum_{j=0..k} n^j * T(n, j) for all n, k in Z. - Michael Somos, Jun 19 2023

Examples

			Triangle begins:
   1;
   1,   1;
  -1,   1,   1;
  -1,  -2,   1,   1;
   1,  -2,  -3,   1,   1;
   1,   3,  -3,  -4,   1,   1;
  -1,   3,   6,  -4,  -5,   1,   1;
  -1,  -4,   6,  10,  -5,  -6,   1,   1;
   1,  -4, -10,  10,  15,  -6,  -7,   1,   1;
   1,   5, -10, -20,  15,  21,  -7,  -8,   1,   1;
		

Crossrefs

Programs

  • Maple
    A267482 := proc (n, k) local y: y := expand(subs(t = 0, diff((1+t)*product(1+t^2+2*t*ChebyshevT(i, x/2), i = 1 .. n),t))): if k = 0 then subs(x = 0, y) else subs(x = 0, diff(y, x$k)/k!) end if: end proc: seq(seq(A267482(n, k), k = 0 .. n), n = 0 .. 20);
  • Mathematica
    row[n_] := D[(1+t)*Product[1+t^2+2*t*ChebyshevT[i, x/2], {i, 1, n}], t] /. t -> 0 // CoefficientList[#, x]&; Table[row[n], {n, 0, 20}] // Flatten (* Jean-François Alcover, Jan 16 2016 *)
  • PARI
    T(n,k) = (-1)^((n-k)\2)*binomial((n+k)\2, k); \\François Marques, Sep 28 2021

Formula

G.f. for row polynomial: G(n,x) = (d^2/dt^2)((1+t)*Product_{i=1..n+1}(1+t^2+2t*T(i,x/2)))|_{t=0}.
From Wolfdieter Lang, Oct 19 2019: (Start)
Row polynomial R(n, x) = S(2*n, sqrt(2+x)) = S(n, x) + S(n-1, x) = Sum_{k=0..n} (-1)^k*binomial(2*n-k, k)*(2 + x)^(n-k), for n >= 0. (See the Thomas Baruchel conjecture and the proof above.) For the S(n, x) coefficients see A049310.
R(n, x) = Sum_{j=0} (-1)^e(n,j)*binomial(e(n,j) + j, j)*x^j*, with e(n,j) := floor((n-j)/2). See eq. (12) of the Brändli and Beyne paper.
G.f. for row polynomials R(n, x) (that is of the triangle): G(x,z) = (1 + z)/(1 - x*z + z^2).
Recurrence for R(n, x): R(-1, x) = -1, R(0, x) = 1, R(n, x) = x*R(n-1, x) - R(n-2, x), for n >= 1. (See the Brändli and Beyne link, polynomials P_m(s) in Definition 6.)
(End)
T(n,k) = (-1)^(floor((n-k)/2))*binomial(floor((n+k)/2), k). - François Marques, Sep 28 2021
Showing 1-10 of 14 results. Next