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

A029635 The (1,2)-Pascal triangle (or Lucas triangle) read by rows.

Original entry on oeis.org

2, 1, 2, 1, 3, 2, 1, 4, 5, 2, 1, 5, 9, 7, 2, 1, 6, 14, 16, 9, 2, 1, 7, 20, 30, 25, 11, 2, 1, 8, 27, 50, 55, 36, 13, 2, 1, 9, 35, 77, 105, 91, 49, 15, 2, 1, 10, 44, 112, 182, 196, 140, 64, 17, 2, 1, 11, 54, 156, 294, 378, 336, 204, 81, 19, 2, 1, 12, 65, 210, 450, 672, 714, 540, 285, 100
Offset: 0

Views

Author

Keywords

Comments

This is also called Vieta's array. - N. J. A. Sloane, Nov 22 2017
Dropping the first term and changing the boundary conditions to T(n,1)=n, T(n,n-1)=2 (n>=2), T(n,n)=1 yields the number of nonterminal symbols (which generate strings of length k) in a certain context-free grammar in Chomsky normal form that generates all permutations of n symbols. Summation over k (1<=k<=n) results in A003945. For the number of productions of this grammar: see A090327. Example: 1; 2, 1; 3, 2, 1; 4, 5, 2, 1; 5, 9, 7, 2, 1; 6, 14, 16, 9, 2, 1; In addition to the example of A090327 we have T(3,3)=#{S}=1, T(3,2)=#{D,E}=2 and T(3,1)=#{A,B,C}=3. - Peter R. J. Asveld, Jan 29 2004
Much as the original Pascal triangle gives the Fibonacci numbers as sums of its diagonals, this triangle gives the Lucas numbers (A000032) as sums of its diagonals; see Posamentier & Lehmann (2007). - Alonso del Arte, Apr 09 2012
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 04 2013
It appears that for the infinite set of (1,N) Pascal's triangles, the binomial transform of the n-th row (n>0), followed by zeros, is equal to the n-th partial sum of (1, N, N, N, ...). Example: for the (1,2) Pascal's triangle, the binomial transform of the second row followed by zeros, i.e., of (1, 3, 2, 0, 0, 0, ...), is equal to the second partial sum of (1, 2, 2, 2, ...) = (1, 4, 9, 16, ...). - Gary W. Adamson, Aug 11 2015
Given any (1,N) Pascal triangle, let the binomial transform of the n-th row (n>1) followed by zeros be Q(x). It appears that the binomial transform of the (n-1)-th row prefaced by a zero is Q(n-1). Example: In the (1,2) Pascal triangle the binomial transform of row 3: (1, 4, 5, 2, 0, 0, 0, ...) is A000330 starting with 1: (1, 5, 14, 30, 55, 91, ...). The binomial transform of row 2 prefaced by a zero and followed by zeros, i.e., of (0, 1, 3, 2, 0, 0, 0, ...) is (0, 1, 5, 14, 30, 55, ...). - Gary W. Adamson, Sep 28 2015
It appears that in the array accompanying each (1,N) Pascal triangle (diagonals of the triangle), the binomial transform of (..., 1, N, 0, 0, 0, ...) preceded by (n-1) zeros generates the n-th row of the array (n>0). Then delete the zeros in the result. Example: in the (1,2) Pascal triangle, row 3 (1, 5, 14, 30, ...) is the binomial transform of (0, 0, 1, 2, 0, 0, 0, ...) with the resulting zeros deleted. - Gary W. Adamson, Oct 11 2015
Read as a square array (similar to the Example section Sq(m,j), but with Sq(0,0)=0 and Sq(m,j)=P(m+1,j) otherwise), P(n,k) are the multiplicities of the eigenvalues, lambda_n = n(n+k-1), of the Laplacians on the unit k-hypersphere, given by Teo (and Choi) as P(n,k) = (2n-k+1)(n+k-2)!/(n!(k-1)!). P(n,k) is also the numerator of a Dirichlet series for the Minakashisundarum-Pleijel zeta function for the sphere. Also P(n,k) is the dimension of the space of homogeneous, harmonic polynomials of degree k in n variables (Shubin, p. 169). For relations to Chebyshev polynomials and simple Lie algebras, see A034807. - Tom Copeland, Jan 10 2016
For a relation to a formulation for a universal Lie Weyl algebra for su(1,1), see page 16 of Durov et al. - Tom Copeland, Jan 15 2016

Examples

			Triangle begins:
  [0] [2]
  [1] [1, 2]
  [2] [1, 3,  2]
  [3] [1, 4,  5,  2]
  [4] [1, 5,  9,  7,   2]
  [5] [1, 6, 14, 16,   9,  2]
  [6] [1, 7, 20, 30,  25, 11,  2]
  [7] [1, 8, 27, 50,  55, 36, 13,  2]
  [8] [1, 9, 35, 77, 105, 91, 49, 15, 2]
.
Read as a square, the array begins:
  n\k| 0  1   2    3    4    5
  --------------------------------------
  0 |  2  2   2    2    2    2   A040000
  1 |  1  3   5    7    9   11   A005408
  2 |  1  4   9   16   25   36   A000290
  3 |  1  5  14   30   55   91   A000330
  4 |  1  6  20   50  105  196   A002415
  5 |  1  7  27   77  182  378   A005585
  6 |  1  8  35  112  294  672   A040977
		

References

  • Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8.
  • Alfred S. Posamentier & Ingmar Lehmann, The (Fabulous) Fibonacci Numbers. New York: Prometheus Books (2007): 97 - 105.
  • M. Shubin and S. Andersson, Pseudodifferential Operators and Spectral Theory, Springer Series in Soviet Mathematics, 1987.

Crossrefs

Cf. A003945 (row sums), A007318, A034807, A061896, A029653 (row-reversed), A157000.
Sums along ascending antidiagonals give Lucas numbers, n>0.

Programs

  • Haskell
    a029635 n k = a029635_tabl !! n !! k
    a029635_row n = a029635_tabl !! n
    a029635_tabl = [2] : iterate
       (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [1,2]
    -- Reinhard Zumkeller, Mar 12 2012, Feb 23 2012
    
  • Maple
    T := proc(n, k) option remember;
    if n = k then 2 elif k = 0 then 1 else T(n-1, k-1) + T(n-1, k) fi end:
    for n from 0 to 8 do seq(T(n, k), k = 0..n) od;  # Peter Luschny, Dec 22 2024
  • Mathematica
    t[0, 0] = 2; t[n_, k_] := If[k < 0 || k > n, 0, Binomial[n, k] + Binomial[n-1, k-1]]; Flatten[Table[t[n, k], {n, 0, 11}, {k, 0, n}]] (* Jean-François Alcover, May 03 2011 *)
    (* The next program cogenerates A029635 and A029638. *)
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := u[n - 1, x] + v[n - 1, x]
    v[n_, x_] := x*u[n - 1, x] + x*v[n - 1, x] + 1
    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[%]   (* A029638  *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]   (* A029635 *)
    (* Clark Kimberling, Feb 20 2012 *)
    Table[Binomial[n,k]+Binomial[n-1,k-1],{n,0,20},{k,0,n}]//Flatten (* Harvey P. Dale, Feb 08 2024 *)
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, (n==0) + binomial(n, k) + binomial(n-1, k-1))}; /* Michael Somos, Jul 15 2003 */
    
  • Sage
    # uses[riordan_array from A256893]
    riordan_array((2-x)/(1-x), x/(1-x), 8) # Peter Luschny, Nov 09 2019

Formula

From Henry Bottomley, Apr 26 2002; (Start)
T(n, k) = T(n-1, k-1) + T(n-1, k).
T(n, k) = C(n, k) + C(n-1, k-1).
T(n, k) = C(n, k)*(n + k)/n.
T(n, k) = A007318(n, k) + A007318(n-1, k-1).
T(n, k) = A061896(n + k, k) but with T(0, 0) = 1 and T(1, 1) = 2.
Row sum is floor(3^2(n-1)) i.e., A003945. (End)
G.f.: 1 + (1 + x*y) / (1 - x - x*y). - Michael Somos, Jul 15 2003
G.f. for n-th row: (x+2*y)*(x+y)^(n-1).
O.g.f. for row n: (1+x)/(1-x)^(n+1). The entries in row n are the nonzero entries in column n of A053120 divided by 2^(n-1). - Peter Bala, Aug 14 2008
T(2n, n) - T(2n, n+1)= Catalan(n)= A000108(n). - Philippe Deléham, Mar 19 2009
With T(0, 0) = 1 : Triangle T(n, k), read by rows, given by [1,0,0,0,0,0,...] DELTA [2,-1,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Oct 10 2011
With T(0, 0) = 1, as in the Example section below, this is known as Vieta's array. The LU factorization of the square array is given by Yang and Leida, equation 20. - Peter Bala, Feb 11 2012
For n > 0: T(n, k) = A097207(n-1, k), 0 <= k < n. - Reinhard Zumkeller, Mar 12 2012
For n > 0: T(n, k) = A029600(n, k) - A007318(n, k), 0 <= k <= n. - Reinhard Zumkeller, Apr 16 2012
Riordan array ((2-x)/(1-x), x/(1-x)). - Philippe Deléham, Mar 15 2013
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(1 + 4*x + 5*x^2/2! + 2*x^3/3!) = 1 + 5*x + 14*x^2/2! + 30*x^3/3! + 55*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ). - Peter Bala, Dec 22 2014
For n>=1: T(n, 0) + T(n, 1) + T(n, 2) = A000217(n+1). T(n, n-2) = (n-1)^2. - Bob Selcoe, Mar 29 2016:

Extensions

More terms from David W. Wilson
a(0) changed to 2 (was 1) by Daniel Forgues, Jul 06 2010

A003558 Least number m > 0 such that 2^m == +-1 (mod 2n + 1).

Original entry on oeis.org

1, 1, 2, 3, 3, 5, 6, 4, 4, 9, 6, 11, 10, 9, 14, 5, 5, 12, 18, 12, 10, 7, 12, 23, 21, 8, 26, 20, 9, 29, 30, 6, 6, 33, 22, 35, 9, 20, 30, 39, 27, 41, 8, 28, 11, 12, 10, 36, 24, 15, 50, 51, 12, 53, 18, 36, 14, 44, 12, 24, 55, 20, 50, 7, 7, 65, 18, 36, 34, 69, 46
Offset: 0

Views

Author

Keywords

Comments

