A250309 a(n) = a(n-1)*(1 + a(n-1)/a(n-3)), with a(0) = a(1) = a(2) = 1.

1, 1, 1, 2, 6, 42, 924, 143220, 488523420, 258285263294520, 465795819523189050504840, 444125576385425970712647062585372630520, 763680920404535561780141108036287312478667174369871222219397040

Offset

0,4

Links

Seiichi Manyama, Table of n, a(n) for n = 0..17

Formula

0 = a(n)*(a(n+2) - a(n+3)) + a(n+2)*a(n+2) for all n>=0.

A007660(n+1) = a(n)/a(n-1).

a(n) ~ b * f^(d^n), where b = 0.270887790039424376..., f = c^(2+sqrt(5)) = 1.574173161904651669837597516422779594... and d = (1+sqrt(5))/2. For the constant c = A258113 = 1.11305797590293193285359770716758491... see A007660. - Vaclav Kotesovec, Jan 18 2015

Mathematica

RecurrenceTable[{a[n]==a[n-1]*(1 + a[n-1]/a[n-3]), a[0]==1, a[1]==1, a[2]==1}, a, {n, 0, 12}] (* Vaclav Kotesovec, Jan 18 2015 *)

Prog

(PARI) {a(n) = if( n<3, n>=0, a(n-1)*(1 + a(n-1)/a(n-3)))};
(Magma) I:=[1, 1, 1]; [n le 3 select I[n] else Self(n-1)*(1 + Self(n-1)/Self(n-3)): n in [1..15]]; // G. C. Greubel, Aug 03 2018

Crossrefs

Cf. A007660, A258113.

Sequence in context: A326343 A276453 A152479 * A115961 A355042 A367132

Adjacent sequences: A250306 A250307 A250308 * A250310 A250311 A250312

Keyword

nonn

Author

Michael Somos, Jan 16 2015

Status

approved