A006721 -id:A006721 - cp's OEIS frontend

A038754 a(2n) = 3^n, a(2n+1) = 2*3^n.

1, 2, 3, 6, 9, 18, 27, 54, 81, 162, 243, 486, 729, 1458, 2187, 4374, 6561, 13122, 19683, 39366, 59049, 118098, 177147, 354294, 531441, 1062882, 1594323, 3188646, 4782969, 9565938, 14348907, 28697814, 43046721, 86093442, 129140163, 258280326, 387420489

Offset: 0

Henry Bottomley, May 03 2000

In general, for the recurrence a(n) = a(n-1)*a(n-2)/a(n-3), all terms are integers iff a(0) divides a(2) and first three terms are positive integers, since a(2n+k) = a(k)*(a(2)/a(0))^n for all nonnegative integers n and k.
Equals eigensequence of triangle A070909; (1, 1, 2, 3, 6, 9, 18, ...) shifts to the left with multiplication by triangle A070909. - Gary W. Adamson, May 15 2010
The a(n) represent all paths of length (n+1), n >= 0, starting at the initial node on the path graph P_5, see the second Maple program. - Johannes W. Meijer, May 29 2010
a(n) is the difference between numbers of multiple of 3 evil (A001969) and odious (A000069) numbers in interval [0, 2^(n+1)). - Vladimir Shevelev, May 16 2012
A "half-geometric progression": to obtain a term (beginning with the third one) we multiply the before previous one by 3. - Vladimir Shevelev, May 21 2012
Pisano periods: 1, 2, 1, 4, 8, 2, 12, 4, 1, 8, 10, 4, 6, 12, 8, 8, 32, 2, 36, 8, ... . - R. J. Mathar, Aug 10 2012
Numbers k such that the k-th cyclotomic polynomial has a root mod 3. - Eric M. Schmidt, Jul 31 2013
Range of row n of the circular Pascal array of order 6. - Shaun V. Ault, Jun 05 2014
Also, the number of walks of length n on the graph 0--1--2--3--4 starting at vertex 1. - Sean A. Irvine, Jun 03 2025

			In the interval [0,2^5) we have 11 multiples of 3 numbers, from which 10 are evil and only one (21) is odious. Thus a(4) = 10 - 1 = 9. - _Vladimir Shevelev_, May 16 2012

Indranil Ghosh, Table of n, a(n) for n = 0..1500 (first 401 terms from T. D. Noe)
S. V. Ault and C. Kicey, Counting paths in corridors using circular Pascal arrays, Discrete Mathematics, Vol. 332 (2014), pp. 45-54.
Paul Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.
Nachum Dershowitz, Between Broadway and the Hudson, arXiv:2006.06516 [math.CO], 2020.
Sean A. Irvine, Walks on Graphs.
Richard L. Ollerton and Anthony G. Shannon, Some properties of generalized Pascal squares and triangles, Fib. Q., 36 (1998), 98-109. See pages 106-7.
Vladimir Shevelev, On monotonic strengthening of Newman-like phenomenon on (2m+1)-multiples in base 2m, arXiv:0710.3177 [math.NT], 2007.
M. B. Wells, Elements of Combinatorial Computing, Pergamon, Oxford, 1971. [Annotated scanned copy of pages 237-240]
Index entries for linear recurrences with constant coefficients, signature (0,3).

Cf. A000045, A000069, A000079, A000120, A000244, A001969, A006720, A006721.
Cf. A006722, A006723, A008776, A028495, A030436, A037124, A048328, A061551.
Cf. A068911, A070909, A087503, A124302, A128588, A133464, A140730, A140740.
Cf. A178381, A182751, A182752, A182753, A182754, A182755, A182756, A182757.
Fifth row of the array A094718.

import Data.List (transpose)
a038754 n = a038754_list !! n
a038754_list = concat $ transpose [a000244_list, a008776_list]
-- Reinhard Zumkeller, Oct 19 2015

[n le 2 select n else 3*Self(n-2): n in [1..40]]; // Vincenzo Librandi, Aug 18 2016

a[0]:=0:a[1]:=1:for n from 2 to 50 do a[n]:=3*a[n-2]+2 od: seq(a[n]+1, n=0..34); # Zerinvary Lajos, Mar 20 2008
with(GraphTheory): P:=5: G:=PathGraph(P): A:= AdjacencyMatrix(G): nmax:=35; for n from 1 to nmax do B(n):=A^n; a(n):=add(B(n)[1,k],k=1..P) od: seq(a(n),n=1..nmax); # Johannes W. Meijer, May 29 2010

LinearRecurrence[{0,3},{1,2},40] (* Harvey P. Dale, Jan 26 2014 *)
CoefficientList[Series[(1+2x)/(1-3x^2), {x, 0, 40}], x] (* Vincenzo Librandi, Aug 18 2016 *)
Module[{nn=20,c},c=3^Range[0,nn];Riffle[c,2c]] (* Harvey P. Dale, Aug 21 2021 *)

PARI
```
a(n)=(1/6)*(5-(-1)^n)*3^floor(n/2)
```
PARI
```
a(n)=3^(n>>1)<
				
```

[2^(n%2)*3^((n-(n%2))/2) for n in range(61)] # G. C. Greubel, Oct 10 2022

