login
Maximum value on Lower Shuffle Part of Kimberling's Expulsion Array (A035486).
5

%I #20 Oct 19 2019 10:50:46

%S 1,3,5,7,10,12,15,17,20,22,25,28,31,33,36,39,42,44,47,50,53,55,58,61,

%T 64,67,70,73,76,78,81,84,87,90,93,96,99,101,104,107,110,113,116,119,

%U 122,124,127,130,133,136,139,142,145,148,151,154,157,160,163,166,169,171

%N Maximum value on Lower Shuffle Part of Kimberling's Expulsion Array (A035486).

%C a(n) is the maximum value on or below diagonal of Kimberling's Expulsion Array; this part could be called the Lower Shuffle.

%D D. Gale, Tracking the Automatic Ant: And Other Mathematical Explorations, ch. 5, p. 27. Springer, 1998

%D R. K. Guy, Unsolved Problems Number Theory, Sect E35.

%H Enrique Pérez Herrero, <a href="/A175312/b175312.txt">Table of n, a(n) for n=1..20000</a>

%H Clark Kimberling, <a href="https://cms.math.ca/crux/backfile/Crux_v17n02_Feb.pdf">Problem 1615</a>, Crux Mathematicorum, Vol. 17 (2) 44 1991; <a href="https://cms.math.ca/crux/backfile/Crux_v18n03_Mar.pdf">Solution to Problem 1615</a>, Crux Mathematicorum, Vol. 18, March 1992, p. 82-83.

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/KimberlingSequence.html">Kimberling Sequence</a>.

%F a(n) = 1 + 3(n-lambda(n)) - floor((n+2)/2^lambda(n)), lambda(n) = floor(log_2((n+2)/3)).

%F a(n) >= A007063(n); a(n) = max(K(n,1),K(n,2),...,K(n,n)), where K(i,j) is an element of Kimberling's Array given by A035486.

%F From _Enrique Pérez Herrero_, Mar 30 2010: (Start)

%F a(theta(k)) = A007063(theta(k)), where theta(k) = Sum_{i=0..k-1} 2^floor(i/3).

%F At these values the maximum in the Lower Shuffle is the diagonal expelled element. (End)

%t (* By direct computation *)

%t K[i_, j_] := i + j - 1 /; (j >= 2 i - 3);

%t K[i_, j_] := K[i - 1, i - (j + 2)/2] /; (EvenQ[j] && (j < 2 i - 3));

%t K[i_, j_] := K[i - 1, i + (j - 1)/2] /; (OddQ[j] && (j < 2 i - 3));

%t K[i_] := K[i] = K[i, i]; SetAttributes[K, Listable];

%t A175312[n_] := Max[Table[K[n, i], {i, 1, n}]] (* _Enrique Pérez Herrero_, Mar 30 2010 *)

%t (* By the Formula *)

%t \[Lambda][n_] := Floor[Log[2, (n + 2)/3]];

%t A175312[n_] := 1 + 3*(n - \[Lambda][n]) - Floor[(n + 2)/(2^\[Lambda][n]) (* _Enrique Pérez Herrero_, Mar 30 2010 *)

%o (PARI) lambda(n)= floor(log((n + 2)/3)/log(2));

%o A175312(n)= 1 + 3*(n - lambda(n)) - floor((n + 2)/(2^lambda(n))); \\ _Enrique Pérez Herrero_, Mar 30 2010

%Y Cf. A006852, A007063, A035486, A035505.

%K nonn

%O 1,2

%A _Enrique Pérez Herrero_, Mar 28 2010