|
%I #31 Jan 06 2019 15:59:17
%S 4,7,12,21,37,65,114,200,351,616,1081,1897,3329,5842,10252,17991,
%T 31572,55405,97229,170625,299426,525456,922111,1618192,2839729,
%U 4983377,8745217,15346786,26931732,47261895,82938844,145547525,255418101,448227521,786584466
%N Pisot sequences E(4,7), P(4,7).
%C Essentially the same as A005251: a(n) = A005251(n+5).
%C See A008776 for definitions of Pisot sequences.
%H Colin Barker, <a href="/A010901/b010901.txt">Table of n, a(n) for n = 0..1000</a>
%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_03">Index entries for linear recurrences with constant coefficients</a>, signature (2, -1, 1).
%F a(n) = 2a(n-1) - a(n-2) + a(n-3) for n>=3. (Proved using the PtoRv program of Ekhad-Sloane-Zeilberger.) - _N. J. A. Sloane_, Sep 09 2016
%t LinearRecurrence[{2, -1, 1}, {4, 7, 12}, 35] (* _Jean-François Alcover_, Oct 05 2018 *)
%o (PARI) pisotE(nmax, a1, a2) = {
%o a=vector(nmax); a[1]=a1; a[2]=a2;
%o for(n=3, nmax, a[n] = floor(a[n-1]^2/a[n-2]+1/2));
%o a
%o }
%o pisotE(50, 4, 7) \\ _Colin Barker_, Jul 27 2016
%Y Cf. A005251, A008776, A010925.
%K nonn
%O 0,1
%A _Simon Plouffe_
%E Edited by _N. J. A. Sloane_, Jul 26 2016
|