a(n) = a(n-1)*a(n-2)/a(n-3) with a(0)=1, a(1)=2, a(2)=3.
a(2*n) = (3/2)*a(2*n-1) = 3^n, a(2*n+1) = 2*a(2*n) = 2*3^n.
From Benoit Cloitre, Apr 27 2003: (Start)
a(1)=1, a(n)= 2*a(n-1) if a(n-1) is odd, or a(n)= (3/2)*a(n-1) if a(n-1) is even.
a(n) = (1/6)*(5-(-1)^n)*3^floor(n/2).
a(2*n) = a(2*n-1) + a(2*n-2) + a(2*n-3).
a(2*n+1) = a(2*n) + a(2*n-1). (End)
G.f.: (1+2*x)/(1-3*x^2). - Paul Barry, Aug 25 2003
From Reinhard Zumkeller, Sep 11 2003: (Start)
a(n) = (1 + n mod 2) * 3^floor(n/2).
a(n) = A087503(n) - A087503(n-1). (End)
a(n) = sqrt(3)*(2+sqrt(3))*(sqrt(3))^n/6 - sqrt(3)*(2-sqrt(3))*(-sqrt(3))^n/6. - Paul Barry, Sep 16 2003
From Reinhard Zumkeller, May 26 2008: (Start)
a(n) = A140740(n+2,2).
a(n+1) = a(n) + a(n - n mod 2). (End)
If p(i) = Fibonacci(i-3) and if A is the Hessenberg matrix of order n defined by A(i,j) = p(j-i+1), (i<=j), A(i,j)=-1, (i=j+1), and A(i,j)=0 otherwise. Then, for n>=1, a(n-1) = (-1)^n det A. - Milan Janjic, May 08 2010
a(n) = A182751(n) for n >= 2. - Jaroslav Krizek, Nov 27 2010
a(n) = Sum_{i=0..2^(n+1), i==0 (mod 3)} (-1)^A000120(i). - Vladimir Shevelev, May 16 2012
a(0)=1, a(1)=2, for n>=3, a(n)=3*a(n-2). - Vladimir Shevelev, May 21 2012
Sum_(n>=0) 1/a(n) = 9/4. - Alexander R. Povolotsky, Aug 24 2012
a(n) = sqrt(3*a(n-1)^2 + (-3)^(n-1)). - Richard R. Forberg, Sep 04 2013
a(n) = 2^((1-(-1)^n)/2)*3^((2*n-1+(-1)^n)/4). - Luce ETIENNE, Aug 11 2014
From Reinhard Zumkeller, Oct 19 2015: (Start)
a(2*n) = A000244(n), a(2*n+1) = A008776(n).
For n > 0: a(n+1) = a(n) + if a(n) odd then min{a(n), a(n-1)} else max{a(n), a(n-1)}, see also A128588. (End)
E.g.f.: (7*cosh(sqrt(3)*x) + 4*sqrt(3)*sinh(sqrt(3)*x) - 4)/3. - Stefano Spezia, Feb 17 2022
Sum_{n>=0} (-1)^n/a(n) = 3/4. - Amiram Eldar, Dec 02 2022

A006720 Somos-4 sequence: a(0)=a(1)=a(2)=a(3)=1; for n >= 4, a(n) = (a(n-1) * a(n-3) + a(n-2)^2) / a(n-4).

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 7, 23, 59, 314, 1529, 8209, 83313, 620297, 7869898, 126742987, 1687054711, 47301104551, 1123424582771, 32606721084786, 1662315215971057, 61958046554226593, 4257998884448335457, 334806306946199122193, 23385756731869683322514, 3416372868727801226636179

Offset: 0

N. J. A. Sloane

