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-23 of 23 results.

A015544 Lucas sequence U(5,-8): a(n+1) = 5*a(n) + 8*a(n-1), a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 5, 33, 205, 1289, 8085, 50737, 318365, 1997721, 12535525, 78659393, 493581165, 3097180969, 19434554165, 121950218577, 765227526205, 4801739379641, 30130517107845, 189066500576353, 1186376639744525, 7444415203333449, 46713089134623445
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [n le 2 select n-1 else 5*Self(n-1) + 8*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 13 2012
    
  • Mathematica
    a[n_]:=(MatrixPower[{{1,2},{1,-6}},n].{{1},{1}})[[2,1]]; Table[Abs[a[n]],{n,-1,40}] (* Vladimir Joseph Stephan Orlovsky, Feb 19 2010 *)
    LinearRecurrence[{5, 8}, {0, 1}, 30] (* Vincenzo Librandi, Nov 13 2012 *)
  • PARI
    A015544(n)=imag((2+quadgen(57))^n) \\ M. F. Hasler, Mar 06 2009
    
  • PARI
    x='x+O('x^30); concat([0], Vec(x/(1 - 5*x - 8*x^2))) \\ G. C. Greubel, Jan 01 2018
  • Sage
    [lucas_number1(n,5,-8) for n in range(0, 21)] # Zerinvary Lajos, Apr 24 2009
    

Formula

a(n) = 5*a(n-1) + 8*a(n-2).
G.f.: x/(1 - 5*x - 8*x^2). - M. F. Hasler, Mar 06 2009

Extensions

More precise definition by M. F. Hasler, Mar 06 2009

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).

A206800 Riordan array (1/(1-3*x+x^2), x*(1-x)/(1-3*x+x^2)).

Original entry on oeis.org

1, 3, 1, 8, 5, 1, 21, 19, 7, 1, 55, 65, 34, 9, 1, 144, 210, 141, 53, 11, 1, 377, 654, 534, 257, 76, 13, 1, 987, 1985, 1905, 1111, 421, 103, 15, 1, 2584, 5911, 6512, 4447, 2041, 641, 134, 17, 1, 6765, 17345, 21557, 16837, 9038, 3440, 925, 169, 19, 1
Offset: 0

Views

Author

Philippe Deléham, Feb 12 2012

Keywords

Examples

			Triangle begins :
1
3, 1
8, 5, 1
21, 19, 7, 1
55, 65, 34, 9, 1
144, 210, 141, 53, 11, 1
377, 654, 534, 257, 76, 13, 1
987, 1985, 1905, 1111, 421, 103, 15, 1
2584, 5911, 6512, 4447, 2041, 641, 134, 17, 1
6765, 17345, 21557, 16837, 9038, 3440, 925, 169, 19, 1
Triangle (0,3,-1/3,1/3,0,0,0,0,0,...) DELTA (1,0,-1/3,1/3,0,0,0,0,...) begins :
1
0, 1
0, 3, 1
0, 8, 5, 1
0, 21, 19, 7, 1
0, 55, 65, 34, 9, 1...
		

References

  • Subtriangle of the triangle given by (0, 3, -1/3, 1/3, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, -1/3, 1/3, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938.
  • Antidiagonal sums are A072264(n).

Crossrefs

Formula

T(n,k) = 3*T(n-1,k) + T(n-1,k-1) - T(n-2,k) - T(n-2,k-1).
G.f.: 1/(1-(y+3)*x+(y+1)*x^2).
Sum_{k, 0<=k<=n} T(n,k)*x^k = (-1)^n* A015587(n+1), (-1)^n*A190953(n+1), (-1)^n*A015566(n+1), (-1)*A189800(n+1), (-1)^n*A015541(n+1), (-1)^n*A085939(n+1), (-1)^n*A015523(n+1), (-1)^n*A063727(n), (-1)^n*A006130(n), A077957(n), A000045(n+1), A000079(n), A001906(n+1), A007070(n), A116415(n), A084326(n+1), A190974(n+1), A190978(n+1), A190984(n+1), A190990(n+1), A190872(n) for x = -12, -11, -10, -9, -8, -7, -6, -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8 respectively.
Previous Showing 21-23 of 23 results.