cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A190958 a(n) = 2*a(n-1) - 10*a(n-2), with a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 2, -6, -32, -4, 312, 664, -1792, -10224, -2528, 97184, 219648, -532544, -3261568, -1197696, 30220288, 72417536, -157367808, -1038910976, -504143872, 9380822016, 23803082752, -46202054656, -330434936832, -198849327104, 2906650714112, 7801794699264
Offset: 0

Views

Author

Keywords

Comments

For the difference equation a(n) = c*a(n-1) - d*a(n-2), with a(0) = 0, a(1) = 1, the solution is a(n) = d^((n-1)/2) * ChebyshevU(n-1, c/(2*sqrt(d))) and has the alternate form a(n) = ( ((c + sqrt(c^2 - 4*d))/2)^n - ((c - sqrt(c^2 - 4*d))/2)^n )/sqrt(c^2 - 4*d). In the case c^2 = 4*d then the solution is a(n) = n*d^((n-1)/2). The generating function is x/(1 - c*x + d^2) and the exponential generating function takes the form (2/sqrt(c^2 - 4*d))*exp(c*x/2)*sinh(sqrt(c^2 - 4*d)*x/2) for c^2 > 4*d, (2/sqrt(4*d - c^2))*exp(c*x/2)*sin(sqrt(4*d - c^2)*x/2) for 4*d > c^2, and x*exp(sqrt(d)*x) if c^2 = 4*d. - G. C. Greubel, Jun 10 2022

Crossrefs

