约瑟夫问题:
据说著名犹太历史学家 Josephus有过以下的故事:在罗马人占领乔塔帕特后,39 个犹太人与Josephus及他的朋友躲到一个洞中,39个犹太人决定宁愿死也不要被敌人抓到,于是决定了一个自杀方式,41个人排成一个圆圈,由第1个人开始报数,每报数到第3人该人就必须自杀,然后再由下一个重新报数,直到所有人都自杀身亡为止。然而Josephus 和他的朋友并不想遵从,Josephus要他的朋友先假装遵从,他将朋友与自己安排在第16个与第31个位置,于是逃过了这场死亡游戏。
编了一段程序,模拟约瑟夫问题,如下:
#define MAX 41#define interval 3
#include <stdio.h>
//#include <stdlib.h>
char A[MAX];
int main()
{
int i=0;
int jishu=0;
int totaldead=0;
int ptr=-1;
while(totaldead<MAX-1)
{
while(jishu<interval)
{
ptr=(++ptr)%MAX;
if(A[ptr]==0)
jishu++;
}
A[ptr]=1;
totaldead++;
jishu=0;
}
for(i=0;i<MAX;i++)
{
if(A[i]==0)
printf("%d\n",i+1);
}
return 0;
}
这个方法仅仅是模拟,还有更加简便的方法。
第一种方法,我们可以通过从后往前看,在整个过程的最后一步会发生什么事情?整个队列只剩下一个人没有被杀,这个人的编号是1,于是不管隔多少个人(interval)杀一次,这个人总会把自己给杀了,当然这是违背约瑟夫的本意,他会跑掉。最后一局的前一局,发生了什么?两个人包括约瑟夫在对峙,结果是约瑟夫胜出,那么约瑟夫在这一局的编号是多少?我约瑟夫最后一局之所以能够编号1是因为倒数第二局,循环数数中我逻辑前面的那个人的数字是interval,因此我理应的序号就是(interval+1)%2。得到了倒数第二局我应该有的序列号。这是这种方法的引导,下面开始详细讨论这种方法。
假设第X局被杀的人是k(第k局还有n个人活着),那么开始对剩下的人的编号更新如下:
k+1-->1
k+2-->2
.......
k-2-->n-2
k-1-->n-1
可见,k局过后剩下了n-1个人且都进行了重新排序。可以通过后一局的编号(a)推导出前一局的编号(b):b=(a+k)%n。
k是什么?k就是说interval的那个人,也就是这个人报的那个数。因此上面的式子还可以写成b=(a+interval)%n。
推导出来了。
下面考虑一种情况就是 剩下的人数n<interval 了,看看红线部分是否成立。
设n=3,interval=4
第一局:1号被杀,编号更新如下:
2-->1
3-->2
但是根据b=(a+interval)%n式子,临界值求值出现了问题:0=(2+4)%3。但是并没有设置0号编号,为什么会出现这种情况呢?
仔细想一下:
这是由于运算符%本身的性质引起的,因为%的结果只会是0--n-1,并不能是n,所以上例中,如果是本应输出3,但是却错误地输出了0,因此可以做如下修改:
假设某一局中有n个人,那么将这些人从0到n-1编号。
会这样修改问题:假设第k-1号人物被杀,那么这个k也就是所谓的interval,因为从0开始嘛。如此更新上面红线部分:
k-->0
k+1-->1
....
k-3-->n-3
k-2-->n-2
得到同样地规律:b=(a+interval)%n。同样考虑n<interval 的情况:
设n=3,interval=4。三个人的编号是{0,1,2}
第一局:1号被杀,编号更新如下:
1-->0
2-->1
根据式子b=(a+interval)%n,可得(1+4)%3=2;(0+4)%3=1;这次就满足情况 了。
可以如此编程了:
number=0;
for (i=2; i<=n; i++)
number=(number+interval)%i;
测试成功。
另一种方法是相对于interval=2实现的二进制左循环算法。
这种算法如果是7个人的话那么最后留下的一定是第7个人(111->111),六个人最后留下的是第5个人(110->101),8个人最后留下的是第一个人(1000->0001).......
很强大啊,想想吧为什么?
http://zh.wikipedia.org/wiki/%E7%BA%A6%E7%91%9F%E5%A4%AB%E6%96%AF%E9%97%AE%E9%A2%98好点看看,这个解释太牛逼了