约瑟夫问题



约瑟夫 问题（有时也称为约瑟夫斯置换，是一个出现在计算机科学和数学中的问题。在计算机编程的算法中，类似问题又称为约瑟夫环。又称“丢手绢问题”.）

一般形式

编辑

约瑟夫问题是个有名的问题：N个人围成一圈，从第一个开始报数，第M个将被杀掉，最后剩下一个，其余人都将被杀掉。例如N=6，M=5，被杀掉的顺序是：5，4，6，2，3，1。

分析：

（1）由于对于每个人只有死和活两种状态，因此可以用布朗型数组标记每个人的状态，可用true表示死，false表示活。

（2）开始时每个人都是活的，所以数组初值全部赋为false。

（3）模拟杀人过程，直到所有人都被杀死为止。

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

为了讨论方便，先把问题稍微改变一下，并不影响原意：

问题描述：n个人（编号0~(n-1）），从0开始报数，报到（m-1）的退出，剩下的人继续从0开始报数。求胜利者的编号。

我们知道第一个人（编号一定是（m-1）) 出列之后，剩下的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

变换后就完完全全成为了（n-1）个人报数的子问题，假如我们知道这个子问题的解：例如x是最终的胜利者，那么根据上面这个表把这个x变回去不刚好就是n个人情况的解吗？！！变回去的公式很简单，相信大家都可以推出来：x'=(x+k) mod n

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

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

递推公式

f[1]=0;

f=(f+m) mod i; (i>1）

有了这个公式，我们要做的就是从1-n顺序算出f的数值，最后结果是f[n]。因为实际生活中编号总是从1开始，我们输出f[n]+1

由于是逐级递推，不需要保存每个f，程序也是异常简单：

#include <iostream>
using namespace std;
const int m = 3;
int main()
{
    int n, f = 0;
    cin >> n;
    for (int i = 1; i <= n; i++) f = (f + m) % i;
    cout << f + 1 << endl;
}

拓展：HDU 5860

给你n个人排成一列编号，每次杀第一个人第i×k+1个人一直杀到没的杀。然后剩下的人重新编号从1～剩余的人数。按照上面的方式杀。问第几次杀的是谁。

解析：这道题首先我们先将编号改成从0开始。这样如果i%k==0那么我们可以知道第i个数字第一轮就被杀死。如果i%k!=0，那么我们知道在这之前一轮杀死了i/k+1个人。这样我们就得到递推式。

定义dp[i]代表第i轮杀死多少人。num[i]表示i在当前轮第几个被杀死。

对于i%k==0我们可以知道第一轮被杀死，所以dp[i]=0,num[i]=i/k+1。

对于i%k!=0通过前面我们知道他们的状态和i-i/k-1有关。故dp[i]=dp[i-i/k-1]+1,num[i]=num[i-i/k-1]。

代码如下：

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<limits.h>
#include<algorithm>
#include<queue>
#include<math.h>
#include<stack>
#include<vector>
using namespace std;
#define maxn 100005
int dp[maxn],ans[maxn];
int s[maxn];
int n,m,k;
void work()
{
	int i,j;
	int cnt=0,t=n;
	int kill[maxn]={0};//kill[i]表示前i轮杀了多少人 
	while(t)
	{
		cnt++;
		kill[cnt]=kill[cnt-1]+(t-1)/k+1;
		t-=(t-1)/k+1;
	}
	for(i=0;i<n;i++)
	{
		if(i%k==0)
		{
		    dp[i]=0;
		    s[i]=i/k+1;
		}
		else
		{
			dp[i]=dp[i-i/k-1]+1;
			s[i]=s[i-i/k-1];
		}
	}
	for(i=0;i<n;i++)
	 ans[kill[dp[i]]+s[i]]=i;
}
int  main()
{
	int T,i,j,x;
	queue<int>q;
	scanf("%d",&T);
	while(T--)
	{
		memset(dp,0,sizeof(dp));
		memset(s,0,sizeof(s));
		scanf("%d%d%d",&n,&k,&m);
		work();
		for(i=1;i<=m;i++)
		{
			scanf("%d",&x);
			printf("%d\n",ans[x]+1);
		}
	}
}