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

A041041 Denominators of continued fraction convergents to sqrt(26).

Original entry on oeis.org

1, 10, 101, 1020, 10301, 104030, 1050601, 10610040, 107151001, 1082120050, 10928351501, 110365635060, 1114584702101, 11256212656070, 113676711262801, 1148023325284080, 11593909964103601, 117087122966320090, 1182465139627304501, 11941738519239365100
Offset: 0

Views

Author

Keywords

Comments

Generalized Fibonacci sequence.
Sqrt(26) = 10/2 + 10/101 + 10/(101*10301) + 10/(10301*1050601) + ... - Gary W. Adamson, Jun 13 2008
For positive n, a(n) equals the permanent of the n X n tridiagonal matrix with 10's along the main diagonal and 1's along the superdiagonal and the subdiagonal. - John M. Campbell, Jul 08 2011
a(n) equals the number of words of length n on alphabet {0, 1, ..., 10} avoiding runs of zeros of odd lengths. - Milan Janjic, Jan 28 2015
From Bruno Berselli, May 03 2018: (Start)
Numbers k for which m*k^2 + (-1)^k is a perfect square:
m = 2: 0, 1, 2, 5, 12, 29, 70, 169, ... (A000129);
m = 3: 0, 4, 56, 780, 10864, 151316, ... (4*A007655);
m = 5: 0, 1, 4, 17, 72, 305, 1292, ... (A001076);
m = 6: 0, 2, 20, 198, 1960, 19402, ... (A001078);
m = 7: 0, 48, 12192, 3096720, ... (2*A175672);
m = 8: 0, 6, 204, 6930, 235416, ... (A082405);
m = 10: 0, 1, 6, 37, 228, 1405, 8658, ... (A005668);
m = 11: 0, 60, 23880, 9504180, ... [°];
m = 12: 0, 2, 28, 390, 5432, 75658, ... (A011944);
m = 13: 0, 5, 180, 6485, 233640, ... (5*A041613);
m = 14: 0, 4, 120, 3596, 107760, ... (A068204);
m = 15: 0, 8, 496, 30744, 1905632, ... [°];
m = 17: 0, 1, 8, 65, 528, 4289, 34840, ... (A041025);
m = 18: 0, 4, 136, 4620, 156944, ... (A202299);
m = 19: 0, 13260, 1532829480, ... [°];
m = 20: 0, 2, 36, 646, 11592, 208010, ... (A207832);
m = 21: 0, 12, 1320, 145188, ... (A174745);
m = 22: 0, 42, 16548, 6519870, ... (A174766);
m = 23: 0, 240, 552480, 1271808720, ... [°];
m = 24: 0, 10, 980, 96030, 9409960, ... (A168520);
m = 26: 0, 1, 10, 101, 1020, 10301, ... (this sequence);
m = 27: 0, 260, 702520, 1898208780, ... [°];
m = 28: 0, 24, 6096, 1548360, ... (A175672);
m = 29: 0, 13, 1820, 254813, 35675640, ... [°];
m = 30: 0, 2, 44, 966, 21208, 465610, ... (2*A077421), etc.
[°] apparently without related sequences in the OEIS.
(End)
From Michael A. Allen, Mar 12 2023: (Start)
Also called the 10-metallonacci sequence; the g.f. 1/(1-k*x-x^2) gives the k-metallonacci sequence.
a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are 10 kinds of squares available. (End)

Crossrefs

Programs

  • Magma
    I:=[1,10]; [n le 2 select I[n] else 10*Self(n-1) + Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 24 2018
  • Maple
    seq(combinat:-fibonacci(n+1, 10), n=0..19); # Peter Luschny, May 04 2018
  • Mathematica
    Denominator[Convergents[Sqrt[26], 30]] (* Vincenzo Librandi, Dec 10 2013 *)
    LinearRecurrence[{10,1}, {1,10}, 30] (* G. C. Greubel, Jan 24 2018 *)
  • PARI
    x='x+O('x^30); Vec(1/(1-10*x-x^2)) \\ G. C. Greubel, Jan 24 2018
    
  • Sage
    [lucas_number1(n,10,-1) for n in range(1, 19)] # Zerinvary Lajos, Apr 26 2009
    

Formula

G.f.: 1/(1 - 10*x - x^2).
a(n) = 10*a(n-1) + a(n-2), n>=1; a(-1):=0, a(0)=1.
a(n) = S(n, 10*i)*(-i)^n where i^2:=-1 and S(n, x)=U(n, x/2) Chebyshev's polynomials of the second kind. See A049310.
a(n) = (ap^(n+1) - am^(n+1))/(ap-am) with ap = 5+sqrt(26), am = -1/ap = 5-sqrt(26).
a(n) = F(n+1, 10), the (n+1)-th Fibonacci polynomial evaluated at x=10. - T. D. Noe, Jan 19 2006
a(n) = Sum_{i=0..floor(n/2)} binomial(n-i,i)*10^(n-2*i). - Sergio Falcon, Sep 24 2007

Extensions

Extended by T. D. Noe, May 23 2011

A207832 Numbers x such that 20*x^2 + 1 is a perfect square.

Original entry on oeis.org

0, 2, 36, 646, 11592, 208010, 3732588, 66978574, 1201881744, 21566892818, 387002188980, 6944472508822, 124613502969816, 2236098580947866, 40125160954091772, 720016798592704030
Offset: 0

