login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A179554 a(n) = ( ( a(n-1) * a(n-5) )^2 + ( a(n-2) * a(n-3)^2 * a(n-4) )^2 ) / a(n-6), a(1) = ... = a(6) = 1. 6
1, 1, 1, 1, 1, 1, 2, 5, 29, 1241, 3642581, 80305336110269, 93845241317908146817718208193, 1831028415236189172750287983272896126569309241020364043033034 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,7

COMMENTS

Similar to Somos-6 sequence but the recurrence is not bilinear.

log( log( a(n))) ~ n * log(c) where c = 2.0810189966245... is largest root of 0 = x^4 - x^3 - 2*x^2 - x + 1 according to Hone.

LINKS

G. C. Greubel, Table of n, a(n) for n = 1..17

A. P. Fordy and R. J. Marsh, Cluster Mutation-Periodic Quivers and Associated Laurent Sequences, arXiv:0904.0200 [math.CO], 2009-2011.

A. N. W. Hone, Laurent Polynomials and Superintegrable Maps, arXiv:math/0702280 [math.NT], 2007.

FORMULA

a(n) = a(7 - n).

MATHEMATICA

a[1] = a[2] = a[3] = a[4] = a[5] = a[6] = 1; a[n_ /; n < 1] := a[7 - n];

a[n_ /; n > 6] := a[n] = ((a[n - 1] a[n - 5])^2 + (a[n - 2] a[n - 3]^2 a[n - 4])^2)/ a[n - 6];

PROG

(MAGMA) I:=[1, 1, 1, 1, 1, 1]; [n le 6 select I[n] else ((Self(n-1)*Self(n-5) )^2 + (Self(n-2)*Self(n-3)^2*Self(n-4))^2 )/Self(n-6): n in [1..20]]; // G. C. Greubel, Aug 14 2018

CROSSREFS

Sequence in context: A121910 A073833 A229918 * A086383 A118612 A187628

Adjacent sequences:  A179551 A179552 A179553 * A179555 A179556 A179557

KEYWORD

nonn

AUTHOR

Michael Somos, Jan 08 2011

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 23 04:54 EDT 2020. Contains 337295 sequences. (Running on oeis4.)