Multiplicative suborder of 2 (mod 2n+1) (or sord(2, 2n+1)).
This is called quasi-order of 2 mod b, with b = 2*n+1, for n >= 1, in the Hilton/Pederson reference.
For the complexity of computing this, see A002326.
Also, the order of the so-called "milk shuffle" of a deck of n cards, which maps cards (1,2,...,n) to (1,n,2,n-1,3,n-2,...). See the paper of Lévy. - Jeffrey Shallit, Jun 09 2019
It appears that under iteration of the base-n Kaprekar map, for even n > 2 (A165012, A165051, A165090, A151949 in bases 4, 6, 8, 10), almost all cycles are of length a(n/2 - 1); proved under the additional constraint that the cycle contains at least one element satisfying "number of digits (n-1) - number of digits 0 = o(total number of digits)". - Joseph Myers, Sep 05 2009
From Gary W. Adamson, Sep 20 2011: (Start)
a(n) can be determined by the cycle lengths of iterates using x^2 - 2, seed 2*cos(2*Pi/N); as shown in the A065941 comment of Sep 06 2011. The iterative map of the logistic equation 4x*(1-x) is likewise chaotic with the same cycle lengths but initiating the trajectory with sin^2(2*Pi/N), N = 2n+1 [Kappraff & Adamson, 2004]. Chaotic terms with the identical cycle lengths can be obtained by applying Newton's method to i = sqrt(-1) [Strang, also Kappraff and Adamson, 2003], resulting in the morphism for the cot(2*Pi/N) trajectory: (x^2-1)/2x. (End)
From Gary W. Adamson, Sep 11 2019: (Start)
Using x^2 - 2 with seed 2*cos(Pi/7), we obtain the period-three trajectory 1.8019377...-> 1.24697...-> -0.445041... For an odd prime N, the trajectory terms represent diagonal lengths of regular star 2N-gons, with edge the shortest value (0.445... in this case.) (Cf. "Polygons and Chaos", p. 9, Fig 4.) We can normalize such lengths by dividing through with the lowest value, giving 3 diagonals of the 14-gon: (1, 2.801937..., 4.048917...). Label the terms ranked in magnitude with odd integers (1, 3, 5), and we find that the diagonal lengths are in agreement with the diagonal formula (sin(j*Pi)/14)/(sin(Pi/14)), with j = (1,3,5). (End)
Roots of signed n-th row A054142 polynomials are chaotic with respect to the operation (-2, x^2), with cycle lengths a(n). Example: starting with a root to x^3 - 5x^2 + 6x - 1 = 0; (2 + 2*cos(2*Pi/N) = 3.24697...); we obtain the trajectory (3.24697...-> 1.55495...-> 0.198062...); the roots to the polynomial with cycle length 3 matching a(3) = 3. - Gary W. Adamson, Sep 21 2011
From Juhani Heino, Oct 26 2015: (Start)
Start a sequence with numbers 1 and n. For next numbers, add previous numbers going backwards until the sum is even. Then the new number is sum/2. I conjecture that the sequence returns to 1,n and a(n) is the cycle length.
For example:
1,7,4,2,1,7,... so a(7) = 4.
1,6,3,5,4,2,1,6,... so a(6) = 6. (End)
From Juhani Heino, Nov 06 2015: (Start)
Proof of the above conjecture: Let n = -1/2; thus 2n + 1 = 0, so operations are performed mod (2n + 1). When the member is even, it is divided by 2. When it is odd, multiply by n, so effectively divide by -2. This is all well-defined in the sense that new members m are 1 <= m <= n. Now see what happens starting from an odd member m. The next member is -m/2. As long as there are even members, divide by 2 and end up with an odd -m/(2^k). Now add all the members starting with m. The sum is m/(2^k). It's divided by 2, so the next member is m/(2^(k+1)). That is the same as (-m/(2^k))/(-2), as with the definition.
So actually start from 1 and always divide by 2, although the sign sometimes changes. Eventually 1 is reached again. The chain can be traversed backwards and then 2^(cycle length) == +-1 (mod 2n + 1).
To conclude, we take care of a(0): sequence 1,0 continues with zeros and never returns to 1. So let us declare that cycle length 0 means unavailable. (End)
From Gary W. Adamson, Aug 20 2019: (Start)
Terms in the sequence can be obtained by applying the doubling sequence mod (2n + 1), then counting the terms until the next term is == +1 (mod 2n + 1). Example: given 25, the trajectory is (1, 2, 4, 8, 16, 7, 14, 3, 6, 12).
The cycle ends since the next term is 24 == -1 (mod 25) and has a period of 10. (End)
From Gary W. Adamson, Sep 04 2019: (Start)
Conjecture of Kappraff and Adamson in "Polygons and Chaos", p. 13 Section 7, "Chaos and Number": Given the cycle length for N = 2n + 1, the same cycle length is present in bases 4, 9, 16, 25, ..., m^2, for the expansion of 1/N.
Examples: The cycle length for 7 is 3, likewise for 1/7 in base 4: 0.021021021.... In base 9 the expansion of 1/7 is 0.125125125... Check: The first few terms are 1/9 + 2/81 + 5/729 = 104/279 = 0.1426611... (close to 1/7 = 0.142857...). (End)
From Gary W. Adamson, Sep 24 2019: (Start)
An exception to the rule for 1/N in bases m^2: (when N divides m^2 as in 1/7 in base 49, = 7/49, rational). When all terms in the cycle are the same, the identity reduces to 1/N in (some bases) = .a, a, a, .... The minimal values of "a" for 1/N are provided as examples, with the generalization 1/N in base (N-1)^2 = .a, a, a, ... for N odd:
1/3 in base 4 = .1, 1, 1, ...
1/5 in base 16 = .3, 3, 3, ...
1/7 in base 36 = .5, 5, 5, ...
1/9 in base 64 = .7, 7, 7, ...
1/11 in base 100 = .9, 9, 9, ... (Check: the first three terms are 9/100 +9/(100^2) + 9/(100^3) = 0.090909 where 1/11 = 0.09090909...). (End)
For N = 2n+1, the corresponding entry is equal to the degree of the polynomial for N shown in (Lang, Table 2, p. 46). As shown, x^3 - 3x - 1 is the minimal polynomial for N = 9, with roots (1.87938..., -1.53208..., 0.347296...); matching the (abs) values of the 2*cos(Pi/9) trajectory using x^2 - 2. Thus, a(4) = 3. If N is prime, the polynomials shown in Table 2 are the same as those for the same N in A065941. If different, the minimal polynomials shown in Table 2 are factors of those in A065941. - Gary W. Adamson, Oct 01 2019
The terms in the 2*cos(Pi/N) trajectory (roots to the minimal polynomials in A187360 and (Lang)), are quickly obtained from the doubling trajectory (mod N) by using the operation L(m) 2*cos(x)--> 2*cos(m*x), where L(2), the second degree Lucas polynomial (A034807) is x^2 - 2. Relating to the heptagon and using seed 2*cos(Pi/7), we obtain the trajectory 1.8019..., 1.24697..., and 0.445041....; cyclic with period 3. All such roots can be derived from the N-th roots of Unity and can be mapped on the Vesica Piscis. Given the roots of Unity (Polar 1Angle(k*2*Pi/N), k = 1, 2, ..., (N-1)/2) the Vesica Piscis maps these points on the left (L) circle to the (R) circle by adding 1A(0) or (a + b*I) = (1 + 0i). But this operation is the same as vector addition in which the resultant vector is 1 + 1A(k*(2*Pi/N)). Example: given the radius at 2*Pi/7 on the left circle, this maps to (1 + 1A(2*Pi/7)) on the right circle; or 1A(2*Pi/7) --> (1.8019377...A(Pi/7). Similarly, 1 + 1A((2)*2*Pi/7)) maps to (1.24697...A (2*Pi/7); and 1 + 1A(3*2*Pi/7) maps to (.0445041...A(3*Pi/7). - Gary W. Adamson, Oct 23 2019
From Gary W. Adamson, Dec 01 2021: (Start)
As to segregating the two sets: (A014659 terms are those N = (2*n+1), N divides (2^m - 1), and (A014657 terms are those N that divide (2^m + 1)); it appears that the following criteria apply: Given IcoS(N, 1) (cf. Lang link "On the Equivalence...", p. 16, Definition 20), if the number of odd terms is odd, then N belongs to A014659, otherwise A014657. In IcoS(11, 1): (1, 2, 4, 3, 5), three odd terms indicate that 11 is a term in A014657. IcoS(15, 1) has the orbit (1, 2, 4, 7) with two odd terms indicating that 15 is a term in A014659.
It appears that if sin(2^m * Pi/N) has a negative sign, then N is in A014659; otherwise N is in A014657. With N = 15, m is 4 and sin(16 * Pi/15) is -0.2079116... If N is 11, m is 5 and sin(32 * Pi/11) is 0.2817325. (End)
On the iterative map using x^2 - 2, (Devaney, p. 126) states that we must find the function that takes 2*cos(Pi) -> 2*cos(2*Pi). "However, we may write 2*cos(2*Pi) = 2*(2*cos^2(Pi) - 1) = (2*cos(Pi))^2 - 2. So the required function is x^2 - 2." On the period 3 implies chaos theorem of James Yorke and T.Y. Li, proved in 1975; Devaney (p. 133) states that if F is continuous and we find a cycle of period 3, there are infinitely many other cycles for this map with every possible period. Check: The x^2 - 2 orbit for 7 has a period of 3, so this entry has periodic points of all other periods. - Gary W. Adamson, Jan 04 2023
It appears that a(n) is the length of the cycle starting at 2/(2*n+1) for the map x->1 - abs(2*x-1). - Michel Marcus, Jul 16 2025

Examples

			a(3) = 3 since f(x) = x^2 - 2 has a period of 3 using seed 2*cos(2*Pi/7), where 7 = 2*3 + 1.
a(15) = 5 since the iterative map of the logistic equation 4x*(1-x) has a period 5 using seed sin^2(2*Pi)/N; N = 31 = 2*15 + 1.
		

References

  • Peter Hilton and Jean Pedersen, A Mathematical Tapestry: Demonstrating the Beautiful Unity of Mathematics, Cambridge University Press, 2010, pp. 261-264.
  • Carl Schick, Trigonometrie und unterhaltsame Zahlentheorie, Bokos Druck, Zürich, 2003 (ISBN 3-9522917-0-6).
  • Robert L. Devaney, A First Course in Chaotic Dynamical Systems, Theory and Experiment; Perseus Books Publishing, 1992, pp. 121-126.

Crossrefs

Cf. A054142, A065941, A085478, A160657, A179480, A135303 (coach numbers), A216371 (odd primes with one coach), A000215 (Fermat numbers).
A216066 is an essentially identical sequence apart from the offset.
Cf. A329593, A332433 (signs).

Programs

  • Maple
    A003558 := proc(n)
        local m,mo ;
        if n = 0 then
            return 0 ;
        end if;
        for m from 1 do
            mo := modp(2^m,2*n+1) ;
            if mo in {1,2*n} then
                return m;
            end if;
        end do:
    end proc:
    seq(A003558(n),n=0..20) ; # R. J. Mathar, Dec 01 2014
    f:= proc(n) local t;
          t:= numtheory:-mlog(-1,2,n);
          if t = FAIL then numtheory:-order(2,n) else t fi
    end proc:
    0, seq(f(2*k+1),k=1..1000); # Robert Israel, Oct 26 2015
  • Mathematica
    Suborder[a_,n_]:=If[n>1&&GCD[a,n]==1,Min[MultiplicativeOrder[a,n,{-1,1}]],0];
    Join[{1},Table[Suborder[2,2n+1],{n,100}]] (* T. D. Noe, Aug 02 2006 *) (* revised by Vincenzo Librandi, Apr 11 2020 *)
  • PARI
    a(n) = {m=1; while(m, if( (2^m) % (2*n+1) == 1 || (2^m) % (2*n+1) == 2*n, return(m)); m++)} \\ Altug Alkan, Nov 06 2015
    
  • PARI
    isok(m, n) = my(md = Mod(2, 2*n+1)^m); (md==1) || (md==-1);
    A003558(n) = my(m=1); while(!isok(m,n) , m++); m; \\ Michel Marcus, May 06 2020
    
  • Python
    def A003558(n):
        m, k = 1, 2 % (c:=(r:=n<<1)+1)
        while not (k==1 or k==r):
            k = 2*k%c
            m += 1
        return m # Chai Wah Wu, Oct 09 2023

Formula

a(n) = log_2(A160657(n) + 2) - 1. - Nathaniel Johnston, May 22 2009
a(n-1) = card {cos((2^k)*Pi/(2*n-1)): k in N} for n >= 1 (see A216066, an essentially identical sequence, for more information). - Roman Witula, Sep 01 2012
a(n) <= n. - Charles R Greathouse IV, Sep 15 2012 [For n >= 1]
a(n) = min{k > 0 | q_k = q_0} where q_0 = 1 and q_k = |2*n+1 - 2*q_{k-1}| (cf. [Schick, p. 4]; q_k=1 for n=1; q_k=A010684(k) for n=2; q_k=A130794(k) for n=3; q_k=|A154870(k-1)| for n=4; q_k=|A135449(k)| for n=5.) - Jonathan Skowera, Jun 29 2013
2^(a(n)) == A332433(n) (mod (2*n+1)), and (2^(a(n)) - A332433(n))/(2*n+1) = A329593(n), for n >= 0. - Wolfdieter Lang, Apr 09 2020

Extensions

More terms from Harry J. Smith, Feb 11 2005
Entry revised by N. J. A. Sloane, Aug 02 2006 and again Dec 10 2017

A003500 a(n) = 4*a(n-1) - a(n-2) with a(0) = 2, a(1) = 4.

Original entry on oeis.org

2, 4, 14, 52, 194, 724, 2702, 10084, 37634, 140452, 524174, 1956244, 7300802, 27246964, 101687054, 379501252, 1416317954, 5285770564, 19726764302, 73621286644, 274758382274, 1025412242452, 3826890587534, 14282150107684, 53301709843202, 198924689265124
Offset: 0

Views

Author

Keywords

Comments

a(n) gives values of x satisfying x^2 - 3*y^2 = 4; corresponding y values are given by 2*A001353(n).
If M is any given term of the sequence, then the next one is 2*M + sqrt(3*M^2 - 12). - Lekraj Beedassy, Feb 18 2002
For n > 0, the three numbers a(n) - 1, a(n), and a(n) + 1 form a Fleenor-Heronian triangle, i.e., a Heronian triangle with consecutive sides, whose area A(n) may be obtained from the relation [4*A(n)]^2 = 3([a(2n)]^2 - 4); or A(n) = 3*A001353(2*n)/2 and whose semiperimeter is 3*a[n]/2. The sequence is symmetrical about a[0], i.e., a[-n] = a[n].
For n > 0, a(n) + 2 is the number of dimer tilings of a 2*n X 2 Klein bottle (cf. A103999).
Tsumura shows that, for prime p, a(p) is composite (contrary to a conjecture of Juricevic). - Charles R Greathouse IV, Apr 13 2010
Except for the first term, positive values of x (or y) satisfying x^2 - 4*x*y + y^2 + 12 = 0. - Colin Barker, Feb 04 2014
Except for the first term, positive values of x (or y) satisfying x^2 - 14*x*y + y^2 + 192 = 0. - Colin Barker, Feb 16 2014
A268281(n) - 1 is a member of this sequence iff A268281(n) is prime. - Frank M Jackson, Feb 27 2016
a(n) gives values of x satisfying 3*x^2 - 4*y^2 = 12; corresponding y values are given by A005320. - Sture Sjöstedt, Dec 19 2017
Middle side lengths of almost-equilateral Heronian triangles. - Wesley Ivan Hurt, May 20 2020
For all elements k of the sequence, 3*(k-2)*(k+2) is a square. - Davide Rotondo, Oct 25 2020

References

  • B. C. Berndt, Ramanujan's Notebooks Part IV, Springer-Verlag, see p. 82.
  • J. M. Borwein and P. B. Borwein, Pi and the AGM, Wiley, 1987, p.91.
  • Michael P. Cohen, Generating Heronian Triangles With Consecutive Integer Sides. Journal of Recreational Mathematics, vol. 30 no. 2 1999-2000 p. 123.
  • L. E. Dickson, History of The Theory of Numbers, Vol. 2 pp. 197;198;200;201. Chelsea NY.
  • Charles R. Fleenor, Heronian Triangles with Consecutive Integer Sides, Journal of Recreational Mathematics, Volume 28, no. 2 (1996-7) 113-115.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley. Enumerative combinatorics. Vol. 2, volume 62 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 1999.
  • V. D. To, "Finding All Fleenor-Heronian Triangles", Journal of Recreational Mathematics vol. 32 no.4 2003-4 pp. 298-301 Baywood NY.

Crossrefs

Cf. A011945 (areas), A334277 (perimeters).
Cf. this sequence (middle side lengths), A016064 (smallest side lengths), A335025 (largest side lengths).

Programs

  • Haskell
    a003500 n = a003500_list !! n
    a003500_list = 2 : 4 : zipWith (-)
       (map (* 4) $ tail a003500_list) a003500_list
    -- Reinhard Zumkeller, Dec 17 2011
    
  • Magma
    I:=[2,4]; [n le 2 select I[n] else 4*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 14 2018
  • Maple
    A003500 := proc(n) option remember; if n <= 1 then 2*n+2 else 4*procname(n-1)-procname(n-2); fi;
    end proc;
  • Mathematica
    a[0]=2; a[1]=4; a[n_]:= a[n]= 4a[n-1] -a[n-2]; Table[a[n], {n, 0, 23}]
    LinearRecurrence[{4,-1},{2,4},30] (* Harvey P. Dale, Aug 20 2011 *)
    Table[Round@LucasL[2n, Sqrt[2]], {n, 0, 20}] (* Vladimir Reshetnikov, Sep 15 2016 *)
  • PARI
    x='x+O('x^99); Vec(-2*(-1+2*x)/(1-4*x+x^2)) \\ Altug Alkan, Apr 04 2016
    
  • Sage
    [lucas_number2(n,4,1) for n in range(0, 24)] # Zerinvary Lajos, May 14 2009
    

Formula

a(n) = ( 2 + sqrt(3) )^n + ( 2 - sqrt(3) )^n.
a(n) = 2*A001075(n).
G.f.: 2*(1 - 2*x)/(1 - 4*x + x^2). Simon Plouffe in his 1992 dissertation.
a(n) = A001835(n) + A001835(n+1).
a(n) = trace of n-th power of the 2 X 2 matrix [1 2 / 1 3]. - Gary W. Adamson, Jun 30 2003 [corrected by Joerg Arndt, Jun 18 2020]
From the addition formula, a(n+m) = a(n)*a(m) - a(m-n), it is easy to derive multiplication formulas, such as: a(2*n) = (a(n))^2 - 2, a(3*n) = (a(n))^3 - 3*(a(n)), a(4*n) = (a(n))^4 - 4*(a(n))^2 + 2, a(5*n) = (a(n))^5 - 5*(a(n))^3 + 5*(a(n)), a(6*n) = (a(n))^6 - 6*(a(n))^4 + 9*(a(n))^2 - 2, etc. The absolute values of the coefficients in the expansions are given by the triangle A034807. - John Blythe Dobson, Nov 04 2007
a(n) = 2*A001353(n+1) - 4*A001353(n). - R. J. Mathar, Nov 16 2007
From Peter Bala, Jan 06 2013: (Start)
Let F(x) = Product_{n=0..infinity} (1 + x^(4*n + 1))/(1 + x^(4*n + 3)). Let alpha = 2 - sqrt(3). This sequence gives the simple continued fraction expansion of 1 + F(alpha) = 2.24561 99455 06551 88869 ... = 2 + 1/(4 + 1/(14 + 1/(52 + ...))). Cf. A174500.
Also F(-alpha) = 0.74544 81786 39692 68884 ... has the continued fraction representation 1 - 1/(4 - 1/(14 - 1/(52 - ...))) and the simple continued fraction expansion 1/(1 + 1/((4 - 2) + 1/(1 + 1/((14 - 2) + 1/(1 + 1/((52 - 2) + 1/(1 + ...))))))).
F(alpha)*F(-alpha) has the simple continued fraction expansion 1/(1 + 1/((4^2 - 4) + 1/(1 + 1/((14^2 - 4) + 1/(1 + 1/((52^2 - 4) + 1/(1 + ...))))))).
(End)
a(2^n) = A003010(n). - John Blythe Dobson, Mar 10 2014
a(n) = [x^n] ( (1 + 4*x + sqrt(1 + 8*x + 12*x^2))/2 )^n for n >= 1. - Peter Bala, Jun 23 2015
E.g.f.: 2*exp(2*x)*cosh(sqrt(3)*x). - Ilya Gutkovskiy, Apr 27 2016
a(n) = Sum_{k=0..floor(n/2)} (-1)^k*n*(n - k - 1)!/(k!*(n - 2*k)!)*4^(n - 2*k) for n >= 1. - Peter Luschny, May 10 2016
From Peter Bala, Oct 15 2019: (Start)
a(n) = trace(M^n), where M is the 2 X 2 matrix [0, 1; -1, 4].
Consequently the Gauss congruences hold: a(n*p^k) = a(n*p^(k-1)) ( mod p^k ) for all prime p and positive integers n and k. See Zarelua and also Stanley (Ch. 5, Ex. 5.2(a) and its solution).
2*Sum_{n >= 1} 1/( a(n) - 6/a(n) ) = 1.
6*Sum_{n >= 1} (-1)^(n+1)/( a(n) + 2/a(n) ) = 1.
8*Sum_{n >= 1} 1/( a(n) + 24/(a(n) - 12/(a(n))) ) = 1.
8*Sum_{n >= 1} (-1)^(n+1)/( a(n) + 8/(a(n) + 4/(a(n))) ) = 1.
Series acceleration formulas for sums of reciprocals:
Sum_{n >= 1} 1/a(n) = 1/2 - 6*Sum_{n >= 1} 1/(a(n)*(a(n)^2 - 6)),
Sum_{n >= 1} 1/a(n) = 1/8 + 24*Sum_{n >= 1} 1/(a(n)*(a(n)^2 + 12)),
Sum_{n >= 1} (-1)^(n+1)/a(n) = 1/6 + 2*Sum_{n >= 1} (-1)^(n+1)/(a(n)*(a(n)^2 + 2)) and
Sum_{n >= 1} (-1)^(n+1)/a(n) = 1/8 + 8*Sum_{n >= 1} (-1)^(n+1)/(a(n)*(a(n)^2 + 12)).
Sum_{n >= 1} 1/a(n) = ( theta_3(2-sqrt(3))^2 - 1 )/4 = 0.34770 07561 66992 06261 .... See Borwein and Borwein, Proposition 3.5 (i), p.91.
Sum_{n >= 1} (-1)^(n+1)/a(n) = ( 1 - theta_3(sqrt(3)-2)^2 )/4. Cf. A003499 and A153415. (End)
a(n) = tan(Pi/12)^n + tan(5*Pi/12)^n. - Greg Dresden, Oct 01 2020
From Wolfdieter Lang, Sep 06 2021: (Start)
a(n) = S(n, 4) - S(n-2, 4) = 2*T(n, 2), for n >= 0, with S and T Chebyshev polynomials, with S(-1, x) = 0 and S(-2, x) = -1. S(n, 4) = A001353(n+1), for n >= -1, and T(n, 2) = A001075(n).
a(2*k) = A067902(k), a(2*k+1) = 4*A001570(k+1), for k >= 0. (End)
a(n) = sqrt(2 + 2*A011943(n+1)) = sqrt(2 + 2*A102344(n+1)), n>0. - Ralf Steiner, Sep 23 2021
Sum_{n>=1} arctan(3/a(n)^2) = Pi/6 - arctan(1/3) = A019673 - A105531 (Ohtskua, 2024). - Amiram Eldar, Aug 29 2024

Extensions

More terms from James Sellers, May 03 2000
Additional comments from Lekraj Beedassy, Feb 14 2002

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

A122367 Dimension of 3-variable non-commutative harmonics (twisted derivative) of order n. The dimension of the space of non-commutative polynomials of degree n in 3 variables which are killed by all symmetric differential operators (where for a monomial w, d_{xi} ( xi w ) = w and d_{xi} ( xj w ) = 0 for i != j).

Original entry on oeis.org

1, 2, 5, 13, 34, 89, 233, 610, 1597, 4181, 10946, 28657, 75025, 196418, 514229, 1346269, 3524578, 9227465, 24157817, 63245986, 165580141, 433494437, 1134903170, 2971215073, 7778742049, 20365011074, 53316291173, 139583862445, 365435296162, 956722026041
Offset: 0

Views

Author

Mike Zabrocki, Aug 30 2006

Keywords

Comments

Essentially identical to A001519.
From Matthew Lehman, Jun 14 2008: (Start)
Number of monotonic rhythms using n time intervals of equal duration (starting with n=0).
Representationally, let 0 be an interval which is "off" (rest),
1 an interval which is "on" (beep),
1 1 two consecutive "on" intervals (beep, beep),
1 0 1 (beep, rest, beep) and
1-1 two connected consecutive "on" intervals (beeeep).
For f(3)=13:
0 0 0, 0 0 1, 0 1 0, 0 1 1, 0 1-1, 1 0 0, 1 0 1,
1 1 0, 1-1 0, 1 1 1, 1 1-1, 1-1 1, 1-1-1.
(End)
Equivalent to the number of one-dimensional graphs of n nodes, subject to the condition that a node is either 'on' or 'off' and that any two neighboring 'on' nodes can be connected. - Matthew Lehman, Nov 22 2008
Sum_{n>=0} arctan(1/a(n)) = Pi/2. - Jaume Oliver Lafont, Feb 27 2009
This is the limit sequence for certain generalized Pell numbers. - Gregory L. Simay, Oct 21 2024

Examples

			a(1) = 2 because x1-x2, x1-x3 are both of degree 1 and are killed by the differential operator d_x1 + d_x2 + d_x3.
a(2) = 5 because x1*x2 - x3*x2, x1*x3 - x2*x3, x2*x1 - x3*x1, x1*x1 - x2*x1 - x2*x2 + x1*x2, x1*x1 - x3*x1 - x3*x3 + x3*x1 are all linearly independent and are killed by d_x1 + d_x2 + d_x3, d_x1 d_x1 + d_x2 d_x2 + d_x3 d_x3 and Sum_{j = 1..3} (d_xi d_xj, i).
		

Crossrefs

Programs

  • Magma
    [Fibonacci(2*n+1): n in [0..40]]; // Vincenzo Librandi, Jul 04 2015
    
  • Maple
    a:=n->if n=0 then 1; elif n=1 then 2 else 3*a(n-1)-a(n-2); fi;
    A122367List := proc(m) local A, P, n; A := [1,2]; P := [2];
    for n from 1 to m - 2 do P := ListTools:-PartialSums([op(A), P[-1]]);
    A := [op(A), P[-1]] od; A end: A122367List(30); # Peter Luschny, Mar 24 2022
  • Mathematica
    Table[Fibonacci[2 n + 1], {n, 0, 30}] (* Vincenzo Librandi, Jul 04 2015 *)
  • PARI
    Vec((1-x)/(1-3*x+x^2) + O(x^50)) \\ Michel Marcus, Jul 04 2015

Formula

G.f.: (1-q)/(1-3*q+q^2). More generally, (Sum_{d=0..n} (n!/(n-d)!*q^d)/Product_{r=1..d} (1 - r*q)) / (Sum_{d=0..n} q^d/Product_{r=1..d} (1 - r*q)) where n=3.
a(n) = 3*a(n-1) - a(n-2) with a(0) = 1, a(1) = 2.
a(n) = Fibonacci(2n+1) = A000045(2n+1). - Philippe Deléham, Feb 11 2009
a(n) = (2^(-1-n)*((3-sqrt(5))^n*(-1+sqrt(5)) + (1+sqrt(5))*(3+sqrt(5))^n)) / sqrt(5). - Colin Barker, Oct 14 2015
a(n) = Sum_{k=0..n} Sum_{i=0..n} binomial(k+i-1, k-i). - Wesley Ivan Hurt, Sep 21 2017
a(n) = A048575(n-1) for n >= 1. - Georg Fischer, Nov 02 2018
a(n) = Fibonacci(n)^2 + Fibonacci(n+1)^2. - Michel Marcus, Mar 18 2019
a(n) = Product_{k=1..n} (1 + 4*cos(2*k*Pi/(2*n+1))^2). - Seiichi Manyama, Apr 30 2021
From J. M. Bergot, May 27 2022: (Start)
a(n) = F(n)*L(n+1) + (-1)^n where L(n)=A000032(n) and F(n)=A000045(n).
a(n) = (L(n)^2 + L(n)*L(n+2))/5 - (-1)^n.
a(n) = 2*(area of a triangle with vertices at (L(n-1), L(n)), (F(n+1), F(n)), (L(n+1), L(n+2))) - 5*(-1)^n for n > 1. (End)
G.f.: (1-x)/(1-3x+x^2) = 1/(1-2x-x^2-x^3-x^4-...) - Gregory L. Simay, Oct 21 2024
E.g.f.: exp(3*x/2)*(5*cosh(sqrt(5)*x/2) + sqrt(5)*sinh(sqrt(5)*x/2))/5. - Stefano Spezia, Nov 07 2024
From Peter Bala, May 04 2025: (Start)
a(n) = sqrt(2/5) * sqrt( 1 - T(2*n+1, - 3/2) ), where T(k, x) denotes the k-th Chebyshev polynomial of the first kind.
a(2*n+1/2) = sqrt(5)*a(n)^2 - 2/sqrt(5).
a(3*n+1) = 5*a(n)^3 - 3*a(n); hence a(n) divides a(3*n+1).
a(4*n+3/2) = 5^(3/2)*a(n)^4 - 4*sqrt(5)*a(n)^2 + 2/sqrt(5).
a(5*n+2) = (5^2)*a(n)^5 - 5*5*a(n)^3 + 5*a(n); hence a(n) divides a(5*n+2).
See A034807 for the unsigned coefficients [1, 2; 1, 3; 1, 4, 2; 1, 5, 5; ...].
In general, for k >= 0, a(k*n + (k-1)/2) = a(-1/2) * T(k, a(n)/a(-1/2)), where a(n) = (2^(-1-n)*((3 - sqrt(5))^n *(-1 + sqrt(5)) + (1 + sqrt(5))*(3 + sqrt(5))^n)) / sqrt(5), as given above, and a(-1/2) = 2/sqrt(5).
The aerated sequence [b(n)]n>=1 = [1, 0, 2, 0, 5, 0, 13, 0, ...] is a fourth-order linear divisibility sequence; that is, if n | m then b(n) | b(m). It is the case P1 = 0, P2 = -5, Q = 1 of the 3-parameter family of divisibility sequences found by Williams and Guy.
Sum_{n >= 1} 1/(a(n) - 1/a(n)) = 1 (telescoping series: for n >= 1, 1/(a(n) - 1/a(n)) = 1/A001906(n) - 1/A001906(n+1).) (End)

A114525 Triangle of coefficients of the Lucas (w-)polynomials.

Original entry on oeis.org

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

Views

Author

Eric W. Weisstein, Dec 06 2005

Keywords

Comments

Unsigned version of A108045.
The row reversed triangle is A162514. - Paolo Bonzini, Jun 23 2016

Examples

			2, x, 2 + x^2, 3*x + x^3, 2 + 4*x^2 + x^4, 5*x + 5*x^3 + x^5, ... give triangle
  n\k   0  1  2  3  4  5  6  7  8  9 10 ...
  0:    2
  1:    0  1
  2:    2  0  1
  3:    0  3  0  1
  4:    2  0  4  0  1
  5:    0  5  0  5  0  1
  6:    2  0  9  0  6  0  1
  7:    0  7  0 14  0  7  0  1
  8:    2  0 16  0 20  0  8  0  1
  9:    0  9  0 30  0 27  0  9  0  1
  10:   2  0 25  0 50  0 35  0 10  0  1
  n\k   0  1  2  3  4  5  6  7  8  9 10 ...
  .... reformatted by _Wolfdieter Lang_, Feb 10 2023
		

Crossrefs

Cf. A108045 (signed version).
Cf. Sequences L(n,x): A000032(x = 1), A002203 (x = 2), A006497 (x = 3), A014448 (x = 4), A087130 (x = 5), A085447 (x = 6), A086902 (x = 7), A086594 (x = 8), A087798 (x = 9), A086927 (x = 10), A001946 (x = 11), A086928 (x = 12), A088316 (x = 13), A090300 (x = 14), A090301 (x = 15), A090305 (x = 16), A090306 (x = 17), A090307 (x = 18), A090308 (x = 19), A090309 (x = 20), A090310 (x = 21), A090313 (x = 22), A090314 (x = 23), A090316 (x = 24), A087281 (x = 29), A087287 (x = 76), A089772 (x = 199).

Programs

  • Maple
    Lucas := proc(n,x)
        option remember;
        if  n=0 then
            2;
        elif n =1 then
            x ;
        else
            x*procname(n-1,x)+procname(n-2,x) ;
        end if;
    end proc:
    A114525 := proc(n,k)
        coeftayl(Lucas(n,x),x=0,k) ;
    end proc:
    seq(seq(A114525(n,k),k=0..n),n=0..12) ; # R. J. Mathar, Aug 16 2019
  • Mathematica
    row[n_] := CoefficientList[LucasL[n, x], x];
    Table[row[n], {n, 0, 12}] // Flatten (* Jean-François Alcover, Aug 11 2018 *)

Formula

From Peter Bala, Mar 18 2015: (Start)
The Lucas polynomials L(n,x) satisfy the recurrence L(n+1,x) = x*L(n,x) + L(n-1,x) with L(0,x) = 2 and L(1,x) = x.
O.g.f.: Sum_{n >= 0} L(n,x)*t^n = (2 - x*t)/(1 - t^2 - x*t) = 2 + x*t + (x^2 + 2)*t^2 + (3*x + x^3)*t^3 + ....
L(n,x) = trace( [ x, 1; 1, 0 ]^n ).
exp( Sum_{n >= 1} L(n,x)*t^n/n ) = Sum_{n >= 0} F(n+1,x)*t^n, where F(n,x) denotes the n-th Fibonacci polynomial. (see Appendix A3 in Johnson).
exp( Sum_{n >= 1} L(n,x)*L(2*n,x)*t^n/n ) = 1/( F(1,x)*F(2*x)*F(3,x) ) * Sum_{n >= 0} F(n+1,x)*F(n+2,x)*F(n+3,x)*t^n.
exp( Sum_{n >= 1} L(3*n,x)/L(n,x)*t^n/n ) = Sum_{n >= 0} L(2*n + 1,x)*t^n.
L(n,1) = Lucas(n) = A000032(n); L(n,4) = Lucas(3*n) = A014448(n); L(n,11) = Lucas(5*n) = A001946(n); L(n,29) = Lucas(7*n) = A087281(n); L(n,76) = Lucas(9*n) = A087287(n); L(n,199) = Lucas(11*n) = A089772(n). The general result is L(n,Lucas(2*k + 1)) = Lucas((2*k + 1)*n). (End)
From Jeremy Dover, Jun 10 2016: (Start)
Read as a triangle T(n,k), n >= 0, n >= k >= 0, T(n,k) = (Binomial((n+k)/2,k) + Binomial((n+k-2)/2,k))*(1+(-1)^(n-k))/2.
T(n,k) = A046854(n-1,k-1) + A046854(n-1,k) + A046854(n-2,k) for even n+k with n+k > 0, assuming A046854(n,k) = 0 for n < 0, k < 0, k > n.
T(n,k) is the number of binary strings of length n with exactly k pairs of consecutive 0's and no pair of consecutive 1's, where the first and last bits are considered consecutive. (End)
From Peter Bala, Sep 03 2022: (Start)
L(n,x) = 2*(i)^n*T(n,-i*x/2), where i = sqrt(-1) and T(n,x) is the n-th Chebyshev polynomial of the first kind.
d/dx(L(n,x)) = n*F(n,x), where F(n,x) denotes the n-th Fibonacci polynomial.
Let P_n(x,y) = (L(n,x) - L(n,y))/(x - y). Then {P_n(x,y): n >= 1} is a fourth-order linear divisibility sequence of polynomials in the ring Z[x,y]: if m divides n in Z then P_m(x,y) divides P_n(x,y) in Z[x,y].
P_n(1,1) = A045925(n); P_n(1,4) = A273622; P_n(2,2) = A093967(n).
L(2*n,x)^2 - L(2*n-1,x)*L(2*n+1,x) = x^2 + 4 for n >= 1.
Sum_{n >= 1} L(2*n,x)/( L(2*n-1,x) * L(2*n+1,x) ) = 1/x^2 and
Sum_{n >= 1} (-1)^(n+1)/( L(2*n,x) + x^2/L(2*n,x) ) = 1/(x^2 + 4), both valid for all nonzero real x. (End)
From Peter Bala, Nov 18 2022: (Start)
L(n,x) = Sum_{k = 0..floor(n/2)} (n/(n-k))*binomial(n-k,k)*x^(n-2*k) for n >= 1.
For odd m, L(n, L(m,x)) = L(n*m, x).
For integral x, the sequence {u(n)} := {L(n,x)} satisfies the Gauss congruences: u(m*p^r) == u(m*p^(r-1)) (mod p^r) for all positive integers m and r and all primes p.
Let p be an odd prime and let 0 <= k <= p - 1. Let alpha_k = the p-adic limit_{n -> oo} L(p^n,k). Then alpha_k is a well-defined p-adic-integer and the polynomial L(p,x) - x of degree p factorizes as L(p,x) - x = Product_{k = 0..p-1} (x - alpha_k). For example, L(5,x) - x = x^5 + 5*x^3 + 4*x = x*(x - A269591)*(x - A210850)*(x - A210851)*(x - A269592) in the ring of 5-adic integers. (End)
The formula for L(n,x) given in the first line of the preceding section, with L(0, x) = 2, is rewritten L(n, x) = Sum_{k = 0..floor(n/2)} A034807(n, k)*x^(n - 2*k). See the formula by Alexander Elkins in A034807. - Wolfdieter Lang, Feb 10 2023

A263916 Coefficients of the Faber partition polynomials.

Original entry on oeis.org

-1, -2, 1, -3, 3, -1, -4, 4, 2, -4, 1, -5, 5, 5, -5, -5, 5, -1, -6, 6, 6, -6, 3, -12, 6, -2, 9, -6, 1, -7, 7, 7, -7, 7, -14, 7, -7, -7, 21, -7, 7, -14, 7, -1, -8, 8, 8, -8, 8, -16, 8, 4, -16, -8, 24, -8, -8, 12, 24, -32, 8, 2, -16, 20, -8, 1
Offset: 1

Views

Author

Tom Copeland, Oct 29 2015

Keywords

Comments

The coefficients of the Faber polynomials F(n,b(1),b(2),...,b(n)) (Bouali, p. 52) in the order of the partitions of Abramowitz and Stegun. Compare with A115131 and A210258.
These polynomials occur in discussions of the Virasoro algebra, univalent function spaces and the Schwarzian derivative, symmetric functions, and free probability theory. They are intimately related to symmetric functions, free probability, and Appell sequences through the raising operator R = x - d log(H(D))/dD for the Appell sequence inverse pair associated to the e.g.f.s H(t)e^(xt) (cf. A094587) and (1/H(t))e^(xt) with H(0)=1.
Instances of the Faber polynomials occur in discussions of modular invariants and modular functions in the papers by Asai, Kaneko, and Ninomiya, by Ono and Rolen, and by Zagier. - Tom Copeland, Aug 13 2019
The Faber polynomials, denoted by s_n(a(t)) where a(t) is a formal power series defined by a product formula, are implicitly defined by equation 13.4 on p. 62 of Hazewinkel so as to extract the power sums of the reciprocals of the zeros of a(t). This is the Newton identity expressing the power sum symmetric polynomials in terms of the elementary symmetric polynomials/functions. - Tom Copeland, Jun 06 2020
From Tom Copeland, Oct 15 2020: (Start)
With a_n = n! * b_n = (n-1)! * c_n for n > 0, represent a function with f(0) = a_0 = b_0 = 1 as an
A) exponential generating function (e.g.f), or formal Taylor series: f(x) = e^{a.x} = 1 + Sum_{n > 0} a_n * x^n/n!
B) ordinary generating function (o.g.f.), or formal power series: f(x) = 1/(1-b.x) = 1 + Sum_{n > 0} b_n * x^n
C) logarithmic generating function (l.g.f): f(x) = 1 - log(1 - c.x) = 1 + Sum_{n > 0} c_n * x^n /n.
Expansions of log(f(x)) are given in
I) A127671 and A263634 for the e.g.f: log[ e^{a.*x} ] = e^{L.(a_1,a_2,...)x} = Sum_{n > 0} L_n(a_1,...,a_n) * x^n/n!, the logarithmic polynomials, cumulant expansion polynomials
II) A263916 for the o.g.f.: log[ 1/(1-b.x) ] = log[ 1 - F.(b_1,b_2,...)x ] = -Sum_{n > 0} F_n(b_1,...,b_n) * x^n/n, the Faber polynomials.
Expansions of exp(f(x)-1) are given in
III) A036040 for an e.g.f: exp[ e^{a.x} - 1 ] = e^{BELL.(a_1,...)x}, the Bell/Touchard/exponential partition polynomials, a.k.a. the Stirling partition polynomials of the second kind
IV) A130561 for an o.g.f.: exp[ b.x/(1-b.x) ] = e^{LAH.(b.,...)x}, the Lah partition polynomials
V) A036039 for an l.g.f.: exp[ -log(1-c.x) ] = e^{CIP.(c_1,...)x}, the cycle index polynomials of the symmetric groups S_n, a.k.a. the Stirling partition polynomials of the first kind.
Since exp and log are a compositional inverse pair, one can extract the indeterminates of the log set of partition polynomials from the exp set and vice versa. For a discussion of the relations among these polynomials and the combinatorics of connected and disconnected graphs/maps, see Novak and LaCroix on classical moments and cumulants and the two books on statistical mechanics referenced in A036040. (End)

