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

A002315 NSW numbers: a(n) = 6*a(n-1) - a(n-2); also a(n)^2 - 2*b(n)^2 = -1 with b(n) = A001653(n+1).

Original entry on oeis.org

1, 7, 41, 239, 1393, 8119, 47321, 275807, 1607521, 9369319, 54608393, 318281039, 1855077841, 10812186007, 63018038201, 367296043199, 2140758220993, 12477253282759, 72722761475561, 423859315570607, 2470433131948081, 14398739476117879, 83922003724759193
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

Named after the Newman-Shanks-Williams reference.
Also numbers k such that A125650(3*k^2) is an odd perfect square. Such numbers 3*k^2 form a bisection of A125651. - Alexander Adamchuk, Nov 30 2006
For positive n, a(n) corresponds to the sum of legs of near-isosceles primitive Pythagorean triangles (with consecutive legs). - Lekraj Beedassy, Feb 06 2007
Also numbers m such that m^2 is a centered 16-gonal number; or a number of the form 8k(k+1)+1, where k = A053141(m) = {0, 2, 14, 84, 492, 2870, ...}. - Alexander Adamchuk, Apr 21 2007
The lower principal convergents to 2^(1/2), beginning with 1/1, 7/5, 41/29, 239/169, comprise a strictly increasing sequence; numerators=A002315 and denominators=A001653. - Clark Kimberling, Aug 27 2008
The upper intermediate convergents to 2^(1/2) beginning with 10/7, 58/41, 338/239, 1970/1393 form a strictly decreasing sequence; essentially, numerators=A075870, denominators=A002315. - Clark Kimberling, Aug 27 2008
General recurrence is a(n) = (a(1)-1)*a(n-1) - a(n-2), a(1) >= 4, lim_{n->oo} a(n) = x*(k*x+1)^n, k = (a(1)-3), x = (1+sqrt((a(1)+1)/(a(1)-3)))/2. Examples in OEIS: a(1)=4 gives A002878. a(1)=5 gives A001834. a(1)=6 gives A030221. a(1)=7 gives A002315. a(1)=8 gives A033890. a(1)=9 gives A057080. a(1)=10 gives A057081. - Ctibor O. Zizka, Sep 02 2008
Numbers k such that (ceiling(sqrt(k*k/2)))^2 = (1+k*k)/2. - Ctibor O. Zizka, Nov 09 2009
A001109(n)/a(n) converges to cos^2(Pi/8) = 1/2 + 2^(1/2)/4. - Gary Detlefs, Nov 25 2009
The values 2(a(n)^2+1) are all perfect squares, whose square root is given by A075870. - Neelesh Bodas (neelesh.bodas(AT)gmail.com), Aug 13 2010
a(n) represents all positive integers K for which 2(K^2+1) is a perfect square. - Neelesh Bodas (neelesh.bodas(AT)gmail.com), Aug 13 2010
For positive n, a(n) equals the permanent of the (2n) X (2n) tridiagonal matrix with sqrt(8)'s along the main diagonal, and i's along the superdiagonal and subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
Integers k such that A000217(k-2) + A000217(k-1) + A000217(k) + A000217(k+1) is a square (cf. A202391). - Max Alekseyev, Dec 19 2011
Integer square roots of floor(k^2/2 - 1) or A047838. - Richard R. Forberg, Aug 01 2013
Remark: x^2 - 2*y^2 = +2*k^2, with positive k, and X^2 - 2*Y^2 = +2 reduce to the present Pell equation a^2 - 2*b^2 = -1 with x = k*X = 2*k*b and y = k*Y = k*a. (After a proposed solution for k = 3 by Alexander Samokrutov.) - Wolfdieter Lang, Aug 21 2015
If p is an odd prime, a((p-1)/2) == 1 (mod p). - Altug Alkan, Mar 17 2016
a(n)^2 + 1 = 2*b(n)^2, with b(n) = A001653(n), is the necessary and sufficient condition for a(n) to be a number k for which the diagonal of a 1 X k rectangle is an integer multiple of the diagonal of a 1 X 1 square. If squares are laid out thus along one diagonal of a horizontal 1 X a(n) rectangle, from the lower left corner to the upper right, the number of squares is b(n), and there will always be a square whose top corner lies exactly within the top edge of the rectangle. Numbering the squares 1 to b(n) from left to right, the number of the one square that has a corner in the top edge of the rectangle is c(n) = (2*b(n) - a(n) + 1)/2, which is A055997(n). The horizontal component of the corner of the square in the edge of the rectangle is also an integer, namely d(n) = a(n) - b(n), which is A001542(n). - David Pasino, Jun 30 2016
(a(n)^2)-th triangular number is a square; a(n)^2 = A008843(n) is a subsequence of A001108. - Jaroslav Krizek, Aug 05 2016
a(n-1)/A001653(n) is the closest rational approximation of sqrt(2) with a numerator not larger than a(n-1). These rational approximations together with those obtained from the sequences A001541 and A001542 give a complete set of closest rational approximations of sqrt(2) with restricted numerator or denominator. a(n-1)/A001653(n) < sqrt(2). - A.H.M. Smeets, May 28 2017
Consider the quadrant of a circle with center (0,0) bounded by the positive x and y axes. Now consider, as the start of a series, the circle contained within this quadrant which kisses both axes and the outer bounding circle. Consider further a succession of circles, each kissing the x-axis, the outer bounding circle, and the previous circle in the series. See Holmes link. The center of the n-th circle in this series is ((A001653(n)*sqrt(2)-1)/a(n-1), (A001653(n)*sqrt(2)-1)/a(n-1)^2), the y-coordinate also being its radius. It follows that a(n-1) is the cotangent of the angle subtended at point (0,0) by the center of the n-th circle in the series with respect to the x-axis. - Graham Holmes, Aug 31 2019
There is a link between the two sequences present at the numerator and at the denominator of the fractions that give the coordinates of the center of the kissing circles. A001653 is the sequence of numbers k such that 2*k^2 - 1 is a square, and here, we have 2*A001653(n)^2 - 1 = a(n-1)^2. - Bernard Schott, Sep 02 2019
Let G be a sequence satisfying G(i) = 2*G(i-1) + G(i-2) for arbitrary integers i and without regard to the initial values of G. Then a(n) = (G(i+4*n+2) - G(i))/(2*G(i+2*n+1)) as long as G(i+2*n+1) != 0. - Klaus Purath, Mar 25 2021
All of the positive integer solutions of a*b+1=x^2, a*c+1=y^2, b*c+1=z^2, x+z=2*y, 0 < a < b < c are given by a=A001542(n), b=A005319(n), c=A001542(n+1), x=A001541(n), y=A001653(n+1), z=A002315(n) with 0 < n. - Michael Somos, Jun 26 2022
3*a(n-1) is the n-th almost Lucas-cobalancing number of second type (see Tekcan and Erdem). - Stefano Spezia, Nov 26 2022
In Moret-Blanc (1881) on page 259 some solution of m^2 - 2n^2 = -1 are listed. The values of m give this sequence, and the values of n give A001653. - Michael Somos, Oct 25 2023
From Klaus Purath, May 11 2024: (Start)
For any two consecutive terms (a(n), a(n+1)) = (x,y): x^2 - 6xy + y^2 = 8 = A028884(1). In general, the following applies to all sequences (t) satisfying t(i) = 6t(i-1) - t(i-2) with t(0) = 1 and two consecutive terms (x,y): x^2 - 6xy + y^2 = A028884(t(1)-6). This includes and interprets the Feb 04 2014 comment on A001541 by Colin Barker as well as the Mar 17 2021 comment on A054489 by John O. Oladokun and the Sep 28 2008 formula on A038723 by Michael Somos. By analogy to this, for three consecutive terms (x,y,z) y^2 - xz = A028884(t(1)-6) always applies.
If (t) is a sequence satisfying t(k) = 7t(k-1) - 7t(k-2) + t(k-3) or t(k) = 6t(k-1) - t(k-2) without regard to initial values and including this sequence itself, then a(n) = (t(k+2n+1) - t(k))/(t(k+n+1) - t(k+n)) always applies, as long as t(k+n+1) - t(k+n) != 0 for integer k and n >= 0. (End)

