约瑟夫问题

一、约瑟夫问题来历

     据说著名犹太历史学家 Josephus有过以下的故事：在罗马人占领乔塔帕特后，39 个犹太人与Josephus及他的朋友躲到一个洞中，39个犹太人决定宁愿死也不要被敌人抓到，于是决定了一个自杀方式，41个人排成一个圆圈，由第1个人开始报数，每报数到第3的人，该人就必须自杀，然后再由下一个重新报数，直到所有人都自杀身亡为止。然而Josephus 和他的朋友并不想遵从，Josephus要他的朋友先假装遵从，他将朋友与自己安排在第16个与第31个位置，于是逃过了这场死亡游戏。

问题实际就是求解最后2个出列的序号。

      17世纪的法国数学家加斯帕在《数目的游戏问题》中讲了这样一个故事：15个教徒和15 个非教徒在深海上遇险，必须将一半的人投入海中，其余的人才能幸免于难，于是想了一个办法：30个人围成一圆圈，从第一个人开始依次报数1，每数到第九个人就将他扔入大海，如此循环进行直到仅余15个人为止。问怎样排法，才能使每次投入大海的都是非教徒。

问题实际求解最后15个出列的序号。

2个历史事件，实际都是求解出列的顺序，然后进行安排，即可达到要求

二、经典约瑟夫问题数学描述

      设有n个人站成一个圈，其编号从1~N。从编号为1的人开始按顺时针方向“1，2，3，4...”循环报数，数到m的人出列，然后从出列者的下一个人重新开始从1报数，数到m的人又出列，如此重复进行，直到n个人都出列为止。要求输出n个人的出列顺序。

     例如，有8个人的初始序列为

     1 2 3 4 5 6 7 8

当m=4时，出列顺序为

     4 8 5 2 1 3 7 6

三、约瑟夫问题解法及时间空间复杂度

/* 经典约瑟夫问题 模拟解法：数组或者循环单链表 */ #include "stdafx.h" #include <stdio.h> #include <malloc.h> #define MAXSIZE 256 typedef struct _List { int num; struct _List *next; }List,*PList; void Array_Josephus(int n, int m) //数组解法 时间复杂度o(n*n) { int p[MAXSIZE]; int i,j,t; //构建初始序列 printf("初始序列:"); for(i=0; i<n; i++) { p[i]=i+1; printf("%d ",p[i]); } printf("/n"); //出列 t = 0; //起始下标 printf("出列顺序:"); for(i=n; i>=1 ;i--) { t = (t+m-1)%i; //出列坐标 printf("%d ",p[t]); for(j=t;j<i;j++) { p[j]=p[j+1]; //移除坐标t } } printf("/n"); } void List_Josephus(int n, int m) //链表解法 时间复杂度o(n*m) { //尾插法建立初始序列 int i; int j=0; PList head = NULL; PList rear = NULL; PList t = NULL; PList q = NULL; printf("初始序列:"); for(i=0;i<n;i++) { PList p=(PList)malloc(sizeof(List)); p->num = i+1; printf("%d ",p->num); if(i==0) { head = p; rear = p; } else { //尾插法 rear->next=p; rear = p; } } rear->next = head; printf("/n"); //出列 printf("出列顺序:"); q=head; t=rear; while(t->next!=t) { j=j+1; if(j==m) { printf("%d ",q->num); t->next=q->next; free(q); q=t->next; j=0; } else { q=q->next; t=t->next; } } printf("%d/n",q->num); free(q); } /*数学解法 问题描述：n个人（编号0~(n-1))，从0开始报数，报到(m-1)的退出，剩下的人继续从0开始报数。求胜利者的编号。 　　 我们知道第一个人(编号一定是(m-1) mod n) 出列之后，剩下的n-1个人组成了一个新的约瑟夫环（以编号为k=m mod 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 很容易得出递推公式： f(1)=0; f(n)=(f(n-1) + m )%n */ void Math_Josephus(int n, int m) //数学解法 时间复杂度o(n) { int f=0; int i; for(i=2; i<=n; i++) { f=(f+m)%i; } printf("最后出列编号：%d/n",f+1); } int main(int argc, char* argv[]) { int n,m; do { printf("N M = "); scanf("%d %d",&n,&m); //数组模拟 printf("-----数组模拟-----/n"); Array_Josephus(n,m); printf("-----循环单模拟-----/n"); List_Josephus(n,m); printf("----数学求最后出列编号----/n"); Math_Josephus(n,m); } while (1); return 0; }

 四、约瑟夫问题求最后出列的人—数学妙解

/*数学解法 问题描述：n个人（编号0~(n-1))，从0开始报数，报到(m-1)的退出，剩下的人继续从0开始报数。求胜利者的编号。 　　 我们知道第一个人(编号一定是(m-1) mod n) 出列之后，剩下的n-1个人组成了一个新的约瑟夫环（以编号为k=m mod 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 很容易得出递推公式： f(1)=0; f(n)=(f(n-1) + m )%n */ void Math_Josephus(int n, int m) //数学解法 时间复杂度o(n) { int f=0; int i; for(i=2; i<=n; i++) { f=(f+m)%i; } printf("最后出列编号：%d/n",f+1); }

                                                                                                              2010-09-28  17:07:09