Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #35 Dec 12 2024 09:25:24
%S 1,2,1,2,1,3,2,4,3,1,2,4,1,5,3,2,4,6,3,1,5,2,4,6,1,5,3,7,2,4,6,8,3,7,
%T 5,1,2,4,6,8,1,5,9,7,3,2,4,6,8,10,3,7,1,9,5,2,4,6,8,10,1,5,9,3,11,7,2,
%U 4,6,8,10,12,3,7,11,5,1,9,2,4,6,8,10,12,1,5,9,13,7,3,11,2,4,6,8,10,12,14
%N Triangle read by rows: T(n,k) is the number of the k-th eliminated person in the Josephus elimination process for n people and a count of 2, 1 <= k <= n.
%C In the Josephus elimination process for n and k, the numbers 1 through n are written in a circle. A pointer starts at position 1. Each turn, the pointer skips (k-1) non-eliminated number(s) going around the circle and eliminates the k-th number, until no numbers remain. This sequence represents the triangle J(n, i), where n is the number of people in the circle, i is the turn number, and k is fixed at 2 (every other number is eliminated).
%H <a href="/index/J#Josephus">Index entries for sequences related to the Josephus Problem</a>
%F From _Pontus von Brömssen_, Sep 18 2022: (Start)
%F The terms are uniquely determined by the following recursive formulas:
%F T(n,k) = 2*k if k <= n/2;
%F T(2*n,k) = 2*T(n,k-n)-1 if k > n;
%F T(2*n+1,k) = 2*T(n,k-n-1)+1 if k > n+1;
%F T(2*n+1,n+1) = 1.
%F (End)
%F From _Pontus von Brömssen_, Dec 11 2024: (Start)
%F The terms are also uniquely determined by the following recursive formulas:
%F T(1,1) = 1;
%F T(n,1) = 2 if n > 1;
%F T(n,k) = T(n-1,k-1)+2 if k > 1 and T(n-1,k-1) != n-1;
%F T(n,k) = 1 if k > 1 and T(n-1,k-1) = n-1.
%F T(n,A225381(n)) = 1.
%F T(n,A225381(n+1)-1) = n.
%F (End)
%e Triangle begins:
%e 1;
%e 2, 1;
%e 2, 1, 3;
%e 2, 4, 3, 1;
%e 2, 4, 1, 5, 3;
%e 2, 4, 6, 3, 1, 5;
%e 2, 4, 6, 1, 5, 3, 7;
%e 2, 4, 6, 8, 3, 7, 5, 1;
%e 2, 4, 6, 8, 1, 5, 9, 7, 3;
%e 2, 4, 6, 8, 10, 3, 7, 1, 9, 5;
%e 2, 4, 6, 8, 10, 1, 5, 9, 3, 11, 7;
%e 2, 4, 6, 8, 10, 12, 3, 7, 11, 5, 1, 9;
%e 2, 4, 6, 8, 10, 12, 1, 5, 9, 13, 7, 3, 11;
%e ...
%e For n = 5, to get the entries in 5th row from left to right, start with (^1, 2, 3, 4, 5) and the pointer at position 1, indicated by the caret. 1 is skipped and 2 is eliminated to get (1, ^3, 4, 5). (The pointer moves ahead to the next "live" number.) On the next turn, 3 is skipped and 4 is eliminated to get (1, 3, ^5). Then 1, 5, and 3 are eliminated in that order (going through (^3, 5) and (^3)). This gives row 5 of the triangle and entries a(11) through a(15) in this sequence.
%t Table[Rest@ Nest[Append[#1, {Delete[#2, #3 + 1], #2[[#3 + 1]], #3}] & @@ {#, #[[-1, 1]], Mod[#[[-1, -1]] + 1, Length@ #[[-1, 1]]]} &, {{Range@ n, 0, 0}}, n][[All, 2]], {n, 14}] // Flatten (* _Michael De Vlieger_, Nov 13 2018 *)
%o (Python)
%o def A321298(n,k):
%o if 2*k<=n: return 2*k
%o n2,r=divmod(n,2)
%o if r==0: return 2*A321298(n2,k-n2)-1
%o if k==n2+1: return 1
%o return 2*A321298(n2,k-n2-1)+1 # _Pontus von Brömssen_, Sep 18 2022
%Y The right border of this triangle is A006257.
%Y Cf. A032434, A054995, A181281, A225381.
%K easy,nonn,tabl,changed
%O 1,2
%A _Zeph L. Turner_, Nov 02 2018
%E Name clarified by _Pontus von Brömssen_, Sep 18 2022