约瑟夫问题

今天闲着无聊整理下《具体数学》中关于约瑟夫问题的讨论。

$n$个人围成一圈，从$1$到$n$标号，每隔一个人处死一个人，求最后幸存下来的人的编号，例如$n=10$，则被杀的人按顺序依次是：

$$2,4,6,8,10,3,7,1,9$$

当然可以直接模拟杀人的过程直到剩下最后一个，但这种方法不够快。

可以考虑递归做法，如果杀了一圈，那么这个时候我们可以把剩下的人从$1$开始重新排列，这样如果能够找到两次排列的编号关系，就可以用递归来解决了。

假设$n$为偶数，那么因为隔一个杀一个，因此，所有标号为偶数的人都会被杀掉，因此，我们如果将剩下幸存的人重新标号可得：

以$n=10$为例，

杀之前： $[1,2,3,4,5,6,7,8,9,10]$

杀之后： $[1,~~,2,~~,3,~~,4,~~,5,~~~~]$

其中空出来的表示在这轮中阵亡的。

可以很轻松的发现这两者之间的关系（$*2-1$），因此我们可得：

假设$J(n)$为$n$个人的约瑟夫环最后剩下的幸存者的编号，因此有

$$J(2n)=2J(n)-1$$

同理如果是奇数呢

以$n=9$为例，

杀之前： $[1,2,3,4,5,6,7,8,9]$

杀之后： $[~~,~~,1,~~,2,~~,3,~~,4]$

其中空出来的表示在这轮中阵亡的。

因此，我们也可以得到他们之间的关系：($*2+1$)，所以有：

$$J(2n+1)=2J(n)+1$$

因此合起来我们可以得到递推式：

$J(2n)=2J(n)-1$

$J(2n+1)=2J(n)+1$

当$n=1$时，则有$J(n)=1$

接下来求$J(n)$的闭形式：

根据$J(2n)=2J(n)-1$，且$J(1)=1$可得，$J(2^m)=1$恒成立，因此假设对于所有的数都能表示为$2^m+l$，$0 \leq l<2^m$，则有

$$J(2^m+l)=2l+1$$
这是因为，如果杀了$l$个人，此时，剩下的人数则为$2^m$个，因此，此时幸存的人应该是第$l$个被杀的人的后一个，因此也就是$2l+1$。

在书中，作者是用了猜想+数学归纳法的方法来证明这个式子的，之后他发现了一个有趣的模式：

如果用二进制表示递推式，可得：

$n=(1b_{m-1}b_{m-2}...b_1b_0)_2$

$l=(0b_{m-1}b_{m-2}...b_1b_0)_2$

$2l+1=(b_{m-1}b_{m-2}...b_1b_01)_2$

因此最后可得：

$J((b_mb_{m-1}...b_1b_0)_2)=(b_{m-1}...b_1b_0b_m)_2$

也就是将人数$n$的二进制形式向左循环移动一位，就可以得到幸存的人的编号，通过这也就可以得出如果当人数为$2^m-1$时，最后一个人将幸存。

接下来拓展一下，可以解决形如：

$f(1)=\alpha$

$f(2n+j)=2f(n)+\beta_j$

其中$j=0,1$，$n\geq1$

的问题：

假设$n=(b_mb_{m-1}...b_1b_0)_2$

则可以不断把最后一位提出来：

$f((b_mb_{m-1}...b_1b_0)_2)=2f((b_mb_{m-1}...b_1)_2)+\beta_{b_0}$
…
…
$=2^m\alpha+2^{m-1}\beta_{b_{m-1}}+...+2\beta_{b_1}+\beta_{b_0}$

第一位是$\alpha$的原因在于$b_m=1$。

于是就有了如下的结论：

$f(b_mb_{m-1}...b_1b_0)_2=(\alpha\beta_{b_{m-1}}\beta_{b_{m-2}}...\beta_{b_1}\beta_{b_0})_2$

回到约瑟夫问题，这是当$\alpha=1, \beta_0=-1, \beta_1=1$时的特例，可以发现，当$n$写成$(\beta_m\beta_{m-1}...\beta_1\beta_0)_2$时，我们可以把它分成若干段$(10...00)_2$，而每段都会变成$(1-1...-1-1)_2=(00...01)_2$，这也就说明了向左循环移位的事实。

例如，当$n=100$时，有$n=(1100100)_2$，因此可以分解成三段：$(1)_2, (100)_2, (100)_2$，而每段变换后分别得到$(1)_2,(001)_2,(001)_2$，合起来后变成$(1001001)_2$，也就相当于进行了向左循环移位。

当然我们可以再弱化一下递归式，变成：

$f(j)=\alpha_j$, $~~~~~~~~~~~~~~~~~~~~~~~1\leq j < d$

$f(dn+j)=cf(n)+\beta_j$, $0 \leq j <d, n \geq 1$

我们同样可以假设此时的$n$可以表示为$n=(b_mb_{m-1}...b_1b_0)_d$，然后用同样方法展开：

$~~~~f((b_mb_{m-1}...b_1b_0)_d)$

$=cf((b_mb_{m-1}...b_1)_d)+\beta_{b_0}$

$=c^2f((b_mb_{m-1}...b_2)_d)+c\beta_{b_1}+\beta_{b_0}$

$...$

$=c^m\alpha_{b_m}+c^{m-1}\beta_{b_{m-1}}+...+c\beta_{b_1}+\beta_{b_0}$

$=(\alpha_{b_m}\beta_{b_{m-1}}\beta_{b_{m-2}}...\beta_{b_1}\beta_{b_0})_c$

即$f((b_mb_{m-1}...b_1b_0)_d) = (\alpha_{b_m}\beta_{b_{m-1}}\beta_{b_{m-2}}...\beta_{b_1}\beta_{b_0})_c$

所以遇到这种类型的递归式，只需要做进制数的改变就行了。
例如，有递归式

$f(1)=34$

$f(2)=5$

$f(3n)=10f(n)+76, (n\geq 1)$

$f(3n+1)=10f(n)-2, (n\geq 1)$

$f(3n+2)=10f(n)+8, (n\geq 1)$

且我们要求解$f(19)$，因为$d=3, c=10$，因此，