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 21-30 of 34 results. Next

A055243 First differences of A001628 (Fibonacci convolution).

Original entry on oeis.org

1, 2, 6, 13, 29, 60, 122, 241, 468, 894, 1686, 3144, 5807, 10636, 19338, 34931, 62731, 112068, 199264, 352787, 622152, 1093260, 1914780, 3343440, 5821645, 10110278, 17515566, 30276073, 52221929, 89896332, 154461110, 264930661, 453654108, 775598634, 1324053522
Offset: 0

Views

Author

Wolfdieter Lang, May 10 2000

Keywords

Comments

2*a(n) = C_{n+3} of Turban reference eq.(2.17), C_{1}= 0 = C_{2}.
Number of binary sequences of length n+3 such that the sequence has exactly two pairs (which may overlap) of consecutive 1's. - George J. Schaeffer (gschaeff(AT)andrew.cmu.edu), Sep 07 2004

Crossrefs

Programs

  • Maple
    a:= n -> (Matrix([[1,0$4,-1]]). Matrix(6, (i,j)-> if (i=j-1) then 1 elif j=1 then [3,0,-5,0,3,1][i] else 0 fi)^(n))[1,1]: seq(a(n), n=0..30); # Alois P. Heinz, Aug 05 2008
  • Mathematica
    Differences[LinearRecurrence[{3,0,-5,0,3,1},{0,1,3,9,22,51,111},40]] (* Harvey P. Dale, Jun 12 2019 *)

Formula

G.f.: (1-x)/(1-x-x^2)^3. (from Turban reference eq.(2.15)).
a(n)= ((5*n^2+37*n+50)*F(n+1)+4*(n+1)*F(n))/50 with F(n)=A000045(n) (Fibonacci numbers) (from Turban reference eq. (2.17)).
From Peter Bala, Oct 25 2007 (Start):
Since F(-n) = (-1)^(n+1)*F(n), we can use the previous formula to extend the sequence to negative values of n; we find a(-n) = (-1)^n* A129707(n-3).
Recurrence relations: a(n+4) = 2*a(n+3) + a(n+2) - 2*a(n+1) - a(n) + F(n+3), with a(0) = 1, a(1) = 2, a(2) = 6 and a(3) = 13;
a(n+2) = a(n+1) + a(n) + A010049(n+3), with a(0) = 1 and a(1) = 2.
a(n-3) = Sum_{k = 2..floor((n+1)/2)} C(k,2)*C(n-k,k-1) = (1/2)*G''(n,1), where the polynomial G(n,x) := Sum_{k = 1..floor((n+1)/2)} C(n-k,k-1)*x^k = x^((n+1)/2) * F(n, 1/sqrt(x)) and where F(n,x) denotes the n-th Fibonacci polynomial. Since G(n,1) yields the Fibonacci numbers A000045 and G'(n,1) yields the second-order Fibonacci numbers A010049, a(n) may be considered as the sequence of third-order Fibonacci numbers.
For n >= 4, the polynomials Sum_{k = 0..n} C(n,k) * G''(n-k,1)*(-x)^k appear to satisfy a Riemann hypothesis; their zeros appear to lie on the vertical line Re x = 1/2 in the complex plane. Compare with the remarks in A094440 and A010049. (End)
a(n) = A076791(n+3, 2). - Michael Somos, Sep 24 2024
E.g.f.: exp(x/2)*(5*(25 + 23*x + 5*x^2)*cosh(sqrt(5)*x/2) + sqrt(5)*(29 + 65*x + 10*x^2)*sinh(sqrt(5)*x/2))/125. - Stefano Spezia, Sep 26 2024

A368150 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 1 + 3*x, p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where u = p(2,x), v = 1 - x^2.

Original entry on oeis.org

1, 1, 3, 2, 6, 8, 3, 15, 25, 21, 5, 30, 76, 90, 55, 8, 60, 188, 324, 300, 144, 13, 114, 439, 948, 1251, 954, 377, 21, 213, 961, 2529, 4207, 4527, 2939, 987, 34, 390, 2026, 6246, 12606, 17154, 15646, 8850, 2584, 55, 705, 4136, 14640, 34590, 56970, 65840
Offset: 1

