A088461 a(1) = 1, a(n) = ceiling((n-a(a(a(n-1))))/2).

1, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8, 8, 8, 9, 9, 10, 10, 11, 11, 12, 12, 12, 13, 13, 14, 14, 15, 15, 16, 16, 16, 17, 17, 18, 18, 19, 19, 20, 20, 21, 21, 21, 22, 22, 23, 23, 24, 24, 25, 25, 25, 26, 26, 27, 27, 28, 28, 29, 29, 30, 30, 30, 31, 31, 32, 32, 33, 33, 34, 34, 34

Offset

1,4

Links

Table of n, a(n) for n=1..76.

Formula

a(n)=c*n+O(1) where c is the positive root of x^3+2*x-1=0, c=0.453397651516403767644746538...

Maple

a:= proc(n) a(n):= ceil((n-a(a(a(n-1))))/2) end: a(1):=1:
seq(a(n), n=1..100);  # Alois P. Heinz, Sep 23 2023

Crossrefs

Cf. A005374.

Sequence in context: A264810 A176841 A176814 * A135020 A242681 A034888

Adjacent sequences: A088458 A088459 A088460 * A088462 A088463 A088464

Keyword

nonn

Author

Benoit Cloitre, Nov 12 2003

Status

approved