Examples

			G.f. = 1 + 7*x + 41*x^2 + 239*x^3 + 1393*x^4 + 8119*x^5 + 17321*x^6 + ... - _Michael Somos_, Jun 26 2022

References

  • Julio R. Bastida, Quadratic properties of a linearly recurrent sequence. Proceedings of the Tenth Southeastern Conference on Combinatorics, Graph Theory and Computing (Florida Atlantic Univ., Boca Raton, Fla., 1979), pp. 163-166, Congress. Numer., XXIII-XXIV, Utilitas Math., Winnipeg, Man., 1979. MR0561042 (81e:10009)
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 256.
  • Paulo Ribenboim, The Book of Prime Number Records. Springer-Verlag, NY, 2nd ed., 1989, p. 288.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 247.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • P.-F. Teilhet, Reply to Query 2094, L'Intermédiaire des Mathématiciens, 10 (1903), 235-238.
  • P.-F. Teilhet, Query 2376, L'Intermédiaire des Mathématiciens, 11 (1904), 138-139. - N. J. A. Sloane, Mar 08 2022

Links

Crossrefs

Bisection of A001333. Cf. A001109, A001653. A065513(n)=a(n)-1.
First differences of A001108 and A055997. Bisection of A084068 and A088014. Cf. A077444.
Cf. A125650, A125651, A125652.
Row sums of unsigned triangle A127675.
Cf. A053141, A075870. Cf. A000045, A002878, A004146, A026003, A100047, A119915, A192425, A088165 (prime subsequence), A057084 (binomial transform), A108051 (inverse binomial transform).
See comments in A301383.
Cf. similar sequences of the type (1/k)*sinh((2*n+1)*arcsinh(k)) listed in A097775.
Cf. A000129, A001653, A003499, A358682, A002203.

