约瑟夫环问题的数学分析解

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

        约瑟夫问题的解法可以用循环链表模拟过程解（《约瑟夫环问题一解》）或者用队列模拟过程求解，不过这个问题在数学上有确定的解，文章《约瑟夫环问题数学解法》中有这个问题的数学分析结果：

      问题描述：n个人（编号0~(n-1))，从0开始报数，报到m-1的退出
　　，剩下的人继续从0开始报数。求胜利者的编号。
　　我们知道第一个人(编号一定是(m-1)%n) 出列之后，剩下的n-1个人组成了一个新的约瑟夫环（以编号为k=m%n的人开始）:
　　k k+1 k+2 ... n-2, n-1, 0, 1, 2, ... k-2
　　并且从k开始报0。
　　现在我们把他们的编号做一下转换：
　　k --> 0
　　k+1 --> 1
　　k+2 --> 2
　　...
　　...
　　k-3 --> n-3
　　k-2 --> n-2
　　序列1： 0, 1, 2, 3 … n-2, n-1
　　序列2： 0, 1, 2, 3 … k-1, k+1, …, n-2, n-1
　　序列3： k, k+1, k+2, k+3, …, n-2, n-1, 1, 2, 3,…, k-2,
　　序列4：0, 1, 2, 3 …, 5, 6, 7, 8, …, n-3, n-2
　　变换后就完完全全成为了(n-1)个人报数的子问题，假如我们知道这个子问题的解：例如x是最终的胜利者，那么根据上面这个表把这个x变回去不刚好就是n个人情况的解               吗？！！变回去的公式很简单，相信大家都可以推出来：
　　∵ k=m%n;
　　∴ x' = x+k = x+ m%n ; 而 x+ m%n 可能大于n
　　∴x'= (x+ m%n)%n = (x+m)%n
　　得到 x‘=(x+m)%n
　　如何知道(n-1)个人报数的问题的解？对，只要知道(n-2)个人的解就行了。(n-2)个人的解呢？当然是先求(n-3)的情况 ---- 这显然就是一个倒推问题！好了，思路出来了，下           面写递推公式：
　　令f表示i个人玩游戏报m退出最后胜利者的编号，最后的结果自然是f[n].
　　递推公式:
　　f[1]=0;
　　f[i]=(f[i-1]+m)%i; (i>1)
　　有了这个公式，我们要做的就是从1-n顺序算出f的数值，最后结果是f[n]。我们输出f[n]由于是逐级递推，不需要保存每个，程序也是异常简单：

int Josephus(int m,int n)
{
	int i,s;
	s = 0;
	for(i = 2;i <= n; i++)
	{ 
		s = (s + m) % i;	
	} 
	return s + 1; //0-(n-1) => 1-n 
} 

测试程序：

#include <iostream>

using namespace std;

int Josephus(int m,int n)
{
	int i,s;
	s = 0;
	for(i = 2;i <= n; i++)
	{ 
		s = (s + m) % i;	
	} 
	return s + 1; //0-(n-1) => 1-n 
} 

int Josephus2(int m,int n)
{
	int i,s;
	s = 1;
	for(i = 2;i <= n; i++)
	{ 
		s = (s + m) % i;
		if (s == 0)  s = i;		
	} 
	return s;  
} 
int main()
{
	int m,n;
	m = 2;
	n = 7;
	cout<<"Josephus m = "<<m<<" n = "<<n<<" :"<<Josephus2(m,n)<<endl; 
}

输出：

Josephus m = 2 n = 7 :7

REF：

1，http://blog.csdn.net/fzydao/article/details/7041541