约瑟夫问题

约瑟夫问题（有时也称为约瑟夫斯置换，是一个出现在计算机科学和数学中的问题。在计算机编程的算法中，类似问题又称为约瑟夫环。又称“丢手绢问题”.）
【问题来历】
    据说著名犹太历史学家Josephus有过以下的故事：在罗马人占领乔塔帕特后，39 个犹太人与Josephus及他的朋友躲到一个洞中，39个犹太人决定宁愿死也不要被敌人抓到，于是决定了一个自杀方式，41个人排成一个圆圈，由第1个人开始报数，每报数到第3人该人就必须自杀，然后再由下一个重新报数，直到所有人都自杀身亡为止。然而Josephus 和他的朋友并不想遵从。首先从一个人开始，越过k-2个人（因为第一个人已经被越过），并杀掉第k个人。接着，再越过k-1个人，并杀掉第k个人。这个过程沿着圆圈一直进行，直到最终只剩下一个人留下，这个人就可以继续活着。问题是，给定了和，一开始要站在什么地方才能避免被处决？Josephus要他的朋友先假装遵从，他将朋友与自己安排在第16个与第31个位置，于是逃过了这场死亡游戏。
    17世纪的法国数学家加斯帕在《数目的游戏问题》中讲了这样一个故事：15个教徒和15 个非教徒在深海上遇险，必须将一半的人投入海中，其余的人才能幸免于难，于是想了一个办法：30个人围成一圆圈，从第一个人开始依次报数，每数到第九个人就将他扔入大海，如此循环进行直到仅余15个人为止。问怎样排法，才能使每次投入大海的都是非教徒。
*问题分析与算法设计
    约瑟夫问题并不难，但求解的方法很多；题目的变化形式也很多。这里给出一种实现方法。
题目中30个人围成一圈，因而启发我们用一个循环的链来表示，可以使用结构数组来构成一个循环链。结构中有两个成员，其一为指向下一个人的指针，以构成环形的链；其二为该人是否被扔下海的标记，为1表示还在船上。从第一个人开始对还未扔下海的人进行计数，每数到9时，将结构中的标记改为0，表示该人已被扔下海了。这样循环计数直到有15个人被扔下海为止。
约瑟夫问题的一般描述为：现在有n个人围成一个圆圈，编号分别为1,2,...,n，现在从第一个人开始报数1，他后面的人接着报数2，...,当出现第一个人报数为m时，该人退出圆圈，该人后面的人从1开始重新报数，一直循环到只剩下一个人。
约瑟夫问题一般有两种解法：
解法一：直接模拟约瑟夫环的操作，建立一个循环单链表，然后模拟操作。
代码（我这里用next[]数组作为单链表，然后get()函数每次返回的是下一次数"1"的人）

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
int next[110] , n , m;
void init() {
    for(int i=1;i<n;i++) next[i] = i + 1;
    next[n] = 1;
}
int get(int p ， int m) {
    int q;
    for(int i=0;i<m;i++) {
        q = p;
        p = next[p];
    }
    next[q] = next[p];
    return next[p];
}
int main() {
    while(scanf("%d%d" , &n , &m) != EOF) {
        init();
        int p = 1;
        for(int i=0;i<n;i++) p = get(p);
        printf("%d\n" , p);
    }
    return 0;
}

这种算法的时间复杂度是O(n*m)。
解法二：无论是用链表实现还是用数组实现都有一个共同点：要模拟整个游戏过程，不仅程序写起来比较烦，而且时间复杂度高达O(n*m)，当n，m非常大(例如上百万，上千万)的时候，几乎是没有办法在短时间内出结果的。我们注意到原问题仅仅是要求出最后的胜利者的序号，而不是要读者模拟整个过程。因此如果要追求效率，就要打破常规，实施一点数学策略。


为了讨论方便，先把问题稍微改变一下，并不影响原意：
问题描述：n个人（编号0~(n-1))，从0开始报数，报到(m-1)的退出，剩下的人继续从0开始报数。求胜利者的编号。


我们知道第一个人(编号一定是m%n-1) 出列之后，剩下的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-2   --> n-2
k-1   --> n-1
变换后就完完全全成为了(n-1)个人报数的子问题，假如我们知道这个子问题的解：例如x是最终的胜利者，那么根据上面这个表把这个x变回去不刚好就是n个人情况的解吗？！！变回去的公式很简单，相信大家都可以推出来：x'=(x+k)%n


如何知道(n-1)个人报数的问题的解？对，只要知道(n-2)个人的解就行了。(n-2)个人的解呢？当然是先求(n-3)的情况 ---- 这显然就是一个倒推问题！好了，思路出来了，下面写递推公式：


令f[i]表示i个人玩游戏报m退出最后胜利者的编号，最后的结果自然是f[n]


递推公式
f[1]=0;
f[i]=(f[i-1]+m)%i;  (i>1)


有了这个公式，我们要做的就是从1-n顺序算出f[i]的数值，最后结果是f[n]。因为实际生活中编号总是从1开始，我们输出f[n]+1
由于是逐级递推，不需要保存每个f[i]，程序也是异常简单：

#include <cstdio>
int main() {
    int n , m;
    while(scanf("%d%d" , &n,&m) != EOF) {
        int a = 0;
        for(int i=2;i<=n;i++) {
            a = (a+m) % i;
        }
        printf("The winner is %d\n" , a);
    }
    return 0;
}