login
a(n) = (1/2)*(s(n) - 1), where s = A026139.
1

%I #8 Aug 28 2021 02:53:12

%S 0,1,3,4,7,9,10,12,13,16,19,21,22,25,27,28,30,31,34,36,37,39,40,43,46,

%T 48,49,52,55,57,58,61,63,64,66,67,70,73,75,76,79,81,82,84,85,88,90,91,

%U 93,94,97,100,102,103,106,108,109,111,112,115,117,118

%N a(n) = (1/2)*(s(n) - 1), where s = A026139.

%H Michael De Vlieger, <a href="/A026140/b026140.txt">Table of n, a(n) for n = 1..5004</a>

%t Block[{a, r, s, nn = 119}, a[1] = 1; Do[If[! IntegerQ[a[#1]], Set[a[#1], i], Set[a[#2], i]] & @@ {i - #, i + #} &@ Floor[i/2], {i, 2 nn}]; s = TakeWhile[Array[a[#] &, 2 nn], IntegerQ]; Map[(1/2) (# - 1) &, Union@ FoldList[Max, s]]] (* _Michael De Vlieger_, Aug 27 2021 *)

%Y Cf. A026136, A026139.

%K nonn

%O 1,3

%A _Clark Kimberling_