A153111 -id:A153111 - cp's OEIS frontend

A238379 Expansion of (1 - x)/(1 - 36*x + x^2).

1, 35, 1259, 45289, 1629145, 58603931, 2108112371, 75833441425, 2727895778929, 98128414600019, 3529895029821755, 126978092658983161, 4567681440693572041, 164309553772309610315, 5910576254362452399299, 212616435603275976764449

Offset: 0

Bruno Berselli, Feb 25 2014

First bisection of A041611.

Bruno Berselli, Table of n, a(n) for n = 0..100
Tanya Khovanova, Recursive Sequences
Index entries for linear recurrences with constant coefficients, signature (36,-1).

Cf. similar sequences with g.f. (1-x)/(1-k*x+x^2): A122367 (k=3), A079935 (k=4), A004253 (k=5), A001653 (k=6), A049685 (k=7), A070997 (k=8), A070998 (k=9), A138288 (k=10), A078922 (k=11), A077417 (k=12), A085260 (k=13), A001570 (k=14), A160682 (k=15), A157456 (k=16), A161595 (k=17). From 18 to 38, even k only, except k=27 and k=31: A007805 (k=18), A075839 (k=20), A157014 (k=22), A159664 (k=24), A153111 (k=26), A097835 (k=27), A159668 (k=28), A157877 (k=30), A111216 (k=31), A159674 (k=32), A077420 (k=34), this sequence (k=36), A097315 (k=38).

[n le 2 select 35^(n-1) else 36*Self(n-1)-Self(n-2): n in [1..20]];

R:=PowerSeriesRing(Integers(), 20); Coefficients(R!( (1 - x)/(1 - 36*x + x^2))); // Marius A. Burtea, Jan 14 2020

CoefficientList[Series[(1 - x)/(1 - 36 x + x^2), {x, 0, 20}], x] (* or *) LinearRecurrence[{36, -1}, {1, 35}, 20]

a(n)=([0,1; -1,36]^n*[1;35])[1,1] \\ Charles R Greathouse IV, May 10 2016

m = 20; L. = PowerSeriesRing(ZZ, m); f = (1-x)/(1-36*x+x^2)
print(f.coefficients())

G.f.: (1 - x)/(1 - 36*x + x^2).
a(n) = a(-n-1) = 36*a(n-1) - a(n-2).
a(n) = ((19-sqrt(323))/38)*(1+(18+sqrt(323))^(2*n+1))/(18+sqrt(323))^n.
a(n+1) - a(n) = 34*A144128(n+1).
323*a(n+1)^2 - ((a(n+2)-a(n))/2)^2 = 34.
Sum_{n>0} 1/(a(n) - 1/a(n)) = 1/34.
See also Tanya Khovanova in Links field:
a(n) = 35*a(n-1) + 34*Sum_{i=0..n-2} a(i).
a(n+2)*a(n) - a(n+1)^2 = 36-2 = 34 = 34*1,
a(n+3)*a(n) - a(n+1)*a(n+2) = 36*(36-2) = 1224 = 34*36.
Generalizing:
a(n+4)*a(n) - a(n+1)*a(n+3) = 44030 = 34*1295,
a(n+5)*a(n) - a(n+1)*a(n+4) = 1583856 = 34*46584,
a(n+6)*a(n) - a(n+1)*a(n+5) = 56974786 = 34*1675729, etc.,
where 1, 36, 1295, 46584, 1675729, ... is the sequence A144128, which is the second bisection of A041611.
a(n)^2 - 36*a(n)*a(n+1) + a(n+1)^2 + 34 = 0 (see comments by Colin Barker in similar sequences).

A157014 Expansion of x(1-x)/(1 - 22x + 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

Paul Weisenhorn, Feb 21 2009

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
  1     A001653       A002315(n-1)   A078522
  2     A072256       A054320(n-1)   A045502(n-1)
  3     A001570       A028230        A059989(n-1)
  4     A007805       A049629(n-1)   A157459
  5  -> A157014 <-    A133283        A157460
  6     A153111       A157461        A157874
  7     A157877       A157878        A157879
  8     A077420(n-1)  A046176        A157880
  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)

Vincenzo Librandi, Table of n, a(n) for n = 1..200
Index entries for linear recurrences with constant coefficients, signature (22,-1).

Cf. similar sequences listed in A238379.

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

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

seq( simplify(ChebyshevU(n-1,11) - ChebyshevU(n-2,11)), n=1..20); # G. C. Greubel, Jan 14 2020

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

Vec((1-x)/(1-22*x+x^2)+O(x^20)) \\ Charles R Greathouse IV, Sep 23 2012

