|
%I #18 Dec 06 2023 13:57:24
%S 4,14,49,171,596,2077,7238,25223,87897,306303,1067403,3719680,
%T 12962320,45171020,157411717,548547468,1911575138,6661446313,
%U 23213770727,80895217952,281903201529,982374694626,3423373822671,11929753885009,41572739387791,144872448909191,504850696923520,1759300875378480
%N Pisot sequence T(4,14).
%H <a href="/index/Ph#Pisot">Index entries for Pisot sequences</a>
%F a(n) = floor(a(n-1)^2/a(n-2)), a(0) = 4, a(1) = 14.
%F Conjectures: (Start)
%F G.f.: (4 - 2*x + x^2 - x^3)/(1 - 4*x + 2*x^2 - x^3 + x^4).
%F a(n) = 4*a(n-1) - 2*a(n-2) + a(n-3) - a(n-4). (End)
%t RecurrenceTable[{a[0] == 4, a[1] == 14, a[n] == Floor[a[n - 1]^2/a[n - 2]]}, a, {n, 27}]
%o (PARI) first(n)=my(v=vector(n+1)); v[1]=4; v[2]=14; for(i=3,#v, v[i]=v[i-1]^2\v[i-2]); v \\ _Charles R Greathouse IV_, Nov 28 2016
%o (Python)
%o from itertools import islice
%o def A278692_gen(): # generator of terms
%o a, b = 4, 14
%o yield from (a,b)
%o while True:
%o a, b = b, b**2//a
%o yield b
%o A278692_list = list(islice(A278692_gen(),30)) # _Chai Wah Wu_, Dec 06 2023
%Y Cf. A008776 for definitions of Pisot sequences.
%Y Cf. A010919, A019495, A022031.
%Y Cf. A010904 (Pisot sequence E(4,14)), A251221 (seems to be Pisot sequence P(4,14)), A277084 (Pisot sequence L(4,14)).
%K nonn,easy
%O 0,1
%A _Ilya Gutkovskiy_, Nov 28 2016
|