Examples

			F(1,b1) = - b1
F(2,b1,b2) = -2 b2 + b1^2
F(3,b1,b2,b3) = -3 b3 + 3 b1 b2 - b1^3
F(4,b1,...) = -4 b4 + 4 b1 b3 + 2 b2^2  - 4 b1^2 b2 + b1^4
F(5,...) = -5 b5 + 5 b1 b4 + 5 b2 b3 - 5 b1^2 b3 - 5 b1 b2^2 + 5 b1^3 b2 - b1^5
------------------------------
IF(1,b1) = -b1
IF(2,b1,,b2) = -b2 + b1^2
IF(3,b1,b2,b3) = -2 b3 + 3 b1 b2 - b1^3
IF(4,b1,...) = -6 b4 + 8 b1 b3 + 3 b2^2  - 6 b1^2 b2 + b1^4
IF(5,...) = -24 b5 + 30 b1 b4 + 20 b2 b3 - 20 b1^2 b3 - 15 b1 b2^2 + 10 b1^3 b2 - b1^5
------------------------------
For 1/(1+x)^2 = 1- 2x + 3x^2 - 4x^3 + 5x^4 - ..., F(n,-2,3,-4,...) = (-1)^(n+1) 2.
------------------------------
F(n,x,2x,...,nx), F(n,-x,2x,-3x,...,(-1)^n n*x), and F(n,(2-x),1,0,0,...) are related to the Chebyshev polynomials through A127677 and A111125. See also A110162, A156308, A208513, A217476, and A220668.
------------------------------
For b1 = p, b2 = q, and all other indeterminates 0, see A113279 and A034807.
For b1 = -y, b2 = 1 and all other indeterminates 0, see A127672.
		

