%I #21 Sep 08 2022 08:45:55
%S 4,7,10,14,17,21,24,28,31,34,38,41,45,48,51,55,58,62,65,68,72,75,79,
%T 82,86,89,92,96,99,103,106,109,113,116,120,123,127,130,133,137,140,
%U 144,147,150,154,157,161,164,168,171,174,178,181,185,188,191,195,198,202,205,208,212,215,219,222,226,229,232,236,239,243,246,249,253,256,260,263,267,270,273,277,280,284,287,290,294,297,301,304,307,311,314,318,321,325,328,331,335,338,342
%N Upper s(n)-Wythoff sequence, where s(n) = 2n + 1.
%C See A184117 (the lower s(n)-Wythoff sequence).
%H Alois P. Heinz, <a href="/A184118/b184118.txt">Table of n, a(n) for n = 1..1000</a>
%F a(n) = A184117(n) + s(n) for all n. - _M. F. Hasler_, Jan 07 2019
%p a:=n->floor((2+sqrt(2))*n+sqrt(2)/2): seq(a(n),n=1..80); # _Muniru A Asiru_, Jan 08 2019
%t k=2; r=-1;
%t mex:=First[Complement[Range[1, Max[#1]+1], #1]]&;
%t S[n_]:=k n-r; A[1]=1; B[n_]:=B[n]=S[n]+A[n];
%t A[n_]:=A[n]=mex[Flatten[Table[{A[i], B[i]}, {i, 1, n-1}]]];
%t Table[S[n], {n, 30}]
%t Table[A[n], {n, 100}]
%t Table[B[n], {n, 100}]
%t Table[Floor[(2 + Sqrt[2]) n + Sqrt[2] / 2], {n, 80}] (* _Vincenzo Librandi_, Jan 07 2019 *)
%o (PARI) A184118_upto(N,s(n)=2*n+1,U=[0],b=[])={until(b[#b]>=N, b=concat(b,s(1+#b)+U[1]+=1); U=setunion(U,[b[#b]]); while(#U>1&&U[2]==U[1]+1,U=U[^1]));b} \\ _M. F. Hasler_, Jan 07 2019
%o (Magma) [Floor((2+Sqrt(2))*n+Sqrt(2)/2): n in [1..100]]; // _Vincenzo Librandi_, Jan 07 2019
%Y Cf. A000201, A001950, A001950, A001951, A136119, A184117.
%K nonn
%O 1,1
%A _Clark Kimberling_, Jan 09 2011
|