程序员十大常用算法
文章目录
二分查找算法(非递归)
二分查找算法(非递归)介绍
- 前面我们讲过了二分查找算法，是使用递归的方式，下面我们讲解二分查找算法的非递归方式
- 二分查找法只适用于从有序的数列中进行查找(比如数字和字母等)，将数列排序后再进行查找
- 二分查找法的运行时间为对数时间O(㏒₂n)，
即查找到需要的目标位置最多只需要㏒₂n步，
假设从[0,99]的队列(100个数，即n=100)中寻到目标数30，
则需要查找步数为㏒₂100, 即最多需要查找7次( 2^6 < 100 < 2^7)
数组 {1,3, 8, 10, 11, 67, 100}, 编程实现二分查找.
要求使用非递归的方式完成.
代码实现：
public class BinarySearchNoRecur {
public static void main(String[] args) {
//测试
int[] arr = {1,3,8,10,11,67,100};
int index = binarySearch(arr, -8);
System.out.println("index=" + index);//1
}
//二分查找的非递归实现
/**
*
* @param arr 待查找的数组,arr时升序排序
* @param target 需要查找的数
* @return 返回对应下标,-1表示没有找到
*/
public static int binarySearch(int arr[], int target) {
int left = 0;
int right = arr.length - 1;
while (left <= right){//说明可以继续查找
int mid = (left + right) / 2;
if (arr[mid] == target){
return mid;
}else if (arr[mid] > target){
right = mid - 1;//需要向左边查找
}else {
left = mid + 1;//需要向右边查找
}
}
return -1;
}
}
分治算法
1.分治算法介绍
分治法是一种很重要的算法。字面上的解释是“分而治之”，
就是把一个复杂的问题分成两个或更多的相同或相似的子问题，
再把子问题分成更小的子问题……直到最后子问题可以简单的直接求解，
原问题的解即子问题的解的合并。这个技巧是很多高效算法的基础，
如排序算法(快速排序，归并排序)，傅立叶变换(快速傅立叶变换)……
2.分治算法可以求解的一些经典问题:
二分搜索
大整数乘法
棋盘覆盖
合并排序
快速排序
线性时间选择
最接近点对问题
循环赛日程表
汉诺塔
分治算法的基本步骤
分治法在每一层递归上都有三个步骤：
- 1.分解：将原问题分解为若干个规模较小，相互独立，与原问题形式相同的子问题
- 2.解决：若子问题规模较小而容易被解决则直接解，否则递归地解各个子问题
- 3.合并：将各个子问题的解合并为原问题的解。
分治(Divide-and-Conquer§)算法设计模式如下：
其中|P|表示问题P的规模；n0为一阈值，表示当问题P的规模不超过n0时，问题已容易直接解出，不必再继续分解。
ADHOC§是该分治法中的基本子算法，用于直接解小规模的问题P。
因此，当P的规模不超过n0时直接用算法ADHOC§求解。
算法MERGE(y1,y2,…,yk)是该分治法中的合并子算法，用于将P的子问题P1 ,P2 ,…,Pk的相应的解y1,y2,…,yk合并为P的解。
分治算法最佳实践——汉诺塔
汉诺塔的传说
汉诺塔：汉诺塔（又称河内塔）问题是源于印度一个古老传说的益智玩具。
大梵天创造世界的时候做了三根金刚石柱子,
在一根柱子上从下往上按照大小顺序摞着64片黄金圆盘。
大梵天命令婆罗门把圆盘从下面开始按大小顺序重新摆放在另一根柱子上。
并规定，在小圆盘上不能放大圆盘，在三根柱子之间一次只能移动一个圆盘。
假如每秒钟一次，共需多长时间呢？移完这些金片需要5845.54亿年以上，
太阳系的预期寿命据说也就是数百亿年。真的过了5845.54亿年，
地球上的一切生命，连同梵塔、庙宇等，都早已经灰飞烟灭。
汉诺塔游戏的演示和思路分析:
1)如果是有一个盘，A->C
2)如果我们有 n >= 2 情况，我们总是可以看做是两个盘:
1.最下边的盘
2.上面的盘
操作步骤：
1)先把最上面的盘 A->B
2)把最下边的盘 A->C
3)把B塔的所有盘 从 B->C
代码实现:
public class Hanoitower {
public static void main(String[] args) {
hanoiTower(3, 'A', 'B', 'C');
}
//汉诺塔移动的方法
//使用分治算法
public static void hanoiTower(int num, char a, char b, char c) {
//如果只有一个盘
if (num == 1) {
System.out.println("第1个盘从 " + a + "->" + c);
}else {
//如果我们有 n>=2 情况,我们总是可以看作是两个盘 1.最下边的一个盘 2.上面的所有盘
//1.先把最上面的盘A->B,移动过程会使用到c
hanoiTower(num-1, a, c, b);
//2.把最下边的盘 A->C
System.out.println("第" + num + "个盘从 " + a + "->" + c);
//3.把B塔的所有盘 从B->C,移动过程使用到a塔
hanoiTower(num-1, b, a, c);
}
}
}
动态规划算法
背包问题：有一个背包，容量为4磅 ， 现有如下物品
| 物品 | 重量 | 价格 |
|---|---|---|
| 吉他(G) | 1 | 1500 |
| 音响(S) | 4 | 3000 |
| 电脑(L) | 3 | 2000 |
- 要求达到的目标为装入的背包的总价值最大，并且重量不超出
- 要求装入的物品不能重复
动态规划算法介绍
1.动态规划(Dynamic Programming)算法的核心思想是： 将大问题划分为小问题进行解决，从而一步步获取最优解的处理算法
2.动态规划算法与分治算法类似，其基本思想也是 将待求解问题分解成若干个子问题，先求解子问题，然后从这些子问题的解得到原问题的解
3.与分治法不同的是，适合于用动态规划求解的问题，经分解得到子问题往往不是互相独立的。 ( 即下一个子阶段的求解是建立在上一个子阶段的解的基础上，进行进一步的求解 )
4.动态规划可以通过填表的方式来逐步推进，得到最优解.
思路分析和图解:
背包问题主要是指一个给定容量的背包、若干具有一定价值和重量的物品，
如何选择物品放入背包使物品的价值最大。
其中又分01背包和完全背包(完全背包指的是：每种物品都有无限件可用)
这里的问题属于01背包，即每个物品最多放一个。
而无限背包可以转化为01背包。
算法的主要思想，利用动态规划来解决。
每次遍历到的第i个物品，根据w[i]和v[i]来确定是否需要将该物品放入背包中。
即对于给定的n个物品，设v[i]、w[i]分别为第i个物品的价值和重量，C为背包的容量。
再令v[i][j]表示在前i个物品中能够装入容量为j的背包中的最大价值。则我们有下面的结果：
(1) v[i][0]=v[0][j]=0;
表示 填入表 第一行和第一列是0。
(2) 当w[i]> j 时：v[i][j]=v[i-1][j]
当准备加入新增的商品的容量 > 当前背包的容量时，
就直接使用上一个单元格的装入策略。
(3) 当j>=w[i]时： v[i][j]=max{v[i-1][j], v[i]+v[i-1][j-w[i]]}
当 准备加入的新增的商品的容量 <= 当前背包的容量 装入的方式: v[i-1][j]： 就是上一个单元格的装入的最大值 v[i] : 表示当前商品的价值 v[i-1][j-w[i]] ： 装入i-1商品，到剩余空间j-w[i]的最大值 当j>=w[i]时: v[i][j]=max{v[i-1][j],v[i]+v[i-1][j-w[i]]} :
对上述公式进行验证：
验证1：v[1][1] = 1500
- i = 1, j = 1
- w[i] = w[1] = 1 w [1] = 1 j = 1 v[i][j]=max{v[i-1][j], v[i]+v[i-1][j-w[i]]} :
- v[1][1] = max {v[0][1], v[1] + v[0][1-1]} = max{0, 1500 + 0} = 1500
验证2：v[3][4] = 2000+1500
- i = 3；j = 4
- w[i] = w[3] =3, j = 4 j = 4 >= w[i] = 3 => 4 >= 3
- v[3][4] = max {v[2][4], v[3] + v[2][1]} = max{3000, 2000+1500}
= 2000+1500
代码实现:
public class KnapsackProblem {
public static void main(String[] args) {
int[] w = {1, 4, 3, 5, 2, 6, 7, 9, 8};//物品的重量
int[] val = {1500, 3000, 2000, 1000, 850, 2700, 1650, 850, 1000};//物品的价值 这是的val[i] 就是前面讲的v[i]
int m = 20;//背包的容量
int n = val.length;
//创建二维数组
//v[i][j] 表示在 前i个物品中 能够装入容量为j的 背包中的最大值。
int[][] v = new int[n + 1][m + 1];
//为了记录放入商品的情况,我们定义一个二维数组
int[][] path = new int[n + 1][m + 1];
//初始化第一行和第一类,这里在本程序中,可以不去处理,因为默认就是0
for (int i = 0; i < v.length; i++) {
v[i][0] = 0;//将第一列设置为0
}
for (int i = 0; i < v[0].length; i++) {
v[0][i] = 0;//将第一行设置为0
}
//根据前面得到的公式来动态规划处理
for (int i = 1; i < v.length; i++) {//不处理第一行,i从1开始
for (int j = 1; j < v[0].length; j++) {//不处理第一列,j从1开始
//公式
if (w[i - 1] > j) {//因为程序i是从1开始的,因此公式中的w[i]->w[i-1]
v[i][j] = v[i - 1][j];
} else {
/*说明:
因为我们的i 从1开始的,因此公式需要调整
v[i][j] = Math.max(v[i - 1][j], val[i] + v[i - 1][j - w[i]]);
v[i][j] = Math.max(v[i - 1][j], val[i - 1] + v[i - 1][j - w[i - 1]]);
为了记录商品存放到背包的情况,我们不能直接的使用上面公式
需要用if-else来体现
*/
if (v[i - 1][j] < val[i - 1] + v[i - 1][j - w[i - 1]]) {
v[i][j] = val[i - 1] + v[i - 1][j - w[i - 1]];
//把当前的情况记录到path
//因为只有这个条件是最优的,我们只要把这个条件写下来就ok了.
path[i][j] = 1;
} else {
v[i][j] = v[i - 1][j];
}
}
}
}
//输出一下v 看看目前的情况
for (int i = 0; i < v.length; i++) {
for (int j = 0; j < v[i].length; j++) {
System.out.print(v[i][j] + " ");
}
System.out.println();
}
System.out.println("===========");
//输出最后我们是放入的哪些商品
//遍历path,这样输出会把所有的放入情况都得到,其实我们只需要最后的放入情况
/*for (int i = 0; i < path.length; i++) {
for (int j = 0; j < path[i].length; j++) {
if (path[i][j] == 1){
System.out.print("第" + i + "个商品放入到背包\n");
}
}
}*/
//
int i = path.length - 1;//行的最大下标
int j = path[0].length - 1;//列的最大下标
while (i > 0 && j > 0) {//逆向遍历,从path最后开始找
if (path[i][j] == 1) {
System.out.print("第" + i + "个商品放入到背包\n");
j -= w[i - 1];
}
i--;
}
}
}
KMP算法
应用场景——字符串匹配问题
- 有一个字符串 str1= “数据 数据结构 数据结构与算法 Java数据结构”，和一个子串 str2=“数据结构Java”
- 现在要判断 str1 是否含有 str2, 如果存在，就返回第一次出现的位置, 如果没有，则返回-1
暴力匹配算法
如果用暴力匹配的思路，并假设现在str1匹配到 i 位置，子串str2匹配到 j 位置，则有:
1.如果当前字符匹配成功（即str1[i] ==str2[j]），则i++，j++，继续匹配下一个字符
2.如果失配（即str1[i]! = str2[j]），令i = i - (j -1)，j = 0。 相当于每次匹配失败时，i 回溯，j 被置为0。
3.用暴力方法解决的话就会有大量的回溯，每次只移动一位，若是不匹配，移动到下一位接着判断，浪费了大量的时间。(不可行!)
暴力匹配算法实现:
public class ViolenceMatch {
public static void main(String[] args) {
//测试暴力匹配法
String str1 = "尚硅谷谷 尚硅谷你尚硅 尚硅谷你尚硅谷你尚硅你好";
String str2 = "尚硅谷你尚硅你";
int index = violenceMatch(str1, str2);
System.out.println("index=" + index);
}
//暴力匹配算法实现
public static int violenceMatch(String str1, String str2) {
char[] s1 = str1.toCharArray();
char[] s2 = str2.toCharArray();
int s1Len = s1.length;
int s2Len = s2.length;
int i = 0;//i索引指向s1
int j = 0;//j索引指向s2
while (i < s1Len && j < s2Len) {//保证匹配时,不越界
if (s1[i] == s2[j]) {//匹配成功
i++;
j++;
} else {//没有匹配成功
//如果失败(即str1[i]!=str2[j]),令i = i - (j - 1),j = 0;
i = i - (j - 1);
j = 0;
}
}
//判断是否匹配成功
if (j == s2Len) {
return i - j;
} else {
return -1;
}
}
}
KMP算法简介
KMP是一个解决模式串在文本串是否出现过，如果出现过，
得到最早出现的位置的经典算法
Knuth-Morris-Pratt 字符串查找算法，简称为 “KMP算法”，
常用于在一个文本串S内查找一个模式串P 的出现位置，
这个算法由Donald Knuth、Vaughan Pratt、James H. Morris
三人于1977年联合发表，故取这3人的姓氏命名此算法.
KMP方法算法就利用之前判断过信息，通过一个next数组，
保存模式串中前后最长公共子序列的长度，每次回溯时，
通过next数组找到，前面匹配过的位置，省去了大量的计算时间
参考资料：KMP算法参考资料[转载]
字符串匹配问题：
- 有一个字符串 str1= “BBC ABCDAB ABCDABCDABDE”，和一个子串 str2=“ABCDABD”
- 现在要判断 str1 是否含有 str2, 如果存在，就返回第一次出现的位置, 如果没有，则返回-1
- 要求：使用KMP算法完成判断，不能使用简单的暴力匹配算法.
思路分析
- 首先，用Str1的第一个字符和Str2的第一个字符去比较，不符合，关键词向后移动一位
- 重复第一步，还是不符合，再后移
-
一直重复，直到Str1有一个字符与Str2的第一个字符符合为止
-
接着比较字符串和搜索词的下一个字符，还是符合。
-
遇到Str1有一个字符与Str2对应的字符不符合。
6. 这时候，想到的是继续遍历Str1的下一个字符，重复第1步。
(其实是很不明智的，因为此时BCD已经比较过了，没有必要再做重复的工作，一个基本事实是，当空格与D不匹配时，你其实知道前面六个字符是”ABCDAB”。KMP 算法的想法是，设法利用这个已知信息，不要把”搜索位置”移回已经比较过的位置，继续把它向后移，这样就提高了效率。)
- 怎么做到把刚刚重复的步骤省略掉？可以对Str2计算出一张《部分匹配表》，这张表的产生在后面介绍
- 已知空格与D不匹配时，前面六个字符”ABCDAB”是匹配的。查表可知，最后一个匹配字符B对应的”部分匹配值”为2，因此按照下面的公式算出向后移动的位数：
-
移动位数 = 已匹配的字符数 - 对应的部分匹配值
-
因为 6 - 2 等于4，所以将搜索词向后移动 4 位。
- 因为空格与Ｃ不匹配，搜索词还要继续往后移。这时，已匹配的字符数为2（”AB”），对应的”部分匹配值”为0。所以，移动位数 = 2 - 0，结果为 2，于是将搜索词向后移 2 位。
- 因为空格与A不匹配，继续后移一位。
- 逐位比较，直到发现C与D不匹配。于是，移动位数 = 6 - 2，继续将搜索词向后移动 4 位。
- 逐位比较，直到搜索词的最后一位，发现完全匹配，于是搜索完成。如果还要继续搜索（即找出全部匹配），移动位数 = 7 - 0，再将搜索词向后移动 7 位，这里就不再重复了。
- 介绍《部分匹配表》怎么产生的?
- 先介绍前缀，后缀是什么
“部分匹配值”就是”前缀”和”后缀”的最长的共有元素的长度。以”ABCDABD”为例
”A”的前缀和后缀都为空集，共有元素的长度为0；
”AB”的前缀为[A]，后缀为[B]，共有元素的长度为0；
”ABC”的前缀为[A, AB]，后缀为[BC, C]，共有元素的长度0；
”ABCD”的前缀为[A, AB, ABC]，后缀为[BCD, CD, D]，共有元素的长度为0；
”ABCDA”的前缀为[A, AB, ABC, ABCD]，后缀为[BCDA, CDA, DA, A]，
共有元素为”A”，长度为1；
”ABCDAB”的前缀为[A, AB, ABC, ABCD, ABCDA]，
后缀为[BCDAB, CDAB, DAB, AB, B]，共有元素为”AB”，长度为2；
”ABCDABD”的前缀为[A, AB, ABC, ABCD, ABCDA, ABCDAB]，
后缀为[BCDABD, CDABD, DABD, ABD, BD, D]，共有元素的长度为0。
- ”部分匹配”的实质是，有时候，字符串头部和尾部会有重复。比如，”ABCDAB”之中有两个”AB”，那么它的”部分匹配值”就是2（”AB”的长度）。搜索词移动的时候，第一个”AB”向后移动 4 位（字符串长度-部分匹配值），就可以来到第二个”AB”的位置。
KMP算法思想分析完毕!
KMP实现
- 先得到字串的部分匹配表
- 使用部分匹配表完成KMP匹配
代码实现:
public class KMPAlgorithm {
public static void main(String[] args) {
String str1 = "BBC ABCDAB ABCDABCDABDE";
String str2 = "ABCDABD";
int[] next = kmpNext(str2);//[0,1,2,0]
int i = kmpSearch(str1, str2, next);
System.out.println(i);
}
//写出我们的kmp搜索算法
/**
* @param str1 原字符串
* @param str2 子串(需要去匹配的字符串)
* @param next 部分匹配表,是子串对应的部分匹配表
* @return 如果是-1就是没有匹配到,否则返回第一个匹配的位置
*/
public static int kmpSearch(String str1, String str2, int[] next) {
//遍历
for (int i = 0, j = 0; i < str1.length(); i++) {
//需要处理str1.charAt(i) != str2.charAt(j),去调整j的大小
//KMP算法核心点,可以验证
while (j > 0 && str1.charAt(i) != str2.charAt(j)){
j = next[j-1];
}
if (str1.charAt(i) == str2.charAt(j)) {
j++;
}
if (j == str2.length()) {//找到了 //j = 3,i =
return i - j + 1;//i长度不是固定的,j是固定的
}
}
return -1;
}
//获取到一个字符串(子串)的部分匹配值表
public static int[] kmpNext(String dest) {
//创建一个next 数组保存部分匹配值
int[] next = new int[dest.length()];
next[0] = 0;//如果字符串长度为1,那么它的部分匹配值就是0
for (int i = 1, j = 0; i < dest.length(); i++) {
//当当(dest.charAt(i) != dest.charAt(j),我们需要从next[j-1]获取新的j
//直到我们发现有dest.charAt(i) == dest.charAt(j)成立时才退出
//这是KMP算法的核心点
while (j > 0 && dest.charAt(i) != dest.charAt(j)) {
j = next[j - 1];
}
//当(dest.charAt(i) == dest.charAt(j) 满足时,部分匹配值就是+1
if (dest.charAt(i) == dest.charAt(j)) {
j++;
}
next[i] = j;
}
return next;
}
}
注意！！
KMP算法核心点:
while (j > 0 && str1.charAt(i) != str2.charAt(j)) {
j = next[j - 1];
}
对应:
贪心算法
应用场景
假设存在下面需要付费的广播台，以及广播台信号可以覆盖的地区。 如何选择最少的广播台，让所有的地区都可以接收到信号。
| 广播台 | 覆盖地区 |
|---|---|
| K1 | “北京”, “上海”, “天津” |
| K2 | “广州”, “北京”, “深圳” |
| K3 | “成都”, “上海”, “杭州” |
| K4 | “上海”, “天津” |
| K5 | “杭州”, “大连” |
-
具体介绍：
- 贪婪算法(贪心算法)是指在对问题进行求解时，在每一步选择中都采取最好或者最优(即最有利)的选择，从而希望能够导致结果是最好或者最优的算法。
- 贪婪算法所得到的结果不一定是最优的结果(有时候会是最优解)，但是都是相对近似(接近)最优解的结果
-
思路分析：
- 如何找出覆盖所有地区的广播台的集合呢，使用穷举法实现,列出每个可能的广播台的集合，这被称为幂集。假设总的有n个广播台，则广播台的组合总共有2ⁿ -1 个,假设每秒可以计算10个子集， 如图:
| 广播台数量n | 子集总数2ⁿ | 需要的时间 |
|---|---|---|
| 5 | 32 | 3.2秒 |
| 10 | 1024 | 102.4秒 |
| 32 | 4294967296 | 13.6年 |
| 100 | 1.26*100³º | 4x10²³年 |
使用贪婪算法，效率高:
目前并没有算法可以快速计算得到准备的值,使用贪婪算法,则可以得到非常
接近的解，并且效率高。选择策略上，因为需要覆盖全部地区的最小集合:
1.遍历所有的广播电台, 找到一个覆盖了最多未覆盖的地区的电台
(此电台可能包含一些已覆盖的地区，但没有关系）
2.将这个电台加入到一个集合中(比如ArrayList),
想办法把该电台覆盖的地区在下次比较时去掉。
3.重复第1步直到覆盖了全部的地区
代码实现:
public class GreedyAlgorithm {
public static void main(String[] args) {
//创建广播电台,放入到Map
HashMap<String, HashSet<String>> broadcasts = new HashMap<>();
//将各个电台放入broadcasts
HashSet<String> hashSet1 = new HashSet<>();
hashSet1.add("北京");
hashSet1.add("上海");
hashSet1.add("天津");
HashSet<String> hashSet2 = new HashSet<>();
hashSet2.add("广州");
hashSet2.add("北京");
hashSet2.add("深圳");
HashSet<String> hashSet3 = new HashSet<>();
hashSet3.add("成都");
hashSet3.add("上海");
hashSet3.add("杭州");
HashSet<String> hashSet4 = new HashSet<>();
hashSet4.add("上海");
hashSet4.add("天津");
HashSet<String> hashSet5 = new HashSet<>();
hashSet5.add("杭州");
hashSet5.add("大连");
//加入map
broadcasts.put("K1", hashSet1);
broadcasts.put("K2", hashSet2);
broadcasts.put("K3", hashSet3);
broadcasts.put("K4", hashSet4);
broadcasts.put("K5", hashSet5);
//allAreas 存放所有的地区
HashSet<String> allAreas = new HashSet<>();
allAreas.add("北京");
allAreas.add("上海");
allAreas.add("天津");
allAreas.add("广州");
allAreas.add("深圳");
allAreas.add("成都");
allAreas.add("杭州");
allAreas.add("大连");
//创建ArrayList,存放选择的电台集合
ArrayList<String> selects = new ArrayList<>();
//定义一个临时的集合,在遍历的过程中,存放遍历过程中的电台覆盖的地区和当前还没有覆盖的地区的交集
HashSet<String> tempSet = new HashSet<>();
//定义maxKey,保存在一次遍历过程中,能够覆盖 最多 未覆盖的地区 对应的电台的key
//如果maxKey不为null,则会加入到selects
String maxKey = null;
while (allAreas.size() != 0){ //如果allAreas不为0,则表示还没有覆盖到所有的地区
//每进行一次while,需要
maxKey = null;
//遍历broadcasts,取出对应key
for (String key : broadcasts.keySet()) {
//每进行一次for
tempSet.clear();
//当前这个key能够覆盖的地区
HashSet<String> ares = broadcasts.get(key);
tempSet.addAll(ares);//key=K1,ares-tempSet=北京上海天津
//allAreas-8个地区
//求出tempSet和allAreas集合的交集,共有的部分,交集会赋给tempSet
tempSet.retainAll(allAreas);
//如果当前这个集合包含的未覆盖地区的数量,比maxKey指向的集合未覆盖的地区还要多
//就需要重置maxKey
//tempSet.size() > broadcasts.get(maxKey).size() 体现出贪心算法的特性,每次都选最优的
if (tempSet.size() > 0 && (maxKey == null || tempSet.size() > broadcasts.get(maxKey).size()) ){
maxKey = key;
}
}
//maxKey != null,就应该将maxKey 加入selects
if (maxKey != null){
selects.add(maxKey);
//将maxKey指向的广播电台覆盖的地区,从allAreas去掉
allAreas.removeAll(broadcasts.get(maxKey));
}
}
System.out.println("得到的选择结果是" + selects);
}
}
贪心算法注意事项和细节:
-
贪婪算法所得到的结果不一定是最优的结果(有时候会是最优解)，但是都是相对近似(接近)最优解的结果
-
比如上题的算法选出的是K1, K2, K3, K5，符合覆盖了全部的地区
-
但是我们发现 K2, K3,K4,K5 也可以覆盖全部地区，如果K2 的使用成本低于K1,那么我们上题的 K1, K2, K3, K5 虽然是满足条件，但是并不是最优的.
普利姆算法
先看一个应用场景:
有胜利乡有7个村庄(A, B, C, D, E, F, G) ，现在需要修路把7个村庄连通
各个村庄的距离用边线表示(权) ，比如 A –B 距离 5公里
问：如何修路保证各个村庄都能连通，并且总的修建公路总里程最短?
思路: 将10条边，连接即可，但是总的里程数不是最小.
正确的思路，就是尽可能的选择少的路线，并且每条路线最小，保证总里程数最少.
最小生成树:
修路问题本质就是就是最小生成树问题， 先介绍一下最小生成树(Minimum Cost Spanning Tree)，简称MST。
给定一个带权的无向连通图,如何选取一棵生成树,使树上所有边上权的总和为最小,这叫最小生成树 N个顶点，一定有N-1条边
包含全部顶点
N-1条边都在图中
举例说明(如图:)
求最小生成树的算法主要是普里姆算法和克鲁斯卡尔算法
普利姆算法介绍：
-
普利姆(Prim)算法求最小生成树，也就是在包含n个顶点的连通图中，找出只有(n-1)条边包含所有n个顶点的连通子图，也就是所谓的极小连通子图
-
普利姆的算法如下:
-
设G=(V,E)是连通网，T=(U,D)是最小生成树，V,U是顶点集合，E,D是边的集合
-
若从顶点u开始构造最小生成树，则从集合V中取出顶点u放入集合U中，标记顶点v的visited[u]=1
-
若集合U中顶点ui与集合V-U中的顶点vj之间存在边，则寻找这些边中权值最小的边，但不能构成回路，将顶点vj加入集合U中，将边（ui,vj）加入集合D中，标记visited[vj]=1
-
重复步骤②，直到U与V相等，即所有顶点都被标记为访问过，此时D中有n-1条边
-
提示: 单独看步骤很难理解，我们通过代码来讲解，比较好理解.
-
从<A>顶点开始处理 => <A,G>
A-C [7] A-G[2] A-B[5] =><A,G> 开始 , 将A 和 G 顶点和他们相邻的还没有访问的顶点进行处理 =><A,G,B> A-C[7] A-B[5] G-B[3] G-E[4] G-F[6]
<A,G,B> 开始，将A,G,B 顶点 和他们相邻的还没有访问的顶点进行处理=><A,G,B,E> A-C[7] G-E[4] G-F[6] B-D[9]
…{A,G,B,E}->F//第4次大循环 ， 对应 边<E,F> 权值：5
{A,G,B,E,F}->D//第5次大循环 ， 对应 边<F,D> 权值：4
{A,G,B,E,F,D}->C//第6次大循环 ， 对应 边<A,C> 权值：7
=> <A,G,B,E,F,D,C>
代码实现:
public class PrimAlgorithm {
public static void main(String[] args) {
//测试看看图是否创建成功
char[] data = new char[]{'A', 'B', 'C', 'D', 'E', 'F', 'G'};
int verxs = data.length;
//邻接矩阵的关系使用二维数组表示,用10000这个较大的数,表示两个点不连通
int[][] weight = new int[][]{
//A B C D E F G
{10000, 5, 7, 10000, 10000, 10000, 2}, //A
{5, 10000, 10000, 9, 10000, 10000, 3}, //B
{7, 10000, 10000, 10000, 8, 10000, 10000}, //C
{10000, 9, 10000, 10000, 10000, 4, 10000}, //D
{10000, 10000, 8, 10000, 10000, 5, 4}, //E
{10000, 10000, 10000, 4, 5, 10000, 6}, //F
{2, 3, 10000, 10000, 4, 6, 10000},}; //G
//创建MGraph对象
MGraph graph = new MGraph(verxs);
MinTree minTree = new MinTree();
minTree.createGraph(graph, verxs, data, weight);
//输出
minTree.showGraph(graph);
//测试普里姆算法
minTree.prim(graph, 0);
}
}
//创建最小生成树->村庄的图
class MinTree {
//创建图的邻接矩阵
/**
* @param graph 图对象
* @param verxs 图对应的顶点个数
* @param data 图的各个顶点的值
* @param weight 图的邻接矩阵
*/
public void createGraph(MGraph graph, int verxs, char data[], int[][] weight) {
int i, j;
for (i = 0; i < verxs; i++) {//顶点
graph.data[i] = data[i];//把传进来顶点的值传给图里面的data
for (j = 0; j < verxs; j++) {
graph.weight[i][j] = weight[i][j];
}
}
}
//显示图的方法
public void showGraph(MGraph graph) {
for (int[] link : graph.weight) {
System.out.println(Arrays.toString(link));
}
}
//编写prim算法,得到最小生成树
/**
*
* @param graph 图
* @param v 表示从图的第几个顶点开始生成'A'->0 'B'->1...
*/
public void prim(MGraph graph,int v){
//visited[] 标记结点(顶点)是否被访问过
int visited[] = new int[graph.verxs];
//visited[] 默认元素的值都是0,表示没有访问过
/*for (int i = 0; i < graph.verxs; i++) {
visited[i] = 0;
}*/
//把当前这结点标记为已访问
visited[v] = 1;
//用h1和h2记录两个顶点的下标
int h1 = -1;
int h2 = -1;
int minWeight = 10000;//将minWeight 初始成一个大数,后面在遍历过程中,会被替换
for (int k = 1; k < graph.verxs; k++) {
//因为有graph.verxs个顶点,普里姆算法结束后,有graph.verxs-1边
for (int i = 0; i < graph.verxs; i++) {//i 结点表示被访问过的结点
for (int j = 0; j < graph.verxs; j++) {//j 结点表示没有被访问过的结点
//这个在确保每一次生成的子图,和哪个结点的距离最近
//确保一个访问过的结点和一个没有访问过的结点探视,同时当前这个结点小于最小值,
//那就把真正最小的权值记录下来,并且记录
if (visited[i] == 1 && visited[j] == 0 && graph.weight[i][j] < minWeight){
//替换minWeight(寻找已经访问过的结点和未访问过的结点间的权值最小的边)
minWeight = graph.weight[i][j];
h1 = i;
h2 = j;
//这层循环出去,这个子图中哪两个结点是最近的,同时还保存了i和j结点。
}
}
}
//找到一条边是最小的
System.out.println("边<" + graph.data[h1] + "," + graph.data[h2] + "> 权值:" + minWeight);
//将当前这个结点标记为已经访问
visited[h2] = 1;
//重置minWeight为最大值
minWeight = 10000;
}
}
}
class MGraph {
int verxs;//表示图的结点个数
char[] data;//存放结点数据
int[][] weight;//存放边,就是我们的邻接矩阵
public MGraph(int verxs) {
this.verxs = verxs;
data = new char[verxs];
weight = new int[verxs][verxs];
}
}
克鲁斯卡尔算法分析
来看一个应用场景
- 某城市新增7个站点(A, B, C, D, E, F, G) ，现在需要修路把7个站点连通
- 各个站点的距离用边线表示(权) ，比如 A – B 距离 12公里
- 问：如何修路保证各个站点都能连通，并且总的修建公路总里程最短?
克鲁斯卡尔算法介绍
- 克鲁斯卡尔(Kruskal)算法，是用来求加权连通图的最小生成树的算法。
- 基本思想：按照权值从小到大的顺序选择n-1条边，并保证这n-1条边不构成回路
- 具体做法：首先构造一个只含n个顶点的森林，然后依权值从小到大从连通网中选择边加入到森林中，并使森林中不产生回路，直至森林变成一棵树为止。
克鲁斯卡尔算法图解说明
在含有n个顶点的连通图中选择n-1条边，构成一棵极小连通子图，并使该连通子图中n-1条边上权值之和达到最小，则称其为连通网的最小生成树。
- 例如，对于如上图G4所示的连通网可以有多棵权值总和不相同的生成树。
- 以上图G4为例，来对克鲁斯卡尔进行演示(假设，用数组R保存最小生成树结果)。
第1步：将边<E,F>加入R中。
边<E,F>的权值最小，因此将它加入到最小生成树结果R中。
第2步：将边<C,D>加入R中。
上一步操作之后，边<C,D>的权值最小，因此将它加入到最小生成树结果R中。
第3步：将边<D,E>加入R中。
上一步操作之后，边<D,E>的权值最小，因此将它加入到最小生成树结果R中。
第4步：将边<B,F>加入R中。
上一步操作之后，边<C,E>的权值最小，但<C,E>会和已有的边构成回路；因此，跳过边<C,E>。同理，跳过边<C,F>。将边<B,F>加入到最小生成树结果R中。
第5步：将边<E,G>加入R中。
上一步操作之后，边<E,G>的权值最小，因此将它加入到最小生成树结果R中。
第6步：将边<A,B>加入R中。
上一步操作之后，边<F,G>的权值最小，但<F,G>会和已有的边构成回路；因此，跳过边<F,G>。同理，跳过边<B,C>。将边<A,B>加入到最小生成树结果R中。
此时，最小生成树构造完成！它包括的边依次是：<E,F> <C,D> <D,E> <B,F> <E,G> <A,B>。
克鲁斯卡尔算法分析
根据前面介绍的克鲁斯卡尔算法的基本思想和做法，我们能够了解到，克鲁斯卡尔算法重点需要解决的以下两个问题：
- 问题一 对图的所有边按照权值大小进行排序。
- 问题二 将边添加到最小生成树中时，怎么样判断是否形成了回路。
问题一很好解决，采用排序算法进行排序即可。
问题二，处理方式是：记录顶点在"最小生成树"中的终点，顶点的终点是"在最小生成树中与它连通的最大顶点"。然后每次需要将一条边添加到最小生存树时，判断该边的两个顶点的终点是否重合，重合的话则会构成回路。
如何判断是否构成回路-举例说明？
在将<E,F> <C,D> <D,E>加入到最小生成树R中之后，这几条边的顶点就都有了终点：
(01) C的终点是F。
(02) D的终点是F。
(03) E的终点是F。
(04) F的终点是F。
关于终点的说明：
1.就是将所有顶点按照从小到大的顺序排列好之后； 某个顶点的终点就是"与它连通的最大顶点"。
2.因此，接下来，虽然<C,E>是权值最小的边。
但是C和E的终点都是F，即它们的终点相同，
因此，将<C,E>加入最小生成树的话，会形成回路。
这就是判断回路的方式。也就是说，
我们加入的边的两个顶点不能都指向同一个终点，否则将构成回路。
代码实现:
public class KruskalCase {
private int edgeNum;//记录边的个数
private char[] vertexs;//顶点数组
private int[][] matrix;//邻接矩阵
//使用 INF 表示两个顶点不能连通
private static final int INF = Integer.MAX_VALUE;
public static void main(String[] args) {
char[] vertexs = {'A', 'B', 'C', 'D', 'E', 'F', 'G'};
int matrix[][] = {
//自己和自己连,是一个对角线用0表示.
/*A*//*B*//*C*//*D*//*E*//*F*//*G*/
/*A*/ {0, 12, INF, INF, INF, 16, 14},
/*B*/ {12, 0, 10, INF, INF, 7, INF},
/*C*/ {INF, 10, 0, 3, 5, 6, INF},
/*D*/ {INF, INF, 3, 0, 4, INF, INF},
/*E*/ {INF, INF, 5, 4, 0, 2, 8},
/*F*/ {16, 7, 6, INF, 2, 0, 9},
/*G*/ {14, INF, INF, INF, 8, 9, 0}};
//创建一个KruskalCase 对象实例
KruskalCase kruskalCase = new KruskalCase(vertexs, matrix);
//输出构建的是否正确
kruskalCase.print();
EData[] edges = kruskalCase.getEdges();
System.out.println("排序前=" + Arrays.toString(edges));
kruskalCase.sortEdge(edges);
System.out.println("排序后=" + Arrays.toString(edges));
kruskalCase.Kruskal();
}
//构造器
public KruskalCase(char[] vertexs, int[][] matrix) {
//初始化顶点数和边的个数
int vlen = vertexs.length;//vlen数组个数
//初始化顶点,复制拷贝的方式
this.vertexs = new char[vlen];
//vertexs顶点数组
for (int i = 0; i < vertexs.length; i++) {
//把传进来的vertexs赋给属性里的vertexs
this.vertexs[i] = vertexs[i];
}
//第二种方式:
//this.vertexs = vertexs;这样一步到位.
//第一种方式就是这个构造器里面的,在里面的改变,不影响外面传进来的数组.
//初始化边,使用的是复制拷贝
this.matrix = new int[vlen][vlen];
for (int i = 0; i < vlen; i++) {
for (int j = 0; j < vlen; j++) {
this.matrix[i][j] = matrix[i][j];
}
}
//统计边的条数
for (int i = 0; i < vlen; i++) {
for (int j = i + 1; j < vlen; j++) {
if (this.matrix[i][j] != INF) {//!=INF说明这条边是有效的
edgeNum++;
}
}
}
}
public void Kruskal() {
int index = 0;//表示最后结果数组的索引
//用于保存"已有最小生成树" 中的每个顶点在最小生成树中的终点
int[] ends = new int[edgeNum];
//创建结果数组,保存最后的最小生成树
EData[] rets = new EData[edgeNum];
//获取图中 所有的边的集合 , 一共有12条
EData[] edges = getEdges();
System.out.println("图的边集合=" + Arrays.toString(edges) + "共" + edges.length + "条边");
//按照边的权值大小进行排序(从小到大)
sortEdge(edges);
//遍历edges数组
//将边添加到最小生成树中,判断准备加入的边是否形成了回路,如果没有,就加入rets,否则不能加入.
for (int i = 0; i < edgeNum; i++) {
//获取到第i条边的第1个顶点(起点)
int p1 = getPosition(edges[i].start);//e-> p1=4
//获取到第i条边的第2个顶点(终点)
int p2 = getPosition(edges[i].end);//f-> p2=5
//获取p1、p2在已有的最小生成树中的终点(逐步形成)
int m1 = getEnd(ends, p1);//m1 = 4,还没加入到最小生成树时,我们认为p1的终点就是它自己
int m2 = getEnd(ends, p2);//m2 = 5.
//是否构成回路
if (m1 != m2){//没有构成回路
ends[m1] = m2;//设置m1 在"已有最小生成树"中的终点<E,F> [0,0,0,0,5,0,0,0,0,0,0,0]
// 下标为4的顶点,在最小生成树中的终点为5.E的终点是F.
rets[index++] = edges[i];//这时,有一条边加入到rets
}
}
//统计并打印"最小生成树",输出rets
System.out.println("最小生成树为");
for (int i = 0; i < index; i++) {
System.out.println(rets[i]);
}
}
//打印邻接矩阵
public void print() {
System.out.println("邻接矩阵为:");
for (int i = 0; i < vertexs.length; i++) {
for (int j = 0; j < vertexs.length; j++) {
System.out.printf("%15d", matrix[i][j]);
}
System.out.println();
}
}
//对边进行排序处理,冒泡排序
/**
* 功能:用冒泡对边进行排序
*
* @param edges 边的集合
*/
private void sortEdge(EData[] edges) {
for (int i = 0; i < edges.length - 1; i++) {
for (int j = 0; j < edges.length - i - 1; j++) {
if (edges[j].weight > edges[j + 1].weight) {//交换
EData tmp = edges[j];
edges[j] = edges[j + 1];
edges[j + 1] = tmp;
}
}
}
}
/**
* @param ch 顶点的值,比如'A','B'
* @return 返回ch顶点对应的下标, 如果找不到, 返回-1
*/
private int getPosition(char ch) {
for (int i = 0; i < vertexs.length; i++) {
if (vertexs[i] == ch) {//找到
return i;
}
}
//找不到,返回-1
return -1;
}
/**
* 功能:获取图中的边,放到EData[] 数组中,后面我们需要遍历该数组
* 通过matrix 邻接矩阵
* EData[] 形式[['A','B',12],['B','F',7],...]
*
* @return 形式[['A','B',12],['B','F',7],...]
*/
private EData[] getEdges() {
int index = 0;
EData[] edges = new EData[edgeNum];//数组长度为边的个数
for (int i = 0; i < vertexs.length; i++) {
for (int j = i + 1; j < vertexs.length; j++) {
if (matrix[i][j] != INF){
edges[index] = new EData(vertexs[i],vertexs[j],matrix[i][j]);
index++;
}
}
}
return edges;
}
/**
* 功能:获取下标为i的顶点的终点(动态加入),用于后面判断两个顶点的终点是否相同
* @param ends : 顶点 对应的终点-数组,顶点下标为i的终点=ends[i]
* ends 数组是在遍历过程中,逐步形成的.
* @param i : 表示传入的顶点对应的下标.
* @return 返回的就是下标为i的顶点 所对应的终点的下标.
*/
private int getEnd(int[] ends,int i) {
while (ends[i] != 0){
i = ends[i];
}
return i;
}
}
//创建一个类EData,它的对象实例就表示一条边
class EData {
char start;//边的起点,边的一个点
char end;//边的终点,边的另外一个点
int weight;//边的权值
//构造器
public EData(char start, char end, int weight) {
this.start = start;
this.end = end;
this.weight = weight;
}
//重写toString,便于输出边的信息
@Override
public String toString() {
return "EData{" +
"<" + start +
", " + end +
">= =" + weight +
'}';
}
}
迪杰斯特拉算法
看一个应用场景和问题：
- 战争时期，胜利乡有7个村庄(A, B, C, D, E, F, G) ，现在有六个邮差，从G点出发，需要分别把邮件分别送到 A, B, C , D, E, F 六个村庄
- 各个村庄的距离用边线表示(权) ，比如 A – B 距离 5公里
- 问：如何计算出G村庄到 其它各个村庄的最短距离?
如果从其它点出发到各个点的最短距离又是多少?
迪杰斯特拉(Dijkstra)算法是典型最短路径算法，用于计算一个结点到其他结点的最短路径。
它的主要特点是以起始点为中心向外层层扩展 (广度优先 搜索思想)，直到扩展到终点为止。
迪杰斯特拉(Dijkstra)算法过程
- 设置出发顶点为v，顶点集合V{v1,v2,vi…}，v到V中各顶点的距离构成距离集合Dis，Dis{d1,d2,di…}，Dis集合记录着v到图中各顶点的距离(到自身可以看作0，v到vi距离对应为di)
- 从Dis中选择值最小的di并移出Dis集合，同时移出V集合中对应的顶点vi，此时的v到vi即为最短路径
- 更新Dis集合，更新规则为：比较v到V集合中顶点的距离值，与v通过vi到V集合中顶点的距离值，保留值较小的一个(同时也应该更新顶点的前驱节点为vi，表明是通过vi到达的)
- 重复执行两步骤，直到最短路径顶点为目标顶点即可结束
代码实现:
public class DijkstraAlgorithm {
public static void main(String[] args) {
char[] vertex = {'A', 'B', 'C', 'D', 'E', 'F', 'G'};
//邻接矩阵
int[][] matrix = new int[vertex.length][vertex.length];
final int N = 61234;//表示不可以连接
matrix[0] = new int[]{N, 5, 7, N, N, N, 2};
matrix[1] = new int[]{5, N, N, 9, N, N, 3};
matrix[2] = new int[]{7, N, N, N, 8, N, N};
matrix[3] = new int[]{N, 9, N, N, N, 4, N};
matrix[4] = new int[]{N, N, 8, N, N, 5, 4};
matrix[5] = new int[]{N, N, N, 4, 5, N, 6};
matrix[6] = new int[]{2, 3, N, N, 4, 6, N};
//创建 Graph对象
Graph graph = new Graph(vertex, matrix);
//测试,看看图的邻接矩阵对不对
//graph.showGraph();
//测试迪杰斯特拉算法
graph.dijkstra(6);
graph.showDijkstra();
}
}
class Graph {
private char[] vertex;//顶点数组
private int[][] matrix;//邻接矩阵
private VisitedVertex vv;
//构造器
public Graph(char[] vertex, int[][] matrix) {
this.vertex = vertex;
this.matrix = matrix;
}
//显示图
public void showGraph() {
for (int[] link : matrix) {
System.out.println(Arrays.toString(link));
}
}
//显示结果
public void showDijkstra() {
vv.show();
}
//迪杰斯特拉算法实现
/**
* @param index 表示出发顶点对应的下标
*/
public void dijkstra(int index) {
vv = new VisitedVertex(vertex.length, index);
update(index);//更新index顶点到周围顶点的距离和前驱顶点
for (int j = 1; j < vertex.length; j++) {
index = vv.updateArr();//选择并返回新的访问顶点
update(index);
}
}
//更新index下标顶点到周围顶点的距离和周围顶点的前驱顶点
private void update(int index) {
int len = 0;
//根据遍历我们的邻接矩阵的matrix[index]行
for (int j = 0; j < matrix[index].length; j++) {
//len : 出发顶点到index顶点的距离 + 从index顶点到j顶点的距离 的和
len = vv.getDis(index) + matrix[index][j]; //获取出发顶点到j的距离
//如果j顶点没有被访问过,并且len 小于出发顶点到j距离的距离
if (!vv.in(j) && len < vv.getDis(j)) {
vv.updatePre(j, index);//更新j顶点的前驱为index顶点
vv.updateDis(j, len);//更新出发顶点到j顶点的距离
}
}
}
}
class VisitedVertex {//已访问顶点集合
//记录各个顶点是否访问过->1表示访问过,0未访问,会动态更新
public int[] already_arr;
//每个下标对应的值为前一个顶点下标,会动态更新
public int[] pre_visited;
//记录出发顶点到其他所有顶点的距离,出发顶点为G,就记录G到其他顶点的距离,会动态更新,求得最短距离就会存放到dis
public int[] dis;
//构造器
/**
* @param length 表示顶点的个数
* @param index 出发的顶点对应的下标,比如G顶点,下标就是6
*/
public VisitedVertex(int length, int index) {
this.already_arr = new int[length];
this.pre_visited = new int[length];
this.dis = new int[length];
//初始化 dis数组
Arrays.fill(dis, 61234);
this.already_arr[index] = 1;//设置出发顶点被访问过
this.dis[index] = 0;//设置出发顶点的访问距离为0
}
/**
* 功能:判断index顶点是否被访问过
*
* @param index
* @return 如果访问过就返回true, 没访问过返回false
*/
public boolean in(int index) {
return already_arr[index] == 1;
}
/**
* 功能:更新出发顶点到index顶点的距离
*
* @param index
* @param len
*/
public void updateDis(int index, int len) {
//前提:前驱结点到你的这个点距离+你这个点到下一个点的距离 < 顶点到你这个点到下一个点的距离
// EG: (G->A->D) < 上一次的(G->D)的距离 时,调用这个方法
dis[index] = len;
}
/**
* 功能:更新pre这个顶点的前驱为index顶点
* <p>
* 下标pre的前驱结点->index
*
* @param pre
* @param index
*/
public void updatePre(int pre, int index) {
pre_visited[pre] = index;
}
/**
* 功能:返回出发顶点为index顶点的距离
*
* @param index
*/
public int getDis(int index) {
return dis[index];
}
/**
* 功能:继续选择并返回新的访问顶点,比如这里的G 完后,就是A点作为新的访问顶点(注意不是出发顶点)
* @return
*/
public int updateArr() {
int min = 61234, index = 0;
for (int i = 0; i < already_arr.length; i++) {
if (already_arr[i] == 0 && dis[i] < min) { //i还没有被访问过
min = dis[i];
index = i;
}
}
//同时更新index顶点被访问过
already_arr[index] = 1;
return index;
}
//显示最后的结果
//即将三个数组的情况输出
public void show() {
System.out.println("========================");
//输出already_arr
for (int i :already_arr) {
System.out.print(i + " ");
}
System.out.println();
//输出pre_visited
for (int i :
pre_visited) {
System.out.print(i + " ");
}
System.out.println();
//输出dis
for (int i :
dis) {
System.out.print( i + " ");
}
System.out.println();
//为了提高代码可读性,处理显示最短距离格式
char[] vertex = {'A', 'B', 'C', 'D', 'E', 'F', 'G'};
int count = 0;
for (int i :
dis) {
if (i != 61234) {
System.out.print(vertex[count] + "[" + i + "]");
}else {
System.out.println("N ");
}
count++;
}
System.out.println();
}
}
弗洛伊德算法
弗洛伊德算法介绍
和Dijkstra算法一样，弗洛伊德(Floyd)算法也是一种用于寻找给定的加权图中顶点间最短路径的算法。该算法名称以创始人之一、1978年图灵奖获得者、斯坦福大学计算机科学系教授罗伯特·弗洛伊德命名
弗洛伊德算法(Floyd)计算图中各个顶点之间的最短路径
迪杰斯特拉算法用于计算图中某一个顶点到其他顶点的最短路径。
弗洛伊德算法 VS 迪杰斯特拉算法：
迪杰斯特拉算法通过选定的被访问顶点，求出从出发访问顶点到其他顶点的最短路径；
弗洛伊德算法中每一个顶点都是出发访问点，所以需要将每一个顶点看做被访问顶点，求出从每一个顶点到其他顶点的最短路径。
弗洛伊德(Floyd)算法图解分析:
- 设置顶点vi到顶点vk的最短路径已知为Lik，顶点vk到vj的最短路径已知为Lkj，顶点vi到vj的路径为Lij，则vi到vj的最短路径为：min((Lik+Lkj),Lij)，vk的取值为图中所有顶点，则可获得vi到vj的最短路径
- 至于vi到vk的最短路径Lik或者vk到vj的最短路径Lkj，是以同样的方式获得
- 弗洛伊德(Floyd)算法图解分析-举例说明:
-
第一轮循环中，以A(下标为：0)作为中间顶点，距离表和前驱关系更新为：
分析如下： -
以A顶点作为中间顶点是，B->A->C的距离由N->9，同理C到B；C->A->G的距离由N->12，同理G到C
-
更换中间顶点，循环执行操作，直到所有顶点都作为中间顶点更新后，计算结束。
离表的数据来源：
A A A F G G A
(A到A的最短路径是0) (A到B的最短路径是5) (A到C的最短路径是7) (A到D的最短路径是12 (A到E的最短路径是6) (A到F的最短路径是8) (A到G的最短路径是2)
B B A B G G B
(B到A的最短路径是5) (B到B的最短路径是0) (B到C的最短路径是12) (B到D的最短路径是9) (B到E的最短路径是7) (B到F的最短路径是9) (B到G的最短路径是3)
C A C F C E A
(C到A的最短路径是7) (C到B的最短路径是12) (C到C的最短路径是0) (C到D的最短路径是17) (C到E的最短路径是8) (C到F的最短路径是13) (C到G的最短路径是9)
G D E D F D F
(D到A的最短路径是12) (D到B的最短路径是9) (D到C的最短路径是17) (D到D的最短路径是0) (D到E的最短路径是9) (D到F的最短路径是4) (D到G的最短路径是10)
G G E F E E E
(E到A的最短路径是6) (E到B的最短路径是7) (E到C的最短路径是8) (E到D的最短路径是9) (E到E的最短路径是0) (E到F的最短路径是5) (E到G的最短路径是4)
G G E F F F F
(F到A的最短路径是8) (F到B的最短路径是9) (F到C的最短路径是13) (F到D的最短路径是4) (F到E的最短路径是5) (F到F的最短路径是0) (F到G的最短路径是6)
G G A F G G G
(G到A的最短路径是2) (G到B的最短路径是3) (G到C的最短路径是9) (G到D的最短路径是10) (G到E的最短路径是4) (G到F的最短路径是6) (G到G的最短路径是0)
代码实现:
public class FloydAlgorithm {
public static void main(String[] args) {
//测试->图是否创建成功
char[] vertex = {'A', 'B', 'C', 'D', 'E', 'F', 'G'};
//创建邻接矩阵
int[][] matrix = new int[vertex.length][vertex.length];
final int N = 61234;
matrix[0] = new int[]{0, 5, 7, N, N, N, 2};
matrix[1] = new int[]{5, 0, N, 9, N, N, 3};
matrix[2] = new int[]{7, N, 0, N, 8, N, N};
matrix[3] = new int[]{N, 9, N, 0, N, 4, N};
matrix[4] = new int[]{N, N, 8, N, 0, 5, 4};
matrix[5] = new int[]{N, N, N, 4, 5, 0, 6};
matrix[6] = new int[]{2, 3, N, N, 4, 6, 0};
//创建 Graph 对象
Graph graph = new Graph(vertex.length, matrix, vertex);
//调用floyd算法
graph.floyd();
graph.show();
}
}
//创建图
class Graph {
private char[] vertex;//存放顶点的数组
private int[][] dis;// 保存,从各个顶点出发到其他顶点的距离,最后的结果,也是保留在该数组
private int[][] pre;// 保存到达目标顶点的前驱顶点
//构造器
/**
* @param length 大小
* @param matrix 邻接矩阵
* @param vertex 顶点数组
*/
public Graph(int length, int[][] matrix, char[] vertex) {
this.vertex = vertex;
this.dis = matrix;
this.pre = new int[length][length];
//对pre数组初始化,注意存放的是前驱顶点的下标(并不是'A','B')
for (int i = 0; i < length; i++) {
Arrays.fill(pre[i], i);
}
}
//显示pre数组和dis数组
public void show() {
//为了便于阅读,优化输出
char[] vertex = {'A', 'B', 'C', 'D', 'E', 'F', 'G'};
for (int k = 0; k < dis.length; k++) {
//先将pre数组输出的一行数据
for (int i = 0; i < dis.length; i++) {
System.out.print(vertex[pre[k][i]] + " ");
}
System.out.println();
//输出dis数组的一行数据
for (int i = 0; i < dis.length; i++) {
System.out.print("(" + vertex[k] + "到" + vertex[i] + "的最短路径是" + dis[k][i] + ")" + " ");
}
System.out.println();
System.out.println();
}
}
//弗洛伊德算法
public void floyd() {
int len = 0;//变量保存距离
//对从中间顶点遍历,k就是中间顶点的下标[图片中红字下方部分]
for (int k = 0; k < dis.length; k++) {//['A', 'B', 'C', 'D', 'E', 'F', 'G']
//从i顶点开始出发,['A', 'B', 'C', 'D', 'E', 'F', 'G']
for (int i = 0; i < dis.length; i++) {
//到达j顶点,['A', 'B', 'C', 'D', 'E', 'F', 'G']
for (int j = 0; j < dis.length; j++) {
len = dis[i][k] + dis[k][j];//计算出i->k的距离+k->j的距离|从i顶点出发,经过k到达j的距离
if (len < dis[i][j]) {//如果len < dis[i][j]直连距离,就更新I->J距离
dis[i][j] = len;
pre[i][j] = pre[k][j];//更新前驱
}
}
}
}
}
}
回溯算法
算法介绍：
- 马踏棋盘算法也被称为骑士周游问题
- 将马随机放在国际象棋的8×8棋盘Board[0～7][0～7]的某个方格中，马按走棋规则(马走日字)进行移动。要求每个方格只进入一次，走遍棋盘上全部64个方格
- 游戏演示: 小游戏演示
- 马踏棋盘问题(骑士周游问题)实际上是图的深度优先搜索(DFS)的应用。
- 如果使用回溯（就是深度优先搜索）来解决，假如马儿踏了53个点，如图：走到了第53个，坐标（1,0），发现已经走到尽头，没办法，那就只能回退了，查看其他的路径，就在棋盘上不停的回溯…
，思路分析+代码实现- 分析第一种方式的问题，并使用贪心算法（greedyalgorithm）进行优化。解决马踏棋盘问题.
- 使用前面的游戏来验证算法是否正确。
骑士周游问题的解决步骤和思路:
1. 创建棋盘 chessBoard , 是一个二维数组
2. 将当前位置设置为已经访问，然后根据当前位置，计算马儿还能走哪些位置，
并放入到一个集合中(ArrayList), 最多有8个位置， 每走一步，就使用step+1
3. 遍历ArrayList中存放的所有位置，看看哪个可以走通,
如果走通，就继续，走不通，就回溯.
4. 判断马儿是否完成了任务，使用 step 和应该走的步数比较 ，
如果没有达到数量，则表示没有完成任务，将整个棋盘置0
注意：马儿不同的走法（策略），会得到不同的结果，效率也会有影响(优化)
//创建一个Point
Point p1 = new Point();
if((p1.x = curPoint.x - 2) >= 0 && (p1.y = curPoint.y -1) >= 0) {
ps.add(new Point(p1));
}
使用贪心算法对原来的算法优化
1。 我们获取当前位置，可以走的下一个位置的集合
//获取当前位置可以走的下一个位置的集合
ArrayList<Point> ps = next(new Point(column, row));
2. 我们需要对 ps 中所有的Point 的下一步的所有集合的数目，进行非递减排序,就ok ,
9， 7， 6， 5， 3， 2 ， 1 //递减排序
1, 2, 3, 4,5,6, 10, //递增排序
1, 2, 2, 2, 3,3, 4, 5, 6 // 非递减
9， 7， 6，6, 6, 5，5, 3， 2 ， 1 //非递增
代码实现:
马儿不同的走法（策略），会得到不同的结果，效率也会有影响(可以优化).
- 我们应该选择,马儿走的下一步的下一步,方案较少的步骤,因为方案越多回溯越多,时间耗费越多.
- 比如:走0有五个选择,走1有两个选择,我们应该走1.
- 我们可以使用贪心算法优化此代码.
（注意!下方代码未优化,优化代码在最后）
public class HorseChessboard {
//1. 创建棋盘 chessBoard , 是一个二维数组
private static int X;//棋盘的列数
private static int Y;//棋盘的行数
//创建一个数组,标记棋盘的各个位置是否被访问过
private static boolean visited[];
//使用一个属性,标记是否棋盘的所有位置都被访问过
private static boolean finished;//如果为true,表示成功
public static void main(String[] args) {
//测试骑士周游算法是否正确
X = 8;
Y = 8;
int row = 1;//马儿初始位置,从1开始编号
int column = 1;//马儿初始位置,从1开始编号
//创建棋盘
int[][] chessboard = new int[X][Y];
visited = new boolean[X * Y];//初始值都是false
//测试一下耗时
long start = System.currentTimeMillis();
traversalChessboard(chessboard, row - 1, column - 1, 1);
long end = System.currentTimeMillis();
System.out.println("程序耗时=" + (end - start) / 1000);
//输出棋盘的最后情况
for (int rows[] : chessboard
) {
for (int step : rows
) {
System.out.print(step + "\t");
}
System.out.println();
}
}
/**
* 完成骑士周游的算法
*
* @param chessboard 棋盘.
* @param row 马儿当前的行位置,从0开始.
* @param column 马儿当前的列位置,从0开始.
* @param step 是第几步,初始位置就是第1步.
*/
public static void traversalChessboard(int[][] chessboard, int row, int column, int step) {//每走一步，就使用step+1
chessboard[row][column] = step;
// row = 4, X = 8, column = 4, =4 * 8 + 4 = 36 + 1 - 1 = 36(因为从1开始,但是一维数组是从0开始,所以是36).
visited[row * X + column] = true;//标记该位置已经访问
//获取当前位置 可以走的下一个位置的集合
ArrayList<Point> ps = next(new Point(column, row));
//3.遍历ps中存放的所有位置，看看哪个可以走通,
//如果走通，就继续，走不通，就回溯.
while (!ps.isEmpty()) {
Point p = ps.remove(0);//把下一个可以走的位置取出来
//判断该点是否已经访问过/导入
if (!visited[p.y * X + p.x]) {//说明还没有访问过
traversalChessboard(chessboard, p.y, p.x, step + 1);
}
}
/*4. 判断马儿是否完成了任务，使用 step 和应该走的步数比较 ，
如果没有达到数量，则表示没有完成任务，将整个棋盘置0
说明: step < X * Y 成立的情况有两种
1.棋盘到目前为止,仍然没有走完
2.棋盘处于处于一个回溯过程
*/
if (step < X * Y && !finished) {//说明还没完成
chessboard[row][column] = 0;
visited[row * X + column] = false;
} else {
finished = true;
}
}
//2. 将当前位置设置为已经访问，然后根据当前位置，计算马儿还能走哪些位置，
//并放入到一个集合中(ArrayList), 最多有8个位置， 每走一步，就使用step+1
/**
* 功能:根据当前位置(Point对象)，计算马儿还能走哪些位置(Point),并放入一个集合中(ArrayList),最多有8个位置.
*
* @param curPoint
* @return
*/
public static ArrayList<Point> next(Point curPoint) {
//创建一个ArrayList
ArrayList<Point> ps = new ArrayList<Point>();
//创建一个Point
Point p1 = new Point();
//这里不能用if else!!
//表示马儿可以走 5(左2上1) 这个位置
if ((p1.x = curPoint.x - 2) >= 0 && (p1.y = curPoint.y - 1) >= 0) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 6(左1上2)这个位置
if ((p1.x = curPoint.x - 1) >= 0 && (p1.y = curPoint.y - 2) >= 0) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 7(右1上2)这个位置
if ((p1.x = curPoint.x + 1) < X && (p1.y = curPoint.y - 2) >= 0) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 0(右2上1)这个位置
if ((p1.x = curPoint.x + 2) < X && (p1.y = curPoint.y - 1) >= 0) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 1(右2下1)这个位置
if ((p1.x = curPoint.x + 2) < X && (p1.y = curPoint.y + 1) < Y) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 2(右1下2)这个位置
if ((p1.x = curPoint.x + 1) < X && (p1.y = curPoint.y + 2) < Y) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 3(左1下2)这个位置
if ((p1.x = curPoint.x - 1) >= 0 && (p1.y = curPoint.y + 2) < Y) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 4(左2下1)这个位置
if ((p1.x = curPoint.x - 2) >= 0 && (p1.y = curPoint.y + 1) < Y) {
ps.add(new Point(p1));
}
return ps;
}
}
代码优化:
使用贪心算法对原来的代码优化
1. 我们获取当前位置,可以走的下一个位置的集合(ps)
ArrayList<Point> ps = next(new Point(column, row));
2. 我们需要对ps中所有的point的下一步的所有集合的数目,进行非递减排序
3. 递减：7，6，5，3，2
递增：1，2，3，5，6
非递减：1，2，2，3，3，5
非递增：9，8，8，7，7，6
public class HorseChessboard {
//1. 创建棋盘 chessBoard , 是一个二维数组
private static int X;//棋盘的列数
private static int Y;//棋盘的行数
//创建一个数组,标记棋盘的各个位置是否被访问过
private static boolean visited[];
//使用一个属性,标记是否棋盘的所有位置都被访问过
private static boolean finished;//如果为true,表示成功
public static void main(String[] args) {
//测试骑士周游算法是否正确
X = 8;
Y = 8;
int row = 1;//马儿初始位置,从1开始编号
int column = 1;//马儿初始位置,从1开始编号
//创建棋盘
int[][] chessboard = new int[X][Y];
visited = new boolean[X * Y];//初始值都是false
//测试一下耗时
long start = System.currentTimeMillis();
traversalChessboard(chessboard, row - 1, column - 1, 1);
long end = System.currentTimeMillis();
System.out.println("程序耗时=" + (end - start) + "ms");
//输出棋盘的最后情况
for (int rows[] : chessboard
) {
for (int step : rows
) {
System.out.print(step + "\t");
}
System.out.println();
}
}
/**
* 完成骑士周游的算法
*
* @param chessboard 棋盘.
* @param row 马儿当前的行位置,从0开始.
* @param column 马儿当前的列位置,从0开始.
* @param step 是第几步,初始位置就是第1步.
*/
public static void traversalChessboard(int[][] chessboard, int row, int column, int step) {//每走一步，就使用step+1
chessboard[row][column] = step;
// row = 4, X = 8, column = 4, =4 * 8 + 4 = 36 + 1 - 1 = 36(因为从1开始,但是一维数组是从0开始,所以是36).
visited[row * X + column] = true;//标记该位置已经访问
//获取当前位置 可以走的下一个位置的集合
ArrayList<Point> ps = next(new Point(column, row));
//对ps进行排序,排序的规则就是对ps的所有的Point对象的下一步的位置的数目,进行非递减排序
sort(ps);
//3.遍历ps中存放的所有位置，看看哪个可以走通,
//如果走通，就继续，走不通，就回溯.
while (!ps.isEmpty()) {
Point p = ps.remove(0);//把下一个可以走的位置取出来
//判断该点是否已经访问过/导入
if (!visited[p.y * X + p.x]) {//说明还没有访问过
traversalChessboard(chessboard, p.y, p.x, step + 1);
}
}
/*4. 判断马儿是否完成了任务，使用 step 和应该走的步数比较 ，
如果没有达到数量，则表示没有完成任务，将整个棋盘置0
说明: step < X * Y 成立的情况有两种
1.棋盘到目前为止,仍然没有走完
2.棋盘处于处于一个回溯过程
*/
if (step < X * Y && !finished) {//说明还没完成
chessboard[row][column] = 0;
visited[row * X + column] = false;
} else {
finished = true;
}
}
//2. 将当前位置设置为已经访问，然后根据当前位置，计算马儿还能走哪些位置，
//并放入到一个集合中(ArrayList), 最多有8个位置， 每走一步，就使用step+1
/**
* 功能:根据当前位置(Point对象)，计算马儿还能走哪些位置(Point),并放入一个集合中(ArrayList),最多有8个位置.
*
* @param curPoint
* @return
*/
public static ArrayList<Point> next(Point curPoint) {
//创建一个ArrayList
ArrayList<Point> ps = new ArrayList<Point>();
//创建一个Point
Point p1 = new Point();
//这里不能用if else!!
//表示马儿可以走 5(左2上1) 这个位置
if ((p1.x = curPoint.x - 2) >= 0 && (p1.y = curPoint.y - 1) >= 0) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 6(左1上2)这个位置
if ((p1.x = curPoint.x - 1) >= 0 && (p1.y = curPoint.y - 2) >= 0) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 7(右1上2)这个位置
if ((p1.x = curPoint.x + 1) < X && (p1.y = curPoint.y - 2) >= 0) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 0(右2上1)这个位置
if ((p1.x = curPoint.x + 2) < X && (p1.y = curPoint.y - 1) >= 0) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 1(右2下1)这个位置
if ((p1.x = curPoint.x + 2) < X && (p1.y = curPoint.y + 1) < Y) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 2(右1下2)这个位置
if ((p1.x = curPoint.x + 1) < X && (p1.y = curPoint.y + 2) < Y) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 3(左1下2)这个位置
if ((p1.x = curPoint.x - 1) >= 0 && (p1.y = curPoint.y + 2) < Y) {
ps.add(new Point(p1));
}
//判断马儿是否可以走 4(左2下1)这个位置
if ((p1.x = curPoint.x - 2) >= 0 && (p1.y = curPoint.y + 1) < Y) {
ps.add(new Point(p1));
}
return ps;
}
//根据当前这一步的所有的下一步的选择位置,进行非递减排序,减少回溯的次数
public static void sort(ArrayList<Point> ps){
ps.sort(new Comparator<Point>(){
@Override
public int compare(Point o1, Point o2) {
//先获取到o1的下一步的所有位置个数
int count1 = next(o1).size();
//先获取到o2的下一步的所有位置个数
int count2 = next(o2).size();
if (count1 < count2) {
return -1;
} else if (count1 == count2) {
return 0;
} else {
return 1;
}
}
});
}
}
Thank you for your time.
扫一扫
8729
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