Views

Author

Clark Kimberling, Dec 25 2023

Keywords

Comments

Because (p(n,x)) is a strong divisibility sequence, for each integer k, the sequence (p(n,k)) is a strong divisibility sequence of integers.

Examples

			First eight rows:
   1
   1    3
   2    6    8
   3   15   25    21
   5   30   76    90     55
   8   60  188   324    300   144
  13  114  439   948   1251   954   377
  21  213  961  2529   4207  4527  2939   987
Row 4 represents the polynomial p(4,x) = 3 + 15*x + 25*x^2 + 21*x^3, so (T(4,k)) = (3,15,25,21), k=0..3.
		

Crossrefs

Cf. A000045 (column 1); A001906 (p(n,n-1)); A000302 (row sums), (p(n,1)); A122803 (alternating row sums), (p(n,-1)); A190972 (p(n,2)), A116415, (p(n,-2)); A190990, (p(n,3)); A057084, (p(n,-3)); A094440, A367208, A367209, A367210, A367211, A367297, A367298, A367299, A367300, A367301, A368151.

Programs

  • Mathematica
    p[1, x_] := 1; p[2, x_] := 1 + 3 x; u[x_] := p[2, x]; v[x_] := 1 - x^2;
    p[n_, x_] := Expand[u[x]*p[n - 1, x] + v[x]*p[n - 2, x]]
    Grid[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]
    Flatten[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]

Formula

p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where p(1,x) = 1, p(2,x) = 1 + 3*x, u = p(2,x), and v = 1 - x^2.
p(n,x) = k*(b^n - c^n), where k = -1/sqrt(5 + 6*x + 5*x^2), b = (1/2)*(3*x + 1 - 1/k), c = (1/2)*(3*x + 1 + 1/k).

A368151 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 1 + 3x, p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >=3, where u = p(2,x), v = 2 - x^2.

Original entry on oeis.org

1, 1, 3, 3, 6, 8, 5, 21, 25, 21, 11, 48, 101, 90, 55, 21, 123, 290, 414, 300, 144, 43, 282, 850, 1416, 1551, 954, 377, 85, 657, 2255, 4671, 6109, 5481, 2939, 987, 171, 1476, 5883, 13986, 22374, 24300, 18585, 8850, 2584, 341, 3303, 14736, 40320, 74295, 97713
Offset: 1

Views

Author

Clark Kimberling, Dec 31 2023

Keywords

Comments

Because (p(n,x)) is a strong divisibility sequence, for each integer k, the sequence (p(n,k)) is a strong divisibility sequence of integers.

Examples

			First eight rows:
   1
   1     3
   3     6    8
   5    21    25    21
  11    48   101    90    55
  21   123   290   414   300  144
  43   282   850  1416  1551  954    377
  85   657  2255  4671  6109  5481  2939  987
Row 4 represents the polynomial p(4,x) = 5 + 21 x + 25 x^2 + 21 x^3, so (T(4,k)) = (5,21,25,21), k=0..3.
		

Crossrefs

Cf. A001045 (column 1); A001906 (p(n,n-1)); A001076 (row sums), (p(n,1)); A077985 (alternating row sums), (p(n,-1)); A186446 (p(n,2)), A107839, (p(n,-2)); A190989, (p(n,3)); A023000, (p(n,-3)); A094440, A367208, A367209, A367210, A367211, A367297, A367298, A367299, A367300, A367301, A368150.

Programs

  • Mathematica
    p[1, x_] := 1; p[2, x_] := 1 + 3 x; u[x_] := p[2, x]; v[x_] := 2 - x^2;
    p[n_, x_] := Expand[u[x]*p[n - 1, x] + v[x]*p[n - 2, x]]
    Grid[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]
    Flatten[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]

Formula

p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >=3, where p(1,x) = 1, p(2,x) = 1 + 3 x, u = p(2,x), and v = 2 - x^2.
p(n,x) = k*(b^n - c^n), where k = -1/sqrt(9 + 6 x + 5 x^2), b = (1/2) (3 x + 1 - 1/k), c = (1/2) (3 x + 1 + 1/k).

