%I #19 Aug 21 2016 11:07:47

%S 1,1,1,1,6,51,3001,9180001,14050074147451,3870680638643416483474006,

%T 4992392071450646411005278674572370014340582601,

%U 2715030052293379508289500941366397276374058263752394148988972928520177978202810359001

%N a(n) = (a(n-1)^2+a(n-2)^2+a(n-3)^2+a(n-1)*a(n-2)+a(n-1)*a(n-3)+a(n-2)*a(n-3))/a(n-4) with four initial ones.

%H Seiichi Manyama, <a href="/A165896/b165896.txt">Table of n, a(n) for n = 0..15</a>

%H S. Fomin, A. Zelevinsky, <a href="http://dx.doi.org/10.1006/aama.2001.0770">The Laurent Phenomenon</a>, Adv. Appl. Math. 28 (2) (2002) 119-144. [_R. J. Mathar_, Oct 23 2009]

%H Sergey Fomin, Andrei Zelevinsky, <a href="http://arxiv.org/abs/math/0104241">The Laurent phenomenon</a>, arXiv:math/0104241 [math.CO], 2001. [_R. J. Mathar_, Oct 23 2009]

%F a(n) ~ 1/sqrt(10) * c^(t^n), where t = A058265 = 1.8392867552141611325518525646532866..., c = 1.2712241060822553131735186905646486868228186258439... . - _Vaclav Kotesovec_, May 06 2015

%F a(n) = 10*a(n-1)*a(n-2)*a(n-3)-a(n-1)-a(n-2)-a(n-3)-a(n-4). - _Bruno Langlois_, Aug 21 2016

%t RecurrenceTable[{a[0]==a[1]==a[2]==a[3]==1,a[n]==(a[n-1]^2+a[n-2]^2+a[n-3]^2+ a[n-1]a[n-2]+ a[n-1]a[n-3]+a[n-2]a[n-3])/a[n-4]},a,{n,13}] (* _Harvey P. Dale_, May 21 2012 *)

%o (PARI) a(n)=if(n<4,1,(a(n-1)^2+a(n-2)^2+a(n-3)^2+a(n-1)*a(n-2)+a(n-1)*a(n-3)+a(n-2)*a(n-3))/a(n-4))

%Y Cf. A006720, A165903.

%K nonn

%O 0,5

%A _Jaume Oliver Lafont_, Sep 29 2009