From the 5th term on, all terms have a primitive divisor; in other words, a prime divisor that divides no earlier term in the sequence. A proof appears in the Everest-McLaren-Ward paper. - Graham Everest (g.everest(AT)uea.ac.uk), Oct 26 2005
Twelve prime terms are known, occurring at indices 4, 5, 6, 7, 8, 11, 13, 16, 43, 52, 206, 647. The last two have been checked for probable primality only. The 647th term has 18498 decimal digits. Possibly these are the only prime terms in the entire sequence. - Graham Everest (g.everest(AT)uea.ac.uk), Nov 28 2006
The density of primes dividing some term in the sequence is 11/21. - Jeremy Rouse, Sep 18 2013
a(n) is a divisor of a(n+k*(2*n-3)) for all integers n and k. - Peter H van der Kamp, May 18 2015
a(n) is a divisor of A051138(k*(2*n-3)) for all integers n and k. - Helmut Ruhland, Jan 26 2024

Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 565.
G. Everest, A. van der Poorten, I. Shparlinski and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; pp. 9, 179.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Robert G. Wilson v, Table of a(n) for n = 0..100.
Paul Barry, Riordan-Bernstein Polynomials, Hankel Transforms and Somos Sequences, Journal of Integer Sequences, Vol. 15 2012, #12.8.2. - From _N. J. A. Sloane_, Dec 29 2012
Paul Barry, On the Hurwitz Transform of Sequences, Journal of Integer Sequences, Vol. 15 (2012), #12.8.7.
Paul Barry, Riordan arrays, the A-matrix, and Somos 4 sequences, arXiv:1912.01126 [math.CO], 2019.
Paul Barry, Integer sequences from elliptic curves, arXiv:2306.05025 [math.NT], 2023.
H. W. Braden, V. Z. Enolskii and A. N. W. Hone, Bilinear recurrences and addition formulas for hyperelliptic sigma functions, arXiv:math/0501162 [math.NT], 2005.
R. H. Buchholz and R. L. Rathbun, An infinite set of Heron triangles with two rational medians, Amer. Math. Monthly, 104 (1997), 107-115.
Xiangke Chang and Xingbiao Hu, A conjecture based on Somos-4 sequence and its extension, Linear Algebra Appl. 436, No. 11, 4285-4295 (2012).
Harini Desiraju and Brady Haran, The Troublemaker Number, Numberphile video (2022).
S. B. Ekhad and D. Zeilberger, How To Generate As Many Somos-Like Miracles as You Wish, arXiv preprint arXiv:1303.5306[math.CO], 2013.
Graham Everest, Gerard Mclaren and Tom Ward, Primitive divisors of elliptic divisibility sequences, arXiv:math/0409540 [math.NT], 2004-2006.
G. Everest, S. Stevens, D. Tamsett and T. Ward, Primitive divisors of quadratic polynomial sequences, arXiv:math/0412079v1 [math.NT], 2004.
G. Everest et al., Primes generated by recurrence sequences, arXiv:math/0412079 [math.NT], 2006.
G. Everest et al., Primes generated by recurrence sequences, Amer. Math. Monthly, 114 (No. 5, 2007), 417-431.
S. Fomin and A. Zelevinsky, The Laurent phenomenon, arXiv:math/0104241 [math.CO], 2001.
Allan Fordy and Andrew Hone, Discrete integrable systems and Poisson algebras from cluster maps, arXiv preprint arXiv:1207.6072 [nlin.SI], 2012.
A. P. Fordy, Periodic Cluster Mutations and Related Integrable Maps, arXiv preprint arXiv:1403.8061 [math-ph], 2014.
A. P. Fordy, Mutation-periodic quivers, integrable maps and associated Poisson algebras, Phil Trans. R. Soc. Lond. Ser A (Math. Phys. Eng. Sci.) 369 (1939) (2011) 1264-1279.
David Gale, The strange and surprising saga of the Somos sequences, in Mathematical Entertainments, Math. Intelligencer 13(1) (1991), pp. 40-42.
R. W. Gosper and Richard C. Schroeppel, Somos Sequence Near-Addition Formulas and Modular Theta Functions, arXiv:math/0703470 [math.NT], 2007.
A. N. W. Hone, Sigma function solution of the initial value problem for Somos 5 sequences, arXiv:math/0501554 [math.NT], 2005-2006.
A. N. W. Hone, Algebraic curves, integer sequences and a discrete Painlevé transcendent, Proceedings of SIDE 6, Helsinki, Finland, 2004; arXiv:0807.2538 [nlin.SI], 2008.
A. N. W. Hone, Elliptic curves and quadratic recurrence sequences, Bull. Lond. Math. Soc. 37 (2005) 161-171.
Andrew N. W. Hone, Growth of Mahler measure and algebraic entropy of dynamics with the Laurent property, arXiv:2109.08217 [math.NT], 2021.
A. N. W. Hone and R. Inoue, Discrete Painlevé equations from Y-systems, arXiv preprint arXiv:1405.5379 [math-ph], 2014
R. Jones and J. Rouse, Galois Theory of Iterated Endomorphisms, arXiv:0706.2384 [math.NT], 2007-2009; Proceedings of the London Mathematical Society, 100, no. 3 (2010), 763-794.
Xinrong Ma, Magic determinants of Somos sequences and theta functions, Discrete Mathematics 310.1 (2010): 1-5.
J. L. Malouf, An integer sequence from a rational recursion, Discr. Math. 110 (1992), 257-261.
Valentin Ovsienko and Serge Tabachnikov, Dual numbers, weighted quivers, and extended Somos and Gale-Robinson sequences, arXiv:1705.01623 [math.CO], 2017. See p. 3.
Kevin I. Piterman and Leandro Vendramin, Computer algebra with GAP, 2023. See p. 39.
J. Propp, The Somos Sequence Site
J. Propp, The 2002 REACH tee-shirt
R. M. Robinson, Periodicity of Somos sequences, Proc. Amer. Math. Soc., 116 (1992), 613-619.
Helmut Ruhland, Somos-4 and a quartic Surface in RP^3, arXiv:2312.02085 [math.AG], 2023.
Matthew Christopher Russell, Using experimental mathematics to conjecture and prove theorems in the theory of partitions and commutative and non-commutative recurrences, PhD Dissertation, Mathematics Department, Rutgers University, May 2016.
Vladimir Shevelev and Peter J. C. Moses, On a sequence of polynomials with hypothetically integer coefficients, arXiv preprint arXiv:1112.5715 [math.NT], 2011.
Michael Somos, Somos 6 Sequence
Michael Somos, Brief history of the Somos sequence problem
Michael Somos, Four polynomial sequences w,x,y,z are discrete versions of the four Jacobi theta functions or the four Weierstrass sigma functions, 2016.
D. E. Speyer, Perfect matchings and the octahedral recurrence, arXiv:math/0402452 [math.CO], 2004.
Alex Stone, The Astonishing Behavior of Recursive Sequences, Quanta Magazine, Nov 16 2023, 13 pages.
Andrei K. Svinin, Somos-4 equation and related equations, arXiv:2307.05866 [math.CA], 2023.
P. H. van der Kamp, Somos-4 and Somos-5 are arithmetic divisibility sequences, arXiv:1505.00194 [math.NT], 2015.
A. J. van der Poorten, Recurrence relations for elliptic sequences: every Somos 4 is a Somos k, arXiv:math/0412293 [math.NT], 2004.
A. J. van der Poorten, Hyperelliptic curves, continued fractions and Somos sequences, arXiv:math/0608247 [math.NT], 2006.
A. J. van der Poorten, Elliptic curves and continued fractions, J. Int. Sequences, Volume 8, no. 2 (2005), article 05.2.5.
Leandro Vendramin, Mini-couse on GAP - Exercises, Universidad de Buenos Aires (Argentina, 2020).
Eric Weisstein's World of Mathematics, Somos Sequence
Index entries for two-way infinite sequences

Cf. A006721, A006722, A006723, A006769, A048736, A028945, A028935, A151502, A165896, A188313, A051138, A006769.
For primes see A129739, A129740, A129741.
Cf. A227199 (primes dividing some term).
Cf. A178384, A247368.

a006720 n = a006720_list !! n
a006720_list = [1,1,1,1] ++
   zipWith div (foldr1 (zipWith (+)) (map b [1..2])) a006720_list
   where b i = zipWith (*) (drop i a006720_list) (drop (4-i) a006720_list)
-- Reinhard Zumkeller, Jan 22 2012

I:=[1,1,1,1]; [n le 4 select I[n] else (Self(n-1)*Self(n-3)+Self(n-2)^2)/Self(n-4): n in [1..30]]; // Vincenzo Librandi, Aug 07 2017

Digits:=11; f(x):=4*x^3-4*x+1;sols:=evalf(solve(f(x),x)); e1:=Re(sols[1]); e3:=Re(sols[2]); w1:=evalf(Int((f(x))^(-0.5),x=e1..infinity)); w3:=I*evalf(Int((-f(x))^(-0.5),x=-infinity..e3)); k:=2*w1-evalf(Int((f(x))^(-0.5),x=1..infinity)); z0:=w3+evalf(Int((f(x))^(-0.5),x=e3..-1)); A:=1/WeierstrassSigma(z0,4.0,-1.0); B:=WeierstrassSigma(k,4.0,-1.0)/WeierstrassSigma(z0+k,4.0,-1.0)/A; for n from 0 to 10 do a[n]:=A*B^n*WeierstrassSigma(z0+n*k,4.0,-1.0)/(WeierstrassSigma(k,4.0,-1.0))^(n^2) od; # Andrew Hone, Oct 12 2005
A006720 := proc(n)
    option remember;
    if n <= 3 then
        1;
    else
        (procname(n-1)*procname(n-3)+procname(n-2)^2)/procname(n-4) ;
    end if;
end proc: # R. J. Mathar, Jul 12 2012

