算法总结——约瑟夫问题

问题描述

约瑟夫问题：有ｎ只猴子，按顺时针方向围成一圈选大王（编号从１到ｎ），从第１号开始报数，一直数到ｍ，数到ｍ的猴子退出圈外，剩下的猴子再接着从1开始报数。就这样，直到圈内只剩下一只猴子时，这个猴子就是猴王，编程求输入ｎ，ｍ后，输出最后猴王的编号。 
输入数据
每行是用空格分开的两个整数，第一个是 n, 第二个是 m ( 0 < m, n < 300)。最后一行是： 
0 0 
输出要求
对于每行输入数据（最后一行除外)，输出数据也是一行，即最后猴王的编号 
输入样例：
6 2
12 4
8 3
0 0
输出样例：
5
1
7

解题思路

初一看，很可能想把这道题目当作数学题来做，即认为结果也许会是以n和m为自变量的某个函数f(n,m)，只要发现这个函数，问题就迎刃而解。实际上，这样的函数很难找，甚至也许根本就不存在。用人工解决的办法就是将 n 个数写在纸上排成一圈，然后从1开始数，每数到第 m 个就划掉一个数，一遍遍做下去，直到剩下最后一个。有了计算机，这项工作做起来就会快多了，我们只要编写一个程序，模拟人工操作的过程就可以了。

用数组anLoop来存放n个数，相当于n个数排成的圈；用整型变量 nPtr 指向当前数到的数组元素，相当于人的手指；划掉一个数的操作，就用将一个数组元素置0的方法来实现。人工数的时候，要跳过已经被划掉的数，那么程序执行的时候，就要跳过为0的数组元素。需要注意的是，当nPtr指向anLoop中最后一个元素（下标n-1）时，再数下一个，则nPtr要指回到数组的头一个元素（下标0），这样anLoop才象一个圈。

参考程序：

#include <stdio.h>
#include <stdlib.h>
#define MAX_NUM 300
int aLoop[MAX_NUM + 10];
main()
{
	int n, m, i;
	while(1) {
		scanf( "%d%d", & n, & m);
		if ( n == 0 ) 
			break;
		for( i = 0; i < n; i ++ ) 
			aLoop[i] = i+1;
		int nPtr = 0;
		for( i = 0; i < n; i ++ ) {  //每次循环将1个猴子全赶出圈子，
									 //最后被赶出的就是猴王
			int nCounted = 0;	
			while( nCounted < m ) {  //数出m个猴子
				while( aLoop[nPtr] == 0)  //跳过已经出圈的猴子
					nPtr = ( nPtr + 1 ) % n;  //到下一个位置
				nCounted ++;        //找到一只猴子
				nPtr = ( nPtr + 1 ) % n; //到下一个位置
			}
			nPtr --;   //要回退一个位置
			if( nPtr < 0 )
				nPtr = n - 1;
			if( i == n-1 )  //最后一只出圈的猴子
				printf("%d\n", aLoop[nPtr]);  
			aLoop[nPtr] = 0; //猴子出圈
		}
	}
}

上面的程序完全模拟了人工操作的过程，但因为要反复跳过为0的数组元素，因此算法的效率不是很高。后文的“链表”一章，采用单链表进行模拟来解决本题，就能省去跳过已出圈的猴子这个操作，大大提高了效率。
实现技巧
n个元素的数组，从下标0的元素开始存放猴子编号，则循环报数的时候，下一个猴子的下标就是 “(当前猴子下标 + 1 )% n”。这种写法比用分支语句来决定下个猴子的下标是多少，更快捷而且写起来更方便。

常见问题

问题一：在数组里循环计数的时候，一定要小心计算其开始的下标和终止的下标。比如，语句15，循环是从0到n-1，而不是从0到n。
问题二：语句24到26回退一个位置，易被忽略或写错。比如只写了语句24，忘了处理nPtr变成小于0的情况。