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

A097726 Pell equation solutions (5*a(n))^2 - 26*b(n)^2 = -1 with b(n):=A097727(n), n >= 0.

Original entry on oeis.org

1, 103, 10505, 1071407, 109273009, 11144775511, 1136657829113, 115927953794015, 11823514629160417, 1205882564220568519, 122988198035868828521, 12543590317094399940623, 1279323224145592925115025, 130478425272533383961791927, 13307520054574259571177661529
Offset: 0

Views

Author

Wolfdieter Lang, Aug 31 2004

Keywords

Comments

a(-1) = -1. - Artur Jasinski, Feb 10 2010
5*a(n) gives the x-values in the solution to the Pell equation x^2 - 26*y^2 = -1. - Colin Barker, Aug 24 2013

Examples

			(x,y) = (5,1), (515,101), (52525,10301), ... give the positive integer solutions to x^2 - 26*y^2 = -1.
		

Crossrefs

Cf. A097725 for S(n, 102).
Cf. similar sequences of the type (1/k)*sinh((2*n+1)*arcsinh(k)) listed in A097775.

Programs

  • Mathematica
    Table[(1/5) Round[N[Sinh[(2 n - 1) ArcSinh[5]], 100]], {n, 1, 50}] (* Artur Jasinski, Feb 10 2010 *)
    CoefficientList[Series[(1 + x)/(1 - 102 x + x^2), {x, 0, 40}], x] (* Vincenzo Librandi, Apr 13 2014 *)
    LinearRecurrence[{102,-1},{1,103},20] (* Harvey P. Dale, Aug 20 2017 *)
  • PARI
    x='x+O('x^99); Vec((1+x)/(1-102*x+x^2)) \\ Altug Alkan, Apr 05 2018

Formula

G.f.: (1 + x)/(1 - 102*x + x^2).
a(n) = S(n, 2*51) + S(n-1, 2*51) = S(2*n, 2*sqrt(26)), with Chebyshev polynomials of the 2nd kind. See A049310 for the triangle of S(n, x)= U(n, x/2) coefficients. S(-1, x) := 0 =: U(-1, x).
a(n) = ((-1)^n)*T(2*n+1, 5*i)/(5*i) with the imaginary unit i and Chebyshev polynomials of the first kind. See the T-triangle A053120.
a(n) = 102*a(n-1) - a(n-2) for n > 1; a(0)=1, a(1)=103. - Philippe Deléham, Nov 18 2008
a(n) = (1/5)*sinh((2*n-1)*arcsinh(5)), n >= 1. - Artur Jasinski, Feb 10 2010

Extensions

More terms from Harvey P. Dale, Aug 20 2017

A195500 Denominators a(n) of Pythagorean approximations b(n)/a(n) to sqrt(2).

Original entry on oeis.org

3, 228, 308, 5289, 543900, 706180, 1244791, 51146940, 76205040, 114835995824, 106293119818725, 222582887719576, 3520995103197240, 17847666535865852, 18611596834765355, 106620725307595884, 269840171418387336, 357849299891217865
Offset: 1

Views

Author

Clark Kimberling, Sep 20 2011

Keywords

Comments