a[0] = a[1] = a[2] = a[3] = 1; a[n_] := a[n] = (a[n - 1] a[n - 3] + a[n - 2]^2)/a[n - 4]; Array[a, 23] (* Robert G. Wilson v, Jul 04 2007 *)
RecurrenceTable[{a[0]==a[1]==a[2]==a[3]==1,a[n]==(a[n-1]a[n-3]+a[n-2]^2)/ a[n-4]},a,{n,30}] (* Harvey P. Dale, Apr 07 2018 *)
b[ n_] := If[-2<=n<=2, {2, 1, 1, 3, 23}[[n+3]], 2*a[n+2]^3*a[n+3] + a[n+1]^2*(a[n+3]*a[n+4] - a[n+2]*a[n+5])]; a[ n_] := If[OddQ[n], b[(n-3)/2], b[-n/2]]; (* Michael Somos, Feb 28 2022 *)

a=vector(99);a[1]=a[2]=a[3]=a[4]=1;for(n=5,#a,a[n]=(a[n-1]*a[n-3]+a[n-2]^2)/a[n-4]); a \\ Charles R Greathouse IV, Jun 16 2011

from gmpy2 import divexact
A006720 = [1, 1, 1, 1]
for n in range(4, 101):
    A006720.append(divexact(A006720[n-1]*A006720[n-3]+A006720[n-2]**2,A006720[n-4]))
# Chai Wah Wu, Sep 01 2014

a(n) = a(3-n) = (-1)^n * A006769(2*n-3) for all n in Z.
a(n+1)/a(n) seems to be asymptotic to C^n with C = 1.226.... - Benoit Cloitre, Aug 07 2002. Confirmed by Hone - see below.
The terms of the sequence have the leading order asymptotics log a(n) ~ D n^2 with D = zeta(w1)*k^2/(2*w1) - log|sigma(k)| = 0.10222281... where zeta and sigma are the Weierstrass functions with invariants g2 = 4, g3 = -1, w1 = 1.496729323 is the real half-period of the corresponding elliptic curve, k = -1.134273216 as above. This agrees with Benoit Cloitre's numerical result with C = exp(2D) = 1.2268447... - Andrew Hone, Feb 09 2005
a(n) = (a(n-1)*a(n-3) + a(n-2)^2)/a(n-4); a(0) = a(1) = a(2) = a(3) = 1; exact formula is a(n) = A*B^n*sigma (z_0+nk)/(sigma (k))^(n^2), where sigma is the Weierstrass sigma function associated to the elliptic curve y^2 = 4*x^3-4*x+1, A = 1/sigma(z_0) = 0.112724016 - 0.824911687*i, B = sigma(k)*sigma (z_0)/sigma (z_0+k) = 0.215971963 + 0.616028193*i, k = 1.859185431, z_0 = 0.204680500 + 1.225694691*i, sigma(k) = 1.555836426, all to 9 decimal places. This is a special case of a general formula for 4th-order bilinear recurrences. The Somos-4 sequence corresponds to the sequence of points (2n-3)P on the curve, where P = (0, 1). - Andrew Hone, Oct 12 2005
a(2*n) = b(-n), a(2*n+1) = b(n-1) where b(n) = A188313(n) for all n in Z. - Michael Somos, Feb 27 2022

A060006 Decimal expansion of real root of x^3 - x - 1 (the plastic constant).

Original entry on oeis.org

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

Offset: 1

Fabian Rothelius, Mar 14 2001

Has been also called the silver number, also the plastic number.
This is the smallest Pisot-Vijayaraghavan number.
The name "plastic number" goes back to the Dutch Benedictine monk and architect Dom Hans van der Laan, who gave this name 4 years after the discovery of the number by the French engineer Gérard Cordonnier in 1924, who used the name "radiant number". - Hugo Pfoertner, Oct 07 2018
Sometimes denoted by the symbol rho. - Ed Pegg Jr, Feb 01 2019
Also the solution of 1/x + 1/(1+x+x^2) = 1. - Clark Kimberling, Jan 02 2020
Given any complex p such that real(p)>-1, this constant is the only real solution of the equation z^p+z^(p+1)=z^(p+3), and the only attractor of the complex mapping z->M(z,p), where M(z,p)=(z^p+z^(p+1))^(1/(p+3)), convergent from any complex plane point. - Stanislav Sykora, Oct 14 2021
The Pisot-Vijayaraghavan numbers were named after the French mathematician Charles Pisot (1910-1984) and the Indian mathematician Tirukkannapuram Vijayaraghavan (1902-1955). - Amiram Eldar, Apr 02 2022
The sequence a(n) = v_3^floor(n^2/4) where v_n is the smallest, positive, real solution to the equation (v_n)^n = v_n + 1 satisfies the Somos-5 recursion a(n+3)*a(n-2) = a(n+2)*a(n-1) + a(n+1)*a(n) for all n in Z. Also true if floor is removed. - Michael Somos, Mar 24 2023

			1.32471795724474602596090885447809734...

Steven R. Finch, Mathematical Constants, Cambridge, 2003, Section 1.2.2.
Midhat J. Gazalé, Gnomon: From Pharaohs to Fractals, Princeton University Press, Princeton, NJ, 1999, see Chap. VII.
Donald E. Knuth, The Art of Computer Programming, Vol. 4A, Section 7.1.4, p. 236.
Ian Stewart, A Guide to Computer Dating (Feedback), Scientific American, Vol. 275 No. 5, November 1996, p. 118.

Harry J. Smith, Table of n, a(n) for n = 1..20000
Alex Bellos, The golden ratio has spawned a beautiful new curve: the Harriss spiral, The Guardian, Jan 13 2015.
Gamaliel Cerda-Morales, New Identities for Padovan Numbers, arXiv:1904.05492 [math.CO], 2019.
Brady Haran and Edmund Harriss, The Plastic Ratio, Numberphile video (2019).
Ed Pegg Jr., Pictures based on the plastic constant
Simon Plouffe, Smallest Pisot-Vijayaraghavan number to 50000 digits
Simon Plouffe, The Smallest Pisot-Vijayaraghavan number
F. Rothelius, Formulae.
Ian Stewart, Tales of a Neglected Number, Mathematical Recreations, Scientific American, Vol. 274, No. 6 (1996), pp. 102-103.
Ian Stewart, Tales of a Neglected Number, Mathematical Recreations, Scientific American, Vol. 274, No. 6 (1996), pp. 102-103.
Dom Hans van der Laan, Le nombre plastique: Quinze leçons sur l’ordonnance architectonique, Brill Academic Pub., Leiden, 1960.
Michel Waldschmidt, Lectures on Multiple Zeta Values, IMSC 2011.
Eric Weisstein's World of Mathematics, Maverick Graph.
Eric Weisstein's World of Mathematics, Pisot Number.
Eric Weisstein's World of Mathematics, Plastic Constant.
Wikipedia, Plastic number.
Index entries for algebraic numbers, degree 3

Cf. A001622. A072117 gives continued fraction.
Cf. A006888, A010527, A051016, A051017, A084252, A075778 (inverse), A126772.
Other Pisot numbers: A086106, A092526, A228777, A293506, A293508, A293509, A293557.
Cf. A002620, A006721, A134816.

SetDefaultRealField(RealField(100)); ((3+Sqrt(23/3))/6)^(1/3) + ((3-Sqrt(23/3))/6)^(1/3); // G. C. Greubel, Mar 15 2019

(1/2 +sqrt(23/3)/6)^(1/3) + (1/2-sqrt(23/3)/6)^(1/3) ; evalf(%,130) ; # R. J. Mathar, Jan 22 2013

RealDigits[ Solve[x^3 - x - 1 == 0, x][[1, 1, 2]], 10, 111][[1]] (* Robert G. Wilson v, Sep 30 2009 *)
s = Sqrt[23/108]; RealDigits[(1/2 + s)^(1/3) + (1/2 - s)^(1/3), 10, 111][[1]] (* Robert G. Wilson v, Dec 12 2017 *)
RealDigits[Root[x^3-x-1,1],10,120][[1]] (* or *) RealDigits[(Surd[9-Sqrt[69],3]+Surd[9+Sqrt[69],3])/(Surd[2,3]Surd[9,3]),10,120][[1]] (* Harvey P. Dale, Sep 04 2018 *)

allocatemem(932245000); default(realprecision, 20080); x=solve(x=1, 2, x^3 - x - 1); for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b060006.txt", n, " ", d)); \\ Harry J. Smith, Jul 01 2009

(1/2 +sqrt(23/3)/6)^(1/3) + (1/2-sqrt(23/3)/6)^(1/3) \\ Altug Alkan, Apr 10 2016

polrootsreal(x^3-x-1)[1] \\ Charles R Greathouse IV, Aug 28 2016

default(realprecision, 110); digits(floor(solve(x=1, 2, x^3 - x - 1)*10^105)) /* Michael Somos, Mar 24 2023 */

numerical_approx(((3+sqrt(23/3))/6)^(1/3) + ((3-sqrt(23/3))/6)^(1/3), digits=100) # G. C. Greubel, Mar 15 2019

Equals (1/2+sqrt(23/108))^(1/3) + (1/2-sqrt(23/108))^(1/3). - Henry Bottomley, May 22 2003
Equals CubeRoot(1 + CubeRoot(1 + CubeRoot(1 + CubeRoot(1 + ...)))). - Gerald McGarvey, Nov 26 2004
Equals sqrt(1+1/sqrt(1+1/sqrt(1+1/sqrt(1+...)))). - Gerald McGarvey, Mar 18 2006
Equals (1/2 +sqrt(23/3)/6)^(1/3) + (1/2-sqrt(23/3)/6)^(1/3). - Eric Desbiaux, Oct 17 2008
Equals Sum_{k >= 0} 27^(-k)/k!*(Gamma(2*k+1/3)/(9*Gamma(k+4/3)) - Gamma(2*k-1/3)/(3*Gamma(k+2/3))). - Robert Israel, Jan 13 2015
Equals sqrt(Phi) = sqrt(1.754877666246...) (see A109134). - Philippe Deléham, Sep 29 2020
Equals cosh(arccosh(3*c)/3)/c, where c = sqrt(3)/2 (A010527). - Amiram Eldar, May 15 2021
Equals 1/hypergeom([1/5, 2/5, 3/5, 4/5], [2/4, 3/4, 5/4], -5^5/4^4). - Gerry Martens, Mar 16 2025

Edited and extended by Robert G. Wilson v, Aug 03 2002

Removed incorrect comments, Joerg Arndt, Apr 10 2016

A048736 Dana Scott's sequence: a(n) = (a(n-2) + a(n-1) * a(n-3)) / a(n-4), a(0) = a(1) = a(2) = a(3) = 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 5, 13, 22, 41, 111, 191, 361, 982, 1693, 3205, 8723, 15042, 28481, 77521, 133681, 253121, 688962, 1188083, 2249605, 6123133, 10559062, 19993321, 54419231, 93843471, 177690281, 483649942, 834032173, 1579219205, 4298430243, 7412446082, 14035282561, 38202222241, 65877982561

Offset: 0

David Johnson-Davies

The recursion has the Laurent property. If a(0), a(1), a(2), a(3) are variables, then a(n) is a Laurent polynomial (a rational function with a monic monomial denominator). - Michael Somos, Feb 05 2012

A generalization is if the recursion is modified to a(n) = (a(n-2) + a(n-1) * b*a(n-3)) / a(n-4) where b is a constant, and with arbitrary nonzero initial values, (a(0), a(1), a(2), a(3)), then a(n) = c*(a(n-3) - a(n-6)) + a(n-9) for all n in Z where c is another constant. - Michael Somos, Oct 28 2021

			G.f. = 1 + x + x^2 + x^3 + 2*x^4 + 3*x^5 + 13*x^6 + 22*x^7 + 41*x^8 + 111*x^9 + ...

Seiichi Manyama, Table of n, a(n) for n = 0..3165 (first 501 terms from T. D. Noe)
Joshua Alman, Cesar Cuenca, and Jiaoyang Huang, Laurent phenomenon sequences, Journal of Algebraic Combinatorics 43(3) (2015), 589-633.
Hal Canary, The Dana Scott Recurrence [From _Jaume Oliver Lafont_, Sep 25 2009]
S. Fomin and A. Zelevinsky, The Laurent phenomenon, arXiv:math/0104241 [math.CO], 2001.
S. Fomin and A. Zelevinsky, The Laurent Phenomenon, Advances in Applied Mathematics, 28 (2002), 119-144.
David Gale, The strange and surprising saga of the Somos sequences, Math. Intelligencer 13(1) (1991), pp. 40-42.
D. Gale, Tracking the Automatic Ant And Other Mathematical Explorations, A Collection of Mathematical Entertainments Columns from The Mathematical Intelligencer, Springer, 1998, p. 4.
Matthew Christopher Russell, Using experimental mathematics to conjecture and prove theorems in the theory of partitions and commutative and non-commutative recurrences, PhD Dissertation, Mathematics Department, Rutgers University, May 2016.
Eric Weisstein's World of Mathematics, Laurent Polynomial
Index entries for two-way infinite sequences
Index entries for linear recurrences with constant coefficients, signature (0,0,10,0,0,-10,0,0,1)

Cf. A006720, A006721, A006722, A006723, A092420, A072881.
Cf. A192241, A192242 (primes and where they occur).
Cf. A276531.

a048736 n = a048736_list !! n
a048736_list = 1 : 1 : 1 : 1 :
   zipWith div
     (zipWith (+)
       (zipWith (*) (drop 3 a048736_list)
                    (drop 1 a048736_list))
       (drop 2 a048736_list))
     a048736_list
-- Reinhard Zumkeller, Jun 26 2011

I:=[1,1,1,1]; [n le 4 select I[n] else (Self(n-2) + Self(n-1)*Self(n-3)) / Self(n-4): n in [1..30]]; // G. C. Greubel, Feb 20 2018

P:=proc(q) local n,v; v:=[1,1,1,1]; for n from 1 to q do
v:=[op(v),(v[-2]+v[-1]*v[-3])/v[-4]] od: op(v); end: P(35); # Paolo P. Lava, Aug 24 2025

RecurrenceTable[{a[0] == a[1] == a[2] == a[3] == 1, a[n] == (a[n - 2] + a[n - 1]a[n - 3])/a[n - 4]}, a[n], {n, 40}] (* or *) LinearRecurrence[{0, 0, 10, 0, 0, -10, 0, 0, 1}, {1, 1, 1, 1, 2, 3, 5, 13, 22}, 41] (* Harvey P. Dale, Oct 22 2011 *)

Vec((1+x+x^2-9*x^3-8*x^4-7*x^5+5*x^6+3*x^7+2*x^8) / (1-10*x^3+10*x^6-x^9)+O(x^99)) \\ Charles R Greathouse IV, Jul 01 2011

a(n) = 9*a(n-3) - a(n-6) - 3 - ( ceiling(n/3) - floor(n/3) ), with a(0) = a(1) = a(2) = a(3) = 1, a(4) = 2, a(5) = 3. - Michael Somos
From Jaume Oliver Lafont, Sep 17 2009: (Start)
a(n) = 10*a(n-3) - 10*a(n-6) + a(n-9).
G.f.: (1 + x + x^2 - 9*x^3 - 8*x^4 - 7*x^5 + 5*x^6 + 3*x^7 + 2*x^8)/(1 - 10*x^3 + 10*x^6 - x^9). (End)
a(n) = a(3-n) for all n in Z. - Michael Somos, Feb 05 2012

More terms from Michael Somos

A006722 Somos-6 sequence: a(n) = (a(n-1) * a(n-5) + a(n-2) * a(n-4) + a(n-3)^2) / a(n-6), a(0) = ... = a(5) = 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 3, 5, 9, 23, 75, 421, 1103, 5047, 41783, 281527, 2534423, 14161887, 232663909, 3988834875, 45788778247, 805144998681, 14980361322965, 620933643034787, 16379818848380849, 369622905371172929, 20278641689337631649, 995586066665500470689

Offset: 0

N. J. A. Sloane

C. Pickover, Mazes for the Mind, St. Martin's Press, NY, 1992, p. 350.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n=0..100
R. H. Buchholz and R. L. Rathbun, An infinite set of Heron triangles with two rational medians, Amer. Math. Monthly, 104 (1997), 107-115.
Chang, Xiangke; Hu, Xingbiao, A conjecture based on Somos-4 sequence and its extension, Linear Algebra Appl. 436, No. 11, 4285-4295 (2012).
Yuri N. Fedorov and Andrew N. W. Hone, Sigma-function solution to the general Somos-6 recurrence via hyperelliptic Prym varieties, arXiv:1512.00056 [nlin.SI], 2015.
S. Fomin and A. Zelevinsky, The Laurent phenomenon, arXiv:math/0104241 [math.CO], 2001.
David Gale, The strange and surprising saga of the Somos sequences, Math. Intelligencer 13(1) (1991), pp. 40-42.
R. W. Gosper and Richard C. Schroeppel, Somos Sequence Near-Addition Formulas and Modular Theta Functions, arXiv:math/0703470 [math.NT], 2007.
A. N. W. Hone, Analytic solutions and integrability for bilinear recurrences of order six, Appl. Anal. 89, no.4 (2010), 473-492.
J. L. Malouf, An integer sequence from a rational recursion, Discr. Math. 110 (1992), 257-261.
R. M. Robinson, Periodicity of Somos sequences, Proc. Amer. Math. Soc., 116 (1992), 613-619.
Vladimir Shevelev and Peter J. C. Moses, On a sequence of polynomials with hypothetically integer coefficients, arXiv preprint arXiv:1112.5715 [math.NT], 2011.
Michael Somos, Somos 6 Sequence
Michael Somos, Brief history of the Somos sequence problem
Alex Stone, The Astonishing Behavior of Recursive Sequences, Quanta Magazine, Nov 16 2023, 13 pages.
A. van der Poorten, Hyperelliptic curves, continued fractions and Somos sequences, arXiv:math/0608247 [math.NT], 2006.
Eric Weisstein's World of Mathematics, Somos Sequence.
Index entries for two-way infinite sequences

Cf. A006720, A006721, A006723, A048736.

a006722 n = a006722_list !! n
a006722_list = [1,1,1,1,1,1] ++
  zipWith div (foldr1 (zipWith (+)) (map b [1..3])) a006722_list
  where b i = zipWith (*) (drop i a006722_list) (drop (6-i) a006722_list)
-- Reinhard Zumkeller, Jan 22 2012

[n le 6 select 1 else (Self(n-1)*Self(n-5)+Self(n-2)*Self(n-4)+ Self(n-3)^2)/Self(n-6): n in [1..30]]; // Vincenzo Librandi, Dec 02 2015

a[n_ /; 0 <= n <= 5] = 1; a[n_] := a[n] = (a[n-1]*a[n-5] + a[n-2]*a[n-4] + a[n-3]^2) / a[n-6]; Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Nov 22 2013 *)
RecurrenceTable[{a[0]==a[1]==a[2]==a[3]==a[4]==a[5]==1,a[n]==(a[n-1]a[n-5]+ a[n-2]a[n-4]+a[n-3]^2)/a[n-6]},a,{n,30}] (* Harvey P. Dale, Dec 20 2014 *)

