A045502 -id:A045502 - cp's OEIS frontend

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

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

A160682 The list of the A values in the common solutions to 13k+1 = A^2 and 17k+1 = B^2.

Original entry on oeis.org

1, 14, 209, 3121, 46606, 695969, 10392929, 155197966, 2317576561, 34608450449, 516809180174, 7717529252161, 115246129602241, 1720974414781454, 25699370092119569, 383769576967012081, 5730844284413061646, 85578894689228912609, 1277952576054020627489

Offset: 1

Paul Weisenhorn, May 23 2009

This summarizes the case C=13 of common solutions to C*k+1=A^2, (C+4)*k+1=B^2.
The 2 equations are equivalent to the Pell equation x^2-C*(C+4)*y^2=1,
with x=(C*(C+4)*k+C+2)/2; y=A*B/2 and with smallest values x(1) = (C+2)/2, y(1)=1/2.
Generic recurrences are:
A(j+2)=(C+2)*A(j+1)-A(j) with A(1)=1; A(2)=C+1.
B(j+2)=(C+2)*B(j+1)-B(j) with B(1)=1; B(2)=C+3.
k(j+3)=(C+1)*(C+3)*( k(j+2)-k(j+1) )+k(j) with k(1)=0; k(2)=C+2; k(3)=(C+1)*(C+2)*(C+3).
x(j+2)=(C^2+4*C+2)*x(j+1)-x(j) with x(1)=(C+2)/2; x(2)=(C^2+4*C+1)*(C+2)/2;
Binet-type of solutions of these 2nd order recurrences are:
R=C^2+4*C; S=C*sqrt(R); T=(C+2); U=sqrt(R); V=(C+4)*sqrt(R);
A(j)=((R+S)*(T+U)^(j-1)+(R-S)*(T-U)^(j-1))/(R*2^j);
B(j)=((R+V)*(T+U)^(j-1)+(R-V)*(T-U)^(j-1))/(R*2^j);
x(j)+sqrt(R)*y(j)=((T+U)*(C^2*4*C+2+(C+2)*sqrt(R))^(j-1))/2^j;
k(j)=(((T+U)*(R+2+T*U)^(j-1)+(T-U)*(R+2-T*U)^(j-1))/2^j-T)/R. [Paul Weisenhorn, May 24 2009]
.C -A----- -B----- -k-----
01 A001519 A002878 A058038
02 A001653 A002315 A045899/2
03 A004253 A030221 A160695
04 A001653 A002315 A078522/4
05 A049685 A033890 A161582
06 A070997 A057080 A159683/2
07 A070998 A057081 A161585
08 A072256 A054320 A045502/4
09 A078922 A097783 A161586
10 A077417 A077416 A159681/2
11 A085260 A126816 A161588
12 A001570 A028230 A059989/4
13 A160682 A161591 A161584
14 A157456 A159678 A159679/2
15 A161595 A161599 A161583
16 A007805 A049629 A157459/4
For n>=2, a(n) equals the permanent of the (2n-2)X(2n-2) tridiagonal matrix with sqrt(13)'s along the main diagonal, and 1's along the superdiagonal and the subdiagonal. [John M. Campbell, Jul 08 2011]
Positive values of x (or y) satisfying x^2 - 15xy + y^2 + 13 = 0. - Colin Barker, Feb 11 2014

Vincenzo Librandi, Table of n, a(n) for n = 1..200
J.-C. Novelli, J.-Y. Thibon, Hopf Algebras of m-permutations,(m+1)-ary trees, and m-parking functions, arXiv preprint arXiv:1403.5962 [math.CO], 2014.
Index entries for the Pell equation
Index entries for linear recurrences with constant coefficients, signature (15,-1).

Cf. similar sequences listed in A238379.

