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 21 results. Next

A257582 Lexicographically largest increasing sequence of primes for which the continued square root map (see A257574) produces Pi.

Original entry on oeis.org

5, 17, 37, 53, 131, 181, 263, 317, 859, 887, 1637, 2837, 3413, 5861, 6491, 10531, 13399, 14083, 14563, 21433, 29717, 30529, 31663, 31771, 32069, 32587, 36559, 36809, 39359, 39461, 45319, 46933, 49801, 52391, 52579, 52889, 55871, 57493, 59107, 59539, 64187, 64633, 75377, 77491, 82351, 86587
Offset: 1

Views

Author

N. J. A. Sloane, May 03 2015

Keywords

Comments

The continued square root map applied to a sequence (x,y,z,...) is CSR(x,y,z,...) = sqrt(x + sqrt(y + sqrt(z + ...))); this is well defined if the logarithm of the terms is O(2^n).

Crossrefs

Cf. A000796 (Pi), A257764 (analog for e = 2.71828... instead of Pi), A257809 (analog for delta = 4.6692...), A257574.

Programs

  • PARI
    (CSR(v, s)=forstep(i=#v, 1, -1, s=sqrt(v[i]+s)); s); a=[5]; for(n=1, 50, print1(a[#a]", "); for(i=primepi(a[#a])+1, oo, CSR(concat(a, vector(9, j, prime(i+j))))>=Pi && (a=concat(a, prime(i))) && break)) \\ The default precision of 38 digits yields correct terms only below 30000. To compute larger values correctly, realprecision must be increased. - M. F. Hasler, May 03 2018

Extensions

a(15)-a(46) from Chai Wah Wu, May 06 2015
Edited by M. F. Hasler, May 03 2018

A257764 Lexicographically largest increasing sequence of primes for which the continued square root map (see A257574) produces the decimal expansion of e (Euler's number).

Original entry on oeis.org

3, 13, 31, 59, 67, 103, 179, 193, 227, 229, 317, 983, 1201, 1213, 1321, 1787, 1811, 2179, 3571, 4817, 5333, 6803, 10433, 12197, 13063, 19391, 21283, 24571, 31817, 42307, 45377, 49957, 61909, 67933, 70573, 74843, 82421, 85909, 91099, 99241, 101293, 109639, 112087
Offset: 1

Views

Author

Chai Wah Wu, May 09 2015

Keywords

Comments

Similar to A257582, but converging to e.

Examples

			sqrt(3) =  1.7320508075688772...
sqrt(3+sqrt(13)) = 2.570126704165378...
sqrt(3+sqrt(13+sqrt(31))) = 2.703522309917472...
sqrt(3+sqrt(13+sqrt(31+sqrt(59)))) = 2.7173508299457327...
sqrt(3+sqrt(13+sqrt(31+sqrt(59+sqrt(67))))) = 2.718217091497069...
sqrt(3+sqrt(13+sqrt(31+sqrt(59+sqrt(67+sqrt(103)))))) = 2.7182780002752187...
		

Crossrefs

Cf. A001113 (e), A257582 (analog for Pi instead of e), A257809 (analog for delta = 4.6692...), A257574.

Programs

  • PARI
    (CSR(v, s)=forstep(i=#v, 1, -1, s=sqrt(v[i]+s)); s); a=[3]; for(n=1, 50, print1(a[#a]", "); for(i=primepi(a[#a])+1, oo, CSR(concat(a, vector(9, j, prime(i+j))))>=exp(1)&& (a=concat(a, prime(i)))&& break)) \\ The standard precision of 38 digits yields incorrect terms beyond 10433. Increase realprecision to compute larger values. - M. F. Hasler, May 03 2018

Extensions

Edited by M. F. Hasler, May 03 2018

A257858 An increasing sequence of integers for which the continued square root map (see A257574) produces the decimal expansion of Pi.

Original entry on oeis.org

6, 10, 19, 27, 29, 32, 42, 45, 56, 67, 75, 94, 109, 122, 138, 144, 151, 152, 172, 181, 194, 204, 205, 232, 256, 290, 316, 325, 346, 380, 412, 446, 478, 511, 520, 533, 580, 584, 617, 623, 654, 658, 661, 682, 734, 773, 823, 836, 865, 903, 954, 979, 997, 1008, 1059
Offset: 1

Views

Author

Chai Wah Wu, May 13 2015

Keywords

Examples

			sqrt(6) = 2.449489742783178
sqrt(6+sqrt(10)) = 3.0269254467476365
sqrt(6+sqrt(10+sqrt(19))) = 3.1287879095060176
sqrt(6+sqrt(10+sqrt(19+sqrt(27)))) = 3.140462825727146
sqrt(6+sqrt(10+sqrt(19+sqrt(27+sqrt(29))))) = 3.1414928066743406
		

Crossrefs

A001622 Decimal expansion of golden ratio phi (or tau) = (1 + sqrt(5))/2.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Also decimal expansion of the positive root of (x+1)^n - x^(2n). (x+1)^n - x^(2n) = 0 has only two real roots x1 = -(sqrt(5)-1)/2 and x2 = (sqrt(5)+1)/2 for all n > 0. - Cino Hilliard, May 27 2004
The golden ratio phi is the most irrational among irrational numbers; its successive continued fraction convergents F(n+1)/F(n) are the slowest to approximate to its actual value (I. Stewart, in "Nature's Numbers", Basic Books, 1997). - Lekraj Beedassy, Jan 21 2005
Let t=golden ratio. The lesser sqrt(5)-contraction rectangle has shape t-1, and the greater sqrt(5)-contraction rectangle has shape t. For definitions of shape and contraction rectangles, see A188739. - Clark Kimberling, Apr 16 2011
The golden ratio (often denoted by phi or tau) is the shape (i.e., length/width) of the golden rectangle, which has the special property that removal of a square from one end leaves a rectangle of the same shape as the original rectangle. Analogously, removals of certain isosceles triangles characterize side-golden and angle-golden triangles. Repeated removals in these configurations result in infinite partitions of golden rectangles and triangles into squares or isosceles triangles so as to match the continued fraction, [1,1,1,1,1,...] of tau. For the special shape of rectangle which partitions into golden rectangles so as to match the continued fraction [tau, tau, tau, ...], see A188635. For other rectangular shapes which depend on tau, see A189970, A190177, A190179, A180182. For triangular shapes which depend on tau, see A152149 and A188594; for tetrahedral, see A178988. - Clark Kimberling, May 06 2011
Given a pentagon ABCDE, 1/(phi)^2 <= (A*C^2 + C*E^2 + E*B^2 + B*D^2 + D*A^2) / (A*B^2 + B*C^2 + C*D^2 + D*E^2 + E*A^2) <= (phi)^2. - Seiichi Kirikami, Aug 18 2011
If a triangle has sides whose lengths form a geometric progression in the ratio of 1:r:r^2 then the triangle inequality condition requires that r be in the range 1/phi < r < phi. - Frank M Jackson, Oct 12 2011
The graphs of x-y=1 and x*y=1 meet at (tau,1/tau). - Clark Kimberling, Oct 19 2011
Also decimal expansion of the first root of x^sqrt(x+1) = sqrt(x+1)^x. - Michel Lagneau, Dec 02 2011
Also decimal expansion of the root of (1/x)^(1/sqrt(x+1)) = (1/sqrt(x+1))^(1/x). - Michel Lagneau, Apr 17 2012
This is the case n=5 of (Gamma(1/n)/Gamma(3/n))*(Gamma((n-1)/n)/Gamma((n-3)/n)): (1+sqrt(5))/2 = (Gamma(1/5)/Gamma(3/5))*(Gamma(4/5)/Gamma(2/5)). - Bruno Berselli, Dec 14 2012
Also decimal expansion of the only number x>1 such that (x^x)^(x^x) = (x^(x^x))^x = x^((x^x)^x). - Jaroslav Krizek, Feb 01 2014
For n >= 1, round(phi^prime(n)) == 1 (mod prime(n)) and, for n >= 3, round(phi^prime(n)) == 1 (mod 2*prime(n)). - Vladimir Shevelev, Mar 21 2014
The continuous radical sqrt(1+sqrt(1+sqrt(1+...))) tends to phi. - Giovanni Zedda, Jun 22 2019
Equals sqrt(2+sqrt(2-sqrt(2+sqrt(2-...)))). - Diego Rattaggi, Apr 17 2021
Given any complex p such that real(p) > -1, phi is the only real solution of the equation z^p+z^(p+1)=z^(p+2), and the only attractor of the complex mapping z->M(z,p), where M(z,p)=(z^p+z^(p+1))^(1/(p+2)), convergent from any complex plane point. - Stanislav Sykora, Oct 14 2021
The only positive number such that its decimal part, its integral part and the number itself (x-[x], [x] and x) form a geometric progression is phi, with respectively (phi -1, 1, phi) and a ratio = phi. This is the answer to the 4th problem of the 7th Canadian Mathematical Olympiad in 1975 (see IMO link and Doob reference). - Bernard Schott, Dec 08 2021
The golden ratio is the unique number x such that f(n*x)*c(n/x) - f(n/x)*c(n*x) = n for all n >= 1, where f = floor and c = ceiling. - Clark Kimberling, Jan 04 2022
In The Second Scientific American Book Of Mathematical Puzzles and Diversions, Martin Gardner wrote that, by 1910, Mark Barr (1871-1950) gave phi as a symbol for the golden ratio. - Bernard Schott, May 01 2022
Phi is the length of the equal legs of an isosceles triangle with side c = phi^2, and internal angles (A,B) = 36 degrees, C = 108 degrees. - Gary W. Adamson, Jun 20 2022
The positive solution to x^2 - x - 1 = 0. - Michal Paulovic, Jan 16 2023
The minimal polynomial of phi^n, for nonvanishing integer n, is P(n, x) = x^2 - L(n)*x + (-1)^n, with the Lucas numbers L = A000032, extended to negative arguments with L(n) = (-1)^n*L(n). P(0, x) = (x - 1)^2 is not minimal. - Wolfdieter Lang, Feb 20 2025
This is the largest real zero x of (x^4 + x^2 + 1)^2 = 2*(x^8 + x^4 + 1). - Thomas Ordowski, May 14 2025

Examples

			1.6180339887498948482045868343656381177203091798057628621...
		

References

  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 24, 112, 123, 184, 190, 203.
  • Michael Doob, The Canadian Mathematical Olympiad & L'Olympiade Mathématique du Canada 1969-1993 - Canadian Mathematical Society & Société Mathématique du Canada, Problem 4, 1975, pages 76-77, 1993.
  • Richard A. Dunlap, The Golden Ratio and Fibonacci Numbers, World Scientific, River Edge, NJ, 1997.
  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, Vol. 94, Cambridge University Press, 2003, Section 1.2.
  • Martin Gardner, The Second Scientific American Book Of Mathematical Puzzles and Diversions, "Phi: The Golden Ratio", Chapter 8, Simon & Schuster, NY, 1961.
  • Martin Gardner, Weird Water and Fuzzy Logic: More Notes of a Fringe Watcher, "The Cult of the Golden Ratio", Chapter 9, Prometheus Books, 1996, pages 90-97.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.5 The Fibonacci and Related Sequences, p. 287.
  • H. E. Huntley, The Divine Proportion, Dover, NY, 1970.
  • Mario Livio, The Golden Ratio, Broadway Books, NY, 2002. [see the review by G. Markowsky in the links field]
  • Gary B. Meisner, The Golden Ratio: The Divine Beauty of Mathematics, Race Point Publishing (The Quarto Group), 2018. German translation: Der Goldene Schnitt, Librero, 2023.
  • Scott Olsen, The Golden Section, Walker & Co., NY, 2006.
  • Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 137-139.
  • 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).
  • Hans Walser, The Golden Section, Math. Assoc. of Amer. Washington DC 2001.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See pp. 36-40.
  • Claude-Jacques Willard, Le nombre d'or, Magnard, Paris, 1987.

Crossrefs

Programs

  • Maple
    Digits:=1000; evalf((1+sqrt(5))/2); # Wesley Ivan Hurt, Nov 01 2013
  • Mathematica
    RealDigits[(1 + Sqrt[5])/2, 10, 130] (* Stefan Steinerberger, Apr 02 2006 *)
    RealDigits[ Exp[ ArcSinh[1/2]], 10, 111][[1]] (* Robert G. Wilson v, Mar 01 2008 *)
    RealDigits[GoldenRatio,10,120][[1]] (* Harvey P. Dale, Oct 28 2015 *)
  • PARI
    default(realprecision, 20080); x=(1+sqrt(5))/2; for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b001622.txt", n, " ", d));  \\ Harry J. Smith, Apr 19 2009
    
  • PARI
    /* Digit-by-digit method: write it as 0.5+sqrt(1.25) and start at hundredths digit */
    r=11; x=400; print(1); print(6);
    for(dig=1, 110, {d=0; while((20*r+d)*d <= x, d++);
    d--; /* while loop overshoots correct digit */
    print(d); x=100*(x-(20*r+d)*d); r=10*r+d})
    \\ Michael B. Porter, Oct 24 2009
    
  • PARI
    a(n) = floor(10^(n-1)*(quadgen(5))%10);
    alist(len) = digits(floor(quadgen(5)*10^(len-1))); \\ Chittaranjan Pardeshi, Jun 22 2022
    
  • Python
    from sympy import S
    def alst(n): # truncate extra last digit to avoid rounding
      return list(map(int, str(S.GoldenRatio.n(n+1)).replace(".", "")))[:-1]
    print(alst(105)) # Michael S. Branicky, Jan 06 2021

Formula

Equals Sum_{n>=2} 1/A064170(n) = 1/1 + 1/2 + 1/(2*5) + 1/(5*13) + 1/(13*34) + ... - Gary W. Adamson, Dec 15 2007
Equals Hypergeometric2F1([1/5, 4/5], [1/2], 3/4) = 2*cos((3/5)*arcsin(sqrt(3/4))). - Artur Jasinski, Oct 26 2008
From Hieronymus Fischer, Jan 02 2009: (Start)
The fractional part of phi^n equals phi^(-n), if n is odd. For even n, the fractional part of phi^n is equal to 1-phi^(-n).
General formula: Provided x>1 satisfies x-x^(-1)=floor(x), where x=phi for this sequence, then:
for odd n: x^n - x^(-n) = floor(x^n), hence fract(x^n) = x^(-n),
for even n: x^n + x^(-n) = ceiling(x^n), hence fract(x^n) = 1 - x^(-n),
for all n>0: x^n + (-x)^(-n) = round(x^n).
x=phi is the minimal solution to x - x^(-1) = floor(x) (where floor(x)=1 in this case).
Other examples of constants x satisfying the relation x - x^(-1) = floor(x) include A014176 (the silver ratio: where floor(x)=2) and A098316 (the "bronze" ratio: where floor(x)=3). (End)
Equals 2*cos(Pi/5) = e^(i*Pi/5) + e^(-i*Pi/5). - Eric Desbiaux, Mar 19 2010
The solutions to x-x^(-1)=floor(x) are determined by x=(1/2)*(m+sqrt(m^2+4)), m>=1; x=phi for m=1. In terms of continued fractions the solutions can be described by x=[m;m,m,m,...], where m=1 for x=phi, and m=2 for the silver ratio A014176, and m=3 for the bronze ratio A098316. - Hieronymus Fischer, Oct 20 2010
Sum_{n>=1} x^n/n^2 = Pi^2/10 - (log(2)*sin(Pi/10))^2 where x = 2*sin(Pi/10) = this constant here. [Jolley, eq 360d]
phi = 1 + Sum_{k>=1} (-1)^(k-1)/(F(k)*F(k+1)), where F(n) is the n-th Fibonacci number (A000045). Proof. By Catalan's identity, F^2(n) - F(n-1)*F(n+1) = (-1)^(n-1). Therefore,(-1)^(n-1)/(F(n)*F(n+1)) = F(n)/F(n+1) - F(n-1)/F(n). Thus Sum_{k=1..n} (-1)^(k-1)/(F(k)*F(k+1)) = F(n)/F(n+1). If n goes to infinity, this tends to 1/phi = phi - 1. - Vladimir Shevelev, Feb 22 2013
phi^n = (A000032(n) + A000045(n)*sqrt(5)) / 2. - Thomas Ordowski, Jun 09 2013
Let P(q) = Product_{k>=1} (1 + q^(2*k-1)) (the g.f. of A000700), then A001622 = exp(Pi/6) * P(exp(-5*Pi)) / P(exp(-Pi)). - Stephen Beathard, Oct 06 2013
phi = i^(2/5) + i^(-2/5) = ((i^(4/5))+1) / (i^(2/5)) = 2*(i^(2/5) - (sin(Pi/5))i) = 2*(i^(-2/5) + (sin(Pi/5))i). - Jaroslav Krizek, Feb 03 2014
phi = sqrt(2/(3 - sqrt(5))) = sqrt(2)/A094883. This follows from the fact that ((1 + sqrt(5))^2)*(3 - sqrt(5)) = 8, so that ((1 + sqrt(5))/2)^2 = 2/(3 - sqrt(5)). - Geoffrey Caveney, Apr 19 2014
exp(arcsinh(cos(Pi/2-log(phi)*i))) = exp(arcsinh(sin(log(phi)*i))) = (sqrt(3) + i) / 2. - Geoffrey Caveney, Apr 23 2014
exp(arcsinh(cos(Pi/3))) = phi. - Geoffrey Caveney, Apr 23 2014
cos(Pi/3) + sqrt(1 + cos(Pi/3)^2). - Geoffrey Caveney, Apr 23 2014
2*phi = z^0 + z^1 - z^2 - z^3 + z^4, where z = exp(2*Pi*i/5). See the Wikipedia Kronecker-Weber theorem link. - Jonathan Sondow, Apr 24 2014
phi = 1/2 + sqrt(1 + (1/2)^2). - Geoffrey Caveney, Apr 25 2014
Phi is the limiting value of the iteration of x -> sqrt(1+x) on initial value a >= -1. - Chayim Lowen, Aug 30 2015
From Isaac Saffold, Feb 28 2018: (Start)
1 = Sum_{k=0..n} binomial(n, k) / phi^(n+k) for all nonnegative integers n.
1 = Sum_{n>=1} 1 / phi^(2n-1).
1 = Sum_{n>=2} 1 / phi^n.
phi = Sum_{n>=1} 1/phi^n. (End)
From Christian Katzmann, Mar 19 2018: (Start)
phi = Sum_{n>=0} (15*(2*n)! + 8*n!^2)/(2*n!^2*3^(2*n+2)).
phi = 1/2 + Sum_{n>=0} 5*(2*n)!/(2*n!^2*3^(2*n+1)). (End)
phi = Product_{k>=1} (1 + 2/(-1 + 2^k*(sqrt(4+(1-2/2^k)^2) + sqrt(4+(1-1/2^k)^2)))). - Gleb Koloskov, Jul 14 2021
Equals Product_{k>=1} (Fibonacci(3*k)^2 + (-1)^(k+1))/(Fibonacci(3*k)^2 + (-1)^k) (Melham and Shannon, 1995). - Amiram Eldar, Jan 15 2022
From Michal Paulovic, Jan 16 2023: (Start)
Equals the real part of 2 * e^(i * Pi / 5).
Equals 2 * sin(3 * Pi / 10) = 2*A019863.
Equals -2 * sin(37 * Pi / 10).
Equals 1 + 1 / (1 + 1 / (1 + 1 / (1 + 1 / (1 + 1 / ...)))).
Equals (2 + 3 * (2 + 3 * (2 + 3 * ...)^(1/4))^(1/4))^(1/4).
Equals (1 + 2 * (1 + 2 * (1 + 2 * ...)^(1/3))^(1/3))^(1/3).
Equals (1 + phi + (1 + phi + (1 + phi + ...)^(1/3))^(1/3))^(1/3).
Equals 13/8 + Sum_{k=0..oo} (-1)^(k+1)*(2*k+1)!/((k+2)!*k!*4^(2*k+3)).
(End)
phi^n = phi * A000045(n) + A000045(n-1). - Gary W. Adamson, Sep 09 2023
The previous formula holds for integer n, with F(-n) = (-1)^(n+1)*F(n), for n >= 0, with F(n) = A000045(n), for n >= 0. phi^n are integers in the quadratic number field Q(sqrt(5)). - Wolfdieter Lang, Sep 16 2023
Equals Product_{k>=0} ((5*k + 2)*(5*k + 3))/((5*k + 1)*(5*k + 4)). - Antonio Graciá Llorente, Feb 24 2024
From Antonio Graciá Llorente, Apr 21 2024: (Start)
Equals Product_{k>=1} phi^(-2^k) + 1, with phi = A001622.
Equals Product_{k>=0} ((5^(k+1) + 1)*(5^(k-1/2) + 1))/((5^k + 1)*(5^(k+1/2) + 1)).
Equals Product_{k>=1} 1 - (4*(-1)^k)/(10*k - 5 + (-1)^k) = Product_{k>=1} A047221(k)/A047209(k).
Equals Product_{k>=0} ((5*k + 7)*(5*k + 1 + (-1)^k))/((5*k + 1)*(5*k + 7 + (-1)^k)).
Equals Product_{k>=0} ((10*k + 3)*(10*k + 5)*(10*k + 8)^2)/((10*k + 2)*(10*k + 4)*(10*k + 9)^2).
Equals Product_{k>=5} 1 + 1/(Fibonacci(k) - (-1)^k).
Equals Product_{k>=2} 1 + 1/Fibonacci(2*k).
Equals Product_{k>=2} (Lucas(k)^2 + (-1)^k)/(Lucas(k)^2 - 4*(-1)^k). (End)

Extensions

Additional links contributed by Lekraj Beedassy, Dec 23 2003
More terms from Gabriel Cunningham (gcasey(AT)mit.edu), Oct 24 2004
More terms from Stefan Steinerberger, Apr 02 2006
Broken URL to Project Gutenberg replaced by Georg Fischer, Jan 03 2009
Edited by M. F. Hasler, Feb 24 2014

A028387 a(n) = n + (n+1)^2.

Original entry on oeis.org

1, 5, 11, 19, 29, 41, 55, 71, 89, 109, 131, 155, 181, 209, 239, 271, 305, 341, 379, 419, 461, 505, 551, 599, 649, 701, 755, 811, 869, 929, 991, 1055, 1121, 1189, 1259, 1331, 1405, 1481, 1559, 1639, 1721, 1805, 1891, 1979, 2069, 2161, 2255, 2351, 2449, 2549, 2651
Offset: 0

Views

Author

Keywords

Comments

a(n+1) is the least k > a(n) + 1 such that A000217(a(n)) + A000217(k) is a square. - David Wasserman, Jun 30 2005
Values of Fibonacci polynomial n^2 - n - 1 for n = 2, 3, 4, 5, ... - Artur Jasinski, Nov 19 2006
A127701 * [1, 2, 3, ...]. - Gary W. Adamson, Jan 24 2007
Row sums of triangle A135223. - Gary W. Adamson, Nov 23 2007
Equals row sums of triangle A143596. - Gary W. Adamson, Aug 26 2008
a(n-1) gives the number of n X k rectangles on an n X n chessboard (for k = 1, 2, 3, ..., n). - Aaron Dunigan AtLee, Feb 13 2009
sqrt(a(0) + sqrt(a(1) + sqrt(a(2) + sqrt(a(3) + ...)))) = sqrt(1 + sqrt(5 + sqrt(11 + sqrt(19 + ...)))) = 2. - Miklos Kristof, Dec 24 2009
When n + 1 is prime, a(n) gives the number of irreducible representations of any nonabelian group of order (n+1)^3. - Andrew Rupinski, Mar 17 2010
a(n) = A176271(n+1, n+1). - Reinhard Zumkeller, Apr 13 2010
The product of any 4 consecutive integers plus 1 is a square (see A062938); the terms of this sequence are the square roots. - Harvey P. Dale, Oct 19 2011
Or numbers not expressed in the form m + floor(sqrt(m)) with integer m. - Vladimir Shevelev, Apr 09 2012
Left edge of the triangle in A214604: a(n) = A214604(n+1,1). - Reinhard Zumkeller, Jul 25 2012
Another expression involving phi = (1 + sqrt(5))/2 is a(n) = (n + phi)(n + 1 - phi). Therefore the numbers in this sequence, even if they are prime in Z, are not prime in Z[phi]. - Alonso del Arte, Aug 03 2013
a(n-1) = n*(n+1) - 1, n>=0, with a(-1) = -1, gives the values for a*c of indefinite binary quadratic forms [a, b, c] of discriminant D = 5 for b = 2*n+1. In general D = b^2 - 4ac > 0 and the form [a, b, c] is a*x^2 + b*x*y + c*y^2. - Wolfdieter Lang, Aug 15 2013
a(n) has prime factors given by A038872. - Richard R. Forberg, Dec 10 2014
A253607(a(n)) = -1. - Reinhard Zumkeller, Jan 05 2015
An example of a quadratic sequence for which the continued square root map (see A257574) produces the number 2. There are infinitely many sequences with this property - another example is A028387. See Popular Computing link. - N. J. A. Sloane, May 03 2015
Left edge of the triangle in A260910: a(n) = A260910(n+2,1). - Reinhard Zumkeller, Aug 04 2015
Numbers m such that 4m+5 is a square. - Bruce J. Nicholson, Jul 19 2017
The numbers represented as 131 in base n: 131_4 = 29, 131_5 = 41, ... . If 'digits' larger than the base are allowed then 131_2 = 11 and 131_1 = 5 also. - Ron Knott, Nov 14 2017
From Klaus Purath, Mar 18 2019: (Start)
Let m be a(n) or a prime factor of a(n). Then, except for 1 and 5, there are, if m is a prime, exactly two squares y^2 such that the difference y^2 - m contains exactly one pair of factors {x,z} such that the following applies: x*z = y^2 - m, x + y = z with
x < y, where {x,y,z} are relatively prime numbers. {x,y,z} are the initial values of a sequence of the Fibonacci type. Thus each a(n) > 5, if it is a prime, and each prime factor p > 5 of an a(n) can be assigned to exactly two sequences of the Fibonacci type. a(0) = 1 belongs to the original Fibonacci sequence and a(1) = 5 to the Lucas sequence.
But also the reverse assignment applies. From any sequence (f(i)) of the Fibonacci type we get from its 3 initial values by f(i)^2 - f(i-1)*f(i+1) with f(i-1) < f(i) a term a(n) or a prime factor p of a(n). This relation is also valid for any i. In this case we get the absolute value |a(n)| or |p|. (End)
a(n-1) = 2*T(n) - 1, for n>=1, with T = A000217, is a proper subsequence of A089270, and the terms are 0,-1,+1 (mod 5). - Wolfdieter Lang, Jul 05 2019
a(n+1) is the number of wedged n-dimensional spheres in the homotopy of the neighborhood complex of Kneser graph KG_{2,n}. Here, KG_{2,n} is a graph whose vertex set is the collection of subsets of cardinality 2 of set {1,2,...,n+3,n+4} and two vertices are adjacent if and only if they are disjoint. - Anurag Singh, Mar 22 2021
Also the number of squares between (n+2)^2 and (n+2)^4. - Karl-Heinz Hofmann, Dec 07 2021
(x, y, z) = (A001105(n+1), -a(n-1), -a(n)) are solutions of the Diophantine equation x^3 + 4*y^3 + 4*z^3 = 8. - XU Pingya, Apr 25 2022
The least significant digit of terms of this sequence cycles through 1, 5, 1, 9, 9. - Torlach Rush, Jun 05 2024

Examples

			From _Ilya Gutkovskiy_, Apr 13 2016: (Start)
Illustration of initial terms:
                                        o               o
                        o           o   o o           o o
            o       o   o o       o o   o o o       o o o
    o   o   o o   o o   o o o   o o o   o o o o   o o o o
o   o o o   o o o o o   o o o o o o o   o o o o o o o o o
n=0  n=1       n=2           n=3               n=4
(End)
From _Klaus Purath_, Mar 18 2019: (Start)
Examples:
a(0) = 1: 1^1-0*1 = 1, 0+1 = 1 (Fibonacci A000045).
a(1) = 5: 3^2-1*4 = 5, 1+3 = 4 (Lucas A000032).
a(2) = 11: 4^2-1*5 = 11, 1+4 = 5 (A000285); 5^2-2*7 = 11, 2+5 = 7 (A001060).
a(3) = 19: 5^2-1*6 = 19, 1+5 = 6 (A022095); 7^2-3*10 = 19, 3+7 = 10 (A022120).
a(4) = 29: 6^2-1*7 = 29, 1+6 = 7 (A022096); 9^2-4*13 = 29, 4+9 = 13 (A022130).
a(11)/5 = 31: 7^2-2*9 = 31, 2+7 = 9 (A022113); 8^2-3*11 = 31, 3+8 = 11 (A022121).
a(24)/11 = 59: 9^2-2*11 = 59, 2+9 = 11 (A022114); 12^2-5*17 = 59, 5+12 = 17 (A022137).
(End)
		

Crossrefs

Complement of A028392. Third column of array A094954.
Cf. A000217, A002522, A062392, A062786, A127701, A135223, A143596, A052905, A162997, A062938 (squares of this sequence).
A110331 and A165900 are signed versions.
Cf. A002327 (primes), A094210.
Frobenius number for k successive numbers: this sequence (k=2), A079326 (k=3), A138984 (k=4), A138985 (k=5), A138986 (k=6), A138987 (k=7), A138988 (k=8).

Programs

Formula

a(n) = sqrt(A062938(n)). - Floor van Lamoen, Oct 08 2001
a(0) = 1, a(1) = 5, a(n) = (n+1)*a(n-1) - (n+2)*a(n-2) for n > 1. - Gerald McGarvey, Sep 24 2004
a(n) = A105728(n+2, n+1). - Reinhard Zumkeller, Apr 18 2005
a(n) = A109128(n+2, 2). - Reinhard Zumkeller, Jun 20 2005
a(n) = 2*T(n+1) - 1, where T(n) = A000217(n). - Gary W. Adamson, Aug 15 2007
a(n) = A005408(n) + A002378(n); A084990(n+1) = Sum_{k=0..n} a(k). - Reinhard Zumkeller, Aug 20 2007
Binomial transform of [1, 4, 2, 0, 0, 0, ...] = (1, 5, 11, 19, ...). - Gary W. Adamson, Sep 20 2007
G.f.: (1+2*x-x^2)/(1-x)^3. a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - R. J. Mathar, Jul 11 2009
a(n) = (n + 2 + 1/phi) * (n + 2 - phi); where phi = 1.618033989... Example: a(3) = 19 = (5 + .6180339...) * (3.381966...). Cf. next to leftmost column in A162997 array. - Gary W. Adamson, Jul 23 2009
a(n) = a(n-1) + 2*(n+1), with n > 0, a(0) = 1. - Vincenzo Librandi, Nov 18 2010
For k < n, a(n) = (k+1)*a(n-k) - k*a(n-k-1) + k*(k+1); e.g., a(5) = 41 = 4*11 - 3*5 + 3*4. - Charlie Marion, Jan 13 2011
a(n) = lower right term in M^2, M = the 2 X 2 matrix [1, n; 1, (n+1)]. - Gary W. Adamson, Jun 29 2011
G.f.: (x^2-2*x-1)/(x-1)^3 = G(0) where G(k) = 1 + x*(k+1)*(k+4)/(1 - 1/(1 + (k+1)*(k+4)/G(k+1))); (continued fraction, 3-step). - Sergei N. Gladkovskii, Oct 16 2012
Sum_{n>0} 1/a(n) = 1 + Pi*tan(sqrt(5)*Pi/2)/sqrt(5). - Enrique Pérez Herrero, Oct 11 2013
E.g.f.: exp(x) (1+4*x+x^2). - Tom Copeland, Dec 02 2013
a(n) = A005408(A000217(n)). - Tony Foster III, May 31 2016
From Amiram Eldar, Jan 29 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = -Pi*sec(sqrt(5)*Pi/2).
Product_{n>=1} (1 - 1/a(n)) = -Pi*sec(sqrt(5)*Pi/2)/6. (End)
a(5*n+1)/5 = A062786(n+1). - Torlach Rush, Jun 05 2024

Extensions

Minor edits by N. J. A. Sloane, Jul 04 2010, following suggestions from the Sequence Fans Mailing List

A072449 Decimal expansion of the limit of the nested radical sqrt(1 + sqrt(2 + sqrt(3 + sqrt(4 + ... )))).

Original entry on oeis.org

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

Views

Author

Robert G. Wilson v, Aug 01 2002

Keywords

Comments

Herschfeld calls this the Kasner number, after Edward Kasner. - Charles R Greathouse IV, Dec 30 2008
No closed-form expression is known for this constant.
"It was discovered by T. Vijayaraghavan that the infinite radical sqrt( a_1 + sqrt( a_2 + sqrt ( a_3 + sqrt( a_4 + ...)))), where a_n >= 0, will converge to a limit if and only if the limit of (log a_n)/2^n exists" - Clawson, p. 229. Obviously if a_n = n, the limit of (log a_n) / 2^n as n -> infinity is 0.
The continued fraction is A072450.
Clawson misstates Vijayaraghavan's theorem. Vijayaraghavan proved that for a_n > 0, the infinite radical sqrt(a_1 + sqrt(a_2 + sqrt(a_3 + ...))) converges if and only if limsup (log a_n)/2^n < infinity. (For example, suppose a_n = 1 if n is odd, and a_n = e^2^n if n is even. Then (log a_n)/2^n = 0, 1, 0, 1, 0, 1, ... for n >= 1, so the limit does not exist. However, limsup (log a_n)/2^n = 1 and the infinite radical converges.) - Jonathan Sondow, Mar 25 2014

Examples

			1.757932756618004532708819638218138527653...
		

References

  • Calvin C. Clawson, "Mathematical Mysteries, the beauty and magic of numbers," Perseus Books, Cambridge, Mass., 1996, pages 142 & 229.
  • S. R. Finch, Mathematical Constants, Cambridge, 2003, Section 1.2.1.
  • David Wells, "The Penguin Dictionary of Curious and Interesting Numbers," Revised Edition, London, England, 1997, page 30.
  • Stephen Wolfram, "A New Kind Of Science," Wolfram Media, 2002, page 915.

Crossrefs

Programs

  • Mathematica
    RealDigits[ Fold[ Sqrt[ #1 + #2] &, 0, Reverse[ Range[100]]], 10, 111][[1]] (* A New Kind Of Science *)
  • PARI
    s=200; for(n=1,199,t=200-n+sqrt(s); s=t);sqrt(s) \\ gives at least 180 correct digits

A105546 Decimal expansion of prime nested radical.

Original entry on oeis.org

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

Views

Author

Jonathan Vos Post, Apr 12 2005

Keywords

Comments

sqrt(1 + sqrt(2 + sqrt(3 + sqrt(4 + ...)))) = 1.75793275661800...
"It was discovered by T. Vijayaraghavan that the infinite radical, sqrt( a_1 + sqrt( a_2 + sqrt ( a_3 + sqrt( a_4 + ... )))) where a_n >= 0, will converge to a limit if and only if the limit of log(a_n)/2^n exists." [Clawson, 229; cf. A072449].
We know the asymptotic limit of primes and hence that the Prime Nested Radical converges.
Clawson misstates Vijayaraghavan's theorem. Vijayaraghavan proved that for a_n > 0, the infinite radical sqrt(a_1 + sqrt(a_2 + sqrt(a_3 + ...))) converges if and only if limsup (log a_n)/2^n < infinity. (For example, suppose a_n = 1 if n is odd, and a_n = e^2^n if n is even. Then (log a_n)/2^n = 0, 1, 0, 1, 0, 1, ... for n >= 1, so the limit does not exist. However, limsup (log a_n)/2^n = 1 and the infinite radical converges.) - Jonathan Sondow, Mar 25 2014

Examples

			2.10359749633989726261993964968532544404216228824001387298687284563...
		

References

  • Calvin C. Clawson, "Mathematical Mysteries, the beauty and magic of numbers," Perseus Books, Cambridge, Mass., 1996, pages 142 and 229.
  • S. R. Finch, Analysis of a Radical Expansion, Section 1.2.1 in Mathematical Constants. Cambridge, England: Cambridge University Press, p. 8, 2003.

Crossrefs

A105548 is the continued fraction representation of this prime nested radical.
A105815 is the similar semiprime nested radical.
A105817 is the Fibonacci nested radical.

Programs

  • Mathematica
    RealDigits[Fold[Sqrt[#1 + #2] &, 0, Reverse[Prime[Range[ 80]]]], 10, 111][[1]] (* Robert G. Wilson v, May 31 2005 *)

Formula

sqrt(2 + sqrt(3 + sqrt(5 + sqrt(7 + sqrt(11 + ... + sqrt(prime(n) + ...)))).

Extensions

Crossrefs corrected by Jaroslav Krizek, Jan 03 2015

A099879 Decimal expansion of a nested radical: sqrt(1^2 + sqrt(2^2 + sqrt(3^2 + ...

Original entry on oeis.org

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

Views

Author

Mark Hudson (mrmarkhudson(AT)hotmail.com), Nov 03 2004

Keywords

Examples

			1.94265542276398732822141329141266723768807363...
		

Crossrefs

Cf. A072449, A099874, A099876 to A099878 for other nested radicals.

Programs

  • Mathematica
    k = 64; r = 65; While[k > 0, r = Sqrt[k^2 + r]; k-- ]; RealDigits[r, 10, 111][[1]] (* Robert G. Wilson v, Nov 04 2004 *)
  • PARI
    t=0; forstep(n=100,1,-1,t=sqrt(t+n^2)); print(t)
    
  • PARI
    \\ We need about b/log(b) steps, where epsilon = 2^-b.
    my(b=bitprecision(1.),t); forstep(n=b\log(b)+9,1,-1, t=sqrt(t+n^2)); t \\ Charles R Greathouse IV, Aug 19 2025

A105817 Decimal expansion of the Fibonacci nested radical.

Original entry on oeis.org

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

Views

Author

Jonathan Vos Post, Apr 21 2005

Keywords

Comments

The continued fraction expression of this is A105818. "It was discovered by T. Vijayaraghavan that the infinite radical, sqrt( a_1 + sqrt( a_2 + sqrt ( a_3 + sqrt( a_4 + ... where a_n => 0, will converge to a limit if and only if the limit of (ln a_n)/2^n exists." [Clawson, 229; Sloane]. We know the asymptotic limit of Fibonacci numbers is Phi^n (Binet expansion) and that Phi^n < 2^n and hence that the Fibonacci Nested Radical converges.
Clawson misstates Vijayaraghavan's theorem. Vijayaraghavan proved that for a_n > 0, the infinite radical sqrt(a_1 + sqrt(a_2 + sqrt(a_3 + ...))) converges if and only if limsup (log a_n)/2^n < infinity. (For example, suppose a_n = 1 if n is odd, and a_n = e^2^n if n is even. Then (log a_n)/2^n = 0, 1, 0, 1, 0, 1, ... for n >= 1, so the limit does not exist. However, limsup (log a_n)/2^n = 1 and the infinite radical converges.) - Jonathan Sondow, Mar 25 2014

Examples

			1.66198246232781155796760608181513129505616756246503500829906806743...
		

References

  • Calvin C. Clawson, "Mathematical Mysteries, the beauty and magic of numbers," Perseus Books, Cambridge, Mass., 1996, pages 142 & 229.
  • S. R. Finch, "Analysis of a Radical Expansion." Section 1.2.1 in Mathematical Constants. Cambridge, England: Cambridge University Press, p. 8, 2003.

Crossrefs

Programs

  • Mathematica
    RealDigits[ Fold[ Sqrt[ #1 + #2] &, 0, Reverse[ Fibonacci[ Range[50]]]], 10, 111][[1]] (* Robert G. Wilson v, Apr 21 2005 *)

Formula

Sqrt(1 + sqrt(1 + sqrt(2 + sqrt(3 + sqrt(5 + ... + sqrt(Fibonacci(n)=A000045)))).

A257352 G.f.: (1-2*x+51*x^2)/(1-x)^3.

Original entry on oeis.org

1, 1, 51, 151, 301, 501, 751, 1051, 1401, 1801, 2251, 2751, 3301, 3901, 4551, 5251, 6001, 6801, 7651, 8551, 9501, 10501, 11551, 12651, 13801, 15001, 16251, 17551, 18901, 20301, 21751, 23251, 24801, 26401, 28051, 29751, 31501, 33301
Offset: 0

Views

Author

N. J. A. Sloane, May 03 2015

Keywords

Comments

An example of a quadratic sequence for which the continued square root map (see A257574) produces the number 2. There are infinitely many sequences with this property - another example is A028387.

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[(1-2x+51x^2)/(1-x)^3,{x,0,40}],x] (* or  *) LinearRecurrence[{3,-3,1},{1,1,51},40] (* Harvey P. Dale, Sep 01 2025 *)
  • PARI
    a(n)=25*n^2-25*n+1 \\ Charles R Greathouse IV, Jun 17 2017
Showing 1-10 of 21 results. Next