约瑟夫环递归算法（C++）（初学者也能看懂逻辑分析）

题目：

n个人围成一圈（编号从1到n），从第1个人开始报数，报到m的人出列，从下一个人再重新报数，报到m的人出列，如此下去，直至所有人都出列。求最后一个出列的人的编号。

先给出核心代码：

#include <iostream>
using namespace std;
int josephus(int n, int m) 
{
	if(n == 1) return 0;
	else return ( josephus(n-1,m)+m ) % n;
}
 
int main() 
{
	int n, m;
	cin >> n >> m;
	int result = josephus(n, m);
	cout << result+1 << endl;
}

举例：n=9，m=4

初始	1	2	3	4	5	6	7	8	9
第一轮	1	2	3		5	6	7	8	9
第二轮	1	2	3		5	6	7		9
第三轮	1	2			5	6	7		9
第四轮	1	2			5	6	7
第五轮	1	2			5		7
第六轮	1	2					7
第七轮	1	2
第八轮	1

思路分析1：

1. 既然是递归，那么就是把复杂的问题一步一步分解为最基本、最简单的问题来解决。即把上例的n=9，逐步分解为n=8，n=7，n=6，……，n=2，n=1来解决。
2. 据我们观察，其中输入的变量有两个：一个是人数n，一个是标号数m。所以，假设我们求的递归函数就是$f(n,m)$。又标号数m是固定的，每次计数的时候都是这个标号。
3. 递归说明，前一个项和后一个项肯定是有关系的，不然递归进行不下去。即$f(n,m)$与$f(n-1,m)$之间有一定关系

那么思路便已经建立，我们需要找到$f(n,m)$与$f(n-1,m)$之间的关系即可，即找到n=9与n=8，n=8与n=7，……，n=3与n=2，n=2与n=1之间的关系。

既然是把复杂的问题一步一步分解为最基本、最简单的问题来解决，那我们便从n=1的问题来一步步逆向推回去。

思路分析2：

当n=1时，$f(1,m)$：

当n=1时，环中只有1个人（编号为1），那么最后一个出列的一定是编号为1的人。为了跟之后的推导统一，我们这里得出的result不记为1，而记为0。所以，当n=1时，result=0。

m	2	3	4	5	6	7	8	9	10	11	12
n=1，最后一个出列的人的编号	1	1	1	1	1	1	1	1	1	1	1

当n=2时，$f(2,m)$：

当n=2时，环中有2个人（编号为1，2），那么最后一个出列的人可能编号为1，也可能编号为2，这要看m的奇偶决定。如果m为奇数，那么最后一个出列的是编号为2的人；如果m为偶数，那么最后一个出列的人是编号为1的人。

m	2	3	4	5	6	7	8	9	10	11	12
n=2，最后一个出列的人的编号	1	2	1	2	1	2	1	2	1	2	1

当n=3时，$f(3,m)$：

当n=3时，环中有3个人（编号为1，2，3），那么最后一个出列的人可能编号为1，可能编号为2，也可能是编号为3。

m	2	3	4	5	6	7	8	9	10	11	12
n=3，最后一个出列的人的编号	3	2	2	1	1	3	3	2	2	1	1

可以这样想，当在m确定的情况下，则当n=3时，环中每个人依次报数，报到m的人就离开，那么报完一轮后，肯定走掉一个人，那么还剩2个人，即n=2。用数学的语言说就是从$f(3,m)$递推为$f(2,m)$了，即每进行一轮，n都减1。

那么问题思路更加清晰了，就是找$f(3,m)$和$f(2,m)$之间的递推关系，我们用最笨的方法（穷举法）来找规律：

先对比两张表格，为了方便看，我把上面两张表格合并了：

m	2	3	4	5	6	7	8	9	10	11	12
n=2，最后一个出列的人的编号	1	2	1	2	1	2	1	2	1	2	1
n=3，最后一个出列的人的编号	3	2	2	1	1	3	3	2	2	1	1

得出一张新的规律表格：(n=3确定)（编号为0即编号为3）

