约瑟夫问题小理解

问题由来

下面这是其中一个版本
传说有41个犹太人战败，他们不想被俘，因此围成一圈，从一个人开始由1 顺时针 依次标号。游戏中，从1号开始依次报号，每次报到2的人被上一个报到1的人杀死 ，直到只剩最后一人，他自杀 （犹太教义禁止自杀，他们也许想用这种方式在一定程度上规避自杀行为）。然而其中有1人不想死，于是他站到了19号位置逃避了死亡。

一般形式

有n个人围成一圈，从1号开始报号，报到m的人退出，直到只剩1人，问最后剩下的人是谁。

不同求法

m=2时的O(log n)算法

我从b站学的¹
手算结果如下（答案记作w(n)）：

n	w(n)
1	1
2	1
3	3
4	1
5	3
6	5
7	7
8	1
9	3
10	5
11	7
12	9
13	11

首先发现，w(n)全为奇数，这是因为所有偶数在第一圈筛除中就被淘汰。
然后发现，所有w(2^a)=1。证明：w(1)=1。若w(2^k)=1，当n=2^k+1时，第一圈淘汰后，留下2^k个数，由于在剩下2^k个人中第一个人会留到最后，因此w(2^k+1)=1。因此，w(2^a)=1。
这样，对任意一个正整数n，将n写成2^a+b的形式，那么我们在淘汰b个人后，接下来报1的那个人就是答案。因此，我们可以找到最大的a，答案就是(2*b+1)。例如如13=2³+5，w(13)=2*5+1=11。
找a时间复杂度O(log n)。

递归O(n)

如51nod1073 约瑟夫环。
我们将n个人，报m退出时最后剩下的那个人记作f(n,m)。
现在我们求f(n,m)，假设我们已经知道了f(n-1,m)，我们打算递归解决问题。
这n个人中，第一个退出的一定是m号。他退出后，剩下n-1个人，将从m+1号开始报数，最终留下的是从m+1号开始数的第f(n-1,m)个，即f(n,m)=m+f(n-1,m)。（略去了膜n的细节）

递归O(m+log n)

如51nod 1074 约瑟夫环 V2。
下面是我在51nod讨论区看到的一篇题解。²
沿用上面f(n,m)的定义，依然用递归，略去%n。
当n=1时，f(n,m)=1；
当n<=m时，f(n,m)=m+(n-1,m)；
当n>m时，令x=f(n-n/m,m)。第一圈筛除了n/m个，当前位置是cur=n/m*m。若x<=n%m，f(n,m)=cur+x。若x>n%m，则答案在下一圈中，令x-=n%m，在经历这x个人中，会遇到(x-1)/(m-1)个已退出的人，因此f(n,m)=x+(x-1)/(m-1)。时间复杂度为O(log n)级别。

下面是代码。不知道怎么mod n比较好，写得比较丑。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;

LL n, k;

LL f(LL n, LL k)
{
    if(n==1) return 1;
    if(n<=k) return (k+f(n-1,k)-1)%n+1;
    LL x=f(n-n/k,k);
    if(x<=n%k) return (n/k*k+x-1)%n+1;
    else
    {
        x-=n%k;
        return (x-1)/(k-1)+x;
    }
}

int main()
{
    cin>>n>>k;
    cout<<f(n,k)<<endl;
    return 0;
}

有好多人用循环写的，我看不懂，有哪位大神能教教我吗qwq

---

1. 互砍惨案！约瑟夫斯问题 飞刀好逗 ↩︎

2. 51nod1074讨论区 ↩︎