Programs

  • Magma
    I:=[0,1]; [n le 2 select I[n] else 2*Self(n-1)-10*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Sep 17 2011
    
  • Mathematica
    LinearRecurrence[{2,-10}, {0,1}, 50]
  • PARI
    a(n)=([0,1; -10,2]^n*[0;1])[1,1] \\ Charles R Greathouse IV, Apr 08 2016
    
  • SageMath
    [lucas_number1(n,2,10) for n in (0..50)] # G. C. Greubel, Jun 10 2022

Formula

G.f.: x / ( 1 - 2*x + 10*x^2 ). - R. J. Mathar, Jun 01 2011
E.g.f.: (1/3)*exp(x)*sin(3*x). - Franck Maminirina Ramaharo, Nov 13 2018
a(n) = 10^((n-1)/2) * ChebyshevU(n-1, 1/sqrt(10)). - G. C. Greubel, Jun 10 2022
a(n) = (1/3)*10^(n/2)*sin(n*arctan(3)) = Sum_{k=0..floor(n/2)} (-1)^k*3^(2*k)*binomial(n,2*k+1). - Gerry Martens, Oct 15 2022

A207538 Triangle of coefficients of polynomials v(n,x) jointly generated with A207537; see Formula section.

Original entry on oeis.org

1, 2, 4, 1, 8, 4, 16, 12, 1, 32, 32, 6, 64, 80, 24, 1, 128, 192, 80, 8, 256, 448, 240, 40, 1, 512, 1024, 672, 160, 10, 1024, 2304, 1792, 560, 60, 1, 2048, 5120, 4608, 1792, 280, 12, 4096, 11264, 11520, 5376, 1120, 84, 1, 8192, 24576, 28160, 15360
Offset: 1

Views

Author

Clark Kimberling, Feb 18 2012

Keywords

Comments

As triangle T(n,k) with 0<=k<=n and with zeros omitted, it is the triangle given by (2, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, 1/2, -1/2, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Mar 04 2012
The numbers in rows of the triangle are along "first layer" skew diagonals pointing top-left in center-justified triangle given in A013609 ((1+2*x)^n) and along (first layer) skew diagonals pointing top-right in center-justified triangle given in A038207 ((2+x)^n), see links. - Zagros Lalo, Jul 31 2018
If s(n) is the row sum at n, then the ratio s(n)/s(n-1) is approximately 2.414213562373095... (A014176: Decimal expansion of the silver mean, 1+sqrt(2)), when n approaches infinity. - Zagros Lalo, Jul 31 2018

Examples

			First seven rows:
1
2
4...1
8...4
16..12..1
32..32..6
64..80..24..1
(2, 0, 0, 0, 0, ...) DELTA (0, 1/2, -1/2, 0, 0, 0, ...) begins:
    1
    2,   0
    4,   1,  0
    8,   4,  0, 0
   16,  12,  1, 0, 0
   32,  32,  6, 0, 0, 0
   64,  80, 24, 1, 0, 0, 0
  128, 192, 80, 8, 0, 0, 0, 0
		

References

  • Shara Lalo and Zagros Lalo, Polynomial Expansion Theorems and Number Triangles, Zana Publishing, 2018, ISBN: 978-1-9995914-0-3, pp. 80-83, 357-358.

Crossrefs

Programs

  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := u[n - 1, x] + (x + 1)*v[n - 1, x]
    v[n_, x_] := u[n - 1, x] + v[n - 1, x]
    Table[Factor[u[n, x]], {n, 1, z}]
    Table[Factor[v[n, x]], {n, 1, z}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%]  (* A207537, |A028297| *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]  (* A207538, |A133156| *)
    t[0, 0] = 1; t[n_, k_] := t[n, k] = If[n < 0 || k < 0, 0, 2 t[n - 1, k] + t[n - 2, k - 1]]; Table[t[n, k], {n, 0, 15}, {k, 0, Floor[n/2]}] // Flatten (* Zagros Lalo, Jul 31 2018 *)
    t[n_, k_] := t[n, k] = 2^(n - 2 k) * (n -  k)!/((n - 2 k)! k!) ; Table[t[n, k], {n, 0, 15}, {k, 0, Floor[n/2]} ]  // Flatten (* Zagros Lalo, Jul 31 2018 *)

Formula

u(n,x) = u(n-1,x)+(x+1)*v(n-1,x), v(n,x) = u(n-1,x)+v(n-1,x), where u(1,x) = 1, v(1,x) = 1. Also, A207538 = |A133156|.
From Philippe Deléham, Mar 04 2012: (Start)
With 0<=k<=n:
Mirror image of triangle in A099089.
Skew version of A038207.
Riordan array (1/(1-2*x), x^2/(1-2*x)).
G.f.: 1/(1-2*x-y*x^2).
Sum_{k, 0<=k<=n} T(n,k)*x^k = A190958(n+1), A127357(n), A090591(n), A089181(n+1), A088139(n+1), A045873(n+1), A088138(n+1), A088137(n+1), A099087(n), A000027(n+1), A000079(n), A000129(n+1), A002605(n+1), A015518(n+1), A063727(n), A002532(n+1), A083099(n+1), A015519(n+1), A003683(n+1), A002534(n+1), A083102(n), A015520(n+1), A091914(n) for x = -10, -9, -8, -7, -6, -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 respectively.
T(n,k) = 2*T(n-1,k) + T(-2,k-1) with T(0,0) = 1, T(1,0) = 2, T(1,1) = 0 and T(n, k) = 0 if k<0 or if k>n. (End)
T(n,k) = A013609(n-k, n-2*k+1). - Johannes W. Meijer, Sep 05 2013
From Tom Copeland, Feb 11 2016: (Start)
A053117 is a reflected, aerated and signed version of this entry. This entry belongs to a family discussed in A097610 with parameters h1 = -2 and h2 = -y.
Shifted o.g.f.: G(x,t) = x / (1 - 2 x - t x^2).
The compositional inverse of G(x,t) is Ginv(x,t) = -[(1 + 2x) - sqrt[(1+2x)^2 + 4t x^2]] / (2tx) = x - 2 x^2 + (4-t) x^3 - (8-6t) x^4 + ..., a shifted o.g.f. for A091894 (mod signs with A091894(0,0) = 0).
(End)

A090590 (1,1) entry of powers of the orthogonal design shown below.

Original entry on oeis.org

1, -6, -20, 8, 176, 288, -832, -3968, -1280, 29184, 68608, -96256, -741376, -712704, 4505600, 14712832, -6619136, -130940928, -208928768, 629669888, 2930769920, 824180736, -21797797888, -50189041664, 74004299776, 549520932864
Offset: 1

Views

Author

Simone Severini, Dec 04 2003

Keywords

Comments

+1 +1 +1 +1 +1 +1 +1 +1
-1 +1 +1 -1 +1 -1 -1 +1
-1 -1 +1 +1 +1 +1 -1 -1
-1 +1 -1 +1 +1 -1 +1 -1
-1 -1 -1 -1 +1 +1 +1 +1
-1 +1 -1 +1 -1 +1 -1 +1
-1 +1 +1 -1 -1 +1 +1 -1
-1 -1 +1 +1 -1 -1 +1 +1
Also real part of (1 +- i*sqrt(7))^n. - Bruno Berselli, Jun 24-25 2011

Crossrefs

Programs

  • Magma
    m:=27; R:=PowerSeriesRing(Integers(), m); Coefficients(R!(x*(1-8*x)/(1-2*x+8*x^2))); // Bruno Berselli, Jun 24-25 2011
    
  • Maple
    a := proc(n) option remember: if(n=1)then return 1:elif(n=2)then return -6:fi: return 2*a(n-1)-8*a(n-2): end: seq(a(n),n=1..26); # Nathaniel Johnston, Jun 25 2011
  • Mathematica
    LinearRecurrence[{2,-8},{1,-6},30] (* Harvey P. Dale, Mar 30 2019 *)
  • Maxima
    makelist(expand(((1+sqrt(-1)*sqrt(7))^n+(1-sqrt(-1)*sqrt(7))^n)/2),n,1,26); /* Bruno Berselli, Jun 24-25 2011 */
    
  • PARI
    a=vector(26); a[1]=1; a[2]=-6; for(i=3, #a, a[i]=2*a[i-1]-8*a[i-2]); a \\ Bruno Berselli, Jun 24-25 2011

Formula

G.f.: x*(1-8*x)/(1-2*x+8*x^2). - T. D. Noe, Dec 11 2006
From Bruno Berselli, Jun 24-25 2011: (Start)
a(n) = (1/2)*((1+i*sqrt(7))^n + (1-i*sqrt(7))^n), where i=sqrt(-1).
a(n) = cos(n*arctan(sqrt(7)))*sqrt(8)^n.
a(n) = 2*a(n-1) - 8*a(n-2) (n > 2). (End)

Extensions

Corrected by T. D. Noe, Dec 11 2006
More terms from Bruno Berselli, Jun 24 2011

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

Original entry on oeis.org

1, 1, 1, 2, 2, 0, 3, 5, 1, -1, 5, 10, 4, -2, -1, 8, 20, 12, -4, -4, 0, 13, 38, 31, -4, -13, -2, 1, 21, 71, 73, 3, -33, -11, 3, 1, 34, 130, 162, 34, -74, -42, 6, 6, 0, 55, 235, 344, 128, -146, -130, 0, 24, 3, -1
Offset: 0

Views

Author

Philippe Deléham, Dec 06 2011

Keywords

Comments

Row-reversed variant of A123585. Row sums: 2^n.

Examples

			Triangle begins:
1
1, 1
2, 2, 0
3, 5, 1, -1
5, 10, 4, -2, -1
8, 20, 12, -4, -4, 0
13, 38, 31, -4, -13, -2, 1
21, 71, 73, 3, -33, -11, 3, 1
34, 130, 162, 34, -74, -42, 6, 6, 0
55, 235, 344, 128, -146, -130, 0, 24, 3, -1
		

Crossrefs

Cf. Columns: A000045, A001629, A129707.
Diagonals: A010892, A099254, Antidiagonal sums: A158943.

Formula

G.f.: 1/(1-(1+y)*x+(y+1)*(y-1)*x^2).
T(n,0) = A000045(n+1).
T(n,k) = T(n-1,k) + T(n-1,k-1) + T(n-2,k) - T(n-2,k-2) with T(0,0)= 1 and T(n,k)= 0 if n
Sum_{k, 0<=k<=n} T(n,k)*x^k = (-1)^n*A090591(n), (-1)^n*A106852(n), A000007(n), A000045(n+1), A000079(n), A057083(n), A190966(n+1) for n = -3, -2, -1, 0, 1, 2, 3 respectively.
Sum_{k, 0<=k<=n} T(n,k)*x^(n-k) = A010892(n), A000079(n), A030195(n+1), A180222(n+2) for x = 0, 1, 2, 3 respectively.
Showing 1-4 of 4 results.