Programs

  • Haskell
    a002315 n = a002315_list !! n
a002315_list = 1 : 7 : zipWith (-) (map (* 6) (tail a002315_list)) a002315_list
-- Reinhard Zumkeller, Jan 10 2012
  • Magma
    I:=[1,7]; [n le 2 select I[n] else 6*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Mar 22 2015
  • Maple
    A002315 := proc(n)
    option remember;
    if n = 0 then
        1 ;
    elif n = 1 then
        7;
    else
        6*procname(n-1)-procname(n-2) ;
    end if;
end proc: # Zerinvary Lajos, Jul 26 2006, modified R. J. Mathar, Apr 30 2017
a:=n->abs(Im(simplify(ChebyshevT(2*n+1,I)))):seq(a(n),n=0..20); # Leonid Bedratyuk, Dec 17 2017
# third Maple program:
a:= n-> (<<0|1>, <-1|6>>^n. <<1, 7>>)[1, 1]:
seq(a(n), n=0..22);  # Alois P. Heinz, Aug 25 2024
  • Mathematica
    a[0] = 1; a[1] = 7; a[n_] := a[n] = 6a[n - 1] - a[n - 2]; Table[ a[n], {n, 0, 20}] (* Robert G. Wilson v, Jun 09 2004 *)
Transpose[NestList[Flatten[{Rest[#],ListCorrelate[{-1,6},#]}]&, {1,7},20]][[1]]  (* Harvey P. Dale, Mar 23 2011 *)
Table[ If[n>0, a=b; b=c; c=6b-a, b=-1; c=1], {n, 0, 20}] (* Jean-François Alcover, Oct 19 2012 *)
LinearRecurrence[{6, -1}, {1, 7}, 20] (* Bruno Berselli, Apr 03 2018 *)
a[ n_] := -I*(-1)^n*ChebyshevT[2*n + 1, I]; (* Michael Somos, Jun 26 2022 *)
  • PARI
    {a(n) = subst(poltchebi(abs(n+1)) - poltchebi(abs(n)), x, 3)/2};
  • PARI
    {a(n) = if(n<0, -a(-1-n), polsym(x^2-2*x-1, 2*n+1)[2*n+2]/2)};
  • PARI
    {a(n) = my(w=3+quadgen(32)); imag((1+w)*w^n)};
  • PARI
    for (i=1,10000,if(Mod(sigma(i^2+1,2),2)==1,print1(i,",")))
  • PARI
    {a(n) = -I*(-1)^n*polchebyshev(2*n+1, 1, I)}; /* Michael Somos, Jun 26 2022 */

Formula

a(n) = (1/2)*((1+sqrt(2))^(2*n+1) + (1-sqrt(2))^(2*n+1)).
a(n) = A001109(n)+A001109(n+1).
a(n) = (1+sqrt(2))/2*(3+sqrt(8))^n+(1-sqrt(2))/2*(3-sqrt(8))^n. - Ralf Stephan, Feb 23 2003
a(n) = sqrt(2*(A001653(n+1))^2-1), n >= 0. [Pell equation a(n)^2 - 2*Pell(2*n+1)^2 = -1. - Wolfdieter Lang, Jul 11 2018]
G.f.: (1 + x)/(1 - 6*x + x^2). - Simon Plouffe in his 1992 dissertation
a(n) = S(n, 6)+S(n-1, 6) = S(2*n, sqrt(8)), S(n, x) = U(n, x/2) are Chebyshev's polynomials of the 2nd kind. Cf. A049310. S(n, 6)= A001109(n+1).
a(n) ~ (1/2)*(sqrt(2) + 1)^(2*n+1). - Joe Keane (jgk(AT)jgk.org), May 15 2002
Limit_{n->oo} a(n)/a(n-1) = 3 + 2*sqrt(2). - Gregory V. Richardson, Oct 06 2002
Let q(n, x) = Sum_{i=0..n} x^(n-i)*binomial(2*n-i, i); then (-1)^n*q(n, -8) = a(n). - Benoit Cloitre, Nov 10 2002
With a=3+2*sqrt(2), b=3-2*sqrt(2): a(n) = (a^((2n+1)/2)-b^((2n+1)/2))/2. a(n) = A077444(n)/2. - Mario Catalani (mario.catalani(AT)unito.it), Mar 31 2003
a(n) = Sum_{k=0..n} 2^k*binomial(2*n+1, 2*k). - Zoltan Zachar (zachar(AT)fellner.sulinet.hu), Oct 08 2003
Same as: i such that sigma(i^2+1, 2) mod 2 = 1. - Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 26 2004
a(n) = L(n, -6)*(-1)^n, where L is defined as in A108299; see also A001653 for L(n, +6). - Reinhard Zumkeller, Jun 01 2005
a(n) = A001652(n)+A046090(n); e.g., 239=119+120. - Charlie Marion, Nov 20 2003
A001541(n)*a(n+k) = A001652(2n+k) + A001652(k)+1; e.g., 3*1393 = 4069 + 119 + 1; for k > 0, A001541(n+k)*a(n) = A001652(2n+k) - A001652(k-1); e.g., 99*7 = 696 - 3. - Charlie Marion, Mar 17 2003
a(n) = Jacobi_P(n,1/2,-1/2,3)/Jacobi_P(n,-1/2,1/2,1). - Paul Barry, Feb 03 2006
P_{2n}+P_{2n+1} where P_i are the Pell numbers (A000129). Also the square root of the partial sums of Pell numbers: P_{2n}+P_{2n+1} = sqrt(Sum_{i=0..4n+1} P_i) (Santana and Diaz-Barrero, 2006). - David Eppstein, Jan 28 2007
a(n) = 2*A001652(n) + 1 = 2*A046729(n) + (-1)^n. - Lekraj Beedassy, Feb 06 2007
a(n) = sqrt(A001108(2*n+1)). - Anton Vrba (antonvrba(AT)yahoo.com), Feb 14 2007
a(n) = sqrt(8*A053141(n)*(A053141(n) + 1) + 1). - Alexander Adamchuk, Apr 21 2007
a(n+1) = 3*a(n) + sqrt(8*a(n)^2 + 8), a(1)=1. - Richard Choulet, Sep 18 2007
a(n) = A001333(2*n+1). - Ctibor O. Zizka, Aug 13 2008
a(n) = third binomial transform of 1, 4, 8, 32, 64, 256, 512, ... . - Al Hakanson (hawkuu(AT)gmail.com), Aug 15 2009
a(n) = (-1)^(n-1)*(1/sqrt(-1))*cos((2*n - 1)*arcsin(sqrt(2)). - Artur Jasinski, Feb 17 2010 *WRONG*
a(n+k) = A001541(k)*a(n) + 4*A001109(k)*A001653(n); e.g., 8119 = 17*239 + 4*6*169. - Charlie Marion, Feb 04 2011
In general, a(n+k) = A001541(k)*a(n)) + sqrt(A001108(2k)*(a(n)^2+1)). See Sep 18 2007 entry above. - Charlie Marion, Dec 07 2011
a(n) = floor((1+sqrt(2))^(2n+1))/2. - Thomas Ordowski, Jun 12 2012
(a(2n-1) + a(2n) + 8)/(8*a(n)) = A001653(n). - Ignacio Larrosa Cañestro, Jan 02 2015
(a(2n) + a(2n-1))/a(n) = 2*sqrt(2)*( (1 + sqrt(2))^(4*n) - (1 - sqrt(2))^(4*n))/((1 + sqrt(2))^(2*n+1) + (1 - sqrt(2))^(2*n+1)). [This was my solution to problem 5325, School Science and Mathematics 114 (No. 8, Dec 2014).] - Henry Ricardo, Feb 05 2015
From Peter Bala, Mar 22 2015: (Start)
The aerated sequence (b(n))n>=1 = [1, 0, 7, 0, 41, 0, 239, 0, ...] is a fourth-order linear divisibility sequence; that is, if n | m then b(n) | b(m). It is the case P1 = 0, P2 = -4, Q = -1 of the 3-parameter family of divisibility sequences found by Williams and Guy. See A100047.
b(n) = 1/2*((-1)^n - 1)*Pell(n) + 1/2*(1 + (-1)^(n+1))*Pell(n+1). The o.g.f. is x*(1 + x^2)/(1 - 6*x^2 + x^4).
Exp( Sum_{n >= 1} 2*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 2*A026003(n-1)*x^n.
Exp( Sum_{n >= 1} (-2)*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 2*A026003(n-1)*(-x)^n.
Exp( Sum_{n >= 1} 4*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 4*Pell(n)*x^n.
Exp( Sum_{n >= 1} (-4)*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 4*Pell(n)*(-x)^n.
Exp( Sum_{n >= 1} 8*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 8*A119915(n)*x^n.
Exp( Sum_{n >= 1} (-8)*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 8*A119915(n)*(-x)^n. Cf. A002878, A004146, A113224, and A192425. (End)
E.g.f.: (sqrt(2)*sinh(2*sqrt(2)*x) + cosh(2*sqrt(2)*x))*exp(3*x). - Ilya Gutkovskiy, Jun 30 2016
a(n) = Sum_{k=0..n} binomial(n,k) * 3^(n-k) * 2^k * 2^ceiling(k/2). - David Pasino, Jul 09 2016
a(n) = A001541(n) + 2*A001542(n). - A.H.M. Smeets, May 28 2017
a(n+1) = 3*a(n) + 4*b(n), b(n+1) = 2*a(n) + 3*b(n), with b(n)=A001653(n). - Zak Seidov, Jul 13 2017
a(n) = |Im(T(2n-1,i))|, i=sqrt(-1), T(n,x) is the Chebyshev polynomial of the first kind, Im is the imaginary part of a complex number, || is the absolute value. - Leonid Bedratyuk, Dec 17 2017
a(n) = sinh((2*n + 1)*arcsinh(1)). - Bruno Berselli, Apr 03 2018
a(n) = 5*a(n-1) + A003499(n-1), a(0) = 1. - Ivan N. Ianakiev, Aug 09 2019
From Klaus Purath, Mar 25 2021: (Start)
a(n) = A046090(2*n)/A001541(n).
a(n+1)*a(n+2) = a(n)*a(n+3) + 48.
a(n)^2 + a(n+1)^2 = 6*a(n)*a(n+1) + 8.
a(n+1)^2 = a(n)*a(n+2) + 8.
a(n+1) = a(n) + 2*A001541(n+1).
a(n) = 2*A046090(n) - 1. (End)
3*a(n-1) = sqrt(8*b(n)^2 + 8*b(n) - 7), where b(n) = A358682(n). - Stefano Spezia, Nov 26 2022
a(n) = -(-1)^n - 2 + Sum_{i=0..n} A002203(i)^2. - Adam Mohamed, Aug 22 2024
From Peter Bala, May 09 2025: (Start)
a(n) = Dir(n, 3), where Dir(n, x) denotes the n-th row polynomial of the triangle A244419.
For arbitrary x, a(n+x)^2 - 6*a(n+x)*a(n+x+1) + a(n+x+1)^2 = 8 with a(n) := (1/2)*((1+sqrt(2))^(2*n+1) + (1-sqrt(2))^(2*n+1)) as above. The particular case x = 0 is noted above,
a(n+1/2) = sqrt(2) * A001542(n+1).
Sum_{n >= 1} (-1)^(n+1)/(a(n) - 1/a(n)) = 1/8 (telescoping series: for n >= 1, 1/(a(n) - 1/a(n)) = 1/A081554(n) + 1/A081554(n+1)).
Product_{n >= 1} (a(n) + 1)/(a(n) - 1) = sqrt(2) (telescoping product: Product_{n = 1..k} ((a(n) + 1)/(a(n) - 1))^2 = 2*(1 - 1/A055997(k+2))). (End)

A001542 a(n) = 6*a(n-1) - a(n-2) for n > 1, a(0)=0 and a(1)=2.

Original entry on oeis.org

0, 2, 12, 70, 408, 2378, 13860, 80782, 470832, 2744210, 15994428, 93222358, 543339720, 3166815962, 18457556052, 107578520350, 627013566048, 3654502875938, 21300003689580, 124145519261542, 723573111879672
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

Consider the equation core(x) = core(2x+1) where core(x) is the smallest number such that x*core(x) is a square: solutions are given by a(n)^2, n > 0. - Benoit Cloitre, Apr 06 2002
Terms > 0 give numbers k which are solutions to the inequality |round(sqrt(2)*k)/k - sqrt(2)| < 1/(2*sqrt(2)*k^2). - Benoit Cloitre, Feb 06 2006
Also numbers m such that A125650(6*m^2) is an even perfect square, where A124650(m) is a numerator of m*(m+3)/(4*(m+1)*(m+2)) = Sum_{k=1..m} 1/(k*(k+1)*(k+2)). Sequence A033581 is a bisection of A125651. - Alexander Adamchuk, Nov 30 2006
The upper principal convergents to 2^(1/2), beginning with 3/2, 17/12, 99/70, 577/408, comprise a strictly decreasing sequence; essentially, numerators = A001541 and denominators = {a(n)}. - Clark Kimberling, Aug 26 2008
Even Pell numbers. - Omar E. Pol, Dec 10 2008
Numbers k such that 2*k^2+1 is a square. - Vladimir Joseph Stephan Orlovsky, Feb 19 2010
These are the integer square roots of the Half-Squares, A007590(k), which occur at values of k given by A001541. Also the numbers produced by adding m + sqrt(floor(m^2/2) + 1) when m is in A002315. See array in A227972. - Richard R. Forberg, Aug 31 2013
A001541(n)/a(n) is the closest rational approximation of sqrt(2) with a denominator not larger than a(n), and 2*a(n)/A001541(n) is the closest rational approximation of sqrt(2) with a numerator not larger than 2*a(n). These rational approximations together with those obtained from the sequences A001653 and A002315 give a complete set of closest rational approximations of sqrt(2) with restricted numerator as well as denominator. - A.H.M. Smeets, May 28 2017
Conjecture: Numbers k such that c/m < k for all natural a^2 + b^2 = c^2 (Pythagorean triples), a < b < c and a+b+c = m. Numbers which correspondingly minimize c/m are A002939. - Lorraine Lee, Jan 31 2020
All of the positive integer solutions of a*b + 1 = x^2, a*c + 1 = y^2, b*c + 1 = z^2, x + z = 2*y, 0 < a < b < c are given by a=a(n), b=A005319(n), c=a(n+1), x=A001541(n), y=A001653(n+1), z=A002315(n) with 0 < n. - Michael Somos, Jun 26 2022

Examples

			G.f. = 2*x + 12*x^2 + 70*x^3 + 408*x^4 + 2378*x^5 + 13860*x^6 + ...

References

  • Jay Kappraff, Beyond Measure, A Guided Tour Through Nature, Myth and Number, World Scientific, 2002; pp. 480-481.
  • Thomas Koshy, Fibonacci and Lucas Numbers with Applications, 2001, Wiley, pp. 77-79.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 257-258.
  • P.-F. Teilhet, Query 2376, L'Intermédiaire des Mathématiciens, 11 (1904), 138-139. - N. J. A. Sloane, Mar 08 2022

Links

Crossrefs

Bisection of Pell numbers A000129: {a(n)} and A001653(n+1), n >= 0.
Cf. A001108, A001353, A001541, A001835, A003499, A007805, A007913, A115598, A125650, A125651, A125652.

Programs

  • GAP
    a:=[0,2];; for n in [3..20] do a[n]:=6*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 23 2019
  • Haskell
    a001542 n = a001542_list !! n
a001542_list =
   0 : 2 : zipWith (-) (map (6 *) $ tail a001542_list) a001542_list
-- Reinhard Zumkeller, Aug 14 2011
  • Magma
    I:=[0,2]; [n le 2 select I[n] else 6*Self(n-1) -Self(n-2): n in [1..20]]; // G. C. Greubel, Dec 23 2019
  • Maple
    A001542:=2*z/(1-6*z+z**2); # conjectured by Simon Plouffe in his 1992 dissertation
seq(combinat:-fibonacci(2*n, 2), n = 0..20); # Peter Luschny, Jun 28 2018
  • Mathematica
    LinearRecurrence[{6, -1}, {0, 2}, 30] (* Harvey P. Dale, Jun 11 2011 *)
Fibonacci[2*Range[0,20], 2] (* G. C. Greubel, Dec 23 2019 *)
Table[2 ChebyshevU[-1 + n, 3], {n, 0, 20}] (* Herbert Kociemba, Jun 05 2022 *)
  • Maxima
    a[0]:0$
a[1]:2$
a[n]:=6*a[n-1]-a[n-2]$
A001542(n):=a[n]$
makelist(A001542(x),x,0,30); /* Martin Ettl, Nov 03 2012 */
  • PARI
    {a(n) = imag( (3 + 2*quadgen(8))^n )}; /* Michael Somos, Jan 20 2017 */
  • PARI
    vector(21, n, 2*polchebyshev(n-1, 2, 33) ) \\ G. C. Greubel, Dec 23 2019
  • Python
    l=[0, 2]
for n in range(2, 51): l+=[6*l[n - 1] - l[n - 2], ]
print(l) # Indranil Ghosh, Jun 06 2017
  • Sage
    [2*chebyshev_U(n-1,3) for n in (0..20)] # G. C. Greubel, Dec 23 2019

Formula

a(n) = 2*A001109(n).
a(n) = ((3+2*sqrt(2))^n - (3-2*sqrt(2))^n) / (2*sqrt(2)).
G.f.: 2*x/(1-6*x+x^2).
a(n) = sqrt(2*(A001541(n))^2 - 2)/2. - Barry E. Williams, May 07 2000
a(n) = (C^(2n) - C^(-2n))/sqrt(8) where C = sqrt(2) + 1. - Gary W. Adamson, May 11 2003
For all terms x of the sequence, 2*x^2 + 1 is a square. Limit_{n->oo} a(n)/a(n-1) = 3 + 2*sqrt(2). - Gregory V. Richardson, Oct 10 2002
For n > 0: a(n) = A001652(n) + A046090(n) - A001653(n); e.g., 70 = 119 + 120 - 169. Also a(n) = A001652(n - 1) + A046090(n - 1) + A001653(n - 1); e.g., 70 = 20 + 21 + 29. Also a(n)^2 + 1 = A001653(n - 1)*A001653(n); e.g., 12^2 + 1 = 145 = 5*29. Also a(n + 1)^2 = A084703(n + 1) = A001652(n)*A001652(n + 1) + A046090(n)*A046090(n + 1). - Charlie Marion, Jul 01 2003
a(n) = ((1+sqrt(2))^(2*n) - (1-sqrt(2))^(2*n))/(2*sqrt(2)). - Antonio Alberto Olivares, Dec 24 2003
2*A001541(k)*A001653(n)*A001653(n+k) = A001653(n)^2 + A001653(n+k)^2 + a(k)^2; e.g., 2*3*5*29 = 5^2 + 29^2 + 2^2; 2*99*29*5741 = 29^2 + 5741^2 + 70^2. - Charlie Marion, Oct 12 2007
a(n) = sinh(2*n*arcsinh(1))/sqrt(2). - Herbert Kociemba, Apr 24 2008
For n > 0, a(n) = A001653(n) + A002315(n-1). - Richard R. Forberg, Aug 31 2013
a(n) = 3*a(n-1) + 2*A001541(n-1); e.g., a(4) = 70 = 3*12 + 2*17. - Zak Seidov, Dec 19 2013
a(n)^2 + 1^2 = A115598(n)^2 + (A115598(n)+1)^2. - Hermann Stamm-Wilbrandt, Jul 27 2014
Sum {n >= 1} 1/( a(n) + 1/a(n) ) = 1/2. - _Peter Bala, Mar 25 2015
E.g.f.: exp(3*x)*sinh(2*sqrt(2)*x)/sqrt(2). - Ilya Gutkovskiy, Dec 07 2016
A007814(a(n)) = A001511(n). See Mathematical Reflections link. - Michel Marcus, Jan 06 2017
a(n) = -a(-n) for all n in Z. - Michael Somos, Jan 20 2017
From A.H.M. Smeets, May 28 2017: (Start)
A051009(n) = a(2^(n-2)).
a(2n) = 2*a(2)*A001541(n).
A001541(n)/a(n) > sqrt(2) > 2*a(n)/A001541(n). (End)
a(A298210(n)) = A002349(2*n^2). - A.H.M. Smeets, Jan 25 2018
a(n) = A000129(n)*A002203(n). - Adam Mohamed, Jul 20 2024

A125650 Numerator of n(n+3)/(4(n+1)(n+2)) = sum(k=1..n, 1/(k(k+1)(k+2)) ).

Original entry on oeis.org

0, 1, 5, 9, 7, 5, 27, 35, 11, 27, 65, 77, 45, 26, 119, 135, 38, 85, 189, 209, 115, 63, 275, 299, 81, 175, 377, 405, 217, 116, 495, 527, 140, 297, 629, 665, 351, 185, 779, 819, 215, 451, 945, 989, 517, 270, 1127, 1175, 306, 637, 1325, 1377, 715, 371, 1539, 1595, 413, 855
Offset: 0

Views

Author

Alexander Adamchuk, Nov 29 2006

Keywords

Comments

3^2 divides a(3k). p divides a(p) for an odd prime p. p divides a(p-3) for prime p>3. p^k divides a(p^k) for an odd prime p. a(n) = m^2 is a perfect square for n = {1,3,24,147,864,5043,29400,171363,...} = A125651(n). Corresponding numbers m such that m^2 = a[ A125651(n) ] are listed in A125652(n) = {1,3,9,105,306,3567,10395,121173,...}.

Examples

			The rationals n(n+3)/(4(n+1)(n+2)) = a(n)/A230328(n) begin:
0, 1/6, 5/24, 9/40, 7/30, 5/21, 27/112, 35/144, 11/45, 27/110, 65/264, 77/312, 45/182, 26/105, 119/480, ... - _Wolfdieter Lang_, Mar 08 2018

References

  • L. B.W. Jolley, Summation of Series, Second revised ed., Dover, 1961, p.38, (201).

Links

Crossrefs

Cf. A125651, A125652. A160050, A230328 (denominators).

Programs

  • Magma
    [Numerator(n*(n+3)/(4*(n+1)*(n+2))): n in [0..60]]; // Vincenzo Librandi, May 21 2012
  • Mathematica
    Table[Numerator[n(n+3)/(4(n+1)(n+2))],{n,0,100}]
  • PARI
    a(n)=n*(n+3)/2^min(3,valuation(n*(n+3),2)); \\ Max Alekseyev, Jan 11 2007

Formula

a(n) = Sum_{k=1..n} 1/(k(k+1)(k+2)).
a(n) = n*(n+3)/2^min(3,valuation(n*(n+3),2)). a(n)=n*(n+3)/4 for n=1 or 4 (mod 8); a(n)=n*(n+3)/8 for n=0 or 5 (mod 8); a(n) = n*(n+3)/2 for n=2, 3, 6, or 7 (mod 8). - Max Alekseyev, Jan 11 2007
a(n) = A106609(n)*A106609(n+3). - Paul Curtz, Jan 13 2011
G.f.: x*(x^19 -2*x^18 +3*x^17 -5*x^16 +3*x^15 -6*x^14 +7*x^13 -11*x^12 -12*x^11 +24*x^10 -36*x^9 +24*x^8 -38*x^7 +28*x^6 -18*x^5 -3*x^4 -2*x -1) / ((x-1)^3*(x^2+1)^3*(x^4+1)^3). - Colin Barker, Feb 21 2013
G.f. for rationals r(n) = a(n)/A230328(n): (1/4)*(1 - hypergeometric([1, 2], [3], -x/(1-x)))/(1-x) = (- 2*x + 3*x^2 + 2*(2*x - (1 + x^2))*log(1-x))/(4*(1-x)*x^2). For the r(n) formula see Jolley's general remark (201) on p.38. Thanks to Gary Detlefs for pointing to this remark. - Wolfdieter Lang, Mar 08 2018

A125651 Numbers k such that A125650(k) is a perfect square.

Original entry on oeis.org

1, 3, 24, 147, 864, 5043, 29400, 171363, 998784, 5821347, 33929304, 197754483, 1152597600, 6717831123, 39154389144, 228208503747, 1330096633344, 7752371296323, 45184131144600, 263352415571283, 1534930362283104, 8946229758127347, 52142448186480984
Offset: 1

Views

Author

Alexander Adamchuk, Nov 29 2006

Keywords

Comments

Corresponding numbers m such that m^2 = A125650(a(n)) are listed in A125652.
3 divides a(n) for n>1. For n>1 a(n) = 3*A001108(n-1), where A001108(k) = {0, 1, 8, 49, 288, 1681, ...}, A001108(k)-th triangular number is a square. - Alexander Adamchuk, Jan 19 2007
Disregarding the term 1, numbers k such that A071910(k) is a nonzero square; i.e., numbers k such that A000096(k) = k*(k+3)/2 is a nonzero square. - Rick L. Shepherd, Jul 13 2012

Examples

			a(2)=3 because A125650(3)=9=3^2; a(3)=24 because A125650(24)=81=9^2.

Links

Crossrefs

Cf. A125650, A125652.
Cf. A001108, A071910, A000096, A000217.

Programs

  • Magma
    I:=[1, 3, 24, 147]; [n le 4 select I[n] else 7*Self(n-1)-7*Self(n-2)+Self(n-3): n in [1..30]]; // Vincenzo Librandi, May 21 2012
  • Mathematica
    Join[{1},LinearRecurrence[{7,-7,1},{3,24,147},35]] (* or *) CoefficientList[Series[(-1+x(4+(-10+x)x))/((-1+x)(1+(-6+x) x)),{x,0,35}],x] (* Harvey P. Dale, May 15 2011 *)

Formula

For n>1, a(n+2) = 6*a(n+1) - a(n) + 6.
For n>1, a(n) = ((3+2*sqrt(2))^(n-1) + (3-2*sqrt(2))^(n-1))*3/4 - 3/2.
For n>0, a(2*n) = 3*A002315(n-1)^2; a(2*n+1) = 6*A001542(n)^2.
a(n) = 3*A001108(n-1) for n>1. - Alexander Adamchuk, Jan 19 2007
From Harvey P. Dale, May 15 2011: (Start)
For n>1, a(2)=3, a(3)=24, a(4)=147, a(n) = 7*a(n-1)-7*a(n-2)+a(n-3).
G.f.: x*(-1+x*(4+(-10+x)*x))/((-1+x)*(1+(-6+x)*x)). (End)

Extensions

Edited by Max Alekseyev, Jan 11 2007
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.