References

  • H. Airault, "Symmetric sums associated to the factorization of Grunsky coefficients," in Groups and Symmetries: From Neolithic Scots to John McKay, CRM Proceedings and Lecture Notes: Vol. 47, edited by J. Harnad and P. Winternitz, American Mathematical Society, 2009.
  • D. Bleeker and B. Booss, Index Theory with Applications to Mathematics and Physics, International Press, 2013, (see section 16.7 Characteristic Classes and Curvature).
  • M. Hazewinkel, Formal Groups and Applications, Academic Press, New York San Francisco London, 1978, p. 120.
  • F. Hirzebruch, Topological methods in algebraic geometry. Second, corrected printing of the third edition. Die Grundlehren der Mathematischen Wissenschaften, Band 131 Springer-Verlag, Berlin Heidelberg New York, 1978, p. 11 and 92.
  • D. Knutson, λ-Rings and the Representation Theory of the Symmetric Group, Lect. Notes in Math. 308, Springer-Verlag, 1973, p. 35.
  • D. Yau, Lambda-Rings, World Scientific Publishing Co., Singapore, 2010, p. 45.

Crossrefs

Programs

  • Mathematica
    F[0] = 1; F[1] = -b[1]; F[2] = b[1]^2 - 2 b[2]; F[n_] := F[n] = -b[1] F[n - 1] - Sum[b[n - k] F[k], {k, 1, n - 2}] - n b[n] // Expand;
    row[n_] := (List @@ F[n]) /. b[_] -> 1 // Reverse;
    Table[row[n], {n, 1, 8}] // Flatten // Rest (* Jean-François Alcover, Jun 12 2017 *)