{a(n) = if( n>-1 && n<6, 1, if( n<0, a(5 - n), (a(n - 1) * a(n - 5) + a(n - 2) * a(n - 4) + a(n-3) * a(n-3)) / a(n - 6)))}; /* Michael Somos, Jan 30 2012 */

from gmpy2 import divexact
A006722 = [1,1,1,1,1,1]
for n in range(6,101):
    A006722.append(divexact(A006722[n-1]*A006722[n-5]+A006722[n-2]*A006722[n-4]+A006722[n-3]**2,A006722[n-6]))
# Chai Wah Wu, Sep 01 2014

a(n) = a(5-n).
Michael Somos found an explicit formula for a(n) in 1993, which is not as widely known as it should be. The following is a quotation from the "Somos 6 sequence" document mentioned in the Links section: (Start)
This sequence is one of a large class of sequences of numbers that satisfy a non-linear recurrence relation depending on previous terms. It is also one of the class of sequences which can be computed from a theta series, hence I call them theta sequences. Here are the details:
Fix the following seven constants:
c1 = 0.875782749065950194217251...,
c2 = 1.084125925473763343779968...,
c3 = 0.114986002186402203509006...,
c4 = 0.077115634258697284328024...,
c5 = 1.180397390176742642553759...,
c6 = 1.508030831265086447098989..., and
c7 = 2.551548771413081602906643... .
Consider the doubly indexed series: f(x,y) = c1*c2^(x*y)*sum(k2, (-1)^k2*sum(k1, g(k1,k2,x,y))) , where g(k1,k2,x,y) = c3^(k1*k1)*c4^(k2*k2)*c5^(k1*k2)*cos(c6*k1*x+c7*k2*y) . Here both sums range over all integers.
Then the sequence defined by a(n) = f(n-2.5,n-2.5) is the Somos 6 sequence. I announced this in 1993. (End) - N. J. A. Sloane, Dec 06 2015
From Andrew Hone and Yuri Fedorov, Nov 27 2015: (Start)
The following is an exact formula for a(n):
a(n+3) = A*B^n*C^(n^2 -1)*sigma(v_0 + n*v) / sigma(v)^(n^2),
where
A = C / sigma(v_0),
B = A^(-1)*sigma(v) / sigma(v_0+v),
C = i/sqrt(20) (with i the imaginary unit),
sigma is the two-variable Kleinian sigma-function associated with the genus two curve X: y^2 = 4*x^5 - 233*x^4 + 1624*x^3 - 422*x^2 + 36*x - 1, and
v and v_0 are two-component vectors in the Jacobian of X, being the images under the Abel map of the divisors P_1+P_2 - 2*infinity,  Q_1 + Q_2 - 2*infinity, respectively, where points P_j and Q_j on X are given by
P_1 = ( -8 + sqrt(65), 20*i*(129 -16*sqrt(65)) ),
P_2 = ( -8 - sqrt(65), 20*i*(129 +16*sqrt(65)) ),
Q_1 = ( 5 + 2*sqrt(6), 4*i*(71 +sqrt(6)) ),
Q_2 = ( 5 - 2*sqrt{6}, 4*i*(71 -sqrt(6)) ).
The Abel map is based at infinity and calculated with respect to the basis of holomorphic differentials dx/y, x dx/y.
Approximate values from Maple are A =  0.0619-0.0317*i, B = -0.0000973-0.0000158*i, v = (-.341*i, .477*i), v_0 = (-.379-.150*i, -.259+.576*i).
(End)