Views

Author

Gary Detlefs, Feb 20 2012

Keywords

Comments

Denote as {a,b,c,d} the second-order linear recurrence a(n) = c*a(n-1) + d*a(n-2) with initial terms a, b. The following sequences and recurrence formulas are related to integer solutions of k*x^2 + 1 = y^2.
.
k x y
- ----------------------- -----------------------
2 A001542 {0,2,6,-1} A001541 {1,3,6,-1}
3 A001353 {0,1,4,-1} A001075 {1,2,4,-1}
5 A060645 {0,4,18,-1} A023039 {1,9,18,-1}
6 A001078 {0,2,10,-1} A001079 {1,5,10,-1}
7 A001080 {0,3,16,-1} A001081 {1,8,16,-1}
8 A001109 {0,1,6,-1} A001541 {1,3,6,-1}
10 A084070 {0,1,38,-1} A078986 {1,19,38,-1}
11 A001084 {0,3,20,-1} A001085 {1,10,20,-1}
12 A011944 {0,2,14,-1} A011943 {1,7,14,-1}
13 A075871 {0,180,1298,-1} A114047 {1,649,1298,-1}
14 A068204 {0,4,30,-1} A069203 {1,15,30,-1}
15 A001090 {0,1,8,-1} A001091 {1,4,8,-1}
17 A121740 {0,8,66,-1} A099370 {1,33,66,-1}
18 A202299 {0,4,34,-1} A056771 {1,17,34,-1}
19 A174765 {0,39,340,-1} A114048 {1,179,340,-1}
20 a(n) {0,2,18,-1} A023039 {1,9,18,-1}
21 A174745 {0,12,110,-1} A114049 {1,55,110,-1}
22 A174766 {0,42,394,-1} A114050 {1,197,394,-1}
23 A174767 {0,5,48,-1} A114051 {1,24,48,-1}
24 A004189 {0,1,10,-1} A001079 {1,5,10,-1}
26 A174768 {0,10,102,-1} A099397 {1,51,102,-1}
The sequence of the c parameter is listed in A180495.

Crossrefs

Programs

  • Magma
    m:=16; R:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!(2*x/(1-18*x+x^2))); // Bruno Berselli, Jun 19 2019
    
  • Maple
    readlib(issqr):for x from 1 to 720016798592704030 do if issqr(20*x^2+1) then print(x) fi od;
  • Mathematica
    LinearRecurrence[{18, -1}, {0, 2}, 16] (* Bruno Berselli, Feb 21 2012 *)
    Table[2 ChebyshevU[-1 + n, 9], {n, 0, 16}]  (* Herbert Kociemba, Jun 05 2022 *)
  • Maxima
    makelist(expand(((2+sqrt(5))^(2*n)-(2-sqrt(5))^(2*n))/(4*sqrt(5))), n, 0, 15); /* Bruno Berselli, Jun 19 2019 */

Formula

a(n) = 18*a(n-1) - a(n-2).
From Bruno Berselli, Feb 21 2012: (Start)
G.f.: 2*x/(1-18*x+x^2).
a(n) = -a(-n) = 2*A049660(n) = ((2 + sqrt(5))^(2*n)-(2 - sqrt(5))^(2*n))/(4*sqrt(5)). (End)
a(n) = Fibonacci(6*n)/4. - Bruno Berselli, Jun 19 2019
For n>=1, a(n) = A079962(6n-3). - Christopher Hohl, Aug 22 2021

A114050 x-values in the solution to x^2 - 22*y^2 = 1.

Original entry on oeis.org

1, 197, 77617, 30580901, 12048797377, 4747195585637, 1870383011943601, 736926159510193157, 290347036464004160257, 114395995440658128948101, 45071731856582838801391537, 17758147955498197829619317477
Offset: 1

Views

Author

Cino Hilliard, Feb 01 2006

Keywords

Comments

A Pellian equation (Pell's equation). - Benoit Cloitre, Feb 03 2006
The corresponding values of y are in A174766. - Vincenzo Librandi, Dec 21 2011

Programs

  • Magma
    I:=[1,197]; [n le 2 select I[n] else 394*Self(n-1)-Self(n-2): n in [1..40]]; // Vincenzo Librandi, Dec 21 2011
  • Mathematica
    LinearRecurrence[{394,-1},{1,197},20] (* Harvey P. Dale, Nov 03 2011 *)
  • PARI
    g(n,k=22) = for(y=0,n,x=k*y^2+1;if(issquare(x),print1(sqrtint(x), ", ")))
    
  • PARI
    a=vector(12); a[1]=1; a[2]=197; for(i=3, #a, a[i]=394*a[i-1]-a[i-2]); a \\ Benoit Cloitre, Feb 03 2006
    

Formula

a(n) = 394*a(n-1) - a(n-2); a(1)=1, a(2)=197. - Benoit Cloitre, Feb 03 2006
G.f.: x*(1-197*x)/(1 - 394*x + x^2). - Philippe Deléham, Nov 18 2008

Extensions

More terms from Benoit Cloitre, Feb 03 2006
Offset changed from 0 to 1, g.f. and formula adapted by Vincenzo Librandi, Dec 21 2011
Showing 1-3 of 3 results.