约瑟夫环的O(n)解决算法

问题 ：已知n个人（以编号0，1，2，3...n-1分别表示）围坐在一张圆桌周围。从编号为k的人开始报数，数到m的那个人出列；他的下一个人又从1开始报数，数到m的那个人又出列；依此规律重复下去，直到圆桌周围的人全部出列，求最后一个出列人的编号。

分析 ：
1.    考虑当队列中还有Z个人，我们按照他们报数的先后顺序来排列，生成如下左序列。
2.    这时，剔除第m个人时，我们按照下一次的报数先后顺序列出每个人员，m%(Z-1)则为下次第一个报名的人，生成序列如下中所示。
3.    接下来，我们重新对他们的报数顺序进行映射，生成如下右序列。
我们可以确信，对于中间序列最终剩下来的人，会和右边序列最终剩下来的人是同一个人。只不过在两个序列中，这个人的下标不一样，假设该人在右边的下标（Z-1个人时的下标）为x，则在中间的下标（Z个人时的下标）为(x+m)%Z，如果右边最后一个剩下的人下标为f(Z-1)，则中间最后剩下的人的下标为f(Z) = (f(Z-1) + m) % Z。

由此，我们得出了递推公式，f(Z) = (f(Z - 1) + m) % Z。

 

---

 

 

第一轮从第0号人开始

所以，如果刚开始，报数从第0个人开始，则程序如下：

/** * @brief Josephus problem * * @param n: how many people * @param m: the m'th person will be checked out * * @return The last person in the circle */ int Josephus(int n, int m) { //when there is only one person, the last person's index is 0 int s = 0; for (int i = 2; i <= n; ++i) { s = (s + m) % i; } return s; }

 

---

 

 

第一轮从第k号人开始

如果一开始从下标为k的人开始报数，那么相当于在初始阶段进行一次下标转换，如下图所示：

假设后者的下坐标为x，则它在左边的下坐标为(x + k) % n。
所以程序变为：
/** * @brief Josephus problem start at k * * @param n: how many people * @param m: the m'th person will be checked out * @param k: start at k in the first round * * @return The last person in the circle */ int JosephusK(int n, int m, int k) { //when there is only one person, the last person's index is 0 int s = 0; for (int i = 2; i <= n; ++i) { s = (s + m) % i; } return (s + k) % n; }