# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/
Search: id:a014004
Showing 1-1 of 1
%I A014004 #13 Sep 22 2017 15:42:41
%S A014004 9,17,32,60,113,213,401,755,1422,2678,5043,9497,17885,33682,63432,
%T A014004 119459,224972,423680,797898,1502646,2829867,5329364,10036557,
%U A014004 18901407,35596190,67036742,126247353,237756097,447755619,843238499,1588034044,2990674795,5632206541
%N A014004 Pisot sequence E(9,17), a(n) = floor( a(n-1)^2/a(n-2) + 1/2 ).
%H A014004 Colin Barker, Table of n, a(n) for n = 0..1000
%H A014004 D. W. Boyd, Some integer sequences related to the Pisot sequences, Acta Arithmetica, 34 (1979), 295-305.
%H A014004 D. W. Boyd, Linear recurrence relations for some generalized Pisot sequences, Advances in Number Theory ( Kingston ON, 1991) 333-340, Oxford Sci. Publ., Oxford Univ. Press, New York, 1993.
%F A014004 Known not to satisfy any linear recurrence.
%t A014004 nxt[{a_,b_}]:={b,Floor[b^2/a+1/2]}; NestList[nxt,{9,17},40][[All,1]] (* _Harvey P. Dale_, Sep 22 2017 *)
%o A014004 (PARI) pisotE(nmax, a1, a2) = {
%o A014004 a=vector(nmax); a[1]=a1; a[2]=a2;
%o A014004 for(n=3, nmax, a[n] = floor(a[n-1]^2/a[n-2]+1/2));
%o A014004 a
%o A014004 }
%o A014004 pisotE(50, 9, 17) \\ _Colin Barker_, Jul 28 2016
%K A014004 nonn
%O A014004 0,1
%A A014004 _Simon Plouffe_
# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE