约瑟夫问题求解

约瑟夫问题：

据说著名犹太历史学家 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

好点看看，这个解释太牛逼了