Formula

-log(1 + b(1) x + b(2) x^2 + ...) = Sum_{n>=1} F(n,b(1),...,b(n)) * x^n/n.
-d(1 + b(1) x + b(2) x^2 + ...)/dx / (1 + b(1) x + b(2) x^2 + ...) = Sum_{n>=1} F(n,b(1),...,b(n)) x^(n-1).
F(n,b(1),...,b(n)) = -n*b(n) - Sum_{k=1..n-1} b(n-k)*F(k,b(1),...,b(k)).
Umbrally, with B(x) = 1 + b(1) x + b(2) x^2 + ..., B(x) = exp[log(1-F.x)] and 1/B(x) = exp[-log(1-F.x)], establishing a connection to the e.g.f. of A036039 and the symmetric polynomials.
The Stirling partition polynomials of the first kind St1(n,b1,b2,...,bn;-1) = IF(n,b1,b2,...,bn) (cf. the Copeland link Lagrange a la Lah, signed A036039, and p. 184 of Airault and Bouali), i.e., the cyclic partition polynomials for the symmetric groups, and the Faber polynomials form an inverse pair for isolating the indeterminates in their definition, for example, F(3,IF(1,b1),IF(2,b1,b2)/2!,IF(3,b1,b2,b3)/3!)= b3, with bk = b(k), and IF(3,F(1,b1),F(2,b1,b2),F(3,b1,b2,b3))/3!= b3.
The polynomials specialize to F(n,t,t,...) = (1-t)^n - 1.
See Newton Identities on Wikipedia on relation between the power sum symmetric polynomials and the complete homogeneous and elementary symmetric polynomials for an expression in multinomials for the coefficients of the Faber polynomials.
(n-1)! F(n,x[1],x[2]/2!,...,x[n]/n!) = - p_n(x[1],...,x[n]), where p_n are the cumulants of A127671 expressed in terms of the moments x[n]. - Tom Copeland, Nov 17 2015
-(n-1)! F(n,B(1,x[1]),B(2,x[1],x[2])/2!,...,B(n,x[1],...,x[n])/n!) = x[n] provides an extraction of the indeterminates of the complete Bell partition polynomials B(n,x[1],...,x[n]) of A036040. Conversely, IF(n,-x[1],-x[2],-x[3]/2!,...,-x[n]/(n-1)!) = B(n,x[1],...,x[n]). - Tom Copeland, Nov 29 2015
For a square matrix M, determinant(I - x M) = exp[-Sum_{k>0} (trace(M^k) x^k / k)] = Sum_{n>0} [ P_n(-trace(M),-trace(M^2),...,-trace(M^n)) x^n/n! ] = 1 + Sum_{n>0} (d[n] x^n), where P_n(x[1],...,x[n]) are the cycle index partition polynomials of A036039 and d[n] = P_n(-trace(M),-trace(M^2),...,-trace(M^n)) / n!. Umbrally, det(I - x M)= exp[log(1 - b. x)] = exp[P.(-b_1,..,-b_n)x] = 1 / (1-d.x), where b_k = tr(M^k). Then F(n,d[1],...,d[n]) = tr[M^n]. - Tom Copeland, Dec 04 2015
Given f(x) = -log(g(x)) = -log(1 + b(1) x + b(2) x^2 + ...) = Sum_{n>=1} F(n,b(1),...,b(n)) * x^n/n, action on u_n = F(n,b(1),...,b(n)) with A133932 gives the compositional inverse finv(x) of f(x), with F(1,b(1)) not equal to zero, and f(g(finv(x))) = f(e^(-x)). Note also that exp(f(x)) = 1 / g(x) = exp[Sum_{n>=1} F(n,b(1),...,b(n)) * x^n/n] implies relations among A036040, A133314, A036039, and the Faber polynomials. - Tom Copeland, Dec 16 2015
The Dress and Siebeneicher paper gives combinatorial interpretations and various relations that the Faber polynomials must satisfy for integral values of its arguments. E.g., Eqn. (1.2) p. 2 implies [2 * F(1,-1) + F(2,-1,b2) + F(4,-1,b2,b3,b4)] mod(4) = 0. This equation implies that [F(n,b1,b2,...,bn)-(-b1)^n] mod(n) = 0 for n prime. - Tom Copeland, Feb 01 2016
With the elementary Schur polynomials S(n,a_1,a_2,...,a_n) = Lah(n,a_1,a_2,...,a_n) / n!, where Lah(n,...) are the refined Lah polynomials of A130561, F(n,S(1,a_1),S(2,a_1,a_2),...,S(n,a_1,...,a_n)) = -n * a_n since sum_{n > 0} a_n x^n = log[sum{n >= 0} S(n,a_1,...,a_n) x^n]. Conversely, S(n,-F(1,a_1),-F(2,a_1,a_2)/2,...,-F(n,a_1,...,a_n)/n) = a_n. - Tom Copeland, Sep 07 2016
See Corollary 3.1.3 on p. 38 of Ardila and Copeland's two MathOverflow links to relate the Faber polynomials, with arguments being the signed elementary symmetric polynomials, to the logarithm of determinants, traces of powers of an adjacency matrix, and number of walks on graphs. - Tom Copeland, Jan 02 2017
The umbral inverse polynomials IF appear on p. 19 of Konopelchenko as partial differential operators. - Tom Copeland, Nov 19 2018

