The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A174168 A (1,3) Somos-4 sequence. 2
 1, 2, 5, 17, 109, 706, 9529, 149057, 3464585, 141172802, 5987285341, 540281496721, 53055650250901, 8453268527614082, 2173247123169247025, 610193552206250342657, 364280506115937098013649, 225791960727416355391439746 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,2 COMMENTS Hankel transform of A128720. LINKS Reinhard Zumkeller, Table of n, a(n) for n = 0..100 Paul Barry, Generalized Catalan Numbers, Hankel Transforms and Somos-4 Sequences , J. Int. Seq. 13 (2010) #10.7.2. Paul Barry, Invariant number triangles, eigentriangles and Somos-4 sequences, arXiv:1107.5490 [math.CO], 2011. Paul Barry, On Motzkin-Schröder Paths, Riordan Arrays, and Somos-4 Sequences, J. Int. Seq. (2023) Vol. 26, Art. 23.4.7. P. Rajkovic, P. Barry, and N. Savic, Number Sequences in an Integral Form with a Generalized Convolution Property and Somos-4 Hankel Determinants, Math. Balkanica, Vol. 26 (2012), Fasc. 1-2. FORMULA a(n) = (a(n-1)*a(n-3) + 3*a(n-2)^2)/a(n-4) [index corrected by Franklin T. Adams-Watters, Oct 26 2011]. a(n) = a(-2-n), a(n)*a(n-5) = -3*a(n-1])*a(n-4) + 16*a(n-2)*a(n-3) for all n in Z. - Michael Somos, Apr 30 2020 MAPLE A174168 := proc(n) if n <=4 then op(n+1, [1, 2, 5, 17, 109]) ; else (procname(n-1)*procname(n-3)+3*procname(n-2)^2)/procname(n-4) ; end if; end proc: # R. J. Mathar, Oct 26 2011 MATHEMATICA RecurrenceTable[{a[0]==1, a[1]==2, a[2]==5, a[3]==17, a[n]==(a[n-1]a[n-3]+ 3a[n-2]^2)/a[n-4]}, a[n], {n, 20}] (* Harvey P. Dale, Oct 25 2011 *) a[ n_] := Which[n < -1, a[-2 - n], n < 2, {1, 1, 2}[[n + 2]], True, (a[n - 1] a[n - 3] + 3 a[n - 2]^2)/a[n - 4]]; (* Michael Somos, Apr 30 2020 *) PROG (Haskell) a174168 n = a174168_list !! n a174168_list = [1, 2, 5, 17] ++ zipWith div (zipWith (+) (zipWith (*) (tail a174168_list) (drop 3 a174168_list)) (map ((* 3) . (^ 2)) (drop 2 a174168_list))) a174168_list -- Reinhard Zumkeller, Feb 13 2012 (PARI) m=20; v=concat([1, 2, 5, 17], vector(m-4)); for(n=5, m, v[n] = ( v[n-1]*v[n-3] +3*v[n-2]^2)/v[n-4]); v \\ G. C. Greubel, Sep 25 2018 (Magma) I:=[1, 2, 5, 17]; [n le 4 select I[n] else (Self(n-1)*Self(n-3) + 3*Self(n-2)^2)/Self(n-4): n in [1..20]]; // G. C. Greubel, Sep 25 2018 CROSSREFS Cf. A128720. Sequence in context: A254132 A161609 A228238 * A097980 A210525 A074046 Adjacent sequences: A174165 A174166 A174167 * A174169 A174170 A174171 KEYWORD easy,nonn AUTHOR Paul Barry, Mar 10 2010 EXTENSIONS More terms from Sean A. Irvine, Feb 12 2012 STATUS approved

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

Last modified February 21 04:30 EST 2024. Contains 370219 sequences. (Running on oeis4.)