《具体数学》1.3递归(三)

文章开头:先吐槽一下长春的天气,冬天是一星期有6天在下雪,现在算是夏天了吧,动不动下雨下冰雹,天气变化的有点快,还有就是学校的妹子,分为女神和女屌丝,女神一般在长春流行穿风衣,女屌丝一般流行穿牛仔系列,最后吐槽一下,我曾经的喜欢过的人变了,唏嘘不已,好了吐槽完毕,这一篇介绍一下递归的最后一个故事Josephus Circle(约瑟夫环问题),这个部分有点多估计得分两次了。。。。

1.3Josephus Circle(约瑟夫环问题)

这本书还是老规矩,先以故事引出问题:这最后的例子是一个古老问题的变型,它是以公元一世纪著名历史学家Flavius Josephs命名的。传说是Josephus的数学才干让他逃过一劫。
在犹太人和罗马人战争期间,Josephs等41个犹太反抗者陷入包围。他们宁死不做俘虏,决定围成一个圆圈,每次由剩下的人中的第三个自杀,直到所有人都死去为止。
Josephus和他的1个朋友认为,自杀是愚蠢的。因此他很快地算出死亡之圈中他俩应该站的位置(当然,是两个应该最后自杀的位置,没有人监督他们是否自杀,所以你懂得)

我们的问题是,从编号为1到n的围成一圈的人开始,依次排除剩下的人中的第2个(注意:是从当前“指针”开始),直到仅剩1个幸存者。下面是n = 10的情况:
145332_KRba_1589875.jpg

排除次序是2,4,6,8,10,3,7,1,9。因此5幸存。
问题目标是:如何确定幸存者的号码J (n)?

我们刚才看到J(10)=5.我们可能猜测当n是偶数时J(n)=n/2;而且情形n=2支持了猜测:J(2)=1

但对于n=4,n=6猜测失败,从上面的表看来,J(n)都是奇数,因此将其作为新的猜想:事实上第1轮筛选就排除了所有的偶数,因此猜想成立。
另外还观察到,若n是偶数,那么在做完第1轮筛选之后,剩下的情形和最初的开始情况是非常类似的,但是人数要少一半,而且编号也有所不同。
假设原先有2n个人,在第1轮筛选之后,剩下的情况是
145814_48dJ_1589875.jpg

此时圆圈剩下了n个人:1、3、5、7、……、2n – 3、2n – 1。而且“指针”又回到了1的位置。这与最初有n个人的情形是相同的,区别仅在于编号有变化。

编号的变化会带来什么影响?
从编号1开始,下1个被排除的人是3,其后的筛选过程与n个人的情形一致:
当前的情形:2、4、6、8、……(差为2的等差数列)
 n个人的情形: 3、7、11、15(差为4的等差数列)
相同点:人数相同
不同点:编号变为2倍减1

注意到,幸存者J(2n)必然从当前剩下的人中出现,因此:
J(2n) = 2J(n) – 1

其实,我们对于人数为奇数页可以得到J(2n+1)=2J(n)+1

到目前为止,我们得到了两个“递归”色彩的式子:
J(2n) = 2J(n) - 1
J(2n + 1) = 2J(n) + 1
上面涵盖了所有“偶数”和“奇数” ,即所有正整数,但还缺少boundary value,作为“递归系统”的第一推动力。
很简单,J(1) = 1(最后剩下1个人,不必再筛选)。这样就得到了完整的J的递归方程。
这个递归方程的效率非常高,因为每次将n缩小至少1/2。例如,我们用19次就能计算出J(1000000):

但如果像上面那样就完 了,这本书真的很厉害,有好多新的思想和套路

但是无论从数学的完美意义,还是实际的计算效率上来考虑,最好找到闭形式解。我们仍然使用“猜想证明”思路来求解封闭形式解。首先根据递归方程,对较小的n建立J(n)表:

150354_sAnv_1589875.jpg

可以发现什么规律?(若把n记为形式n=2^m+l,其中2^m是不超过n的2的最大幂,而l则是留下的量)
1、每隔若干个位置就从1开始
2、每次从1开始的子序列都是级差为2的等差数列
3、每个子序列的长度是递增的
继续观察,可以看到:
1、第1子序列长度为1,此后长度为2、4、8、……
(谨慎怀疑子序列长度构成比值为2的等比数列!!)
2、如果上述怀疑成立,那么我们只要知道当前n处于第几个子序列的第几个位置,就能够从1逐渐累加得到J(n)了!!!
继续找n和J(n)的数量关系:
假设n之前已经有m个子序列,而n是第(m+1)个子序列的第l个位置,那么显然有
J(n) = 1 + 2(l - 1)

n与m、l之间有什么关系呢?
n之前的整数数量为(1 + 2 + 4 + …… + 2^m-1) + (l - 1),因此有
n = (2^m  – 1) + l = 2^m + (l – 1)

最后我们可以猜测得到封闭形式为:J(2^m+l) = 2l + 1(m>=0,0<=l<2^m)

(注意,若2^m<=n<2^m+1,则剩余部分l=n-2^m满足0<=l<2^m+1-2^m=2^m)

要证明:J(2^m+l) = 2l + 1(m>=0,0<=l<2^m)这肯定要用到数学归纳法了,但是这里是对m归纳。当m=0时,我们一定有l=0;

因此递归式基础化为J(1)=1,归纳步还分为两部分,依赖于l是偶数还是奇数,若m>0和2^m+l=2n,则 l是偶数并由式J(2n) = 2J(n) - 1和归纳假设可得

J(2^m+l)=2J(2^m-1+l/2)-1=2(2l/2+1)-1=2l+1,这恰好是我们要的,同理可得,奇数情形,完成归纳,这次我们计算J(100),此时我们得到100=2^6+36,所以J(100)=2*36+1=73(其实,在这里我已经觉得很厉害了,下面的更厉害)

大师发现解中2的幂起了重要作用,所以自然要看二进制表示n和J(n),假设n的二进制展开为

n=(bmbm-1...b1b0)2;

也就是说,n=bm2^m+bm-12^m-1+...+b12+b0,其中每个bi是0或1,且其中的第一位bm是1.想到n=2^m+1,我们相继得到

n=(1bm-1bm-2...b1b0)2,l=(0bm-1bm-2...b1b0)2,2l=(bm-1bm-2...b1b01)2,2l+1=(bm-1bm-2...b1b01)2,J(n)=(bm-1bm-2...b1b0bm)2

最后一步是由J(n)=2l+1和bm=1而得出的所以我们证明了J((bmbm-1...b1b0)2)=(bm-1...b1b0bm)2;

也就是说,以计算机程序设计的用语,循环左移一位我们从n取得J(n)不可思议,例如若n=100=(1100100)2,则J(n)=J((1100100)2)=(1001001)2,它是64+8+1=73。

好了先就写到这了,感觉这本书在数学上真的是下了很多功夫,一步一步的引导我们去阅读与思考,有点上瘾了,我怎么了,等过几天,再把另一部分补上,因为有点多,讲约瑟夫这块。。。





转载于:https://my.oschina.net/u/1589875/blog/265631

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值