Extensions

More terms from Jean-François Alcover, Jun 12 2017

A028297 Coefficients of Chebyshev polynomials of the first kind: triangle of coefficients in expansion of cos(n*x) in descending powers of cos(x).

Original entry on oeis.org

1, 1, 2, -1, 4, -3, 8, -8, 1, 16, -20, 5, 32, -48, 18, -1, 64, -112, 56, -7, 128, -256, 160, -32, 1, 256, -576, 432, -120, 9, 512, -1280, 1120, -400, 50, -1, 1024, -2816, 2816, -1232, 220, -11, 2048, -6144, 6912, -3584, 840, -72, 1, 4096, -13312, 16640, -9984
Offset: 0

Views

Author

Keywords

Comments

Rows are of lengths 1, 1, 2, 2, 3, 3, ... (A008619).
This triangle is generated from A118800 by shifting down columns to allow for (1, 1, 2, 2, 3, 3, ...) terms in each row. - Gary W. Adamson, Dec 16 2007
Unsigned triangle = A034839 * A007318. - Gary W. Adamson, Nov 28 2008
Triangle, with zeros omitted, given by (1, 1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, -1, 1, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 16 2011
From Wolfdieter Lang, Aug 02 2014: (Start)
This irregular triangle is the row reversed version of the Chebyshev T-triangle A053120 given by A039991 with vanishing odd-indexed columns removed.
If zeros are appended in each row n >= 1, in order to obtain a regular triangle (see the Philippe Deléham comment, g.f. and example) this becomes the Riordan triangle (1-x)/(1-2*x), -x^2/(1-2*x). See also the unsigned version A201701 of this regular triangle.
(End)
Apparently, unsigned diagonals of this array are rows of A200139. - Tom Copeland, Oct 11 2014
It appears that the coefficients are generated by the following: Let SM_k = Sum( d_(t_1, t_2)* eM_1^t_1 * eM_2^t_2) summed over all length 2 integer partitions of k, i.e., 1*t_1 + 2*t_2 = k, where SM_k are the averaged k-th power sum symmetric polynomials in 2 data (i.e., SM_k = S_k/2 where S_k are the k-th power sum symmetric polynomials, and where eM_k are the averaged k-th elementary symmetric polynomials, eM_k = e_k/binomial(2,k) with e_k being the k-th elementary symmetric polynomials. The data d_(t_1, t_2) form an irregular triangle, with one row for each k value starting with k=1. Thus this procedure and associated OEIS sequences A287768, A288199, A288207, A288211, A288245, A288188 are generalizations of Chebyshev polynomials of the first kind. - Gregory Gerard Wojnar, Jul 01 2017

Examples

			Letting c = cos x, we have: cos 0x = 1, cos 1x = 1c; cos 2x = 2c^2-1; cos 3x = 4c^3-3c, cos 4x = 8c^4-8c^2+1, etc.
T4 = 8x^4 - 8x^2 + 1 = 8, -8, +1 = 2^(3) - (4)(2) + [2^(-1)](4)/2.
From _Wolfdieter Lang_, Aug 02 2014: (Start)
The irregular triangle T(n,k) begins:
n\k     1      2     3      4     5     6   7   8 ....
0:      1
1:      1
2:      2     -1
3:      4     -3
4:      8     -8     1
5:     16    -20     5
6:     32    -48    18     -1
7:     64   -112    56     -7
8:    128   -256   160    -32     1
9:    256   -576   432   -120     9
10:   512  -1280  1120   -400    50    -1
11:  1024  -2816  2816  -1232   220   -11
12:  2048  -6144  6912  -3584   840   -72   1
13:  4096 -13312 16640  -9984  2912  -364  13
14:  8192 -28672 39424 -26880  9408 -1568  98  -1
15: 16384 -61440 92160 -70400 28800 -6048 560 -15
...
T(4,x) = 8*x^4 -8*x^2 + 1*x^0, T(5,x) = 16*x^5 - 20*x^3 + 5*x^1, with Chebyshev's T-polynomials (A053120). (End)
From _Philippe Deléham_, Dec 16 2011: (Start)
The triangle (1,1,0,0,0,0,...) DELTA (0,-1,1,0,0,0,0,...) includes zeros and begins:
   1;
   1,   0;
   2,  -1,  0;
   4,  -3,  0,  0;
   8,  -8,  1,  0, 0;
  16, -20,  5,  0, 0, 0;
  32, -48, 18, -1, 0, 0, 0; (End)
		

References

  • I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series and Products, 5th ed., Section 1.335, p. 35.
  • S. Selby, editor, CRC Basic Mathematical Tables, CRC Press, 1970, p. 106. [From Rick L. Shepherd, Jul 06 2010]

Crossrefs

Cf. A028298.
Reflection of A008310, the main entry. With zeros: A039991.
Cf. A053120 (row reversed table including zeros).
Cf. A001333 (row sums 1), A001333 (alternating row sums). - Wolfdieter Lang, Aug 02 2014

Programs

  • Maple
    b:= proc(n) b(n):= `if`(n<2, 1, expand(2*b(n-1)-x*b(n-2))) end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(n)):
    seq(T(n), n=0..15);  # Alois P. Heinz, Sep 04 2019
  • Mathematica
    t[n_] := (Cos[n x] // TrigExpand) /. Sin[x]^m_ /; EvenQ[m] -> (1 - Cos[x]^2)^(m/2) // Expand; Flatten[Table[ r = Reverse @ CoefficientList[t[n], Cos[x]]; If[OddQ[Length[r]], AppendTo[r,0]]; Partition[r,2][[All, 1]],{n, 0, 13}] ][[1 ;; 53]] (* Jean-François Alcover, May 06 2011 *)
    Tpoly[n_] := HypergeometricPFQ[{(1 - n)/2, -n/2}, {1/2}, 1 - x];
    Table[CoefficientList[Tpoly[n], x], {n, 0, 12}] // Flatten (* Peter Luschny, Feb 03 2021 *)

Formula

cos(n*x) = 2 * cos((n-1)*x) * cos(x) - cos((n-2)*x) (from CRC's Multiple-angle relations). - Rick L. Shepherd, Jul 06 2010
G.f.: (1-x) / (1-2x+y*x^2). - Philippe Deléham, Dec 16 2011
Sum_{k=0..n} T(n,k)*x^k = A011782(n), A000012(n), A146559(n), A087455(n), A138230(n), A006495(n), A138229(n) for x = 0, 1, 2, 3, 4, 5, 6, respectively. - Philippe Deléham, Dec 16 2011
T(n,k) = [x^k] hypergeom([1/2 - n/2, -n/2], [1/2], 1 - x). - Peter Luschny, Feb 03 2021
T(n,k) = (-1)^k * 2^(n-1-2*k) * A034807(n,k). - Hoang Xuan Thanh, Jun 21 2025

Extensions

More terms from David W. Wilson
Row length sequence and link to Abramowitz-Stegun added by Wolfdieter Lang, Aug 02 2014

A102426 Triangle read by rows giving coefficients of polynomials defined by F(0,x)=0, F(1,x)=1, F(n,x) = F(n-1,x) + x*F(n-2,x).

Original entry on oeis.org

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

Views

Author

Russell Walsmith, Jan 08 2005

Keywords

Comments

Essentially the same as A098925: a(0)=0 followed by A098925. - R. J. Mathar, Aug 30 2008
F(n) + 2x * F(n-1) gives Lucas polynomials (cf. A034807). - Maxim Krikun (krikun(AT)iecn.u-nancy.fr), Jun 24 2007
After the initial 0, these are the nonzero coefficients of the Fibonacci polynomials; see the Mathematica section. - Clark Kimberling, Oct 10 2013
Aside from signs and index shift, the coefficients of the characteristic polynomial of the Coxeter adjacency matrix for the Coxeter group A_n related to the Chebyshev polynomial of the second kind (cf. Damianou link p. 19). - Tom Copeland, Oct 11 2014
Aside from the initial zeros, these are the antidiagonals read from bottom to top of the numerical coefficients of the Maurer-Cartan form matrix of the Leibniz group L^(n)(1,1) presented on p. 9 of the Olver paper, which is generated as exp[c. * M] with (c.)^n = c_n and M the Lie infinitesimal generator A218272. Reverse of A011973. - Tom Copeland, Jul 02 2018

Examples

			The first few polynomials are:
  0
  1
  1
  x + 1
  2*x + 1
  x^2 + 3*x + 1
  3*x^2 + 4*x + 1
------------------
From _Tom Copeland_, Jan 19 2016: (Start)
[n]:
0:  0
1:  1
2:  1
3:  1  1
4:  2  1
5:  1  3  1
6:  3  4  1
7:  1  6  5   1
8:  4 10  6   1
9:  1 10 15   7   1
10: 5 20 21   8   1
11: 1 15 35  28   9  1
12: 6 35 56  36  10  1
13: 1 21 70  84  45 11 1
(End)
		

References

  • Dominique Foata and Guo-Niu Han, Multivariable tangent and secant q-derivative polynomials, Manuscript, Mar 21 2012.

Crossrefs

Upward diagonals sums are A062200. Downward rows are A102427. Row sums are A000045. Row terms reversed = A011973. Also A102428, A102429.
All of A011973, A092865, A098925, A102426, A169803 describe essentially the same triangle in different ways.

Programs

  • Magma
    [0] cat [Binomial(Floor(n/2)+k, Floor((n-1)/2-k) ): k in [0..Floor((n-1)/2)], n in [0..17]]; // G. C. Greubel, Oct 13 2019
    
  • Mathematica
    Join[{0}, Table[ Select[ CoefficientList[ Fibonacci[n, x], x], 0 < # &], {n, 0, 17}]//Flatten] (* Clark Kimberling, Oct 10 2013 and slightly modified by Robert G. Wilson v, May 03 2017 *)
  • PARI
    F(n) = if (n==0, 0, if (n==1, 1, F(n-1) + x*F(n-2)));
    tabf(nn) = for (n=0, nn, print(Vec(F(n)))); \\ Michel Marcus, Feb 10 2020

Formula

Alternatively, as n is even or odd: T(n-2, k) + T(n-1, k-1) = T(n, k), T(n-2, k) + T(n-1, k) = T(n, k)
T(n, k) = binomial(floor(n/2)+k, floor((n-1)/2-k) ). - Paul Barry, Jun 22 2005
Beginning with the second polynomial in the example and offset=0, P(n,t)= Sum_{j=0..n}, binomial(n-j,j)*x^j with the convention that 1/k! is zero for k=-1,-2,..., i.e., 1/k! = lim_{c->0} 1/(k+c)!. - Tom Copeland, Oct 11 2014
From Tom Copeland, Jan 19 2016: (Start)
O.g.f.: (x + x^2 - x^3) / (1 - (2+t)*x^2 + x^4) = (x^2 (even part) + x*(1-x^2) (odd)) / (1 - (2+t)*x^2 + x^4).
Recursion relations:
A) p(n,t) = p(n-1,t) + p(n-2,t) for n=2,4,6,8,...
B) p(n,t) = t*p(n-1,t) + p(n-2,t) for n=3,5,7,...
C) a(n,k) = a(n-2,k) + a(n-1,k) for n=4,6,8,...
D) a(n,k) = a(n-2,k) + a(n-1,k-1) for n=3,5,7,...
Relation A generalized to MV(n,t;r) = P(2n+1,t) + r R(2n,t) for n=1,2,3,... (cf. A078812 and A085478) is the generating relation on p. 229 of Andre-Jeannine for the generalized Morgan-Voyce polynomials, e.g., MV(2,t;r) = p(5,t) + r*p(4,t) = (1 + 3t + t^2) + r*(2 + t) = (1 + 2r) + (3 + r)*t + t^2, so P(n,t) = MV(n-4,t;1) for n=4,6,8,... .
The even and odd polynomials are also presented in Trzaska and Ferri.
Dropping the initial 0 and re-indexing with initial m=0 gives the row polynomials Fb(m,t) = p(n+1,t) below with o.g.f. G(t,x)/x, starting with Fb(0,t) = 1, Fb(1,t) = 1, Fb(2,t) = 1 + t, and Fb(3,t) = 2 + t.
The o.g.f. x/G(x,t) = (1 - (2+t)*x^2 + x^4) / (1 + x - x^2) then generates a sequence of polynomials IFb(t) such that the convolution Sum_{k=0..n} IFb(n-k,t) Fb(k,t) vanishes for n>1 and is one for n=0. These linear polynomials have the basic Fibonacci numbers A000045 as an overall factor:
IFb(0,t) = 1
IFb(1,t) = -1
IFb(2,t) = -t
IFb(3,t) = -1 (1-t)
IFb(4,t) = 2 (1-t)
IFb(5,t) = -3 (1-t)
IFb(6,t) = 5 (1-t)
IFb(7,t) = -8 (1-t)
IFb(8,t) = 13 (1-t)
... .
(End)

Extensions

Name corrected by John K. Sikora, Feb 10 2020

A127677 Scaled coefficient table for Chebyshev polynomials 2*T(2*n, sqrt(x)/2) (increasing even scaled powers, without zero entries).

Original entry on oeis.org

2, -2, 1, 2, -4, 1, -2, 9, -6, 1, 2, -16, 20, -8, 1, -2, 25, -50, 35, -10, 1, 2, -36, 105, -112, 54, -12, 1, -2, 49, -196, 294, -210, 77, -14, 1, 2, -64, 336, -672, 660, -352, 104, -16, 1, -2, 81, -540, 1386, -1782, 1287, -546, 135, -18, 1, 2, -100, 825, -2640, 4290, -4004, 2275, -800, 170, -20, 1
Offset: 0

Views

Author

Wolfdieter Lang, Mar 07 2007

Keywords

Comments

2*T(2*n,x) = Sum_{m=0..n} a(n,m)*(2*x)^(2*m).
Closely related to A284982, which has opposite signs and rows begin with 0 of alternating signs instead of +/2. - Eric W. Weisstein, Apr 07 2017
Bisection triangle of A127672 (without zero entries, even part). The odd part is ((-1)^(n-m))*A111125(n,m).
If the leading 2 is replaced by a 1 we get the essentially identical sequence A110162. - N. J. A. Sloane, Jun 09 2007
Also row n gives coefficients of characteristic polynomial of the Cartan matrix for the root system B_n (or, equally, C_n). - Roger L. Bagula, May 23 2007
From Wolfdieter Lang, Oct 04 2013: (Start)
This triangle a(n,m) is used to express the length ratio side/R given by s(4*n+2) = 2*sin(Pi/(4*n+2)) = 2*cos(2*n*Pi/(4*n+2)) in a regular (4*n+2)-gon, inscribed in a circle with radius R, in terms of rho(4*n+2) = 2*cos(Pi/4*n+2), the length ratio of (the smallest diagonal)/side (for n=2 there is no such diagonal).
s(4*n+2) = Sum_{m=0..n}a(n,m)*rho(4*n+2)^(2*m). This formula is needed to show that the total sum of all length ratios in a (4*n+2)-gon is an integer in the algebraic number field Q(rho(4*n+2)). Note that rho(4*n+2) has degree delta(4*n+2) = A055034(4*n+2). Therefore one has to take s(4*n+2) modulo C(4*n+2, x=rho(4*n+2)), the minimal polynomial of rho(4*n+2) (see A187360). Thanks go to Seppo Mustonen for asking me to look into this problem. See ((-1)^(n-m))*A111125(n,m) for the (4*n)-gon situation. (End)

Examples

			The triangle a(n,m) starts:
n\m  0    1    2     3     4     5     6     7    8   9  10 ...
0:   2
1:  -2    1
2:   2   -4    1
3:  -2    9   -6     1
4:   2  -16   20    -8     1
5:  -2   25  -50    35   -10     1
6:   2  -36  105  -112    54   -12     1
7:  -2   49 -196   294  -210    77   -14     1
8:   2  -64  336  -672   660  -352   104   -16    1
9:  -2   81 -540  1386 -1782  1287  -546   135  -18   1
10:  2 -100  825 -2640  4290 -4004  2275  -800  170 -20  1
... Reformatted and extended by _Wolfdieter Lang_, Nov 21 2012.
n=3: [-2,9,-6,1] stands for -2*1 + 9*(2*x)^2 -6*(2*x)^4 +1*(2*x)^6 = 2*(1+18*x^2-48*x^4+32*x^6) = 2*T(6,x).
(4*n+2)-gon side/radius s(4*n+2) as polynomial in rho(4*n+2) = smallest diagonal/side: n=0: s(2) = 2 (rho(2)=0); n=1: s(6) = -2 + rho(6)^2 = -2 + 3 = 1, (C(6,x) = x^2 - 3); n=2: s(10) = 2 - 4*rho(10)^2 + 1*rho(10)^4 = 2 - 4*rho(10)^2 + (5*rho(10)^2 - 5) = -3 + rho(10)^2, (C(10,x) = x^4 - 5*x^2 + 5). - _Wolfdieter Lang_, Oct 04 2013
		

References

  • R. N. Cahn, Semi-Simple Lie Algebras and Their Representations, Dover, NY, 2006, ISBN 0-486-44999-8, p. 62
  • Sigurdur Helgasson,Differential Geometry, Lie Groups and Symmetric Spaces, Graduate Studies in Mathematics, volume 34. A. M. S. :ISBN 0-8218-2848-7, 1978,p. 463.

Crossrefs

Cf. A284982 (opposite signs and rows begin with 0).
Row sums (signed): -A061347(n+3) for n>=0.
Row sums (unsigned): A005248(n) = L(2*n), where L=Lucas.

Programs

  • Mathematica
    T[n_, m_, d_] := If[ n == m, 2, If[n == d && m == d - 1, -2, If[(n == m - 1 || n == m + 1), -1, 0]]] M[d_] := Table[T[n, m, d], {n, 1, d}, {m, 1, d}] a = Join[M[1], Table[CoefficientList[CharacteristicPolynomial[M[d], x], x], {d, 1, 10} ]] (* Roger L. Bagula, May 23 2007 *)
    CoefficientList[2 ChebyshevT[2 Range[0, 10], Sqrt[x]/2], x] // Flatten (* Eric W. Weisstein, Apr 06 2017 *)
    CoefficientList[Table[(-1)^n LucasL[2 n, Sqrt[-x]], {n, 0, 10}], x] // Flatten (* Eric W. Weisstein, Apr 06 2017 *)
  • PARI
    a(n,m) = {if(n>=2, -2*a(n-1,m)+a(n-1,m-1)-a(n-2,m), if(n==0, if(m!=0,0,2), if(m==0,-2, if(m==1,1,0))))};
    for(n=0,10,for(m=0,n,print1(a(n,m),", "))) \\ Hugo Pfoertner, Jul 19 2020

Formula

a(n,m) = 0 if n < m; a(n,0) = 2*(-1)^n; a(n,m) = ((-1)^(n+m))*n*binomial(n+m-1, 2*m-1)/m.
a(n,m) = 0 if n < m, a(0,0) = 2, a(n,m) = (-1)^(n-m)*(2*n/(n+m))*binomial(n+m, n-m), n >= 1. From Waring's formula applied to Chebyshev's T-polynomials. See also A110162. - Wolfdieter Lang, Nov 21 2012
The o.g.f. Sum_{n>=0} p(n,x)*z^n, n>=0, for the row polynomials p(n,x) := Sum_{m=0..n} a(n,m)*x^m is (2 + z*(2-x))/((z+1)^2 - z*x). Here p(n,x) = R(2*n,sqrt(x)) := 2*T(2*n,sqrt(x)/2) with Chebyshev's T-polynomials. For the R-polynomials see A127672. - Wolfdieter Lang, Nov 28 2012
From Tom Copeland, Nov 07 2015: (Start)
A logarithmic generator is 2*(1-log(1+x))-log(1-t*x/(1+x)^2) = 2 - log(1+(2-t)*x+x^2) = 2 + (-2 + t)*x + (2 - 4*t + t^2) x^2/2 + (-2 + 9*t - 6*t^2 + t^3) x^3/3 + ..., so a number of relations to the Faber polynomials of A263916 hold with p(0,x) = 2:
1) p(n,x) = F(n,(2-x),1,0,0,..)
2) p(n,x) = (-1)^n 2 + F(n,-x,2x,-3x,...,(-1)^n n*x)
3) p(n,x) = (-1)^n [2 + F(n,x,2x,3x,...,n*x)].
The unsigned array contains the partial sums of A111125 modified by appending a column of zeros, except for an initial two, to A111125. Then the difference of consecutive rows of unsigned A127677, further modified by appending an initial rows of zeros, generates the modified A111125. Cf. A208513 and A034807.
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. 12, 20, and 21) and Damianou and Evripidou (p. 7).
See A111125 for a relation to the squares of the odd row polynomials here with the constant removed.
p(n,x)^2 = 2 + p(2*n,x). See also A127672. (End)
a(n,m) = -2*a(n-1,m) + a(n-1,m-1) - a(n-2,m) for n >= 2 with initial conditions a(0,0) = 2, a(1,0) = -2, a(1,1) = 1, a(0,m) = 0 for m != 0, a(1,m) = 0 for m != 0,1. - William P. Orrick, Jun 09 2020
p(n,x) = (x-2)*p(n-1,x) - p(n-2,x) for n >= 2. - William P. Orrick, Jun 09 2020

Extensions

Definition corrected by Eric W. Weisstein, Apr 06 2017
Previous Showing 11-20 of 57 results. Next