m	$f(n,m)$	$f(n-1,m)$	规律
2	$f(3,2)=3$	$f(2,2)=1$	$f(3,2)=[f(2,2)+2]\%n$
3	$f(3,3)=2$	$f(2,3)=2$	$f(3,3)=[f(2,3)+3]\%n$
4	$f(3,4)=2$	$f(2,4)=1$	$f(3,4)=[f(2,4)+4]\%n$
5	$f(3,5)=1$	$f(2,5)=2$	$f(3,5)=[f(2,5)+5]\%n$
6	$f(3,6)=1$	$f(2,6)=1$	$f(3,6)=[f(2,6)+6]\%n$
7	$f(3,7)=3$	$f(2,7)=2$	$f(3,7)=[f(2,7)+7]\%n$
8	$f(3,8)=3$	$f(2,8)=1$	$f(3,8)=[f(2,8)+8]\%n$
9	$f(3,9)=2$	$f(2,9)=2$	$f(3,9)=[f(2,9)+9]\%n$
10	$f(3,10)=2$	$f(2,10)=1$	$f(3,10)=[f(2,10)+10]\%n$
11	$f(3,11)=1$	$f(2,11)=2$	$f(3,11)=[f(2,11)+11]\%n$
12	$f(3,12)=1$	$f(2,12)=1$	$f(3,12)=[f(2,12)+12]\%n$

非常容易的观察得到，递推规律为$f(n,m)=[f(n-1,m)+m]\%n$
再把找到的这个规律代入$f(4,m)$和$f(3,m)$之间验证，发现也是正确的。即得出递推关系式。

为什么？

看到这里，你或许想问为什么递推关系式是$f(n,m)=[f(n-1,m)+m]\%n$，我连找规律都没有找出来，又怎么归纳总结递推关系式？

其实，这背后有这样一种思想：
还是拿刚才那个例子说事：
总人数n=9人，从编号为1的人开始报数，每报到4就把一人踢出去（m=4）。

初始	1	2	3	4	5	6	7	8	9
第一轮	1	2	3		5	6	7	8	9

此时，这些编号已经不能组成一个环，因为编号为4的地方产生了一个空位。之后的报数将总要考虑原编号4处的空位问题。

如何才能避免已经产生的空位对报数所造成的影响呢？

不过没有关系，因为下一次报数将从编号为5的人开始，我们可以将剩下的8个人组成一个新的环（5，6，7，8，9，1，2，3）。即，将3和5首尾相连，这样报数的时候就不用在意4的空位了。
但是新产生的环的数字并非连续的，报数时不像之前那样好处理了。

怎么处理新环数字非连续，报数无法简单用[(当前编号)%n]这个式子递推的问题？

所以现在我们必须借助存储结构得知下一个应该报数的现存人员编号。
接下来我们的目的是：使新环上的编号能够递推来简化我们之后的处理。
意思就是我们要改变新环上每个人的编号，从而来建立一种有确定规则的映射，达到映射之后数字可以递推的目的。且可以将在新环中继续按原规则报数得到的结果逆推出在旧环中的对应数字。

方法：将长度为n-1的新环与n-1个人组成的编号为1~n-1的环一 一映射。

之前的例子，将剩余的 8人与 8 人环（编号为0 ~ 8）一 一映射。

初始	1	2	3	4	5	6	7	8	9
第一轮	1	2	3		5	6	7	8	9
旧环	1	2	3		5	6	7	8	9
新环	6	7	8		1	2	3	4	5

注意，这里的旧环就是进行一轮（踢掉编号为4的人）之后的环，新环是按照上面的方法把剩下的8人进行重新编号（1~8），且让原来编号为5的人现在新的编号为1，原来编号为6的人现在编号为2（依此类推）。

这样新环的编号既解决了旧环编号不连续的问题，又解决了新编号与旧编号之间一一映射的关系，即$(旧的编号)=[(新的编号)+m]\%(旧的人数n)$

咦！这个式子好像有点眼熟？
对，就是刚才找规律总结出来的$f(n,m)=[f(n-1,m)+m]\%n$

那再看看这个规律适不适用于第二轮、第三轮……

初始	1	2	3	4	5	6	7	8	9
第一轮	1	2	3		5	6	7	8	9
新环	6	7	8		1	2	3	4	5
第二轮	6	7	8		1	2	3		5
新环	2	3	4		5	6	7		1
第三轮	2	3			5	6	7		1
新环	5	6			1	2	3		4
第四轮	5	6			1	2	3
新环	1	2			3	4	5
第五轮	1	2			3		5
新环	2	3			4		1
第六轮	2	3					1
新环	2	3					1
第七轮	2	3
新环	1	2
第八轮	1
新环	1

验证成功，每一轮与上一轮之间的关系都满足$f(n,m)=[f(n-1,m)+m]\%n$，况且我们之前找规律的那张表也可以充分说明这个关系式的可靠性。
综上所述，约瑟夫递归算法全部讲解完成。
再贴一次代码方便查看：

#include <iostream>
using namespace std;
int josephus(int n, int m) 
{
	if(n == 1) return 0;
	else return ( josephus(n-1,m)+m ) % n;
}
 
int main() 
{
	int n, m;
	cin >> n >> m;
	int result = josephus(n, m);
	cout << result+1 << endl;
}