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-10 of 24 results. Next

A001109 a(n)^2 is a triangular number: a(n) = 6*a(n-1) - a(n-2) with a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 6, 35, 204, 1189, 6930, 40391, 235416, 1372105, 7997214, 46611179, 271669860, 1583407981, 9228778026, 53789260175, 313506783024, 1827251437969, 10650001844790, 62072759630771, 361786555939836, 2108646576008245, 12290092900109634, 71631910824649559, 417501372047787720
Offset: 0

Views

Author

Keywords

Comments

8*a(n)^2 + 1 = 8*A001110(n) + 1 = A055792(n+1) is a perfect square. - Gregory V. Richardson, Oct 05 2002
For n >= 2, A001108(n) gives exactly the positive integers m such that 1,2,...,m has a perfect median. The sequence of associated perfect medians is the present sequence. Let a_1,...,a_m be an (ordered) sequence of real numbers, then a term a_k is a perfect median if Sum_{j=1..k-1} a_j = Sum_{j=k+1..m} a_j. See Puzzle 1 in MSRI Emissary, Fall 2005. - Asher Auel, Jan 12 2006
(a(n), b(n)) where b(n) = A082291(n) are the integer solutions of the equation 2*binomial(b,a) = binomial(b+2,a). - Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de); comment revised by Michael Somos, Apr 07 2003
This sequence gives the values of y in solutions of the Diophantine equation x^2 - 8y^2 = 1. It also gives the values of the product xy where (x,y) satisfies x^2 - 2y^2 = +-1, i.e., a(n) = A001333(n)*A000129(n). a(n) also gives the inradius r of primitive Pythagorean triangles having legs whose lengths are consecutive integers, with corresponding semiperimeter s = a(n+1) = {A001652(n) + A046090(n) + A001653(n)}/2 and area rs = A029549(n) = 6*A029546(n). - Lekraj Beedassy, Apr 23 2003 [edited by Jon E. Schoenfield, May 04 2014]
n such that 8*n^2 = floor(sqrt(8)*n*ceiling(sqrt(8)*n)). - Benoit Cloitre, May 10 2003
For n > 0, ratios a(n+1)/a(n) may be obtained as convergents to continued fraction expansion of 3+sqrt(8): either successive convergents of [6;-6] or odd convergents of [5;1, 4]. - Lekraj Beedassy, Sep 09 2003
a(n+1) + A053141(n) = A001108(n+1). Generating floretion: - 2'i + 2'j - 'k + i' + j' - k' + 2'ii' - 'jj' - 2'kk' + 'ij' + 'ik' + 'ji' + 'jk' - 2'kj' + 2e ("jes" series). - Creighton Dement, Dec 16 2004
Kekulé numbers for certain benzenoids (see the Cyvin-Gutman reference). - Emeric Deutsch, Jun 19 2005
Number of D steps on the line y=x in all Delannoy paths of length n (a Delannoy path of length n is a path from (0,0) to (n,n), consisting of steps E=(1,0), N=(0,1) and D=(1,1)). Example: a(2)=6 because in the 13 (=A001850(2)) Delannoy paths of length 2, namely (DD), (D)NE, (D)EN, NE(D), NENE, NEEN, NDE, NNEE, EN(D), ENNE, ENEN, EDN and EENN, we have altogether six D steps on the line y=x (shown between parentheses). - Emeric Deutsch, Jul 07 2005
Define a T-circle to be a first-quadrant circle with integral radius that is tangent to the x- and y-axes. Such a circle has coordinates equal to its radius. Let C(0) be the T-circle with radius 1. Then for n > 0, define C(n) to be the smallest T-circle that does not intersect C(n-1). C(n) has radius a(n+1). Cf. A001653. - Charlie Marion, Sep 14 2005
Numbers such that there is an m with t(n+m)=2t(m), where t(n) are the triangular numbers A000217. For instance, t(20)=2*t(14)=210, so 6 is in the sequence. - Floor van Lamoen, Oct 13 2005
One half the bisection of the Pell numbers (A000129). - Franklin T. Adams-Watters, Jan 08 2006
Pell trapezoids: for n > 0, a(n) = (A000129(n-1)+A000129(n+1))*A000129(n)/2; see also A084158. - Charlie Marion, Apr 01 2006
Tested for 2 < p < 27: If and only if 2^p - 1 (the Mersenne number M(p)) is prime then M(p) divides a(2^(p-1)). - Kenneth J Ramsey, May 16 2006
If 2^p - 1 is prime then M(p) divides a(2^(p-1)-1). - Kenneth J Ramsey, Jun 08 2006; comment corrected by Robert Israel, Mar 18 2007
If 8*n+5 and 8*n+7 are twin primes then their product divides a(4*n+3). - Kenneth J Ramsey, Jun 08 2006
If p is an odd prime, then if p == 1 or 7 (mod 8), then a((p-1)/2) == 0 (mod p) and a((p+1)/2) == 1 (mod p); if p == 3 or 5 (mod 8), then a((p-1)/2) == 1 (mod p) and a((p+1)/2) == 0 (mod p). Kenneth J Ramsey's comment about twin primes follows from this. - Robert Israel, Mar 18 2007
a(n)*(a(n+b) - a(b-2)) = (a(n+1)+1)*(a(n+b-1) - a(b-1)). This identity also applies to any series a(0) = 0 a(1) = 1 a(n) = b*a(n-1) - a(n-2). - Kenneth J Ramsey, Oct 17 2007
For n < 0, let a(n) = -a(-n). Then (a(n+j) + a(k+j)) * (a(n+b+k+j) - a(b-j-2)) = (a(n+j+1) + a(k+j+1)) * (a(n+b+k+j-1) - a(b-j-1)). - Charlie Marion, Mar 04 2011
Sequence gives y values of the Diophantine equation: 0+1+2+...+x = y^2. If (a,b) and (c,d) are two consecutive solutions of the Diophantine equation: 0+1+2+...+x = y^2 with aMohamed Bouhamida, Aug 29 2009
If (p,q) and (r,s) are two consecutive solutions of the Diophantine equation: 0+1+2+...+x = y^2 with p < r then r = 3*p+4*q+1 and s = 2*p+3*q+1. - Mohamed Bouhamida, Sep 02 2009
a(n)/A002315(n) converges to cos^2(Pi/8) (see A201488). - Gary Detlefs, Nov 25 2009
Binomial transform of A086347. - Johannes W. Meijer, Aug 01 2010
If x=a(n), y=A055997(n+1) and z = x^2+y, then x^4 + y^3 = z^2. - Bruno Berselli, Aug 24 2010
In general, if b(0)=1, b(1)=k and for n > 1, b(n) = 6*b(n-1) - b(n-2), then
for n > 0, b(n) = a(n)*k-a(n-1); e.g.,
for k=2, when b(n) = A038725(n), 2 = 1*2 - 0, 11 = 6*2 - 1, 64 = 35*2 - 6, 373 = 204*2 - 35;
for k=3, when b(n) = A001541(n), 3 = 1*3 - 0, 17 = 6*3 - 1; 99 = 35*3 - 6; 577 = 204*3 - 35;
for k=4, when b(n) = A038723(n), 4 = 1*4 - 0, 23 = 6*4 - 1; 134 = 35*4 - 6; 781 = 204*4 - 35;
for k=5, when b(n) = A001653(n), 5 = 1*5 - 0, 29 = 6*5 - 1; 169 = 35*5 - 6; 985 = 204*5 - 35.
- Charlie Marion, Dec 08 2010
See a Wolfdieter Lang comment on A001653 on a sequence of (u,v) values for Pythagorean triples (x,y,z) with x=|u^2-v^2|, y=2*u*v and z=u^2+v^2, with u odd and v even, generated from (u(0)=1,v(0)=2), the triple (3,4,5), by a substitution rule given there. The present a(n) appears there as b(n). The corresponding generated triangles have catheti differing by one length unit. - Wolfdieter Lang, Mar 06 2012
a(n)*a(n+2k) + a(k)^2 and a(n)*a(n+2k+1) + a(k)*a(k+1) are triangular numbers. Generalizes description of sequence. - Charlie Marion, Dec 03 2012
a(n)*a(n+2k) + a(k)^2 is the triangular square A001110(n+k). a(n)*a(n+2k+1) + a(k)*a(k+1) is the triangular oblong A029549(n+k). - Charlie Marion, Dec 05 2012
From Richard R. Forberg, Aug 30 2013: (Start)
The squares of a(n) are the result of applying triangular arithmetic to the squares, using A001333 as the "guide" on what integers to square, as follows:
a(2n)^2 = A001333(2n)^2 * (A001333(2n)^2 - 1)/2;
a(2n+1)^2 = A001333(2n+1)^2 * (A001333(2n+1)^2 + 1)/2. (End)
For n >= 1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,5}. - Milan Janjic, Jan 25 2015
Panda and Rout call these "balancing numbers" and note that the period of the sequence modulo a prime p is the same as that modulo p^2 when p = 13, 31, 1546463. But these are precisely the p in A238736 such that p^2 divides A000129(p - (2/p)), where (2/p) is a Jacobi symbol. In light of the above observation by Franklin T. Adams-Watters that the present sequence is one half the bisection of the Pell numbers, i.e., a(n) = A000129(2*n)/2, it follows immediately that modulo a fixed prime p, or any power thereof, the period of a(n) is half that of A000129(n). - John Blythe Dobson, Mar 06 2015
The triangular number = square number identity Tri((T(n, 3) - 1)/2) = S(n-1, 6)^2 with Tri, T, and S given in A000217, A053120 and A049310, is the special case k = 1 of the k-family of identities Tri((T(n, 2*k+1) - 1)/2) = Tri(k)*S(n-1, 2*(2*k+1))^2, k >= 0, n >= 0, with S(-1, x) = 0. For k=2 see A108741(n) for S(n-1, 10)^2. This identity boils down to the identities S(n-1, 2*x)^2 = (T(2*n, x) - 1)/(2*(x^2-1)) and 2*T(n, x)^2 - 1 = T(2*n, x) with x = 2*k+1. - Wolfdieter Lang, Feb 01 2016
a(2)=6 is perfect. For n=2*k, k > 0, k not equal to 1, a(n) is a multiple of a(2) and since every multiple (beyond 1) of a perfect number is abundant, then a(n) is abundant. sigma(a(4)) = 504 > 408 = 2*a(4). For n=2*k+1, k > 0, a(n) mod 10 = A000012(n), so a(n) is odd. If a(n) is a prime number, it is deficient; otherwise a(n) has one or two distinct prime factors and is therefore deficient again. So for n=2k+1, k > 0, a(n) is deficient. sigma(a(5)) = 1260 < 2378 = 2*a(5). - Muniru A Asiru, Apr 14 2016
Behera & Panda call these the balancing numbers, and A001541 are the balancers. - Michel Marcus, Nov 07 2017
In general, a second-order linear recurrence with constant coefficients having a signature of (c,d) will be duplicated by a third-order recurrence having a signature of (x,c^2-c*x+d,-d*x+c*d). The formulas of Olivares and Bouhamida in the formula section which have signatures of (7,-7,1) and (5,5,-1), respectively, are specific instances of this general rule for x = 7 and x = 5. - Gary Detlefs, Jan 29 2021
Note that 6 is the largest triangular number in the sequence, because it is proved that 8 and 9 are the largest perfect powers which are consecutive (Catalan's conjecture). 0 and 1 are also in the sequence because they are also perfect powers and 0*1/2 = 0^2 and 8*9/2 = (2*3)^2. - Metin Sariyar, Jul 15 2021

Examples

			G.f. = x + 6*x^2 + 35*x^3 + 204*x^4 + 1189*x^5 + 6930*x^6 + 40391*x^7 + ...
6 is in the sequence since 6^2 = 36 is a triangular number: 36 = 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8. - _Michael B. Porter_, Jul 02 2016
		

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) - From N. J. A. Sloane, May 30 2012
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, pp. 193, 197.
  • D. M. Burton, The History of Mathematics, McGraw Hill, (1991), p. 213.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 10.
  • P. Franklin, E. F. Beckenbach, H. S. M Coxeter, N. H. McCoy, K. Menger, and J. L. Synge, Rings And Ideals, No 8, The Carus Mathematical Monographs, The Mathematical Association of America, (1967), pp. 144-146.
  • A. Patra, G. K. Panda, and T. Khemaratchatakumthorn. "Exact divisibility by powers of the balancing and Lucas-balancing numbers." Fibonacci Quart., 59:1 (2021), 57-64; see B(n).
  • 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