[chebyshev_U(n-1,11) - chebyshev_U(n-2,11) for n in (1..20)] # G. C. Greubel, Jan 14 2020

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

Edited by Alois P. Heinz, Sep 09 2011

A157461 Expansion of x(x+1) / (x^2-26x+1).

Original entry on oeis.org

1, 27, 701, 18199, 472473, 12266099, 318446101, 8267332527, 214632199601, 5572169857099, 144661784084973, 3755634216352199, 97501827841072201, 2531291889651525027, 65716087303098578501, 1706086977990911515999, 44292545340460600837473

Offset: 1

Paul Weisenhorn, Mar 01 2009

This sequence is part of a solution of a more general problem involving two 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, for details see comment in A157014.
A157461 is the b(n) sequence for A=6.
Numbers k such that 42*k^2 + 7 is a square. - Klaus Purath, Jun 12 2021

Colin Barker, Table of n, a(n) for n = 1..700
Index entries for linear recurrences with constant coefficients, signature (26,-1).

6*A157874(n)+1 = A153111(n)^2.

7*A157874(n)+1 = A157461(n)^2.

Vec(x*(x+1)/(x^2-26*x+1)+O(x^20)) \\ Charles R Greathouse IV, Sep 26 2012

a(n) = round((13+2*sqrt(42))^(-n)*(-6-sqrt(42)+(-6+sqrt(42))*(13+2*sqrt(42))^(2*n))/12) \\ Colin Barker, Jul 25 2016

G.f.: x*(x+1) / (x^2-26*x+1).
a(1) = 1, a(2) = 27, a(n) = 26*a(n-1)-a(n-2) for n>2.
a(n) = (13+2*sqrt(42))^(-n)*(-6-sqrt(42)+(-6+sqrt(42))*(13+2*sqrt(42))^(2*n))/12. - Colin Barker, Jul 25 2016
a(n+1) = (a(n)^2 - 28)/a(n-1), n > 1. - Klaus Purath, Jun 12 2021

Edited by Alois P. Heinz, Sep 09 2011

A157874 Expansion of 104x^2 / (-x^3+675x^2-675*x+1).

Original entry on oeis.org

0, 104, 70200, 47314800, 31890105104, 21493883525400, 14486845606014600, 9764112444570315104, 6580997300794786365600, 4435582416623241440099400, 2989575967806763935840630104, 2014969766719342269515144590800, 1358086633192868882889271613569200

Offset: 1

Paul Weisenhorn, Mar 08 2009

This sequence is part of a solution of a more general problem involving two 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, for details see comment in A157014.
A157874 is the c(n) sequence for A=6.

Colin Barker, Table of n, a(n) for n = 1..350
Index entries for linear recurrences with constant coefficients, signature (675,-675,1).

6*A157874(n)+1 = A153111(n)^2.

7*A157874(n)+1 = A157461(n)^2.

Rest[CoefficientList[Series[104x^2/(-x^3+675x^2-675x+1),{x,0,20}],x]] (* or *) LinearRecurrence[{675,-675,1},{0,104,70200},20] (* Harvey P. Dale, Oct 04 2015 *)

concat(0, Vec(104*x^2/(-x^3+675*x^2-675*x+1) + O(x^20))) \\ Charles R Greathouse IV, Sep 26 2012

a(n) = -round((337+52*sqrt(42))^(-n)*(-1+(337+52*sqrt(42))^n)*(13+2*sqrt(42)+(-13+2*sqrt(42))*(337+52*sqrt(42))^n))/168 \\ Colin Barker, Jul 25 2016

G.f.: 104*x^2 / (-x^3+675*x^2-675*x+1).
c(1) = 0, c(2) = 104, c(3) = 675*c(2), c(n) = 675 * (c(n-1)-c(n-2)) + c(n-3) for n>3.
a(n) = -((337+52*sqrt(42))^(-n)*(-1+(337+52*sqrt(42))^n)*(13+2*sqrt(42)+(-13+2*sqrt(42))*(337+52*sqrt(42))^n))/168. - Colin Barker, Jul 25 2016

Edited by Alois P. Heinz, Sep 09 2011

cp's OEIS Frontend

A238379 Expansion of (1 - x)/(1 - 36*x + x^2).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Magma

Mathematica

PARI

Sage

Formula

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

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A157461 Expansion of x*(x+1) / (x^2-26*x+1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

PARI

Formula

Extensions

A157874 Expansion of 104*x^2 / (-x^3+675*x^2-675*x+1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

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

A157461 Expansion of x(x+1) / (x^2-26x+1).

A157874 Expansion of 104x^2 / (-x^3+675x^2-675*x+1).