More terms from James A. Sellers, Aug 22 2000

A006723 Somos-7 sequence: a(n) = (a(n-1) * a(n-6) + a(n-2) * a(n-5) + a(n-3) * a(n-4)) / a(n-7), a(0) = ... = a(6) = 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 3, 5, 9, 17, 41, 137, 769, 1925, 7203, 34081, 227321, 1737001, 14736001, 63232441, 702617001, 8873580481, 122337693603, 1705473647525, 22511386506929, 251582370867257, 9254211194697641, 215321535159114017

Offset: 0

N. J. A. Sloane

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n=0..100
R. H. Buchholz and R. L. Rathbun, An infinite set of Heron triangles with two rational medians, Amer. Math. Monthly, 104 (1997), 107-115.
Chang, Xiangke; Hu, Xingbiao, A conjecture based on Somos-4 sequence and its extension, Linear Algebra Appl. 436, No. 11, 4285-4295 (2012).
S. Fomin and A. Zelevinsky, The Laurent phenomenon, arXiv:math/0104241 [math.CO], 2001.
David Gale, The strange and surprising saga of the Somos sequences, Math. Intelligencer 13(1) (1991), pp. 40-42.
R. W. Gosper and Richard C. Schroeppel, Somos Sequence Near-Addition Formulas and Modular Theta Functions, arXiv:math/0703470 [math.NT], 2007.
J. L. Malouf, An integer sequence from a rational recursion, Discr. Math. 110 (1992), 257-261.
J. Propp, The Somos Sequence Site
R. M. Robinson, Periodicity of Somos sequences, Proc. Amer. Math. Soc., 116 (1992), 613-619.
Vladimir Shevelev and Peter J. C. Moses, On a sequence of polynomials with hypothetically integer coefficients, arXiv preprint arXiv:1112.5715 [math.NT], 2011.
Michael Somos, Somos 7 Sequence
Michael Somos, Brief history of the Somos sequence problem
Alex Stone, The Astonishing Behavior of Recursive Sequences, Quanta Magazine, Nov 16 2023, 13 pages.
A. van der Poorten, Hyperelliptic curves, continued fractions and Somos sequences, arXiv:math/0608247 [math.NT], 2006.
Eric Weisstein's World of Mathematics, Somos Sequence.
Index entries for two-way infinite sequences