For each positive real number r, there is a sequence (a(n),b(n),c(n)) of primitive Pythagorean triples such that the limit of b(n)/a(n) is r and
|r-b(n+1)/a(n+1)| < |r-b(n)/a(n)|. Peter Shiu showed how to find (a(n),b(n)) from the continued fraction for r, and Peter J. C. Moses incorporated Shiu's method in the Mathematica program shown below.
Examples:
r...........a(n)..........b(n)..........c(n)
sqrt(2).....A195500.......A195501.......A195502
sqrt(3).....A195499.......A195503.......A195531
sqrt(5).....A195532.......A195533.......A195534
sqrt(6).....A195535.......A195536.......A195537
sqrt(8).....A195538.......A195539.......A195540
sqrt(12)....A195680.......A195681.......A195682
e...........A195541.......A195542.......A195543
pi..........A195544.......A195545.......A195546
tau.........A195687.......A195688.......A195689
1...........A046727.......A084159.......A001653
2...........A195614.......A195615.......A007805
3...........A195616.......A195617.......A097315
4...........A195619.......A195620.......A078988
5...........A195622.......A195623.......A097727
1/2.........A195547.......A195548.......A195549
3/2.........A195550.......A195551.......A195552
5/2.........A195553.......A195554.......A195555
1/3.........A195556.......A195557.......A195558
2/3.........A195559.......A195560.......A195561
1/4.........A195562.......A195563.......A195564
5/4.........A195565.......A195566.......A195567
7/4.........A195568.......A195569.......A195570
1/5.........A195571.......A195572.......A195573
2/5.........A195574.......A195575.......A195576
3/5.........A195577.......A195578.......A195579
4/5.........A195580.......A195611.......A195612
sqrt(1/2)...A195625.......A195626.......A195627
sqrt(1/3)...{1}+A195503...{0}+A195499...{1}+A195531
sqrt(2/3)...A195631.......A195632.......A195633
sqrt(3/4)...A195634.......A195635.......A195636

Examples

			For r=sqrt(2), the first five fractions b(n)/a(n) can be read from the following five primitive Pythagorean triples (a(n), b(n), c(n)) = (A195500, A195501, A195502):
(3,4,5); |r - b(1)/a(1)| = 0.08...
(228,325,397); |r - b(2)/a(2)| = 0.011...
(308,435,533); |r - b(3)/a(3)| = 0.0018...
(5289,7480,9161); |r - b(4)/a(4)| = 0.000042...
(543900,769189,942061); |r - b(5)/a(5)| = 0.0000003...
		

Crossrefs

