|
|
A274951
|
|
Pisot sequence E(8,12), a(n) = floor( a(n-1)^2/a(n-2) + 1/2 ).
|
|
1
|
|
|
8, 12, 18, 27, 41, 62, 94, 143, 218, 332, 506, 771, 1175, 1791, 2730, 4161, 6342, 9666, 14732, 22453, 34221, 52157, 79494, 121159, 184662, 281449, 428965, 653799, 996476, 1518761, 2314792, 3528048, 5377210, 8195577, 12491140, 19038144, 29016641, 44225186, 67405013, 102734125, 156580349
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,1
|
|
LINKS
|
|
|
MAPLE
|
f:=proc(n, x, y, r) option remember;
if n=0 then x
elif n=1 then y
else floor(f(n-1, x, y, r)^2/f(n-2, x, y, r) + r); fi; end;
[seq(f(n, 8, 12, 1/2), n=0..60)];
|
|
PROG
|
(Python)
a, b = 8, 12
for i in range(1000):
c, d = divmod(b**2, a)
a, b = b, c + (0 if 2*d < a else 1)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|