Cf. A006720, A006721, A006722, A048736.

Cf. A078495.

a006723 n = a006723_list !! n
a006723_list = [1,1,1,1,1,1,1] ++
  zipWith div (foldr1 (zipWith (+)) (map b [1..3])) a006723_list
  where b i = zipWith (*) (drop i a006723_list) (drop (7-i) a006723_list)
-- Reinhard Zumkeller, Jan 22 2012

I:=[1,1,1,1,1,1,1]; [n le 7 select I[n] else (Self(n-1)*Self(n-6) + Self(n-2)*Self(n-5) + Self(n-3)*Self(n-4))/Self(n-7): n in [1..30]]; // G. C. Greubel, Feb 21 2018

RecurrenceTable[{a[0]==a[1]==a[2]==a[3]==a[4]==a[5]==a[6]==1,a[n] == (a[n-1]a[n-6]+a[n-2]a[n-5]+a[n-3]a[n-4])/a[n-7]},a,{n,30}] (* Harvey P. Dale, Jan 19 2012 *)

{a(n) = my(v); if( n<0, n = 6-n); if( n<7, 1, n++; v = vector(n, k, 1); for( k=8, n, v[k] = (v[k-1] * v[k-6] + v[k-2] * v[k-5] + v[k-3] * v[k-4]) / v[k-7]); v[n])};