Programs

  • Maple
    Shiu := proc(r,n)
            t := r+sqrt(1+r^2) ;
            cf := numtheory[cfrac](t,n+1) ;
            mn := numtheory[nthconver](cf,n) ;
            (mn-1/mn)/2 ;
    end proc:
    A195500 := proc(n)
            Shiu(sqrt(2),n) ;
            denom(%) ;
    end proc: # R. J. Mathar, Sep 21 2011
  • Mathematica
    r = Sqrt[2]; z = 18;
    p[{f_, n_}] := (#1[[2]]/#1[[
          1]] &)[({2 #1[[1]] #1[[2]], #1[[1]]^2 - #1[[
             2]]^2} &)[({Numerator[#1], Denominator[#1]} &)[
         Array[FromContinuedFraction[
            ContinuedFraction[(#1 + Sqrt[1 + #1^2] &)[f], #1]] &, {n}]]]];
    {a, b} = ({Denominator[#1], Numerator[#1]} &)[
      p[{r, z}]]  (* A195500, A195501 *)
    Sqrt[a^2 + b^2] (* A195502 *)

A085478 Triangle read by rows: T(n, k) = binomial(n + k, 2*k).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 6, 5, 1, 1, 10, 15, 7, 1, 1, 15, 35, 28, 9, 1, 1, 21, 70, 84, 45, 11, 1, 1, 28, 126, 210, 165, 66, 13, 1, 1, 36, 210, 462, 495, 286, 91, 15, 1, 1, 45, 330, 924, 1287, 1001, 455, 120, 17, 1, 1, 55, 495, 1716, 3003, 3003, 1820, 680, 153, 19, 1
Offset: 0

Views

Author

Philippe Deléham, Aug 14 2003

Keywords

Comments

Coefficient array for Morgan-Voyce polynomial b(n,x). A053122 (unsigned) is the coefficient array for B(n,x). Reversal of A054142. - Paul Barry, Jan 19 2004
This triangle is formed from even-numbered rows of triangle A011973 read in reverse order. - Philippe Deléham, Feb 16 2004
T(n,k) is the number of nondecreasing Dyck paths of semilength n+1, having k+1 peaks. T(n,k) is the number of nondecreasing Dyck paths of semilength n+1, having k peaks at height >= 2. T(n,k) is the number of directed column-convex polyominoes of area n+1, having k+1 columns. - Emeric Deutsch, May 31 2004
Riordan array (1/(1-x), x/(1-x)^2). - Paul Barry, May 09 2005
The triangular matrix a(n,k) = (-1)^(n+k)*T(n,k) is the matrix inverse of A039599. - Philippe Deléham, May 26 2005
The n-th row gives absolute values of coefficients of reciprocal of g.f. of bottom-line of n-wave sequence. - Floor van Lamoen (fvlamoen(AT)planet.nl), Sep 24 2006
Unsigned version of A129818. - Philippe Deléham, Oct 25 2007
T(n, k) is also the number of idempotent order-preserving full transformations (of an n-chain) of height k >=1 (height(alpha) = |Im(alpha)|) and of waist n (waist(alpha) = max(Im(alpha))). - Abdullahi Umar, Oct 02 2008
A085478 is jointly generated with A078812 as a triangular array of coefficients of polynomials u(n,x): initially, u(1,x) = v(1,x) = 1; for n>1, u(n,x) = u(n-1,x)+x*v(n-1)x and v(n,x) = u(n-1,x)+(x+1)*v(n-1,x). See the Mathematica section. - Clark Kimberling, Feb 25 2012
Per Kimberling's recursion relations, see A102426. - Tom Copeland, Jan 19 2016
Subtriangle of the triangle given by (0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Mar 26 2012
T(n,k) is also the number of compositions (ordered partitions) of 2*n+1 into 2*k+1 parts which are all odd. Proof: The o.g.f. of column k, x^k/(1-x)^(2*k+1) for k >= 0, is the o.g.f. of the odd-indexed members of the sequence with o.g.f. (x/(1-x^2))^(2*k+1) (bisection, odd part). Thus T(n,k) is obtained from the sum of the multinomial numbers A048996 for the partitions of 2*n+1 into 2*k+1 parts, all of which are odd. E.g., T(3,1) = 3 + 3 from the numbers for the partitions [1,1,5] and [1,3,3], namely 3!/(2!*1!) and 3!/(1!*2!), respectively. The number triangle with the number of these partitions as entries is A152157. - Wolfdieter Lang, Jul 09 2012
The matrix elements of the inverse are T^(-1)(n,k) = (-1)^(n+k)*A039599(n,k). - R. J. Mathar, Mar 12 2013
T(n,k) = A258993(n+1,k) for k = 0..n-1. - Reinhard Zumkeller, Jun 22 2015
The n-th row polynomial in descending powers of x is the n-th Taylor polynomial of the algebraic function F(x)*G(x)^n about 0, where F(x) = (1 + sqrt(1 + 4*x))/(2*sqrt(1 + 4*x)) and G(x) = ((1 + sqrt(1 + 4*x))/2)^2. For example, for n = 4, (1 + sqrt(1 + 4*x))/(2*sqrt(1 + 4*x)) * ((1 + sqrt(1 + 4*x))/2)^8 = (x^4 + 10*x^3 + 15*x^2 + 7*x + 1) + O(x^5). - Peter Bala, Feb 23 2018
Row n also gives the coefficients of the characteristc polynomial of the tridiagonal n X n matrix M_n given in A332602: Phi(n, x) := Det(M_n - x*1_n) = Sum_{k=0..n} T(n, k)*(-x)^k, for n >= 0, with Phi(0, x) := 1. - Wolfdieter Lang, Mar 25 2020
It appears that the largest root of the n-th degree polynomial is equal to the sum of the distinct diagonals of a (2*n+1)-gon including the edge, 1. The largest root of x^3 - 6*x^2 + 5*x - 1 is 5.048917... = the sum of (1 + 1.80193... + 2.24697...). Alternatively, the largest root of the n-th degree polynomial is equal to the square of sigma(2*n+1). Check: 5.048917... is the square of sigma(7), 2.24697.... Given N = 2*n+1, sigma(N) (N odd) can be defined as 1/(2*sin(Pi/(2*N))). Relating to the 9-gon, the largest root of x^4 - 10*x^3 + 15*x^2 - 7*x + 1 is 8.290859..., = the sum of (1 + 1.879385... + 2.532088... + 2.879385...), and is the square of sigma(9), 2.879385... Refer to A231187 for a further clarification of sigma(7). - Gary W. Adamson, Jun 28 2022
For n >=1, the n-th row is given by the coefficients of the minimal polynomial of -4*sin(Pi/(4*n + 2))^2. - Eric W. Weisstein, Jul 12 2023
Denoting this lower triangular array by L, then L * diag(binomial(2*k,k)^2) * transpose(L) is the LDU factorization of A143007, the square array of crystal ball sequences for the A_n X A_n lattices. - Peter Bala, Feb 06 2024
T(n, k) is the number of occurrences of the periodic substring (01)^k in the periodic string (01)^n (see Proposition 4.7 at page 7 in Fang). - Stefano Spezia, Jun 09 2024

Examples

			Triangle begins as:
  1;
  1    1;
  1    3    1;
  1    6    5    1;
  1   10   15    7    1;
  1   15   35   28    9    1;
  1   21   70   84   45   11    1;
  1   28  126  210  165   66   13    1;
  1   36  210  462  495  286   91   15    1;
  1   45  330  924 1287 1001  455  120   17    1;
  1   55  495 1716 3003 3003 1820  680  153   19    1;
...
From _Philippe Deléham_, Mar 26 2012: (Start)
(0, 1, 0, 1, 0, 0, 0, ...) DELTA (1, 0, 1, -1, 0, 0, 0, ...) begins:
  1
  0, 1
  0, 1,  1
  0, 1,  3,   1
  0, 1,  6,   5,   1
  0, 1, 10,  15,   7,   1
  0, 1, 15,  35,  28,   9,  1
  0, 1, 21,  70,  84,  45, 11,  1
  0, 1, 28, 126, 210, 165, 66, 13, 1. (End)
		

Crossrefs

Programs

  • GAP
    Flat(List([0..12], n-> List([0..n], k-> Binomial(n+k, 2*k) ))); # G. C. Greubel, Aug 01 2019
  • Haskell
    a085478 n k = a085478_tabl !! n !! k
    a085478_row n = a085478_tabl !! n
    a085478_tabl = zipWith (zipWith a007318) a051162_tabl a025581_tabl
    -- Reinhard Zumkeller, Jun 22 2015
    
  • Magma
    [Binomial(n+k, 2*k): k in [0..n], n in [0..12]]; // G. C. Greubel, Aug 01 2019
    
  • Maple
    T := (n,k) -> binomial(n+k,2*k): seq(seq(T(n,k), k=0..n), n=0..11);
  • Mathematica
    (* First program *)
    u[1, x_]:= 1; v[1, x_]:= 1; z = 13;
    u[n_, x_]:= u[n-1, x] + x*v[n-1, x];
    v[n_, x_]:= u[n-1, x] + (x+1)*v[n-1, x];
    Table[Expand[u[n, x]], {n, 1, z/2}]
    Table[Expand[v[n, x]], {n, 1, z/2}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%]   (* A085478 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]   (* A078812 *) (*Clark Kimberling, Feb 25 2012 *)
    (* Second program *)
    Table[Binomial[n + k, 2 k], {n, 0, 12}, {k, 0, n}] // Flatten (* G. C. Greubel, Aug 01 2019 *)
    CoefficientList[Table[Fibonacci[2 n + 1, Sqrt[x]], {n, 0, 10}], x] // Flatten (* Eric W. Weisstein, Jul 03 2023 *)
    Join[{{1}}, CoefficientList[Table[MinimalPolynomial[-4 Sin[Pi/(4 n + 2)]^2, x], {n, 20}], x]] (* Eric W. Weisstein, Jul 12 2023 *)
  • PARI
    T(n,k) = binomial(n+k,n-k)
    
  • Sage
    [[binomial(n+k,2*k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Aug 01 2019
    

Formula

T(n, k) = (n+k)!/((n-k)!*(2*k)!).
G.f.: (1-z)/((1-z)^2-tz). - Emeric Deutsch, May 31 2004
Row sums are A001519 (Fibonacci(2n+1)). Diagonal sums are A011782. Binomial transform of A026729 (product of lower triangular matrices). - Paul Barry, Jun 21 2004
T(n, 0) = 1, T(n, k) = 0 if n=0} T(n-1-j, k-1)*(j+1). T(0, 0) = 1, T(0, k) = 0 if k>0; T(n, k) = T(n-1, k-1) + T(n-1, k) + Sum_{j>=0} (-1)^j*T(n-1, k+j)*A000108(j). For the column k, g.f.: Sum_{n>=0} T(n, k)*x^n = (x^k) / (1-x)^(2*k+1). - Philippe Deléham, Feb 15 2004
Sum_{k=0..n} T(n,k)*x^(2*k) = A000012(n), A001519(n+1), A001653(n), A078922(n+1), A007805(n), A097835(n), A097315(n), A097838(n), A078988(n), A097841(n), A097727(n), A097843(n), A097730(n), A098244(n), A097733(n), A098247(n), A097736(n), A098250(n), A097739(n), A098253(n), A097742(n), A098256(n), A097767(n), A098259(n), A097770(n), A098262(n), A097773(n), A098292(n), A097776(n) for x=0,1,2,...,27,28 respectively. - Philippe Deléham, Dec 31 2007
T(2*n,n) = A005809(n). - Philippe Deléham, Sep 17 2009
A183160(n) = Sum_{k=0..n} T(n,k)*T(n,n-k). - Paul D. Hanna, Dec 27 2010
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k). - Philippe Deléham, Feb 06 2012
O.g.f. for column k: x^k/(1-x)^(2*k+1), k >= 0. [See the o.g.f. of the triangle above, and a comment on compositions. - Wolfdieter Lang, Jul 09 2012]
E.g.f.: (2/sqrt(x + 4))*sinh((1/2)*t*sqrt(x + 4))*cosh((1/2)*t*sqrt(x)) = t + (1 + x)*t^3/3! + (1 + 3*x + x^2)*t^5/5! + (1 + 6*x + 5*x^2 + x^3)*t^7/7! + .... Cf. A091042. - Peter Bala, Jul 29 2013
T(n, k) = A065941(n+3*k, 4*k) = A108299(n+3*k, 4*k) = A194005(n+3*k, 4*k). - Johannes W. Meijer, Sep 05 2013
Sum_{k=0..n} (-1)^k*T(n,k)*A000108(k) = A000007(n) for n >= 0. - Werner Schulte, Jul 12 2017
Sum_{k=0..floor(n/2)} T(n-k,k)*A000108(k) = A001006(n) for n >= 0. - Werner Schulte, Jul 12 2017
From Peter Bala, Jun 26 2025: (Start)
The n-th row polynomial b(n, x) = (-1)^n * U(2*n, (i/2)*sqrt(x)), where U(n,x) is the n-th Chebyshev polynomial of the second kind.
b(n, x) = (-1)^n * Dir(n, -1 - x/2), where Dir(n, x) is the n-th row polynomial of the triangle A244419.
b(n, -1 - x) is the n-th row polynomial of A098493. (End)

A188647 Array read by antidiagonals of a(n) = a(n-1)*k-((k-1)/(k^n)) where a(0)=1 and k=(sqrt(x^2+1)+x)^2 for integers x>=1.

Original entry on oeis.org

1, 5, 1, 29, 17, 1, 169, 305, 37, 1, 985, 5473, 1405, 65, 1, 5741, 98209, 53353, 4289, 101, 1, 33461, 1762289, 2026009, 283009, 10301, 145, 1, 195025, 31622993, 76934989, 18674305, 1050601, 21169, 197, 1, 1136689, 567451585, 2921503573, 1232221121, 107151001, 3090529, 39005, 257, 1
Offset: 0

Views

Author

Charles L. Hohn, Apr 06 2011

Keywords

Comments

Conjecture: Given function f(x, y)=(sqrt(x^2+y)+x)^2; constant k=f(x, y); and initial term a(0)=1; then for all integers x>=1 and y=[+-]1, k may be irrational, but sequence a(n)=a(n-1)*k-((k-1)/(k^n)) always produces integer sequences; y=1 results shown here; y=-1 results are A188646.
Also square array A(n,k), n >= 1, k >= 0, read by antidiagonals, where A(n,k) is (1/sqrt(n^2+1)) * T_{2*k+1}(sqrt(n^2+1)), with T the Chebyshev polynomial of the first kind. - Seiichi Manyama, Jan 02 2019

Examples

			Square array begins:
     | 0    1       2          3             4
-----+---------------------------------------------
   1 | 1,   5,     29,       169,          985, ...
   2 | 1,  17,    305,      5473,        98209, ...
   3 | 1,  37,   1405,     53353,      2026009, ...
   4 | 1,  65,   4289,    283009,     18674305, ...
   5 | 1, 101,  10301,   1050601,    107151001, ...
   6 | 1, 145,  21169,   3090529,    451196065, ...
   7 | 1, 197,  39005,   7722793,   1529074009, ...
   8 | 1, 257,  66305,  17106433,   4413393409, ...
   9 | 1, 325, 105949,  34539049,  11259624025, ...
  10 | 1, 401, 161201,  64802401,  26050404001, ...
  11 | 1, 485, 235709, 114554089,  55673051545, ...
  12 | 1, 577, 333505, 192765313, 111418017409, ...
  13 | 1, 677, 459005, 311204713, 210996336409, ...
  14 | 1, 785, 617009, 484968289, 381184458145, ...
  15 | 1, 901, 812701, 733055401, 661215159001, ...
  ...
		

Crossrefs

Row 1 is A001653, row 2 is A007805, row 3 is A097315, row 4 is A078988, row 5 is A097727, row 6 is A097730, row 7 is A097733, row 8 is A097736, row 9 is A097739, row 10 is A097742, row 11 is A097767, row 12 is A097770, row 13 is A097773.
Column 1 is A053755.
A(n,n) gives A323012.
Cf. A188645, A188646 (f(x, y) as above with y=-1).

Formula

A(n,k) = 2 * A188645(n,k) - A(n,k-1).
A(n,k) = Sum_{j=0..k} binomial(2*k+1,2*j)*(n^2+1)^(k-j)*n^(2*j). - Seiichi Manyama, Jan 02 2019

Extensions

Edited and extended by Seiichi Manyama, Jan 02 2019

A097725 Chebyshev U(n,x) polynomial evaluated at x=51.

Original entry on oeis.org

1, 102, 10403, 1061004, 108212005, 11036563506, 1125621265607, 114802332528408, 11708712296632009, 1194173851923936510, 121794024183944892011, 12421796292910455048612, 1266901427852682470066413, 129211523844680701491725514, 13178308530729578869685936015
Offset: 0

Views

Author

Wolfdieter Lang, Aug 31 2004

Keywords

Comments

Used to form integer solutions of Pell equation a^2 - 26*b^2 =-1. See A097726 with A097727.

Programs

  • Mathematica
    ChebyshevU[Range[0,20],51] (* Harvey P. Dale, Oct 08 2012 *)
    LinearRecurrence[{102, -1},{1, 102},15] (* Ray Chandler, Aug 11 2015 *)

Formula

a(n) = 102*a(n-1) - a(n-2), n>=1, a(0)=1, a(-1):=0.
a(n) = S(n, 2*51)= U(n, 51), Chebyshev's polynomials of the second kind. See A049310.
G.f.: 1/(1-102*x+x^2).
a(n)= sum((-1)^k*binomial(n-k, k)*102^(n-2*k), k=0..floor(n/2)), n>=0.
a(n) = ((51+10*sqrt(26))^(n+1) - (51-10*sqrt(26))^(n+1))/(20*sqrt(26)).

Extensions

More terms from Harvey P. Dale, Oct 08 2012

A195622 Denominators of Pythagorean approximations to 5.

Original entry on oeis.org

20, 2020, 206040, 21014040, 2143226060, 218588044060, 22293837268080, 2273752813300080, 231900493119340100, 23651576545359390100, 2412228907133538450120, 246023696951075562522120, 25092004860102573838806140, 2559138472033511455995704140
Offset: 1

Views

Author

Clark Kimberling, Sep 22 2011

Keywords

Comments

See A195500 for a discussion and references.

Crossrefs

Programs

  • Magma
    I:=[20,2020,206040]; [n le 3 select I[n] else 101*Self(n-1) +101*Self(n-2) -Self(n-3): n in [1..40]]; // G. C. Greubel, Feb 15 2023
    
  • Mathematica
    r = 5; z = 20;
    p[{f_, n_}] := (#1[[2]]/#1[[
          1]] &)[({2 #1[[1]] #1[[2]], #1[[1]]^2 - #1[[
             2]]^2} &)[({Numerator[#1], Denominator[#1]} &)[
         Array[FromContinuedFraction[
            ContinuedFraction[(#1 + Sqrt[1 + #1^2] &)[f], #1]] &, {n}]]]];
    {a, b} = ({Denominator[#1], Numerator[#1]} &)[
      p[{r, z}]]  (* A195622, A195623 *)
    Sqrt[a^2 + b^2] (* A097727 *)
    (* by Peter J. C. Moses, Sep 02 2011 *)
    LinearRecurrence[{101,101,-1},{20,2020,206040},20] (* Harvey P. Dale, Oct 17 2021 *)
  • PARI
    Vec(20*x/((x+1)*(x^2-102*x+1)) + O(x^20)) \\ Colin Barker, Jun 03 2015
    
  • SageMath
    A097726=BinaryRecurrenceSequence(102, -1, 1, 103)
    [(5/26)*(A097726(n) - (-1)^n) for n in range(1, 41)] # G. C. Greubel, Feb 15 2023

Formula

From Colin Barker, Jun 03 2015: (Start)
a(n) = 101*a(n-1) + 101*a(n-2) - a(n-3).
G.f.: 20*x/((1+x)*(1-102*x+x^2)). (End)
a(n) = (5/26)*(A097726(n) - (-1)^n). - G. C. Greubel, Feb 15 2023

A195623 Numerators of Pythagorean approximations to 5.

Original entry on oeis.org

99, 10101, 1030199, 105070201, 10716130299, 1092940220301, 111469186340399, 11368764066500401, 1159502465596700499, 118257882726796950501, 12061144535667692250599, 1230118484755377812610601, 125460024300512869194030699, 12795692360167557279978520701, 1305035160712790329688615080799
Offset: 1

Views

Author

Clark Kimberling, Sep 22 2011

Keywords

Comments

See A195500 for discussion and list of related sequences; see A195622 for Mathematica program.

Crossrefs

Programs

  • Magma
    I:=[99,10101,1030199]; [n le 3 select I[n] else 101*Self(n-1) +101*Self(n-2) -Self(n-3): n in [1..40]]; // G. C. Greubel, Feb 16 2023
    
  • Mathematica
    Table[(5*LucasL[2*n+1,10] +2*(-1)^n)/52, {n,40}] (* G. C. Greubel, Feb 16 2023 *)
  • PARI
    Vec(-x*(x^2-102*x-99) / ((x+1)*(x^2-102*x+1)) + O(x^20)) \\ Colin Barker, Jun 03 2015
    
  • SageMath
    A097726=BinaryRecurrenceSequence(102, -1, 1, 103)
    [(1/26)*(25*A097726(n) + (-1)^n) for n in range(1, 41)] # G. C. Greubel, Feb 16 2023

Formula

From Colin Barker, Jun 03 2015: (Start)
a(n) = 101*a(n-1) + 101*a(n-2) - a(n-3).
G.f.: x*(99+102*x-x^2)/((1+x)*(1-102*x+x^2)). (End)
a(n) = (1/26)*(25*A097726(n) + (-1)^n). - G. C. Greubel, Feb 16 2023
E.g.f.: (5*exp(51*x)*(5*cosh(10*sqrt(26)*x) + sqrt(26)*sinh(10*sqrt(26)*x)) + exp(-x) - 26)/26. - Stefano Spezia, Aug 05 2024
Showing 1-7 of 7 results.