经典问题之约瑟夫问题的快速解决

此问题非常经典，在网上即可找到原题，在此不做描述。

对于原问题模型，一有链表法解决问题，效率极低，在此描述一种用树状数组完成问题的超快速做法。

首先，我们可以有这样递推的思路：不断加k模n，并减去其数字前走了的人即为当前人的真实编号（即是这一轮应踢走的人的编号），如何快速维护每个人其前走了的人的和，答案为树状数组。

现在模拟一下过程，假设有6个人，k=3（每报3个，走一个人）。

初始状态：1 2 3 4 5 6

用树状数组在每个人的位置加一，可得前缀和：1 2 3 4 5 6

现在1+2(其实是k-1)=3走了:1 2 4 5 6

用树状数组在3的位置减1，可得前缀和：1 2 2 3 4 5

再走了3+2=5，5%（6-1）=5（走了一个人，故6需减1）等等，此时并不是走了5，而是在树状数组中前缀和为5的数字，由上可知走了6：1 2 4 5

用树状数组在6的位置减1，可得前缀和：1 2 2 3 4 4

又按5+2=7，7%（6-2）=3，但这时在树状数组的前缀和中查找3，可见是第四个人的状态为3，故此回合走了4：1 2 5

用树状数组在4的位置减1，可得前缀和：1 2 2 2 3 3

又按3+2=5，5%（6-3）=2，在树状数组中查找二，可见即为2，

故此回合走了2：1 5，

前缀和改为：1 1 1 1 2 2

最后2+2=4，4%（6-4）=2，在树状数组中查找2，可见为5，

故最后只剩1了。

因为前缀和是单调的，所以查找可以用二分。

可得序列为：3 6 4 2 5 1

//不信的话可以和链表验证一下。

参考程序：

#include<cstdio>#include<algorithm>#include<cstring>#include<ctime>#define maxn 1000000using namespace std;int bit[maxn];int n,k;int sum(int i){int s=0;while (i>0){s+=bit[i];i-=(i&-i);}return s;}void add(int i,int x){while (i<=n){bit[i]+=x;i+=(i&-i);}}int binary_search(int id){int l=0,r=n;while (l<r){int mid=(l+r)>>1;if (sum(mid)<id)l=mid+1;else r=mid;}return l;}int main(){freopen("input.in","r",stdin);freopen("output.out","w",stdout);while (scanf("%d%d",&n,&k)==2){int id=1;memset(bit,0,sizeof(bit));for (int i=1;i<=n;i++)add(i,1);for (int i=1;i<=n;i++){id=(id+k-2)%(n-i+1)+1;int newid=binary_search(id);printf("%d ",id);add(newid,-1);}}printf("\nUsed Time=%.2f",(double)clock()/CLOCKS_PER_SEC);return 0;}

为什么说这种方法快呢？

1.检验：在n=100000的情况下，仅费时0.24秒。

2.可以明显得O(n*(logn)^2)，应该算是很快了。

对于只需要最后幸存者的问题，可参见递推解决