from gmpy2 import divexact
A006723 = [1,1,1,1,1,1,1]
for n in range(7,101):
    A006723.append(divexact(A006723[n-1]*A006723[n-6]+A006723[n-2]*A006723[n-5]+A006723[n-3]*A006723[n-4],A006723[n-7]))
# Chai Wah Wu, Sep 01 2014

a(6 - n) = a(n) for all n in Z.

a(n) = ((8-2*(-1)^n)*a(n-5)*a(n-3)-a(n-4)^2)/a(n-8). - Bruno Langlois, Aug 09 2016

More terms from James A. Sellers, Aug 22 2000

A078495 a(n) = (a(n-1) * a(n-6) + a(n-3) * a(n-4)) / a(n-7) (a variant of Somos-7).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 3, 4, 6, 12, 24, 72, 144, 288, 864, 3456, 10368, 41472, 124416, 497664, 2985984, 17915904, 71663616, 429981696, 2579890176, 20639121408, 185752092672, 1486016741376, 8916100448256, 106993205379072

Offset: 0

Michael Somos, Nov 26 2002

From Vladimir Shevelev, Apr 16 2016: (Start)
For k >= 0, an infinite sequence {b(k,n)} of Somos's sequences (n >= 0) is:
b(k,0) = b(k,1) = ... = b(k,2*k+2) = 1;
and then for n >= 2*k+3,
b(k,n) = (b(k,n-1)*b(k,n-2*k-2) + b(k,n-k-1)*b(k,n-k-2))/b(k,n-2*k-3).
In particular, {b(0,n)} is essentially A060656, {b(1,n)}=A006721, {a(2,n)}=A078495.
One can prove that the sequence {b(k,n)} has the first 4*(k+1) simple differences: 2k+2 zeros, after that k+1 1's and after that k+1 consecutive doubled triangular numbers (A000217), beginning with 2.
Further we have nontrivial differences. The first of them for k=0,1,2,... are 12, 26, 48, 80, 124, 182, 256, 348, 460, 594, ..., that is, {k^3/3 + 3*k^2 + 32*k/3 + 12}.
(End)

G. Everest, A. van der Poorten, I. Shparlinski and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; see esp. p. 255.

Cf. A000217, A006721, A006723, A025795, A060656, A078495.

a078495 n = a078495_list !! n
a078495_list = [1, 1, 1, 1, 1, 1, 1] ++
  zipWith div (foldr1 (zipWith (+)) (map b [1,3])) a078495_list
  where b i = zipWith (*) (drop i a078495_list) (drop (7-i) a078495_list)
-- Reinhard Zumkeller, May 05 2013

I:=[1,1,1,1,1,1,1]; [n le 7 select I[n] else (Self(n-1)*Self(n-6) + Self(n-3)*Self(n-4))/Self(n-7): n in [1..30]]; // G. C. Greubel, Feb 21 2018

RecurrenceTable[{a[0]==a[1]==a[2]==a[3]==a[4]==a[5]==a[6]==1,a[n] == (a[n-1]*a[n-6]+a[n-3]*a[n-4])/a[n-7]},a,{n,40}] (* Harvey P. Dale, Apr 20 2012 *)

{a(n) = if( n<0, a(6-n), if( n<7, 1, (a(n-1) * a(n-6) + a(n-3) * a(n-4)) / a(n-7)))};

{a(n) = 2^(b(n-9) + b(n-7)) * 3^b(n-8)}; {b(n) = (n^2 + 10*n + 1 - n%2*13) \ 60 + 1}; /* b(n) = A025795(n) */

a(n) = 144 * a(n-6) * a(n-10) / a(n-16), a(n) = a(6-n) for all n in Z.

cp's OEIS Frontend

A333251 Tropical version of Somos-5 sequence A006721.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

A038754 a(2n) = 3^n, a(2n+1) = 2*3^n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

SageMath

Formula

A006720 Somos-4 sequence: a(0)=a(1)=a(2)=a(3)=1; for n >= 4, a(n) = (a(n-1) * a(n-3) + a(n-2)^2) / a(n-4).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Formula

A060006 Decimal expansion of real root of x^3 - x - 1 (the plastic constant).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

Sage

Formula

Extensions

A048736 Dana Scott's sequence: a(n) = (a(n-2) + a(n-1) * a(n-3)) / a(n-4), a(0) = a(1) = a(2) = a(3) = 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

A268199 Somos's sequence {b(3,n)} defined in comment in A078495: a(0)=a(1)=...=a(8)=1; for n>=9, a(n)=(a(n-1)a(n-8)+a(n-4)a(n-5))/a(n-9).

A271949 Somos's sequence {b(4,n)} defined in comment in A078495: a(0)=a(1)=...=a(10)=1; for n>=11, a(n)=(a(n-1)a(n-10)+a(n-5)a(n-6))/a(n-11).