login
Pisot sequence E(4,19), a(n) = floor( a(n-1)^2/a(n-2)+1/2 ).
2

%I #32 Jul 14 2023 14:48:19

%S 4,19,90,426,2016,9541,45154,213697,1011348,4786332,22651920,

%T 107203069,507352048,2401107571,11363544486,53779407822,254517831936,

%U 1204537747753,5700626846950,26978935702753,127681216679304,604267465267128,2859772009358880,13534231802298265,64052459384483260,303136344428812723,1434630991482656082,6789572149788327282

%N Pisot sequence E(4,19), a(n) = floor( a(n-1)^2/a(n-2)+1/2 ).

%D Shalosh B. Ekhad, N. J. A. Sloane and Doron Zeilberger, Automated Proof (or Disproof) of Linear Recurrences Satisfied by Pisot Sequences, Preprint, 2016.

%H Vincenzo Librandi, <a href="/A010907/b010907.txt">Table of n, a(n) for n = 0..200</a>

%H D. W. Boyd, <a href="http://matwbn.icm.edu.pl/ksiazki/aa/aa34/aa3444.pdf">Some integer sequences related to the Pisot sequences</a>, Acta Arithmetica, 34 (1979), 295-305

%H D. W. Boyd, <a href="https://www.researchgate.net/profile/David_Boyd7/publication/262181133_Linear_recurrence_relations_for_some_generalized_Pisot_sequences_-_annotated_with_corrections_and_additions/links/00b7d536d49781037f000000.pdf">Linear recurrence relations for some generalized Pisot sequences</a>, Advances in Number Theory ( Kingston ON, 1991) 333-340, Oxford Sci. Publ., Oxford Univ. Press, New York, 1993.

%H S. B. Ekhad, N. J. A. Sloane, D. Zeilberger, <a href="http://arxiv.org/abs/1609.05570">Automated proofs (or disproofs) of linear recurrences satisfied by Pisot Sequences</a>, arXiv:1609.05570 [math.NT] (2016)

%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (4, 3, 2, 1).

%F Theorem: a(n) = 4 a(n - 1) + 3 a(n - 2) + 2 a(n - 3) + a(n - 4). (Proved using the PtoRv program of Ekhad-Sloane-Zeilberger.) - _N. J. A. Sloane_, Sep 09 2016

%F G.f.: -(x^3+2*x^2+3*x+4)/(x^4+2*x^3+3*x^2+4*x-1). [_Colin Barker_, Nov 29 2012] (This follows from the above recurrence. - _N. J. A. Sloane_, Sep 09 2016)

%t PSE[a_,b_,n_]:=Join[{x=a,y=b}, Table[z=Floor[y^2/x+1/2]; x=y; y=z, {n}]]; A010907=PSE[4,19,20] (* _Zak Seidov_, Mar 24 2011 *)

%t nxt[{a_,b_}]:={b,Floor[b^2/a+1/2]}; Transpose[NestList[nxt,{4,19},20]] [[1]] (* _Harvey P. Dale_, Mar 13 2016 *)

%Y Cf. A077922.

%K nonn

%O 0,1

%A _Simon Plouffe_