A326925 Irregular triangular array read by rows: row n shows the coefficients of this polynomial of degree n: (1/n!)*(numerator of n-th derivative of (1-x)/(1-x-x^2)).

Original entry on oeis.org

1, -1, 0, 2, -1, 1, 0, 3, -1, 1, 4, 0, 4, -1, 2, 5, 10, 0, 5, -1, 3, 12, 15, 20, 0, 6, -1, 5, 21, 42, 35, 35, 0, 7, -1, 8, 40, 84, 112, 70, 56, 0, 8, -1, 13, 72, 180, 252, 252, 126, 84, 0, 9, -1, 21, 130, 360, 600, 630, 504, 210, 120, 0, 10, -1, 34, 231, 715
Offset: 1

Views

Author

Clark Kimberling, Oct 22 2019

Keywords

Comments

Column 1: Fibonacci numbers, F(m), for m >= -1, as in A000045. For n >= 0, the n-th row sum = F(2n), as in A001906.
Conjecture: The odd degree polynomials are irreducible; the even degree (= 2k) polynomials have exactly two irreducible factors, each of degree k.

Examples

			First 7 rows:
1    -1
0     2   -1
1     0    3   -1
1     4    0    4   -1
2     5    0   10    5   -1
3    12   15   20    0    6   -1
5    21   42   35   35    0    7   -1
First 7 polynomials:
1 - x
2 x - x^2
1 + 3 x^2 - x^3
1 + 4 x + 4 x^3 - x^4
2 + 5 x + 10 x^2 + 5 x^4 - x^5
3 + 12 x + 15 x^2 + 20 x^3 + 6 x^5 - x^6
5 + 21 x + 42 x^2 + 35 x^3 + 35 x^4 + 7 x^6 - x^7
Factorizations of even-degree polynomials:
degree 2:  (2 - x)*x
degree 4:  (1 + x^2)*(1 + 4x - x^2)
degree 6:  (1 + 3x + x^3)*(3 + 3x + 6x^2 - x^3)
degree 8:  (2 + 4x + 6x^2 + x^4)*(4 + 12 x + 6x^2 + 8x^3 - x^4)
degree 10: (3 + 10 x + 10 x^2 + 10 x^3 + x^5)*(7 + 20 x + 30 x^2 + 10 x^3 + 10 x^4 - x^5)
It appears that the constant terms of the factors are Fibonacci numbers (A000045) and Lucas numbers (A000032).
		

Crossrefs

Programs

  • Mathematica
    g[x_, n_] := Numerator[(-1)^(n + 1) Factor[D[(1 - x)/(1 - x - x^2), {x, n}]]]
    Column[Expand[Table[g[x, n]/n!, {n, 0, 12}]]] (* polynomials *)
    h[n_] := CoefficientList[g[x, n]/n!, x] (* A326925 *)
    Table[h[n], {n, 0, 10}]
    Column[%]

Formula

G.f. as array: ((y^2 + y - 1)*x - y + 1)/(1 + (y^2 + y - 1)*x^2 + (-2*y - 1)*x). - Robert Israel, Oct 31 2019

A368152 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 1 + 3*x, p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where u = p(2,x), v = 3 - x^2.

Original entry on oeis.org

1, 1, 3, 4, 6, 8, 7, 27, 25, 21, 19, 66, 126, 90, 55, 40, 204, 392, 504, 300, 144, 97, 522, 1363, 1884, 1851, 954, 377, 217, 1425, 4065, 7281, 8011, 6435, 2939, 987, 508, 3642, 12332, 24606, 34044, 31446, 21524, 8850, 2584, 1159, 9441, 35236, 82020, 127830
Offset: 1

Views

Author

Clark Kimberling, Jan 20 2024

Keywords

Comments

Because (p(n,x)) is a strong divisibility sequence, for each integer k, the sequence (p(n,k)) is a strong divisibility sequence of integers.

Examples

			First eight rows:
    1
    1    3
    4    6    8
    7   27   25   21
   19   66  126   90   55
   40  204  392  504  300  144
   97  522 1363 1884 1851  954  377
  217 1425 4065 7281 8011 6435 2939 987
Row 4 represents the polynomial p(4,x) = 7 + 27*x + 25*x^2 + 21*x^3, so (T(4,k)) = (7,27,25,21), k=0..3.
		

Crossrefs

Cf. A006130 (column 1); A001906 (p(n,n-1)); A090017 (row sums), (p(n,1)); A002605 (alternating row sums), (p(n,-1)); A004187, (p(n,2)); A004254, (p(n,-2)); A190988, (p(n,3)); A190978 (unsigned), (p(n,-3)); A094440, A367208, A367209, A367210, A367211, A367297, A367298, A367299, A367300, A367301, A368150, A368151.

Programs

  • Mathematica
    p[1, x_] := 1; p[2, x_] := 1 + 3 x; u[x_] := p[2, x]; v[x_] := 3 - x^2;
    p[n_, x_] := Expand[u[x]*p[n - 1, x] + v[x]*p[n - 2, x]]
    Grid[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]
    Flatten[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]

Formula

p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where p(1,x) = 1, p(2,x) = 1 + 3*x, u = p(2,x), and v = 3 - x^2.
p(n,x) = k*(b^n - c^n), where k = -1/sqrt(13 + 6*x + 5*x^2), b = (1/2)*(3*x + 1 - 1/k), c = (1/2)*(3*x + 1 + 1/k).

A131048 (1/3) * (A007318^2 - A007318^(-1)).

Original entry on oeis.org

1, 1, 2, 3, 3, 3, 5, 12, 6, 4, 11, 25, 30, 10, 5, 21, 66, 75, 60, 15, 6, 43, 147, 231, 175, 105, 21, 7, 85, 344, 588, 616, 350, 168, 28, 8, 171, 765, 1548, 1764, 1386, 630, 252, 36, 9
Offset: 1

Views

Author

Gary W. Adamson, Jun 12 2007

Keywords

Comments

Left border = A001045: (1, 1, 3, 5, 11, 21, 43, 85, ...).
Row sums = (1, 3, 9, 27, ...).
Analogous triangles for other powers of P are: A131047, A131049, A131050 and A131051.

Examples

			First few rows of the triangle:
   1;
   1,   2;
   3,   3,   3;
   5,  12,   6,   4;
  11,  25,  30,  10,   5;
  21,  66,  75,  60,  15,  6;
  43, 147, 231, 175, 105, 21, 7;
  ...
		

Crossrefs

Formula