Crossrefs

Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), this sequence (m=3), A001090 (m=4), A004189 (m=5), A004191 (m=6), A007655 (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.

Programs

  • GAP
    a:=[0,1];; for n in [3..25] do a[n]:=6*a[n-1]-a[n-2]; od; a; # Muniru A Asiru, Dec 18 2018
  • Haskell
    a001109 n = a001109_list !! n :: Integer
    a001109_list = 0 : 1 : zipWith (-)
       (map (* 6) $ tail a001109_list) a001109_list
    -- Reinhard Zumkeller, Dec 17 2011
    
  • Magma
    [n le 2 select n-1 else 6*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Jul 25 2015
    
  • Maple
    a[0]:=1: a[1]:=6: for n from 2 to 26 do a[n]:=6*a[n-1]-a[n-2] od: seq(a[n],n=0..26); # Emeric Deutsch
    with (combinat):seq(fibonacci(2*n,2)/2, n=0..20); # Zerinvary Lajos, Apr 20 2008
  • Mathematica
    Transpose[NestList[Flatten[{Rest[#],ListCorrelate[{-1,6},#]}]&, {0,1}, 30]][[1]]  (* Harvey P. Dale, Mar 23 2011 *)
    CoefficientList[Series[x/(1-6x+x^2),{x,0,30}],x]  (* Harvey P. Dale, Mar 23 2011 *)
    LinearRecurrence[{6, -1}, {0, 1}, 50] (* Vladimir Joseph Stephan Orlovsky, Feb 12 2012 *)
    a[ n_]:= ChebyshevU[n-1, 3]; (* Michael Somos, Sep 02 2012 *)
    Table[Fibonacci[2n, 2]/2, {n, 0, 20}] (* Vladimir Reshetnikov, Sep 16 2016 *)
    TrigExpand@Table[Sinh[2 n ArcCsch[1]]/(2 Sqrt[2]), {n, 0, 10}] (* Federico Provvedi, Feb 01 2021 *)
  • PARI
    {a(n) = imag((3 + quadgen(32))^n)}; /* Michael Somos, Apr 07 2003 */
    
  • PARI
    {a(n) = subst( poltchebi( abs(n+1)) - 3 * poltchebi( abs(n)), x, 3) / 8}; /* Michael Somos, Apr 07 2003 */
    
  • PARI
    {a(n) = polchebyshev( n-1, 2, 3)}; /* Michael Somos, Sep 02 2012 */
    
  • PARI
    is(n)=ispolygonal(n^2,3) \\ Charles R Greathouse IV, Nov 03 2016
    
  • Sage
    [lucas_number1(n,6,1) for n in range(27)] # Zerinvary Lajos, Jun 25 2008
    
  • Sage
    [chebyshev_U(n-1,3) for n in (0..20)] # G. C. Greubel, Dec 23 2019
    

Formula

G.f.: x / (1 - 6*x + x^2). - Simon Plouffe in his 1992 dissertation.
a(n) = S(n-1, 6) = U(n-1, 3) with U(n, x) Chebyshev's polynomials of the second kind. S(-1, x) := 0. Cf. triangle A049310 for S(n, x).
a(n) = sqrt(A001110(n)).
a(n) = A001542(n)/2.
a(n) = sqrt((A001541(n)^2-1)/8) (cf. Richardson comment).
a(n) = 3*a(n-1) + sqrt(8*a(n-1)^2+1). - R. J. Mathar, Oct 09 2000
a(n) = A000129(n)*A001333(n) = A000129(n)*(A000129(n)+A000129(n-1)) = ceiling(A001108(n)/sqrt(2)). - Henry Bottomley, Apr 19 2000
a(n) ~ (1/8)*sqrt(2)*(sqrt(2) + 1)^(2*n). - 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 05 2002
a(n) = ((3 + 2*sqrt(2))^n - (3 - 2*sqrt(2))^n) / (4*sqrt(2)). - Gregory V. Richardson, Oct 13 2002. Corrected for offset 0, and rewritten. - Wolfdieter Lang, Feb 10 2015
a(2*n) = a(n)*A003499(n). 4*a(n) = A005319(n). - Mario Catalani (mario.catalani(AT)unito.it), Mar 21 2003
a(n) = floor((3+2*sqrt(2))^n/(4*sqrt(2))). - Lekraj Beedassy, Apr 23 2003
a(-n) = -a(n). - Michael Somos, Apr 07 2003
For n >= 1, a(n) = Sum_{k=0..n-1} A001653(k). - Charlie Marion, Jul 01 2003
For n > 0, 4*a(2*n) = A001653(n)^2 - A001653(n-1)^2. - Charlie Marion, Jul 16 2003
For n > 0, a(n) = Sum_{k = 0..n-1}((2*k+1)*A001652(n-1-k)) + A000217(n). - Charlie Marion, Jul 18 2003
a(2*n+1) = a(n+1)^2 - a(n)^2. - Charlie Marion, Jan 12 2004
a(k)*a(2*n+k) = a(n+k)^2 - a(n)^2; e.g., 204*7997214 = 40391^2 - 35^2. - Charlie Marion, Jan 15 2004
For j < n+1, a(k+j)*a(2*n+k-j) - Sum_{i = 0..j-1} a(2*n-(2*i+1)) = a(n+k)^2 - a(n)^2. - Charlie Marion, Jan 18 2004
From Paul Barry, Feb 06 2004: (Start)
a(n) = A000129(2*n)/2;
a(n) = ((1+sqrt(2))^(2*n) - (1-sqrt(2))^(2*n))*sqrt(2)/8;
a(n) = Sum_{i=0..n} Sum_{j=0..n} A000129(i+j)*n!/(i!*j!*(n-i-j)!)/2. (End)
E.g.f.: exp(3*x)*sinh(2*sqrt(2)*x)/(2*sqrt(2)). - Paul Barry, Apr 21 2004
A053141(n+1) + A055997(n+1) = A001541(n+1) + a(n+1). - Creighton Dement, Sep 16 2004
a(n) = Sum_{k=0..n} binomial(2*n, 2*k+1)*2^(k-1). - Paul Barry, Oct 01 2004
a(n) = A001653(n+1) - A038723(n); (a(n)) = chuseq[J]( 'ii' + 'jj' + .5'kk' + 'ij' - 'ji' + 2.5e ), apart from initial term. - Creighton Dement, Nov 19 2004, modified by Davide Colazingari, Jun 24 2016
a(n+1) = Sum_{k=0..n} A001850(k)*A001850(n-k), self convolution of central Delannoy numbers. - Benoit Cloitre, Sep 28 2005
a(n) = 7*(a(n-1) - a(n-2)) + a(n-3), a(1) = 0, a(2) = 1, a(3) = 6, n > 3. Also a(n) = ( (1 + sqrt(2) )^(2*n) - (1 - sqrt(2) )^(2*n) ) / (4*sqrt(2)). - Antonio Alberto Olivares, Oct 23 2003
a(n) = 5*(a(n-1) + a(n-2)) - a(n-3). - Mohamed Bouhamida, Sep 20 2006
Define f(x,s) = s*x + sqrt((s^2-1)*x^2+1); f(0,s)=0. a(n) = f(a(n-1),3), see second formula. - Marcos Carreira, Dec 27 2006
The perfect median m(n) can be expressed in terms of the Pell numbers P() = A000129() by m(n) = P(n + 2) * (P(n + 2) + P(n + 1)) for n >= 0. - Winston A. Richards (ugu(AT)psu.edu), Jun 11 2007
For k = 0..n, a(2*n-k) - a(k) = 2*a(n-k)*A001541(n). Also, a(2*n+1-k) - a(k) = A002315(n-k)*A001653(n). - Charlie Marion, Jul 18 2007
[A001653(n), a(n)] = [1,4; 1,5]^n * [1,0]. - Gary W. Adamson, Mar 21 2008
a(n) = Sum_{k=0..n-1} 4^k*binomial(n+k,2*k+1). - Paul Barry, Apr 20 2009
a(n+1)^2 - 6*a(n+1)*a(n) + a(n)^2 = 1. - Charlie Marion, Dec 14 2010
a(n) = A002315(m)*A011900(n-m-1) + A001653(m)*A001652(n-m-1) - a(m) = A002315(m)*A053141(n-m-1) + A001653(m)*A046090(n-m-1) + a(m) with m < n; otherwise a(n) = A002315(m)*A053141(m-n) - A001653(m)*A011900(m-n) + a(m) = A002315(m)*A053141(m-n) - A001653(m)*A046090(m-n) - a(m) = (A002315(n) - A001653(n))/2. - Kenneth J Ramsey, Oct 12 2011
16*a(n)^2 + 1 = A056771(n). - James R. Buddenhagen, Dec 09 2011
A010054(A000290(a(n))) = 1. - Reinhard Zumkeller, Dec 17 2011
In general, a(n+k)^2 - A003499(k)*a(n+k)*a(n) + a(n)^2 = a(k)^2. - Charlie Marion, Jan 11 2012
a(n+1) = Sum_{k=0..n} A101950(n,k)*5^k. - Philippe Deléham, Feb 10 2012
PSUM transform of a(n+1) is A053142. PSUMSIGN transform of a(n+1) is A084158. BINOMIAL transform of a(n+1) is A164591. BINOMIAL transform of A086347 is a(n+1). BINOMIAL transform of A057087(n-1). - Michael Somos, May 11 2012
a(n+k) = A001541(k)*a(n) + sqrt(A132592(k)*a(n)^2 + a(k)^2). Generalizes formula dated Oct 09 2000. - Charlie Marion, Nov 27 2012
a(n) + a(n+2*k) = A003499(k)*a(n+k); a(n) + a(n+2*k+1) = A001653(k+1)*A002315(n+k). - Charlie Marion, Nov 29 2012
From Peter Bala, Dec 23 2012: (Start)
Product_{n >= 1} (1 + 1/a(n)) = 1 + sqrt(2).
Product_{n >= 2} (1 - 1/a(n)) = (1/3)*(1 + sqrt(2)). (End)
G.f.: G(0)*x/(2-6*x), where G(k) = 1 + 1/(1 - x*(8*k-9)/( x*(8*k-1) - 3/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 12 2013
G.f.: H(0)*x/2, where H(k) = 1 + 1/( 1 - x*(6-x)/(x*(6-x) + 1/H(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Feb 18 2014
a(n) = (a(n-1)^2 - a(n-3)^2)/a(n-2) + a(n-4) for n > 3. - Patrick J. McNab, Jul 24 2015
a(n-k)*a(n+k) + a(k)^2 = a(n)^2, a(n+k) + a(n-k) = A003499(k)*a(n), for n >= k >= 0. - Alexander Samokrutov, Sep 30 2015
Dirichlet g.f.: (PolyLog(s,3+2*sqrt(2)) - PolyLog(s,3-2*sqrt(2)))/(4*sqrt(2)). - Ilya Gutkovskiy, Jun 27 2016
4*a(n)^2 - 1 = A278310(n) for n > 0. - Bruno Berselli, Nov 24 2016
From Klaus Purath, Jan 18 2020: (Start)
a(n) = (a(n-3) + a(n+3))/198.
a(n) = Sum_{i=1..n} A001653(i), n>=1.
a(n) = sinh( 2 * n * arccsch(1) ) / ( 2 * sqrt(2) ). - Federico Provvedi, Feb 01 2021
(End)
a(n) = A002965(2*n)*A002965(2*n+1). - Jon E. Schoenfield, Jan 08 2022
a(n) = A002965(4*n)/2. - Gerry Martens, Jul 14 2023
a(n) = Sum_{k = 0..n-1} (-1)^(n+k+1)*binomial(n+k, 2*k+1)*8^k. - Peter Bala, Jul 17 2023

Extensions

Additional comments from Wolfdieter Lang, Feb 10 2000
Duplication of a formula removed by Wolfdieter Lang, Feb 10 2015

A019863 Decimal expansion of sin(3*Pi/10) (sine of 54 degrees, or cosine of 36 degrees).

Original entry on oeis.org

8, 0, 9, 0, 1, 6, 9, 9, 4, 3, 7, 4, 9, 4, 7, 4, 2, 4, 1, 0, 2, 2, 9, 3, 4, 1, 7, 1, 8, 2, 8, 1, 9, 0, 5, 8, 8, 6, 0, 1, 5, 4, 5, 8, 9, 9, 0, 2, 8, 8, 1, 4, 3, 1, 0, 6, 7, 7, 2, 4, 3, 1, 1, 3, 5, 2, 6, 3, 0, 2, 3, 1, 4, 0, 9, 4, 5, 1, 2, 2, 4, 8, 5, 3, 6, 0, 3, 6, 0, 2, 0, 9, 4, 6, 9, 5, 5, 6, 8
Offset: 0

Views

Author

Keywords

Comments

Midsphere radius of regular icosahedron with unit edges.
Also half of the golden ratio (A001622). - Stanislav Sykora, Jan 30 2014
Andris Ambainis (see Aaronson link) observes that combining the results of Barak-Hardt-Haviv-Rao with Dinur-Steurer yields the maximal probability of winning n parallel repetitions of a classical CHSH game (see A201488) asymptotic to this constant to the power of n, an improvement on the naive probability of (3/4)^n. (All the random bits are received upfront but the players cannot communicate or share an entangled state.) - Charles R Greathouse IV, May 15 2014
This is the height h of the isosceles triangle in a regular pentagon, in length units of the circumscribing radius, formed by a side as base and two adjacent radii. h = sin(3*Pi/10) = cos(Pi/5) (radius 1 unit). - Wolfdieter Lang, Jan 08 2018
Also the limiting value(L) of "r" which is abscissa of the vertex of the parabola F(n)*x^2 - F(n+1)*x + F(n + 2)(where F(n)=A000045(n) are the Fibonacci numbers and n>0). - Burak Muslu, Feb 24 2021

Examples

			0.80901699437494742410229341718281905886015458990288143106772431135263...
		

Crossrefs

Platonic solids midradii: A020765 (tetrahedron), A020761 (octahedron), A010503 (cube), A239798 (dodecahedron).

Programs

Formula

Equals (1+sqrt(5))/4 = cos(Pi/5) = sin(3*Pi/10). - R. J. Mathar, Jun 18 2006
Equals 2F1(4/5,1/5;1/2;3/4) / 2 = A019827 + 1/2. - R. J. Mathar, Oct 27 2008
Equals A001622 / 2. - Stanislav Sykora, Jan 30 2014
phi / 2 = (i^(2/5) + i^(-2/5)) / 2 = i^(2/5) - (sin(Pi/5))*i = i^(-2/5) + (sin(Pi/5))*i = i^(2/5) - (cos(3*Pi/10))*i = i^(-2/5) + (cos(3*Pi/10))*i. - Jaroslav Krizek, Feb 03 2014
Equals 1/A134972. - R. J. Mathar, Jan 17 2021
Equals 2*A019836*A019872. - R. J. Mathar, Jan 17 2021
Equals (A094214 + 1)/2 or 1/(2*A094214). - Burak Muslu, Feb 24 2021
Equals hypergeom([-2/5, -3/5], [6/5], -1) = hypergeom([-1/5, 3/5], [6/5], 1) = hypergeom([1/5, -3/5], [4/5], 1). - Peter Bala, Mar 04 2022
Equals Product_{k>=1} (1 - (-1)^k/A001611(k)). - Amiram Eldar, Nov 28 2024
Equals 2*A134944 = 3*A134946 = A187426-11/10 = A296182-1. - Hugo Pfoertner, Nov 28 2024
Equals A134945/4. Root of 4*x^2-2*x-1=0. - R. J. Mathar, Aug 29 2025

A020765 Decimal expansion of 1/sqrt(8).

Original entry on oeis.org

3, 5, 3, 5, 5, 3, 3, 9, 0, 5, 9, 3, 2, 7, 3, 7, 6, 2, 2, 0, 0, 4, 2, 2, 1, 8, 1, 0, 5, 2, 4, 2, 4, 5, 1, 9, 6, 4, 2, 4, 1, 7, 9, 6, 8, 8, 4, 4, 2, 3, 7, 0, 1, 8, 2, 9, 4, 1, 6, 9, 9, 3, 4, 4, 9, 7, 6, 8, 3, 1, 1, 9, 6, 1, 5, 5, 2, 6, 7, 5, 9, 7, 1, 2, 5, 9, 6, 8, 8, 3, 5, 8, 1, 9, 1, 0, 3, 9, 3
Offset: 0

Views

Author

Keywords

Comments

Multiplied by 10, this is the real and the imaginary part of sqrt(25i). - Alonso del Arte, Jan 11 2013
Radius of the midsphere (tangent to the edges) in a regular tetrahedron with unit edges. - Stanislav Sykora, Nov 20 2013
The side of the largest cubical present that can be wrapped (with cutting) by a unit square of wrapping paper. See Problem 10716 link. - Michel Marcus, Jul 24 2018
The ratio between the thickness and diameter of a geometrically fair coin having an equal probability, 1/3, of landing on each of its two faces and on its side after being tossed in the air. The calculation is based on comparing the areal projections of the faces and sides of the coin on a circumscribing sphere. (Mosteller, 1965). See A020760 for a physical solution. - Amiram Eldar, Sep 01 2020

Examples

			1/sqrt(8) = 0.353553390593273762200422181052424519642417968844237018294...
		

References

  • Frederick Mosteller, Fifty challenging problems of probability, Dover, New York, 1965. See problem 38, pp. 10 and 58-60.

Crossrefs

Cf. Midsphere radii in Platonic solids:
A020761 (octahedron),
A010503 (cube),
A019863 (icosahedron),
A239798 (dodecahedron).

Programs

Formula

A010503 divided by 2.
Equals A201488 minus 1/2. Equals 1/(A010487-4) minus 1/4. - Jon E. Schoenfield, Jan 09 2017
Equals Integral_{x=0..oo} x*exp(-x)*BesselJ(0,x) dx. - Kritsada Moomuang, Jun 03 2025

A384277 Decimal expansion of the smallest zero of the Laguerre polynomial of degree 3.

Original entry on oeis.org

4, 1, 5, 7, 7, 4, 5, 5, 6, 7, 8, 3, 4, 7, 9, 0, 8, 3, 3, 1, 1, 5, 3, 3, 8, 7, 3, 1, 2, 8, 2, 7, 4, 4, 7, 3, 5, 4, 6, 6, 1, 7, 4, 1, 2, 6, 9, 3, 1, 1, 8, 4, 6, 5, 0, 9, 3, 9, 6, 5, 9, 5, 4, 3, 2, 2, 3, 2, 5, 0, 1, 9, 9, 3, 6, 9, 1, 3, 3, 1, 4, 9, 5, 7, 1, 9, 6
Offset: 0

Views

Author

A.H.M. Smeets, May 24 2025

Keywords

Examples

			0.41577455678347908331153387312827447354661741269311...
		

Crossrefs

There are k positive real zeros of the Laguerre polynomial of degree k:
k | zeros | corresponding weights for Laguerre-Gauss quadrature
---+------------------------------------------+-----------------------------------------------------
3 | this sequence, A384278, A384279 | A384463, A384464, A384465

Programs

  • Mathematica
    First[RealDigits[Root[LaguerreL[3, #] &, 1], 10, 100]] (* Paolo Xausa, Jun 05 2025 *)

Formula

Smallest root of x^3 - 9 x^2 + 18 x - 6 = 0.

A384279 Decimal expansion of the largest zero of the Laguerre polynomial of degree 3.

Original entry on oeis.org

6, 2, 8, 9, 9, 4, 5, 0, 8, 2, 9, 3, 7, 4, 7, 9, 1, 9, 6, 8, 6, 6, 4, 1, 5, 7, 6, 5, 5, 1, 2, 1, 3, 1, 6, 5, 7, 4, 9, 3, 5, 2, 0, 8, 6, 6, 2, 4, 6, 6, 0, 0, 7, 0, 0, 8, 7, 0, 8, 3, 2, 7, 9, 7, 5, 9, 3, 6, 4, 4, 5, 2, 8, 7, 2, 5, 9, 2, 0, 2, 3, 8, 4, 7, 9, 6, 1
Offset: 1

Views

Author

A.H.M. Smeets, May 26 2025

Keywords

Examples

			6.28994508293747919686641576551213165749352086624660...
		

Crossrefs

Cf. A384590.
There are k positive real zeros of the Laguerre polynomial of degree k:
k | zeros | corresponding weights for Laguerre-Gauss quadrature
---+------------------------------------------+-----------------------------------------------------
3 | A384277, A384278, this sequence | A384463, A384464, A384465

Programs

  • Mathematica
    First[RealDigits[Root[LaguerreL[3, #] &, 3], 10, 100]] (* Paolo Xausa, Jun 05 2025 *)

Formula

largest root of x^3 - 9 x^2 + 18 x - 6 = 0.

A384278 Decimal expansion of the second smallest zero of the Laguerre polynomial of degree 3.

Original entry on oeis.org

2, 2, 9, 4, 2, 8, 0, 3, 6, 0, 2, 7, 9, 0, 4, 1, 7, 1, 9, 8, 2, 2, 0, 5, 0, 3, 6, 1, 3, 5, 9, 5, 9, 3, 8, 6, 8, 9, 5, 9, 8, 6, 1, 7, 2, 1, 0, 6, 0, 2, 8, 0, 8, 3, 4, 0, 3, 5, 2, 0, 1, 2, 4, 8, 0, 8, 4, 0, 3, 0, 4, 5, 1, 3, 3, 7, 1, 6, 6, 4, 4, 6, 5, 6, 3, 1, 8
Offset: 1

Views

Author

A.H.M. Smeets, May 24 2025

Keywords

Examples

			2.29428036027904171982205036135959386895986172106028...
		

Crossrefs

There are k positive real zeros of the Laguerre polynomial of degree k:
k | zeros | corresponding weights for Laguerre-Gauss quadrature
---+------------------------------------------+-----------------------------------------------------
3 | A384277, this sequence, A384279 | A384463, A384464, A384465

Programs

  • Mathematica
    First[RealDigits[Root[LaguerreL[3, #] &, 2], 10, 100]] (* Paolo Xausa, Jun 05 2025 *)

Formula

Second smallest root of x^3 - 9 x^2 + 18 x - 6 = 0.

A384280 Decimal expansion of the smallest zero of the Laguerre polynomial of degree 4.

Original entry on oeis.org

3, 2, 2, 5, 4, 7, 6, 8, 9, 6, 1, 9, 3, 9, 2, 3, 1, 1, 8, 0, 0, 3, 6, 1, 4, 5, 9, 1, 0, 4, 3, 6, 7, 4, 7, 9, 7, 4, 3, 7, 5, 7, 2, 2, 4, 4, 7, 4, 2, 9, 5, 7, 6, 7, 1, 8, 8, 4, 5, 1, 8, 5, 3, 8, 0, 6, 9, 6, 8, 6, 7, 8, 7, 0, 7, 7, 0, 4, 0, 0, 9, 8, 6, 8, 5, 8, 5
Offset: 0

Views

Author

A.H.M. Smeets, May 26 2025

Keywords

Examples

			0.32254768961939231180036145910436747974375722447429...
		

Crossrefs

There are k positive real zeros of the Laguerre polynomial of degree k:
k | zeros | corresponding weights for Laguerre-Gauss quadrature
---+------------------------------------------+-----------------------------------------------------
4 | this sequence, A384281

Programs

  • Mathematica
    First[RealDigits[Root[LaguerreL[4, #] &, 1], 10, 100]] (* Paolo Xausa, Jun 05 2025 *)

Formula

Smallest root of x^4 - 16 x^3 + 72 x^2 - 96 x + 24 = 0.

A384281 Decimal expansion of the second smallest zero of the Laguerre polynomial of degree 4.

Original entry on oeis.org

1, 7, 4, 5, 7, 6, 1, 1, 0, 1, 1, 5, 8, 3, 4, 6, 5, 7, 5, 6, 8, 6, 8, 1, 6, 7, 1, 2, 5, 1, 7, 9, 4, 7, 0, 2, 3, 6, 7, 3, 8, 7, 4, 5, 1, 5, 5, 3, 1, 0, 7, 2, 5, 0, 1, 7, 8, 2, 7, 8, 2, 6, 6, 0, 9, 9, 8, 4, 5, 6, 0, 5, 7, 4, 4, 2, 1, 9, 7, 1, 6, 4, 1, 4, 0, 1, 3
Offset: 1

Views

Author

A.H.M. Smeets, May 26 2025

Keywords

Examples

			1.74576110115834657568681671251794702367387451553107...
		

Crossrefs

There are k positive real zeros of the Laguerre polynomial of degree k:
k | zeros | corresponding weights for Laguerre-Gauss quadrature
---+------------------------------------------+-----------------------------------------------------
4 | A384280, this sequence

Programs

  • Mathematica
    First[RealDigits[Root[LaguerreL[4, #] &, 2], 10, 100]] (* Paolo Xausa, Jun 05 2025 *)

Formula

Second smallest root of x^4 - 16 x^3 + 72 x^2 - 96 x + 24 = 0.

A384463 Decimal expansion of the weight factor for Laguerre-Gauss quadrature corresponding to abscissa A384277.

Original entry on oeis.org

7, 1, 1, 0, 9, 3, 0, 0, 9, 9, 2, 9, 1, 7, 3, 0, 1, 5, 4, 4, 9, 5, 9, 0, 1, 9, 1, 1, 4, 2, 5, 9, 4, 4, 3, 1, 3, 0, 9, 3, 9, 3, 7, 9, 6, 2, 8, 9, 5, 5, 3, 4, 4, 5, 1, 3, 1, 7, 1, 7, 2, 4, 4, 3, 6, 1, 9, 0, 2, 1, 5, 5, 1, 2, 2, 1, 3, 2, 2, 3, 5, 8, 2, 0, 3, 7, 2
Offset: 0

Views

Author

A.H.M. Smeets, May 30 2025

Keywords

Examples

			0.71109300992917301544959019114259443130939379628955...
		

Crossrefs

There are k positive real zeros of the Laguerre polynomial of degree k:
k | zeros | corresponding weights for Laguerre-Gauss quadrature
---+------------------------------------------+-----------------------------------------------------
3 | A384277, A384278, A384279 | this sequence, A384464, A384465

Programs

  • Mathematica
    First[RealDigits[Root[1944*#^3 - 1944*#^2 + 405*# - 4 &, 3], 10, 100]] (* Paolo Xausa, Jun 26 2025 *)
  • PARI
    solve(x = 0.7, 0.8, 1944*x^3 - 1944*x^2 + 405*x - 4)

Formula

Largest root of 1944*x^3 - 1944*x^2 + 405*x - 4 = 0.

A384464 Decimal expansion of the weight factor for Laguerre-Gauss quadrature corresponding to abscissa A384278.

Original entry on oeis.org

2, 7, 8, 5, 1, 7, 7, 3, 3, 5, 6, 9, 2, 4, 0, 8, 4, 8, 8, 0, 1, 4, 4, 4, 8, 8, 8, 4, 5, 6, 7, 2, 6, 4, 8, 1, 0, 3, 4, 8, 9, 0, 0, 3, 0, 9, 8, 6, 3, 8, 8, 6, 7, 1, 8, 5, 6, 7, 3, 4, 9, 4, 8, 4, 3, 4, 4, 9, 4, 0, 9, 6, 5, 7, 9, 3, 6, 5, 7, 5, 3, 0, 3, 5, 7, 4, 2
Offset: 0

Views

Author

A.H.M. Smeets, May 30 2025

Keywords

Examples

			0.27851773356924084880144488845672648103489003098638...
		

Crossrefs

There are k positive real zeros of the Laguerre polynomial of degree k:
k | zeros | corresponding weights for Laguerre-Gauss quadrature
---+------------------------------------------+-----------------------------------------------------
3 | A384277, A384278, A384279 | A384463, this sequence, A384465

Programs

  • Mathematica
    First[RealDigits[Root[1944*#^3 - 1944*#^2 + 405*# - 4 &, 2], 10, 100]] (* Paolo Xausa, Jun 26 2025 *)
  • PARI
    solve(x = 0.2, 0.3, 1944*x^3 - 1944*x^2 + 405*x - 4)

Formula

Second largest root of 1944*x^3 - 1944*x^2 + 405*x - 4 = 0.
Showing 1-10 of 24 results. Next