OFFSET
1,5
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..17
FORMULA
a(n+3) = a(n+2) * A245384(n) for n>0.
MATHEMATICA
a[ n_] := If[ n<5, 1, a[n] = a[n-1] (a[n-1] + a[n-2]) / a[n-4]];
RecurrenceTable[{a[n] == a[n-1]*(a[n-1] + a[n-2])/a[n-4], a[1] == a[2] == a[3] == a[4] == 1}, a, {n, 1, 10}] (* G. C. Greubel, Aug 06 2018 *)
nxt[{a_, b_, c_, d_}]:={b, c, d, d (d+c)/a}; NestList[nxt, {1, 1, 1, 1}, 15][[All, 1]] (* Harvey P. Dale, Jun 25 2021 *)
PROG
(PARI) {a(n) = if( n<5, 1, a(n-1) * (a(n-1) + a(n-2)) / a(n-4))};
(Magma) I:=[1, 1, 1, 1]; [n le 4 select I[n] else Self(n-1)*(Self(n-1) + Self(n-2))/Self(n-4): n in [1..10]]; // G. C. Greubel, Aug 06 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael Somos, Aug 21 2014
STATUS
approved