Josephus Problem的详细算法及其Python、Java实现
笔者昨天看电视,偶尔看到一集讲述古罗马人与犹太人的战争——马萨达战争,深为震撼,有兴趣的同学可以移步:http://finance.ifeng.com/a/20170627/15491157_0.shtml .
这不仅让笔者想起以前在学数据结构时碰到的Josephus问题:
据说著名犹太历史学家 Josephus有过以下的故事:在罗马人占领乔塔帕特后,39 个犹太人与Josephus及他的朋友躲到一个洞中,39个犹太人决定宁愿死也不要被敌人找到,于是决定了一个自杀方式,41个人排成一个圆圈,由第1个人开始报数,每报数到第3人该人就必须自杀,然后再由下一个重新报数,直到所有人都自杀身亡为止。然而Josephus 和他的朋友并不想遵从,Josephus要他的朋友先假装遵从,他将朋友与自己安排在第16个与第31个位置,于是逃过了这场死亡游戏。
以前我们都是用链表的方法编程来解决这个问题的,这次笔者将会讲述两个不同的方法,一个是笔者自己的朴素想法,一个是数学方法。
- 朴素方法
- 数学方法
首先我们先将Josephus问题描述出来,即: 共有N个人围成一圈,编号分别为1,2,...,N,从第一个人开始从1到m报数,报到m的人退出,如此循环下去,直至最后一个人。问最后一个人的最开始的编号是几?
先是笔者的朴素想法。
将N个人储存在列表(list)中,每次报到m的元素剔除,并记录下最后一个人报的数i,然后将缩短后的数组从i+1报数,如此循环下去,直至列表的长度为1,这样剩下来的元素就是我们要求的答案。
这种想法虽然素朴,比较容易实现,但是时间复杂度为O(Nm).
接着是数学方法。
假设一开始的Josephus环编号为0,1,2,...,N-1.我们知道第一个人(编号一定是m%N-1) 出列之后,剩下的N-1个人组成了一个新的Josephus环(以编号为k=m%n的人开始):
\]
并且从k开始报0.
现在我们把他们的编号做一下转换:
k+1 --> 1\\
k+2 --> 2\\
...\\
k-2 --> n-2\\
k-1 --> n-1\]
变换后就成为了(N-1)个人报数的子问题,这启示我们可以用归纳法来解决这个问题。假如我们知道这个子问题的解为\(x\),原来问题的答案为\(x^{'}\),则\(x^{'}=(x+k)\%n.\)因此,递推公式就有了!令\(f(i)\)表示\(i\)个人玩游戏报\(m\)退出最后胜利者的编号,我们要求的结果是\(f(N)\),其递推公式如下:
f(1)=(f(i-1)+m)\% i \qquad (i>1)\]
数学方法简单明了,计算效率高,但是推导比较困难。
最后,我们给出以下两种方法的Python代码和朴素方法的Java代码,希望能给大家一点思考。
完整的Python代码如下:
# -*- coding: utf-8 -*-
# This code is devoted to solve the Josephus Problem by Python.
# N: numper of people
# m: cycle number
def solve1(N, m):
a = list(range(1, N+1)) # sequence
end = 0 # the number of last man in sequence
while len(a) > 1:
b = []
for i in a:
if not (end+a.index(i)+1)%m:
b.append(i)
# print(i, end=' ') # print the order of removing
if a.index(i) == len(a)-1: # last man of sequence
end = (end+a.index(i)+1)%m
# remove elements in b from a
for i in b:
a.remove(i)
return a[0]
# solve the problem by math method
def solve2(N, m):
return 0 if N == 1 else (solve2(N-1, m)+m)%N
# main function for execution
def main():
N, m = 41, 3
left1 = solve1(N, m)
print('\nThe man left: %d' %left1)
left2 = solve2(N, m)+1
print('\nThe man left: %d' % left2)
main()
完整的Java代码如下:
import java.util.ArrayList;
public class Josephus {
public static void main(String[] args) {
int N = 41;
int m = 3;
int left = solve(N, m);
System.out.println("\nThe man left is "+left+".");
}
public static int solve(int N, int m) {
ArrayList<Integer> a = new ArrayList<Integer>();
int end = 0;
for(int i=0; i < N; i++)
a.add(i+1);
while(a.size() > 1) {
ArrayList<Integer> b = new ArrayList<Integer>();
for(int i: a) {
if ((end+a.indexOf(i)+1)%m == 0)
b.add(i);
// System.out.print(i+"-->");
if(a.indexOf(i) == a.size()-1)
end = (end+a.indexOf(i)+1)%m;
}
for(Object i: b) {
a.remove(i);
}
}
return a.get(0);
}
}
本次分享到此结束,欢迎大家交流~~
Josephus Problem的详细算法及其Python、Java实现的更多相关文章
- 算法Sedgewick第四版-第1章基础-017一约瑟夫问题(Josephus Problem)
/************************************************************************* * * Josephus problem * * ...
- 十大经典排序算法(Python,Java实现)
参照:https://www.cnblogs.com/wuxinyan/p/8615127.html https://www.cnblogs.com/onepixel/articles/7674659 ...
- 冒泡排序算法的C++,Java和Python实现和冒泡排序算法三种语言效率的比较
冒泡排序原理: 这一篇百度经验讲得很好,我不多说了 https://jingyan.baidu.com/article/6525d4b13f920bac7d2e9484.html 他讲的是C语言,没有 ...
- 八大排序算法总结与java实现(转)
八大排序算法总结与Java实现 原文链接: 八大排序算法总结与java实现 - iTimeTraveler 概述 直接插入排序 希尔排序 简单选择排序 堆排序 冒泡排序 快速排序 归并排序 基数排序 ...
- josephus Problem 中级(使用数组模拟链表,提升效率)
问题描写叙述: 在<josephus Problem 0基础(使用数组)>中.我们提出了一种最简单直接的解决方式. 可是,细致审视代码之后.发现此种方案的效率并不高,详细体如今.当有人出局 ...
- 机器学习算法与Python实践之(四)支持向量机(SVM)实现
机器学习算法与Python实践之(四)支持向量机(SVM)实现 机器学习算法与Python实践之(四)支持向量机(SVM)实现 zouxy09@qq.com http://blog.csdn.net/ ...
- 机器学习算法与Python实践之(二)支持向量机(SVM)初级
机器学习算法与Python实践之(二)支持向量机(SVM)初级 机器学习算法与Python实践之(二)支持向量机(SVM)初级 zouxy09@qq.com http://blog.csdn.net/ ...
- 约瑟夫问题(Josephus Problem)的两种快速递归算法
博文链接:http://haoyuanliu.github.io/2016/04/18/Josephus/ 对,我是来骗访问量的!O(∩_∩)O~~ 约瑟夫问题(Josephus Problem)也称 ...
- 机器学习算法与Python实践之(五)k均值聚类(k-means)
机器学习算法与Python实践这个系列主要是参考<机器学习实战>这本书.因为自己想学习Python,然后也想对一些机器学习算法加深下了解,所以就想通过Python来实现几个比较常用的机器学 ...
随机推荐
- Mac下Tomcat安装与Intellij IDEA配置Tomcat
Mac下Tomcat安装与Intellij IDEA配置Tomcat 一 安装 1 下载地址:https://tomcat.apache.org/download-90.cgi 2 将压缩包解压后移至 ...
- C#部分试题实例
1.在C#中,下列选项中自定义方法的语句错误的是().(选择一项) 正确答案:AD 解析:本题考查自定义方法的定义及调用.A项void是无返回值类型,D项定义方法的时候没有写返回值类型:故选AD. 2 ...
- Centos7通过yum安装最新MySQL
一:去官网查看最新安装包 https://dev.mysql.com/downloads/repo/yum/ 二:下载MySQL源安装包 wget http://dev.mysql.com/get/m ...
- VUE 动态给对象增加属性,并触发视图更新。
在开发过程中,我们时常会遇到这样一种情况:当vue的data里边声明或者已经赋值过的对象或者数组(数组里边的值是对象)时,向对象中添加新的属性,如果更新此属性的值,是不会更新视图的. 根据官方文档定义 ...
- 在Unity5.6.5f1中使用C#7语法
备忘,记忆力越来越差了,必需把这种琐碎的东西记下来,以防1年后想再用完全没头绪. 之前试过用C#6语法,但是怎么配置操作的完全没印象了. 首先去这下载扩展 https://bitbucket.org/ ...
- canvas 旋转
<!DOCTYPE html><html lang="en"><head> <meta charset="UTF-8" ...
- lgwr的两种模式(post/wait和polling)
11.2之前,oracle的lgwr写入模式为post/wait 11.2之后新增了polling模式,可以与post/wait模式自动切换 通过隐藏参数 _use_adaptive_log_file ...
- 背水一战 Windows 10 (77) - 控件(控件基类): ContentControl, UserControl, Page
[源码下载] 背水一战 Windows 10 (77) - 控件(控件基类): ContentControl, UserControl, Page 作者:webabcd 介绍背水一战 Windows ...
- OpenStack-Ocata版+CentOS7.6 云平台环境搭建 — 2.安装配置OpenStack基础服务
节点配置情况说明: 控制节点:controller: IP:192.168.164.128 hostname&hosts:likeadmin 计算加点:Nova: IP:192.168.164 ...
- zuul熔断代码
package com.sun.fallback; import java.io.ByteArrayInputStream; import java.io.IOException; import ja ...