程序员常用十种算法（7）— 克鲁斯卡尔算法（公交站问题）

程序员常用十种算法（7）— 克鲁斯卡尔算法（公交站问题）

1.克鲁斯卡尔算法介绍
1）克鲁斯卡尔(Kruskal)算法，是用来求加权连通图的最小生成树的算法。
2）基本思想：按照权值从小到大的顺序选择 n-1 条边,并保证这 n-1 条边不构成回路。
3）具体做法：首先构造一个只含 n 个顶点的森林，然后依权值从小到大从连通网中选择边加入到森林中，并使森林中不产生回路，直至森林变成一棵树为止。

*连通网的最小生成树： 在含有 n 个顶点的连通图中选择 n-1 条边，构成一棵极小连通子图，并使该连通子图中 n-1 条边上权值之和达到 最小，则称其为连通网的最小生成树。

2.公交站问题描述
3.克鲁斯卡尔算法图解
假设用数组 R 保存最小生成树结果


这个过程中存在两个不太好解决的问题：
（1）问题一：对图的所有边按照权值大小进行排序。
（2）问题二：将边添加到最小生成树中时，怎样判断是都形成了回路。
问题一：采用冒泡排序的方法解决（其他排序方法也行）。
问题二：记录顶点在“最小生成树中”的终点，顶点的终点是"在最小生成树中与它连通的最大顶点"。 然后每次需要将一条边添加到最小生存树时，判断该边的两个顶点的终点是否重合，重合的话则会构成回路。

举例说明判断是否构成回路
关于终点的说明：
1）就是将所有顶点按照从小到大的顺序排列好之后；某个顶点的终点就是"与它连通的最大顶点"。
2）因此，接下来，虽然<C,E>是权值最小的边，但是 C 和 E 的终点都是 F，即它们的终点相同，因此，将<C,E> 加入最小生成树的话会形成回路。这就是判断回路的方式，也就是说，我们加入的边的两个顶点不能都指向同一 个终点，否则将构成回路。

代码（结合代码理解更容易）：

import java.util.Arrays;

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[][] = {
                      /*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();
        kruskalCase.kruskal();
    }

    //构造器
    public KruskalCase(char[] vertexs, int[][] matrix) {
        //初始化顶点数和边的个数
        int vlen = vertexs.length;

        //初始化顶点，复制拷贝的方式
        this.vertexs = new char[vlen];
        for (int i = 0; i < vertexs.length; i++) {
            this.vertexs[i] = vertexs[i];
        }

        //初始化边，使用的是复制拷贝的方式
        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){
                    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); //12

        //按照边的权值大小进行排序（从小到大）
        sortEdges(edges);

        //遍历edges数组，将边添加到最小生成树中时，判断是准备加入的边是否形成了回路，如果没有就加入rets，否则不能加入
        for (int i = 0; i < edgeNum; i++) {
            //获取到第i条边的第一个顶点（起点）
            int p1 = getPosition(edges[i].start); //p1=4
            //获取到第i条边的第二个顶点
            int p2 = getPosition(edges[i].end); //p2=5

            //获取p1这个顶点在已有最小生成树中的终点
            int m = getEnd(ends, p1); //m=4
            //获取p2这个顶点在已有最小生成树中的终点
            int n = getEnd(ends, p2); //n=5
            //是否构成回路
            if (m != n){ //没有构成回路
                ends[m] = n; // 设置 m 在“已有最小生成树”中的终点 <E,F>[0,0,0,0,5,0,0,0,0,0,0,0]
                rets[index++] = edges[i]; //有一条边加入到rets数组
            }
        }
        //<E,F><C,D><D,E><B,F><E,G><A,B>
        //统计并打印最小生成树，输出 rets
        System.out.println("最小生成树为：");
        for (int i = 0; i < index; i++) {
            System.out.println(rets[i]);
        }
    }

    //打印邻接矩阵
    public void print(){
        System.out.println("邻接矩阵为：\n");
        for (int i = 0; i < vertexs.length; i++) {
            for (int j = 0; j < vertexs.length; j++) {
                System.out.printf("%12d", matrix[i][j]);
            }
            System.out.println(); //换行
        }
    }

    /**
     * 功能：对边进行排序处理，冒泡排序
     * @param edges 边的集合
     */
    private void sortEdges(EData[] edges){
        for (int i = 0; i < edges.length-1; i++) {
            for (int j = 0; j < edges.length-1-i; 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
     */
    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]);
                }
            }
        }
        return edges;
    }

    /**
     * 功能：获取下标为i的顶点的终点（），用于后面判断两个顶点的终点是否相同
     * @param ends：数组就是记录了各个顶点对应的终点是哪个，ends数组是在遍历过程中逐步形成的
     * @param i：表示传入的顶点对应的下标
     * @return 返回的就是下标为i的这个顶点对应的终点的下标
     */
    private int getEnd(int[] ends, int i){  //i=4  [0,0,0,0,5,0,0,0,0,0,0,0]
        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 +"]";
    }
}