OFFSET
1,4
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..18
FORMULA
a(n+1) / a(n) = A006277(n-1). - Michael Somos, Dec 29 2012
MATHEMATICA
a[n_]:= a[n]= If[n<4, 1, (a[n-1]*a[n-2] +a[n-1]^2)/a[n-3]]; Table[a[n], {n, 15}]
RecurrenceTable[{a[1]==a[2]==a[3]==1, a[n]==(a[n-1]a[n-2]+a[n-1]^2)/a[n-3]}, a, {n, 14}] (* Harvey P. Dale, Oct 01 2017 *)
PROG
(Magma) [n le 3 select 1 else (Self(n-1)*Self(n-2) +Self(n-1)^2)/Self(n-3): n in [1..15]]; // G. C. Greubel, Sep 21 2024
(SageMath)
def a(n): # a = A141609
if n<3: return 1
else: return (a(n-1)*a(n-2) +a(n-1)^2)/a(n-3)
[a(n) for n in range(1, 16)] # G. C. Greubel, Sep 21 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Aug 22 2008
EXTENSIONS
Edited by N. J. A. Sloane, Aug 24 2008
STATUS
approved