%I #34 Sep 08 2022 08:44:45
%S 2,7,24,82,280,956,3264,11144,38048,129904,443520,1514272,5170048,
%T 17651648,60266496,205762688,702517760,2398545664,8189147136,
%U 27959497216,95459694592,325919783936,1112759746560,3799199418368,12971278180352,44286713884672
%N Pisot sequence P(2,7): a(0)=2, a(1)=7, thereafter a(n+1) is the nearest integer to a(n)^2/a(n-1).
%C Also Pisot sequence T(2,7). - _R. K. Guy_
%C It appears that a(n) = 4*a(n-1) - 2*a(n-2) (holds at least up to n = 1000 but is not known to hold in general).
%C The recurrence holds up to n = 10^5. - _Ralf Stephan_, Sep 03 2013
%C Empirical g.f.: (2-x)/(1-4*x+2*x^2). - _Colin Barker_, Feb 21 2012
%H Colin Barker, <a href="/A020727/b020727.txt">Table of n, a(n) for n = 0..1000</a>
%t RecurrenceTable[{a[0] == 2, a[1] == 7, a[n] == Floor[a[n - 1]^2/a[n - 2]]}, a, {n, 0, 30}] (* _Bruno Berselli_, Feb 04 2016 *)
%o (Magma) Iv:=[2,7]; [n le 2 select Iv[n] else Floor(Self(n-1)^2/Self(n-2)): n in [1..30]]; // _Bruno Berselli_, Feb 04 2016
%o (PARI) pisotP(nmax, a1, a2) = {
%o a=vector(nmax); a[1]=a1; a[2]=a2;
%o for(n=3, nmax, a[n] = ceil(a[n-1]^2/a[n-2]-1/2));
%o a
%o }
%o pisotP(50, 2, 7) \\ _Colin Barker_, Aug 08 2016
%Y It appears that this is a subsequence of A003480.
%Y See A008776 for definitions of Pisot sequences.
%K nonn
%O 0,1
%A _David W. Wilson_
%E Edited by _N. J. A. Sloane_, Aug 17 2009 at the suggestion of _R. J. Mathar_.