I:=[1,14]; [n le 2 select I[n] else 15*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Feb 12 2014

LinearRecurrence[{15,-1},{1,14},20] (* Harvey P. Dale, Oct 08 2012 *)
CoefficientList[Series[(1 - x)/(1 - 15 x + x^2), {x, 0, 40}], x] (* Vincenzo Librandi, Feb 12 2014 *)

a(n) = round((2^(-1-n)*((15-sqrt(221))^n*(13+sqrt(221))+(-13+sqrt(221))*(15+sqrt(221))^n))/sqrt(221)) \\ Colin Barker, Jul 25 2016

a(n) = 15*a(n-1)-a(n-2).
G.f.: (1-x)*x/(1-15*x+x^2).
a(n) = (2^(-1-n)*((15-sqrt(221))^n*(13+sqrt(221))+(-13+sqrt(221))*(15+sqrt(221))^n))/sqrt(221). - Colin Barker, Jul 25 2016

Edited, extended by R. J. Mathar, Sep 02 2009

First formula corrected by Harvey P. Dale, Oct 08 2012

A245031 Numbers m such that 3m+1 and 8m+1 are both squares.

Original entry on oeis.org

0, 1, 21, 120, 2080, 11781, 203841, 1154440, 19974360, 113123361, 1957283461, 11084934960, 191793804840, 1086210502741, 18793835590881, 106437544333680, 1841604094101520, 10429793134197921, 180458407386358101, 1022013289607062600

Offset: 1

Bruno Berselli, Jul 15 2014

Naturally, all terms are triangular numbers.
Numbers m such that k*m+1 and 8*m+1 are both squares:
k=1: A006454;
k=3: this sequence;
k=4: A029549;
k=5: 0, 3, 231, 4560, 333336, 6575751, ...
k=6: A200999;
k=7: A157879.
Numbers m such that 3*m+1 and k*m+1 are both squares:
k=1: A045899;
k=2: A045502;
k=4: A059989;
k=5: A159683;
k=6: 8*A029546;
k=7: A160695;
k=8: this sequence.

Bruno Berselli, Table of n, a(n) for n = 1..500
Index entries for linear recurrences with constant coefficients, signature (1,98,-98,-1,1).

Cf. A000217.

Cf. A006454, A029546, A029549, A045502, A045899, A059989, A157879, A159683, A160695, A200999.

I:=[0,1,21,120,2080]; [n le 5 select I[n] else Self(n-1)+98*Self(n-2)-98*Self(n-3)-Self(n-4)+Self(n-5): n in [1..20]];

LinearRecurrence[{1, 98, -98, -1, 1}, {0, 1, 21, 120, 2080}, 20] (* or *) CoefficientList[Series[x (1 + 20 x + x^2)/((1 - x) (1 - 10 x + x^2) (1 + 10 x + x^2)), {x, 0, 20}], x]

a[1]:0$ a[2]:1$ a[3]:21$ a[4]:120$ a[5]:2080$ a[n]:=a[n-1]+98*a[n-2]-98*a[n-3]-a[n-4]+a[n-5]$ makelist(a[n], n, 1, 20);

a=vector(20); a[1]=0; a[2]=1; a[3]=21; a[4]=120; a[5]=2080; for(i=6, #a, a[i]=a[i-1]+98*a[i-2]-98*a[i-3]-a[i-4]+a[i-5]); a

G.f.: x^2*(1 + 20*x + x^2)/((1 - x)*(1 - 10*x + x^2)*(1 + 10*x + x^2)).
a(n) = a(n-1) + 98*a(n-2) - 98*a(n-3) - a(n-4) + a(n-5).
G.f. of the quadrisections:
a(4k+1):   40*x*(52 + 3*x)/((1 - x)*(1 - 9602*x + x^2));
a(4k+2): (1 + 2178*x + 21*x^2)/((1 - x)*(1 - 9602*x + x^2));
a(4k+3): (21 + 2178*x + x^2)/((1 - x)*(1 - 9602*x + x^2));
a(4k+4): 40*(3 + 52*x)/((1 - x)*(1 - 9602*x + x^2)).

Changed offset from 0 to 1 and adapted formulas by Bruno Berselli, Mar 03 2016

A278620 Expansion of x/(1 - 99x + 99x^2 - x^3).

Original entry on oeis.org

0, 1, 99, 9702, 950698, 93158703, 9128602197, 894509856604, 87652837344996, 8589083549953005, 841642535058049495, 82472379352138897506, 8081451533974553906094, 791899777950154143899707, 77598096787581131548265193, 7603821585405000737586089208, 745096917272902491151888477192

Offset: 0

Bruno Berselli, Nov 24 2016

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

First differences: A173205.

Cf. A045502, A108741, A123479, A200993, A200994.

P:=proc(q) local a,b,c,n; a:=0; b:=1; print(a); print(b);for n from 1 to q do
c:=98*b-a+1; a:=b; b:=c; print(b); od; end: P(100); # Paolo P. Lava, Nov 30 2016

CoefficientList[x/(1 - 99 x + 99 x^2 - x^3) + O[x]^20, x]
LinearRecurrence[{99,-99,1},{0,1,99},20] (* Harvey P. Dale, Aug 22 2020 *)

makelist(coeff(taylor(x/((1-x)*(1-98*x+x^2)), x, 0, n), x, n), n, 0, 20);

concat(0, Vec(1/(1-99*x+99*x^2-x^3) + O(x^20)))

gf = x/((1-x)*(1-98*x+x^2)); print(taylor(gf, x, 0, 20).list())

O.g.f.: x/((1 - x)*(1 - 98*x + x^2)).
E.g.f.: ((5-2*sqrt(6))*exp((5-2*sqrt(6))^2*x) + (5+2*sqrt(6))*exp((5+2*sqrt(6))^2*x) - 10*exp(x))/960.
a(n) = 99*a(n-1) - 99*a(n-2) + a(n-3) for n>2.
a(n) = 98*a(n-1) - a(n-2) + 1 for n>1.
a(n) = a(-n-1) = ((5+2*sqrt(6))^(2*n+1) + (5-2*sqrt(6))^(2*n+1))/960 - 1/96.
a(n) = floor((5+2*sqrt(6))^(2*n+1)/960).
a(n)*a(n-2) = a(n-1)*(a(n-1)-1) for n>1.
Lim_{i -> infinity} a(i)/a(i-1) = (5 + 2*sqrt(6))^2.
From the closed form: a(n) + a(-n) = A108741(n).
a(n) = A200993(n)/10 = A200994(n)/15.
a(n) = A123479(n)/20 for n>0.
a(n) = A045502(n)/40.

cp's OEIS Frontend

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

A160682 The list of the A values in the common solutions to 13*k+1 = A^2 and 17*k+1 = B^2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

Extensions

A245031 Numbers m such that 3*m+1 and 8*m+1 are both squares.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

Maxima

PARI

Formula

Extensions

A278620 Expansion of x/(1 - 99*x + 99*x^2 - x^3).

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

Maxima

PARI

Sage

Formula

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

A160682 The list of the A values in the common solutions to 13k+1 = A^2 and 17k+1 = B^2.

A245031 Numbers m such that 3m+1 and 8m+1 are both squares.

A278620 Expansion of x/(1 - 99x + 99x^2 - x^3).