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

A157014 Expansion of x*(1-x)/(1 - 22*x + x^2).

Original entry on oeis.org

1, 21, 461, 10121, 222201, 4878301, 107100421, 2351330961, 51622180721, 1133336644901, 24881784007101, 546265911511321, 11992968269241961, 263299036011811821, 5780585823990618101, 126909589091781786401, 2786230374195208682721, 61170158643202809233461
Offset: 1

Views

Author

Paul Weisenhorn, Feb 21 2009

Keywords

Comments

This sequence is part of a solution of a general problem involving 2 equations, three sequences a(n), b(n), c(n) and a constant A:
A * c(n)+1 = a(n)^2,
(A+1) * c(n)+1 = b(n)^2, where solutions are given by the recurrences:
a(1) = 1, a(2) = 4*A+1, a(n) = (4*A+2)*a(n-1)-a(n-2) for n>2, resulting in a(n) terms 1, 4*A+1, 16*A^2+12*A+1, 64*A^3+80*A^2+24*A+1, ...;
b(1) = 1, b(2) = 4*A+3, b(n) = (4*A+2)*b(n-1)-b(n-2) for n>2, resulting in b(n) terms 1, 4*A+3, 16*A^2+20*A+5, 64*A^3+112*A^2+56*A+7, ...;
c(1) = 0, c(2) = 16*A+8, c(3) = (16*A^2+16*A+3)*c(2), c(n) = (16*A^2+16*A+3) * (c(n-1)-c(n-2)) + c(n-3) for n>3, resulting in c(n) terms 0, 16*A+8, 256*A^3+384*A^2+176*A+24, 4096*A^5 + 10240*A^4 + 9472*A^3 + 3968*A^2 + 736*A + 48, ... .
A157014 is the a(n) sequence for A=5.
For other A values the a(n), b(n) and c(n) sequences are in the OEIS:
A a-sequence b-sequence c-sequence
2 A072256 A054320(n-1) A045502(n-1)
9 A097315(n-1) A097314(n-1) A157881
Positive values of x (or y) satisfying x^2 - 22xy + y^2 + 20 = 0. - Colin Barker, Feb 19 2014
From Klaus Purath, Apr 22 2025: (Start)
Nonnegative solutions to the Diophantine equation 5*b(n)^2 - 6*a(n)^2 = -1. The corresponding b(n) are A133283(n). Note that (b(n+1)^2 - b(n)*b(n+2))/4 = 6 and (a(n)*a(n+2) - a(n+1)^2)/4 = 5.
(a(n) + b(n))/2 = (b(n+1) - a(n+1))/2 = A077421(n-1) = Lucas U(22,1). Also b(n)*a(n+1) - b(n+1)*a(n) = -2.
a(n)=(t(i+2*n-1) + t(i))/(t(i+n) + t(i+n-1)) as long as t(i+n) + t(i+n-1) != 0 for any integer i and n >= 1 where (t) is a sequence satisfying t(i+3) = 21*t(i+2) - 21*t(i+1) + t(i) or t(i+2) = 22*t(i+1) - t(i) without regard to initial values and including this sequence itself. (End)

Crossrefs

Cf. similar sequences listed in A238379.

Programs

  • GAP
    a:=[1,21];; for n in [3..20] do a[n]:=22*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Jan 14 2020
  • Magma
    I:=[1,21]; [n le 2 select I[n] else 22*Self(n-1)-Self(n-2): n in [1..20]]; // Vincenzo Librandi, Feb 21 2014
    
  • Maple
    seq( simplify(ChebyshevU(n-1,11) - ChebyshevU(n-2,11)), n=1..20); # G. C. Greubel, Jan 14 2020
  • Mathematica
    CoefficientList[Series[(1-x)/(1-22x+x^2), {x,0,20}], x] (* Vincenzo Librandi, Feb 21 2014 *)
    a[c_, n_] := Module[{},
       p := Length[ContinuedFraction[ Sqrt[ c]][[2]]];
       d := Denominator[Convergents[Sqrt[c], n p]];
       t := Table[d[[1 + i]], {i, 0, Length[d] - 1, p}];
       Return[t];
    ] (* Complement of A041049 *)
    a[30, 20] (* Gerry Martens, Jun 07 2015 *)
    Table[ChebyshevU[n-1, 11] - ChebyshevU[n-2, 11], {n,20}] (* G. C. Greubel, Jan 14 2020 *)
  • PARI
    Vec((1-x)/(1-22*x+x^2)+O(x^20)) \\ Charles R Greathouse IV, Sep 23 2012
    
  • Sage
    [chebyshev_U(n-1,11) - chebyshev_U(n-2,11) for n in (1..20)] # G. C. Greubel, Jan 14 2020
    

Formula

G.f.: x*(1-x)/(1-22*x+x^2).
a(1) = 1, a(2) = 21, a(n) = 22*a(n-1) - a(n-2) for n>2.
5*A157460(n)+1 = a(n)^2 for n>=1.
6*A157460(n)+1 = A133283(n)^2 for n>=1.
a(n) = (6+sqrt(30)-(-6+sqrt(30))*(11+2*sqrt(30))^(2*n))/(12*(11+2*sqrt(30))^n). - Gerry Martens, Jun 07 2015
a(n) = ChebyshevU(n-1, 11) - ChebyshevU(n-2, 11). - G. C. Greubel, Jan 14 2020

Extensions

Edited by Alois P. Heinz, Sep 09 2011

A041048 Numerators of continued fraction convergents to sqrt(30).

Original entry on oeis.org

5, 11, 115, 241, 2525, 5291, 55435, 116161, 1217045, 2550251, 26719555, 55989361, 586613165, 1229215691, 12878770075, 26986755841, 282746328485, 592479412811, 6207540456595, 13007560326001, 136283143716605, 285573847759211, 2992021621308715, 6269617090376641
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

Formula

G.f.: (5+11*x+5*x^2-x^3)/(1-22*x^2+x^4). - Colin Barker, Jan 02 2012
a(n) = 22*a(n-2) - a(n-4). - Wesley Ivan Hurt, Jun 16 2023

Extensions

More terms from Harvey P. Dale, Apr 04 2013
Showing 1-2 of 2 results.