Let A007318 (Pascal's triangle) = P. then A131048 = (1/3) * (P^2 - 1/P). Delete right border of zeros.
From Peter Bala, Oct 24 2007: (Start)
O.g.f.: 1/(1 - (2*x + 1)*t + (x^2 + x - 2)*t^2) = 1 + (1 + 2*x)*t + (3 + 3*x + 3*x^2)*t^2 + ....
T(n,n-k) = (1/3)*C(n,k)*(2^k - (-1)^k) = C(n,k)*A001045(k).
The row polynomials R(n,x) := Sum_{k = 0..n} T(n,n-k)*x^(n-k) = (1/3)*((x + 2)^n - (x - 1)^n) and have the divisibility property R(n,x) divides R(m,x) in the polynomial ring Z[x] if n divides m.
The polynomials R(n,-x), n >= 2, satisfy a Riemann hypothesis: their zeros lie on the vertical line Re x = 1/2 in the complex plane. Compare with A094440. (End)

A368153 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 1 + 2*x, p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where u = p(2,x), v = 1 - 3*x - x^2.

Original entry on oeis.org

1, 1, 2, 2, 1, 3, 3, 4, -2, 4, 5, 5, 4, -10, 5, 8, 10, -3, 4, -25, 6, 13, 16, 1, -29, 14, -49, 7, 21, 28, -8, -24, -78, 56, -84, 8, 34, 47, -12, -88, -26, -162, 168, -132, 9, 55, 80, -31, -140, -200, 100, -330, 408, -195, 10, 89, 135, -58, -301, -230, -296
Offset: 1

Views

Author

Clark Kimberling, Jan 20 2024

Keywords

Comments

Because (p(n,x)) is a strong divisibility sequence, for each integer k, the sequence (p(n,k)) is a strong divisibility sequence of integers.

Examples

			First eight rows:
   1
   1   2
   2   1   3
   3   4  -2    4
   5   5   4  -10    5
   8  10  -3    4  -25    6
  13  16   1  -29   14  -49    7
  21  28  -8  -24  -78   56  -84   8
Row 4 represents the polynomial p(4,x) = 3 + 4*x - 2*x^2 + 4*x^3, so (T(4,k)) = (3,4,-2,4), k=0..3.
		

Crossrefs

Cf. A000045 (column 1); A000027 (p(n,n-1)); A057083 (row sums), (p(n,1)); A182228 (alternating row sums), (p(n,-1)); A190970, (p(n,2)); A030195, (p(n,-2)); A052918, (p(n,-3)); A190972, (p(n,-4)); A057085, (p(n,-5)); A094440, A367208, A367209, A367210, A367211, A367297, A367298, A367299, A367300, A367301, A368150, A368151, A368152.

Programs

  • Mathematica
    p[1, x_] := 1; p[2, x_] := 1 + 2 x; u[x_] := p[2, x]; v[x_] := 1 - 3x - x^2;
    p[n_, x_] := Expand[u[x]*p[n - 1, x] + v[x]*p[n - 2, x]]
    Grid[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]
    Flatten[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]

Formula

p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where p(1,x) = 1, p(2,x) = 1 + 2*x, u = p(2,x), and v = 1 - 3*x - x^2.
p(n,x) = k*(b^n - c^n), where k = -1/sqrt(5 - 8*x), b = (1/2)*(2*x + 1 - 1/k), c = (1/2)*(2*x + 1 + 1/k).

A132148 Triangular array T(n,k) = C(n,k)*Lucas(n-k), 0 <= k <= n.

Original entry on oeis.org

2, 1, 2, 3, 2, 2, 4, 9, 3, 2, 7, 16, 18, 4, 2, 11, 35, 40, 30, 5, 2, 18, 66, 105, 80, 45, 6, 2, 29, 126, 231, 245, 140, 63, 7, 2, 47, 232, 504, 616, 490, 224, 84, 8, 2, 76, 423, 1044, 1512, 1386, 882, 336, 108, 9, 2, 123, 760, 2115, 3480, 3780, 2772, 1470, 480, 135, 10, 2
Offset: 0

Views

Author

Peter Bala, Aug 17 2007

Keywords

Examples

			Triangle starts
   2;
   1,   2;
   3,   2,    2;
   4,   9,    3,   2;
   7,  16,   18,   4,   2;
  11,  35,   40,  30,   5,  2;
  18,  66,  105,  80,  45,  6,  2;
  ...
		

Crossrefs

Cf. A000032 (T(n,0)), A000045, A005248 (row sums), A061084 (alter. row sums), A094440.

Programs

  • Maple
    with(combinat): lucas := n -> fibonacci(n-1) + fibonacci(n+1): T := (n,k) -> binomial(n,k)*lucas(n-k): for n from 0 to 10 do seq( T(n,k), k = 0..n) od;
  • Mathematica
    Flatten[Table[Binomial[n,k]LucasL[n-k],{n,0,10},{k,0,n}]] (* Harvey P. Dale, Nov 06 2011 *)

Formula

G.f.: (2 - (2x + 1)*t)/(1 - (2x + 1)*t + (x^2 + x - 1)*t^2) = 2 + (1 + 2*x)*t + (3 + 2*x + 2*x^2)*t^2 + (4 + 9*x + 3*x^2 + 2*x^3)*t^3 + ... .
The row polynomials L(n,x) = Sum_{k = 0..n} C(n,k)*Lucas(n-k)*x^k satisfy L(n,x)*F(n,x) = F(2*n,x), where F(n,x) = Sum_{k = 0..n} C(n,k)*Fibonacci(n-k)*x^k.
Other identities and formulas include:
L(n+1,x)^2 - L(n,x)*L(n+2,x) = -5*(x^2 + x - 1)^n;
L(n+1,x) - (x^2 + x - 1)*L(n-1,x) = 5*F(n,x) for n >= 1;
L(2*n,x) - 2*(x^2 + x - 1)^n = 5*F(n,x)^2;
L(n,2*x) = Sum_{ k = 0..n} C(n,k)*L(n-k,x)*x^k;
L(n,3*x) = Sum_{ k = 0..n} C(n,k)*L(n-k,2*x)*x^k etc.
Sum_{k = 0..n} C(n,k)*L(k,x)*F(n-k,x) = 2^n*F(n,x).
Row sums: L(n,1) = Lucas(2*n); Alternating row sums: L(n,-1) = (-1)^n*Lucas(n); L(n,1/phi) = (-1)^n*L(n,-phi) = sqrt(5)^n for n >= 1, where phi = (1 + sqrt(5))/2.
The polynomials L(n,-x) satisfy a Riemann hypothesis: the zeros of L(n,-x) lie on the vertical line Re x = 1/2 in the complex plane.
From Peter Bala, Jun 29 2016: (Start)
L(n,x) = (x + phi)^n + (x - 1/phi)^n, where phi = (1 + sqrt(5))/2. The zeros of L(n,x) are given by -1/2 - i*sqrt(5)/2*cot( (2*k - 1)*Pi/(2*n) ) for k = 1..n.
d/dx(L(n,x)) = n*L(n-1,x).
L(-n,x) = L(n,x)/(x^2 + x - 1)^n.
L(n,x - 1) = (-1)^n*L(n,-x).
L(n,x)^2 - L(2*n,x) = 2*(x^2 + x - 1)^n.
L(n,x)^3 - L(3*n,x) = 3*(x^2 + x - 1)^n*L(n,x).
L(n,x)^4 - L(4*n,x) = 4*(x^2 + x - 1)^n*L(n,x)^2 - 2*(x^2 + x - 1)^(2*n).
If n divides m and m/n is odd then L(n,x) divides L(m,x) in the polynomial ring Z[x].
L(n,x) = F(n+1,x) - (x^2 + x - 1)*F(n-1,x) = 2*F(n+1,x) - (2*x + 1)*F(n,x), where F(n,x) = Sum_{k = 0..n} binomial(n,k)*Fibonacci(n-k)*x^k denotes the n-th row polynomial of A094440 (taken with an offset of 0).
exp( Sum_{n >= 1} L(n,x)*z^n/n ) = Sum_{n >= 0} F(n+1,x)*z^n.
exp( Sum_{n >= 1} L(n,x)*L(2*n,x)*z^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)*z^n. (End)
From Peter Bala, Dec 31 2023: (Start)
For n >= 1, the n-th row polynomial L(n,x) is the numerator of 1/(n-1)! * (d/dx)^(n-1) (2*x + 1)/(1 - x - x^2).
Recurrence: for n >= 1, n*L(n+1,x) = n*(2*x + 1)*L(n,x) + (1 - x - x^2)* d/dx(L(n,x)) with L(1,x) = 2*x + 1. (End)

A368155 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 1 + 3*x, p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where u = p(2,x), v = 1 - 3*x - 2*x^2.

Original entry on oeis.org

1, 1, 3, 2, 3, 7, 3, 9, 5, 15, 5, 15, 26, 3, 31, 8, 30, 43, 63, -15, 63, 13, 54, 104, 87, 144, -81, 127, 21, 99, 203, 273, 115, 333, -275, 255, 34, 177, 416, 549, 609, -9, 806, -789, 511, 55, 315, 811, 1263, 1146, 1260, -725, 2043, -2071, 1023, 89, 555, 1573
Offset: 1

Views

Author

Clark Kimberling, Jan 20 2024

Keywords

Comments

Because (p(n,x)) is a strong divisibility sequence, for each integer k, the sequence (p(n,k)) is a strong divisibility sequence of integers.

Examples

			First eight rows:
   1
   1    3
   2    3     7
   3    9     5    15
   5   15    26     3    31
   8   30    43    63   -15    63
  13   54   104    87   144   -81    127
  21   99   203   273   115   333   -275   255
Row 4 represents the polynomial p(4,x) = 3 + 9*x + 5*x^2 + 15*x^3, so (T(4,k)) = (3,9,5,15), k=0..3.
		

Crossrefs

Cf. A000045 (column 1); A000225, (p(n,n-1)); A001787 (row sums), (p(n,1)); A002605 (alternating row sums), (p(n,-1)); A004254, (p(n,-2)); A057084, (p(n,-3)); A094440, A367208, A367209, A367210, A367211, A367297, A367298, A367299, A367300, A367301, A368150, A368151, A368152, A368153, A368154, A368156.

Programs

  • Mathematica
    p[1, x_] := 1; p[2, x_] := 1 + 3 x; u[x_] := p[2, x]; v[x_] := 1 - 3x - 2x^2;
    p[n_, x_] := Expand[u[x]*p[n - 1, x] + v[x]*p[n - 2, x]]
    Grid[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]
    Flatten[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]

Formula

p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where p(1,x) = 1, p(2,x) = 1 + 3*x, u = p(2,x), and v = 1 - 3*x - 2*x^2.
p(n,x) = k*(b^n - c^n), where k = -1/sqrt(5 - 6*x + x^2), b = (1/2)*(3*x + 1 - 1/k), c = (1/2)*(2*x + 1 + 1/k).

A368156 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 1 + 2*x, p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where u = p(2,x), v = 1 + x^2.

Original entry on oeis.org

1, 1, 2, 2, 4, 5, 3, 10, 14, 12, 5, 20, 41, 44, 29, 8, 40, 98, 148, 131, 70, 13, 76, 224, 408, 497, 376, 169, 21, 142, 482, 1044, 1542, 1588, 1052, 408, 34, 260, 1003, 2492, 4351, 5456, 4894, 2888, 985, 55, 470, 2026, 5684, 11359, 16790, 18400, 14672, 7813
Offset: 1

Views

Author

Clark Kimberling, Jan 20 2024

Keywords

Comments

Because (p(n,x)) is a strong divisibility sequence, for each integer k, the sequence (p(n,k)) is a strong divisibility sequence of integers.

Examples

			First eight rows:
   1
   1    2
   2    4    5
   3   10   14    12
   5   20   41    44    29
   8   40   98   148   131    70
  13   76  224   408   497   376   169
  21  142  482  1044  1542  1588  1052  408
Row 4 represents the polynomial p(4,x) = 3 + 10*x + 14*x^2 + 12*x^3, so (T(4,k)) = (3,10,14,12), k=0..3.
		

Crossrefs

Cf. A000045 (column 1); A000129, (p(n,n-1)); A007482 (row sums), (p(n,1)); A077925 (alternating row sums), (p(n,-1)); A057088, (p(n,2)); A015523, (p(n,-2)); A015568, (p(n,3)); A180250, (p(n,-3)); A094440, A367208, A367209, A367210, A367211, A367297, A367298, A367299, A367300, A367301, A368150, A368151, A368152, A368153, A368154, A368155.

Programs

  • Mathematica
    p[1, x_] := 1; p[2, x_] := 1 + 2 x; u[x_] := p[2, x]; v[x_] := 1 + 2x^2;
    p[n_, x_] := Expand[u[x]*p[n - 1, x] + v[x]*p[n - 2, x]]
    Grid[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]
    Flatten[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]

Formula

p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where p(1,x) = 1, p(2,x) = 1 + 2*x, u = p(2,x), and v = 1 + x^2.
p(n,x) = k*(b^n - c^n), where k = -1/sqrt(5 + 4*x + 8*x^2), b = (1/2)*(2*x + 1 - 1/k), c = (1/2)*(2*x + 1 + 1/k).
Previous Showing 21-30 of 34 results. Next