Bioinspired Bare Bones Mayfly Algorithm for Large-Scale Spherical Minimum Spanning Tree

Abstract

Mayfly algorithm (MA) is a bioinspired algorithm based on population proposed in recent years and has been applied to many engineering problems successfully. However, it has too many parameters, which makes it difficult to set and adjust a set of appropriate parameters for different problems. In order to avoid adjusting parameters, a bioinspired bare bones mayfly algorithm (BBMA) is proposed. The BBMA adopts Gaussian distribution and Lévy flight, which improves the convergence speed and accuracy of the algorithm and makes better exploration and exploitation of the search region. The minimum spanning tree (MST) problem is a classic combinatorial optimization problem. This study provides a mathematical model for solving a variant of the MST problem, in which all points and solutions are on a sphere. Finally, the BBMA is used to solve the large-scale spherical MST problems. By comparing and analyzing the results of BBMA and other swarm intelligence algorithms in sixteen scales, the experimental results illustrate that the proposed algorithm is superior to other algorithms for the MST problems on a sphere.

Introduction

Tree is a connected graph with simple structure which contains no loops and widely applied in graph theory (Diestel, 2000). The minimum spanning tree (MST) problem is a practical, well-known, and widely studied problem in the field of combinatorial optimization (Graham and Hell, 1985). This problem has a long history, which was first put forward by Borüvka in 1926. Many engineering problems are solved based on MST (Bo Jiang, 2009), such as communications network design (Hsinghua et al., 2001), the construction of urban roads, the shortest path (Beardwood et al., 1959), distribution network planning, and pavement crack detection. There are some classical algorithms for solving MST, such as the Prim algorithm (Bo Jiang, 2009) and Kruskal algorithm (Joseph, 1956). They all belong to greedy algorithms, and generally, only one minimum spanning tree can be obtained. However, in practical application, it is usually necessary to find a group of minimum or subminimum spanning trees as the basis for scheme evaluation or selection. Therefore, finding an effective algorithm to solve MST problems is still a frontier topic. In recent years, a large number of bioinspired algorithms have been proposed, such as the marine predator algorithm (Faramarzi et al., 2020), chimp optimization algorithm (Khishe and Mosavi, 2020), arithmetic optimization algorithm (Abualigah et al., 2021), bald eagle search algorithm (Alsattar et al., 2020), Harris hawks optimization algorithm (Heidari et al., 2019), squirrel search algorithm (Jain et al., 2018), pathfinder algorithm (Yapici and Cetinkaya, 2019), equilibrium optimizer (Faramarzi et al., 2019). The swarm intelligence algorithm has been widely used in various optimization problems and achieved good results, for example, path planning problems solved by the central force optimization algorithm (Chen et al., 2016), teaching–learning-based optimization algorithm (Majumder et al., 2021), water wave optimization algorithm (Yan et al., 2021), chicken swarm optimization algorithm (Liang et al., 2020), etc. Location problems are solved by the genetic algorithm (Li et al., 2021), particle swarm optimization (Yue et al., 2019), flower pollination algorithm (Singh and Mittal, 2021), etc. Also, the design of a reconfigurable antenna array is solved by the differential evolution algorithm (Li and Yin, 2011a), biogeography-based optimization (Li and Yin, 2011b), etc. In fact, the meta-heuristic algorithm can generate a set of minimum or subminimum spanning trees rather than one minimum spanning tree. The genetic algorithm (Zhou et al., 1996), artificial bee colony algorithm (Singh, 2009), ant colony optimization (Neumann and Witt, 2007), tabu search algorithm (Katagiri et al., 2012), and simulated annealing algorithm have been used for solving the MST problem.

For the MST problem, we usually calculate it in two-dimensional space, but it is of practical significance to study MST in three-dimensional space. For example, sockets are connected with wires in cuboid rooms, and roads on hills and mountains are planned. Also, as we all know, the surface of the Earth where we live is very close to a sphere. In many research fields, atoms, molecules, and proteins are represented as spheres, and foods in life, such as eggs, seeds, onions, and pumpkins, are close to spheres. Some buildings, glass, and plastics are made into spheres. Similar to the traveling salesman problem (TSP), it is also an NP-hard problem. Now scholars have applied the cuckoo search algorithm (Ouyang et al., 2013), glowworm swarm optimization (Chen et al., 2017), and flower pollination algorithm (Zhou et al., 2019) to solve the spherical TSP. Thus, it is of essence crucial to study the MST on a three-dimensional sphere. Bi and Zhou have applied the improved artificial electric field algorithm to the spherical MST problem (Bi et al., 2021). In this article, we will further study the cases of more nodes on the sphere.

The mayfly algorithm (MA) proposed by Konstantinos Zervoudakis and Stelios Tsafarakis (2020) is a population-based intelligent optimization bioinspired algorithm inspired by the flight and mating behavior of adult mayflies. Due to its high calculation accuracy and simple structure, researchers employed it to address problems of numerous disciplines. Guo and Kittisak Jermsittiparsert used improved MA to optimize the component size of high-temperature PEMFC-powered CCHP (Guo et al., 2021). Liu and Jiang proposed a multiobjective MA for a short-term wind speed forecasting system based on optimal sub-model selection (Liu et al., 2021a). Trinav Bhattacharyya and Bitanu Chatterjee combined MA with harmony search algorithm to solve the feature selection problem (Bhattacharyya et al., 2020). Liu and Chai used energy spectrum statistics and improved MA for bearing fault diagnosis (Liu et al., 2021b). Chen and Song proposed the balanced MA to optimize the configuration of electric vehicle charging stations on the distribution system (Chen et al., 2021). MohamedAbd and ElazizaS. Senthilraja used MA to predict the performance of a solar photovoltaic collector and electrolytic hydrogen production system (AbdElaziz et al., 2021). To obtain a group of more perfect minimum spanning trees or subminimum spanning trees on a sphere in finite time, a bare bones mayfly algorithm (BBMA) is proposed to solve spherical MST problems. By simplifying the algorithm parameters and using the statistical update method, the fast convergence and solution accuracy of the proposed algorithm are better than before, and it shows superior ability in solving large-scale problems.

The rest of this article is organized as follows: Related Work describes the related work and basic mayfly algorithm. The Proposed BBMA for Large-Scale Spherical MST introduces the proposed bare bones mayfly algorithm for spherical MST. Comparison and analysis of results evaluated by BBMA and other algorithms are given in Experimental Results and Discussion. This article is concluded in Conclusion and Future Work.

Related Work

Spherical Minimum Spanning Tree Mathematical Model

A semicircle takes its diameter as its axis of rotation, and the surface formed by rotation is called a sphere. The radius of the semicircle is the radius of the sphere. In this study, the coordinate origin (Figures 1A,B ) is set as the center of the sphere. The equation of a sphere with radius iswhere is the coordinate of each point on the sphere.

FIGURE 1

Representation of Points on a Sphere

The coordinate position on the sphere can be expressed by the following formula (Hearn and Pauline Baker, 2004):

Each coordinate is represented by a function of the surface parameters and . Usually, we normalize the three coordinate functions and make and in the range of 0–1. Eqs. 3–5 show a sphere with radius , and the center is at the coordinate origin (Eldem and Ülker, 2017).where parameters and determine a position by representing lines of constant longitude and lines of constant latitude, respectively. To simplify calculations, a sphere with is used in this study. When the parameters and take different values, the coordinate position on the sphere is as shown in Figure 1C (Uğur et al., 2009).

Geodesics Between Point Pairs on a Unit Sphere

The circle of a sphere cut by the plane passing through the center of the sphere is called a great circle (Wikipedia, 2012). On the sphere, the length of the shortest connecting line between two points is the length of an inferior arc between the two points of the great circle passing through the two points. We call this arc length the geodesic (Lomnitz, 1995).

The geodesic between two points and on a sphere is shown in Figure 1D. These two points can be represented by two vectors and . The scalar product of the two vectors iswhere is the angle between two vectors. The scalar product is calculated as

Also, the shortest distance formula is

From Eqs 6–8, we get

The distance from point to point is the same as the distance from to . If there are points on the sphere, an symmetric distance matrix will be obtained by calculating the distance between each two points. The matrix is as follows:where represents the length of the geodesic between point and point . In particular, means that point cannot reach itself.

Spherical Minimum Spanning Tree Mathematical Model

For the two-dimensional MST problem, represents an undirected graph, where is a finite set of nodes and is a finite set of edges. Each has its corresponding weight . is set as 0 or 1. If , is selected; if , is not selected. The variable is the number of nodes of the graph contained in the set . The mathematical model of the minimum spanning tree is as follows:where the constraint condition Eq. 12 ensures that the last generated graph is a spanning tree. Also, constraint condition Eq. 13 ensures that it is not a circle in the process of solving the minimum spanning tree problem.

As for a 3D spherical minimum spanning tree problem, a finite set of nodes are on a sphere. Each node is represented by . is a finite set of geodesics. Each geodesic has its corresponding weight which is calculated by Eq. 9. An symmetric distance matrix can be constructed as shown in Eq. 10. is set as 0 or 1. If , is selected; if , is not selected. The variable is the number of nodes on a sphere. The mathematical model of the spherical MST is as follows:

Similarly, constraint condition Eqs 12–14 are applied to Eq. 15.

Mayfly Algorithm

Mayfly algorithm is a new swarm intelligence bioinspired algorithm proposed in 2020. Its inspiration comes from the flying and mating behavior of male and female mayflies in nature. The algorithm can be considered as a modification of particle swarm optimization (PSO) (Kennedy and Eberhart, 1995), genetic algorithm (GA) (Goldberg, 1989), and firefly algorithm (FA) (Yang, 2009). At present, researchers have applied MA to many engineering problems.

Mayflies are insects that live in water when they are young. The feeding ability will be lost, and they only mate and reproduce when they grow up. In order to attract females, most adult male mayflies gather a few meters above the water to perform a nuptial dance. Then, female mayflies fly into these swarms to mate with male mayflies. After mating, the females lay their eggs on the water, and the mated mayflies will die.

In MA, the two idealized rules should be followed. First, after mayflies are born, they are regarded as adults. Second, the mated mayflies which have stronger ability to adapt to the environment can continue to survive. The algorithm works as follows. First, male and female populations are randomly generated. Each mayfly in the search space is regarded as a candidate solution represented by a -dimensional vector for male and for female. Its performance is evaluated according to the objective function shown in Eq. 15. The velocity of each mayfly is expressed by . The flying direction of each male mayfly is guided by its best location in history and the global optimal position in the population. Meanwhile, the female mayflies fly to the corresponding male mayflies. The main steps of mayfly algorithm are described as follows.

Movement of Male Mayflies

The gathering of male mayflies in a swarm is always a few meters above water for performing the nuptial dance. The position of a male mayfly is updated as follows:where is the position of mayfly at time and is the position at time and is the velocity of mayfly at time . The velocity is adjusted by its own velocity and individual and social experiences at time . However, the best male mayfly in the population is not affected by other mayflies, which helps the algorithm escape the local optimal. The velocity of a male mayfly is calculated aswhere represents the velocity of male mayfly at time in dimension , represents the position of dimension of mayfly at time , and represent positive attraction constants used to scale the contribution of the cognitive and social component, respectively, and is a fixed visibility coefficient used to limit a mayfly’s visibility to others. Furthermore, the best individual historical position of mayfly is represented by and is the global best position at time step , while is the Cartesian distance between and and is the Cartesian distance between and . These distances are calculated according to Eq. 18, 19. Finally, is the nuptial dance coefficient and is a random value.

Movement of Female Mayflies

The female mayflies move toward the males for breeding. The position of a female mayfly is updated as follows:where is the position of female mayfly at time and is the position at time step and represents the velocity of female mayfly at time . Its velocity is affected by its own velocity and the corresponding male mayfly’s position. It means that according to their fitness function, the best female should be attracted by the best male, the second best female by the second best male, and so on. However, the female mayfly which is better than the corresponding male mayfly is not affected by a male, it flies randomly. Consequently, considering minimization problems, their velocities are calculated aswhere is a velocity of female mayfly at time , is the position in dimension at time , is the positive attraction constant, and represents an unchanged visibility coefficient. represents the distance between and calculated according to Eq. 22. Finally, is the random fly coefficient, and is a random value.

Mating of Mayflies

The mating rules are the same as the way females are attracted by males. The best female breeds with the best male, the second best female with the second best male, and so on. The positions of two offspring are generated by the arithmetic weighted sum of the positions of parents as follows:where is the male mayfly’s position, is the female mayfly’s position, and is a random value. Offspring’s initial velocities are set to be zero, which helps the convergence of the algorithm.

After mating, the offspring are mixed with male and female parents. Then, the fitness values are sorted. The mayflies with low adaptability will die, and those with high adaptability will live for the next iteration. Algorithm 1 shows the pseudocode of MA.

Algorithm 1

The Proposed BBMA for Large-Scale Spherical MST

MST Based on Prüfer Coding

Through the abovementioned steps, we can get a Prüfer sequence of tree which is permutations of the numbers between 1 and .

Code Design

We assume that there are points on a sphere, and these points are represented by different integers between 1 and . The dimension of the position of each individual is , and the value in each dimension is a real number between 1 and .

Suppose an individual is represented by

The Prüfer sequence obtained by rounding is as follows:

According to , the spanning tree shown in Figure 2 is obtained.

FIGURE 2

The pseudo code of decoding the Prüfer sequence into a tree is shown in Algorithm 2.

Algorithm 2

The BBMA Algorithm

The basic MA has the problem of many initial parameters which have a great impact on the results. Besides, the accuracy of MA is not high enough due to lack of exploitation ability. Bare bones mayfly algorithm avoids the influence of parameters by cancelling the velocity (Ning and Wang, 2020; Song et al., 2020), and individual position is directly obtained by random sampling obeying Gaussian distribution like bare bones PSO (Kennedy, 2003). In order to enhance the exploitation ability and help the algorithm escape from the local optimal solution, BBMA uses Lévy flight to perform the nuptial dance of the optimal male and the random flight of the excellent female (Nezamivand et al., 2018). In addition, individuals crossing the border are pulled back into the search space instead of the method of placing cross-border individuals on the boundary so that it reduces the waste of search space (Wang et al., 2016). The main steps of BBMA are described as follows.

Movement of Male Mayflies

Male mayflies can be renewed in two ways as before. First, for individuals who are not the best, the Gaussian distribution based on the global optimal position and individual historical optimal position is used to calculate the position. In order to keep a balance between the diversity and convergence of algorithm, a disturbance which changes adaptively based on the diversity of the population and the convergence degree of the current individual is added (Zhang et al., 2014). The new update strategy is described as follows:where is the Gaussian distribution with mean and standard deviation , represents the global optimal individual, is the historical optimal solution of individual , is a random value, and and are two solutions selected from other male mayflies at random.

We know that the population is scattered in the early stage of evolution, so is large, and the Gaussian distribution is scattered, which is conducive for global search. In the later stage of evolution, the population is relatively concentrated, and individuals search carefully around . However, if the of a individual happens to be close or equal to in the evolution process, this individual will stop updating because the variance of Gaussian distribution becomes 0. Also, if most individuals among the swarm stop updating prematurely, the algorithm will converge to a false global optimum with high probability. Thus, assigning a disturbance on the variance of Gaussian distribution is a good way. As shown in Eq. 28, on the assumption that remains constant, the smaller the differential fitness value between and , the higher the disturbance . When the individual has the same fitness as , this individual will be affected by a disturbance with the maximal magnitude. In this case, this disturbance may prevent the algorithm from trapping into a local optimal solution. Furthermore, with the iteration of the algorithm, individuals get denser and denser. The smaller the value of , the smaller the and , which ensures the convergence of the algorithm.

As for the second individual update method, if the individual is the global optimal solution, Lévy flight is adopted. The small step size of Lévy flight improves the exploitation ability of the algorithm, and the less long step increases the ability of avoiding getting stuck in a local optimal value (Dinkar and Deep, 2018; Ren et al., 2021). By using Lévy flight, the overall performance of BBMA in solving large-scale problems has been greatly enhanced. In fact, Lévy flight is a random walk, which follows the Lévy distribution of the following formulas:where represents the step size and is an index by which the peak sharpness of the Lévy distribution can be adjusted. In this work, we set . and follow the Gaussian distribution, and stands for the gamma function which is obtained by Eq. 32. For the best individual, the update formula is as follows:

By using Lévy flight to search the solution space, the global exploration ability and local exploitation ability of the algorithm are better balanced.

Movement of Female Mayflies

Female mayflies can be renewed in two ways as before. Firstly, for individuals who are worse than their corresponding male mayflies, the Gaussian distribution based on the current female mayfly’s position and its corresponding male mayfly’s position is used to calculate the position. The new update strategy is described as follows:where is the Gaussian distribution, is the position of the female mayfly, and is the position of its corresponding male mayfly. The root sign makes the Gaussian distribution relatively concentrated so that it ensures female mayflies approach male mayflies faster, which accelerates the convergence.

As for the second individual update method, if the female mayfly is better than its corresponding male mayfly, the excellent female mayfly, like the best male mayfly, should use the strategy of Lévy flight which will make the algorithm get rid of the local optimum (Barshandeh and Haghzadeh, 2020). For excellent female mayflies, the update formula is as follows:

Both Gaussian distribution and Lévy distribution are statistical random distribution. The distribution of the former is regular, and the distribution of the latter is irregular. Their cooperation can prevent the lack of diversity of the algorithm and improve the convergence speed.

Mating of Mayflies

The mating process is the same as the basic MA as shown in Eq. 23. After mating, the offspring are mixed with parents. Then, the mayflies with low adaptability will die, and those with high adaptability will live for the next iteration.

Handling Cross-Border Mayflies

In the early stage of population evolution, the distance between the historical optimal position and the global optimal position of different individuals is far away, and the standard deviation of Gaussian distribution used for updating positions is relatively large, resulting in a greater opportunity for the new position to cross the boundary of the search space. In basic MA, the position of the cross-border individual is directly placed on the boundary, which will result in a waste of resources. In this study, according to the degree of individuals crossing the boundary, with the expectation of Gaussian distribution as the center, the cross-border individual is pulled back to the search space to obtain , and the cross-border individual is treated according to the following equation:where is the boundary, we assume is the upper bound and is the lower bound, if , and if . According to Eq. 37, when crosses the upper bound, it is pulled back to the interval . The less the crosses , the closer it is pulled back to ; the more the crosses , the closer it is pulled back to the center . When crosses the lower bound, it is pulled back to the interval . Similarly, the degree to which individuals are pulled back into the search space is proportional to the degree of individuals crossing the boundary.

The concrete implementation steps of the bare bones mayfly algorithm for spherical MST are as follows.

Algorithm 3

Experimental Results and Discussion

A large number of cases with different number of points are used to test the ability of BBMA in solving MST problems. All experiments are carried out on a sphere with , and the number of nodes the sphere is n = 25, 50, 75, 100, 150, 200, 250, 300, 350, 400, 500, 600, 700, 800, 900, and 1,000, The data and results of 400 nodes or less come from the literature (Bi et al., 2021) and the node data in higher dimensions are randomly generated. Due to the randomness of meta-heuristic algorithm, each case is run 30 times independently. The structure of this section is as follows: in Experimental Setup, the experimental setup is given; Comparison of Algorithms in Low-Dimensional Cases shows the comparison and analysis of experimental results between BBMA and other algorithms in the cases of low dimension; the comparison for medium-dimensional cases is shown in Comparison of Algorithms for Medium-Dimensional Cases; and the high-dimensional cases are shown in Comparison of Algorithms for High-Dimensional Cases.

Experimental Setup

Comparison of Algorithms in Low-Dimensional Cases

Cases with 25, 50, 75, and 100 points are used to compare the performance of algorithms mentioned above, and the results of 30 runs are obtained. Table 1 gives the best value, worst value, mean value, standard deviation, and the ranking of mean value. The bold data indicate that it is the best value of the twelve algorithms. Figure 3 shows the convergence curves in these four situations, Figure 4 shows the ANOVA test results, Figure 5 shows the fitness values for 30 runs, and Figure 12A–D show the minimum spanning tree for four low-dimensional cases, where “ROOT” is the root of the minimum spanning tree. Figure 13A shows the average running time of 30 runs of 12 algorithms in four dimensions. Finally, the Wilcoxon rank-sum non-parametric test results in low-dimensional cases are shown in Table 2.

FIGURE 3

FIGURE 4

FIGURE 5

The comparison results for 25 points are shown in Table 1. BBMA performs best in the best, worst, and mean value, but its standard deviation is 1.0202 which is worse than that of AMO, while the algorithm with the worst performance is GA. Figure 3A is the convergence curves of all algorithms; obviously, BBMA has the fastest convergence speed of twelve mentioned algorithms. As can be seen from Figure 4A, BBMA has the highest stability, and GWO is the worst. Figure 5A shows that, among the fitness values of 30 runs, BBMA is better than other algorithms in most cases, but AMO is three times better than BBMA, SMA is three times better than it, and MA is six times better than it. Figure 12A shows the minimum spanning tree for 25 points.

The comparison results of twelve algorithms at 50 points are shown in Table 1. It can be seen that the best value, worst value, and mean value of BBMA are the best, but the standard deviation ranks fifth, behind SOA, AMO, GWO, and DE. Figure 3B and Figure 4B show the convergence curve and analysis of variance results, respectively. By observing the convergence curve, we can clearly see that BBMA has the highest accuracy and the fastest convergence speed. Also, the result of variance analysis shows that BBMA is stable for solving this problem. The fitness values for 30 runs is shown in Figure 5B, and it can be seen that BBMA outperforms all other algorithms in 30 runs. Figure 12B shows the minimum spanning tree for 50 points.

The comparison of 75 points is shown in Table 1. BBMA is better than others in the best value, worst value, and mean value, but its standard deviation is 3.5118 which is worse than the standard deviation of PSO, GA, GOA, GWO, and SOA. Besides, GA has the worst accuracy. By observing Figure 3C, it can be noticed that both convergence accuracy and convergence speed are the highest for BBMA. As can be seen from Figure 4C, BBMA is still stable. The fitness values of 30 runs and the minimum spanning tree of 75 points can be seen in Figure 5C and Figure 12C.

The experience results for 100 points are shown in Table 1. The performance of BBMA is superior to that of others in the best, worst, and mean value. The best value of BBMA is 69.2455, and the best value of MA is 93.0942. BBMA is 25.83% better than the original algorithm. Figure 3D shows that BBMA has the highest accuracy and it still has excellent exploitation ability when other algorithms are stuck in a local optimal value. Figure 4D shows that BBMA has high stability. Figure 5D and Figure 12D show the fitness values for 30 runs and the minimum spanning tree path optimized for 100 points.

It can be seen from Figure 13A that BBMA runs the longest at each case, and GA and GWO are the two fastest algorithms. In addition to the running time, by comparing with other eleven algorithms, BBMA has the best performance in low dimensions. In addition, this study statistically tests the proposed algorithm. The Wilcoxon rank-sum non-parametric test results are shown in Table 2. BBMA is tested with others at the significance level. If p values in the table are all less than 0.05, it will prove that BBMA is obviously better than others. Statistically, the experimental results are significant.

Comparison of Algorithms for Medium-Dimensional Cases

In this section, BBMA and other eleven algorithms are tested in six medium-dimensional cases from 150 points, 200 points, 250 points, 300 points, 350 points, and 400 points. Table 3 records the best, worst, and mean value, standard deviation, and the ranking of the mean value. The bold data indicates that it is the best value of the twelve algorithms. Figure 6 shows the convergence curves of these six cases, and Figure 7 shows the ANOVA test results for each case. The fitness values for 30 runs are shown in Figure 8. Also, the minimum spanning tree for these cases is listed in Figure 12E–J, where “ROOT” is the root of the minimum spanning tree. In addition, Figure 13B shows the average running time of the twelve algorithms in different dimensions. Finally, the Wilcoxon rank-sum non-parametric test results in medium-dimensional cases are shown in Table 4.

FIGURE 6

FIGURE 7

FIGURE 8

Table 3 displays the experience results of 150 points and 200 points. Also, the statistical data shown in these tables reflect the great difference between different algorithms in searching ability. We can discover that except standard deviation, the best value, worst value, and mean value of BBMA are all the optimal. Also, the performance of GA is the worst. Figures 6A,B are the convergence curves for the two cases, and it can be seen that BBMA has a faster convergence speed and accuracy and strong exploration ability. Figures 7A,B show the analysis of variance results for the two cases, and we can see that the stability of BBMA is at a relatively high level. Figures 8A,B are the fitness values in 30 runs for 150 and 200 points. The MST for the two cases can be found in Figures 12E,F.

Table 3 also shows the comparison results of different algorithms at 250 points and 300 points. BBMA is the best in the best, worst, and mean value, and GA is the worst. However, as for the standard deviation, GWO is the best at 250 and 300 points. Figures 6C,D show the convergence curves in these two cases. The convergence speed and accuracy of BBMA are much superior to others. When other algorithms fall into local optimization, it still has good performance. The results of analysis of variance can be seen in Figures 7C,D, and BBMA has high stability. Figures 8C,D show the curves of the fitness values of 12 algorithms running independently for 30 times in these two cases. The search accuracy of BBMA is much better than that of the other 11 algorithms. Figures 12G,H show the MST at 250 points and 300 points, respectively.

The situation at 350 points and 400 points is shown in Table 3. BBMA performs best in the best, worst, and mean value. Compared with the best value of MA, the accuracy of BBMA is improved by 26.88% at 350 points and 24.85% at 400 points. As shown in Figures 6E,F, with the growth of dimension, the performance of BBMA is getting better and better. Most algorithms fall into local optimal solution at generation 100, but BBMA always has strong search ability. Figures 7E,F show the analysis of variance results in two cases, Figures 8E,F show the fitness values of 30 runs, and Figures 12I,F show the MST. It can be seen that BBMA has high stability and has better ability to solve spherical MST problems in medium-dimension cases than in low-dimension cases.

In addition, Figure 13B shows that the average running time of BBMA is the longest in the six cases. Compared with other algorithms, MA, DE, and AMO also run longer. Through the abovementioned analysis, we have noticed that BBMA has the outstanding performance in the medium-dimensional cases. The Wilcoxon rank-sum test results are shown in Table 4. Similarly, the p values in the table are all less than 0.05, which proves that BBMA algorithm is better than others in medium dimensions.

Comparison of Algorithms for High-Dimensional Cases

In Comparison of Algorithms in Low-Dimensional Cases and Comparison of Algorithms for Medium-Dimensional Cases, BBMA has been compared with other 11 algorithms in low and medium dimensions. BBMA shows very superior performance. Most of the problems encountered in real life are complex and high-dimensional problems, so in this section, BBMA and MA are tested in higher dimensions where n = 500, 600, 700, 800, 900, and 1,000 (see Table 5).

Table 5 shows the comparison results of BBMA and MA in six high-dimensional cases and also compares the best, worst, and mean value and standard deviation. The bold data indicate that it is the optimal result of the two. It can be seen that BBMA is superior to its original algorithm in the first three items in each case. The convergence curves of the two algorithms are shown in Figures 9A–F; obviously, the convergence speed and convergence accuracy of BBMA are better, and both exploration and exploitation capabilities have been greatly improved. The results of analysis of variance are shown in Figures 10A–F, and MA is more stable than BBMA. The fitness values of 30 independent operations are shown in Figures 11A–F. Also, Figures 12A–P show the MST of BBMA in different cases, where “ROOT” is the root of the MST, and the minimum spanning tree produced by BBMA is of high quality. Figure 13C is a histogram of the average running time of BBMA and MA, and we find that BBMA runs longer. Finally, Table 6 shows the results of the Wilcoxon rank-sum test in high-dimensional cases. The p values are so small that we can know that BBMA is significantly better than MA.

FIGURE 9

FIGURE 10

FIGURE 11

FIGURE 12

FIGURE 13

Conclusion and Future Work

Mayfly algorithm is a new population-based bioinspired algorithm, which has strong ability to solve continuous problems. It combines the advantages of PSO, FA, and GA and has superior exploration ability, high solution accuracy, and fast convergence. It improves the shortcoming that MA has many initial parameters and the parameters have a large impact on the results. Furthermore, Lévy flight is used for updating the position of the optimal male and excellent female to help the algorithm escape from local optimal solution. In addition, in order to make effective use of the search space, a cross-border punishment mechanism similar to “mirror wall” is used to deal with cross-border individuals. In order to demonstrate the effectiveness of BBMA, the MST problems are solved on a sphere. Compared with MA, AEFA, GA, PSO, ICA, SOA, GOA, GWO, SMA, DE, and AMO in 16 different cases, the test results show that BBMA has superior solving ability, and the higher the dimension is, the more obvious the superiority of BBMA will be. Therefore, BBMA is a good method for large-scale problems in real life. According to the NFL theorem, there is no algorithm that has superior performance for any problem. BBMA has some limitations in solving the spherical MST problems: its running time is relatively long, and its stability needs to be improved. In the future, BBMA will be further applied for solving the spherical MST problems in real life, such as removing noise on the femoral surface and directional location estimators (Kirschstein et al., 2013; Kirschstein et al., 2019). Also, it will be applied to other practical applications, such as logistics center location, path planning, weather forecast, and charging station address selection.

References

• 1

  Abd ElazizM.SenthilrajaS.ZayedM. E.ElsheikhA. H.MostafaR. R.LuS. (2021). A New Random Vector Functional Link Integrated with Mayfly Optimization Algorithm for Performance Prediction of Solar Photovoltaic thermal Collector Combined with Electrolytic Hydrogen Production System. Appl. Therm. Eng.193, 117055. 10.1016/j.applthermaleng.2021.117055

  • CrossRef
  • Google Scholar

• 2

  AbualigahL.DiabatA.MirjaliliS.Abd ElazizM. Elaziz.GandomiA. H. (2021). The Arithmetic Optimization Algorithm. Comp. Methods Appl. Mech. Eng.376, 113609. 10.1016/j.cma.2020.113609

  • CrossRef
  • Google Scholar

• 3

  AlsattarH. A.ZaidanA. A.ZaidanB. B. (2020). Novel Meta-Heuristic Bald eagle Search Optimisation Algorithm. Artif. Intell. Rev.53, 2237–2264. 10.1007/s10462-019-09732-5

  • CrossRef
  • Google Scholar

• 4

  AnitaA. Y.YadavA. (2019). AEFA: Artificial Electric Field Algorithm for Global Optimization. Swarm Evol. Comput.48, 93–108. 10.1016/j.swevo.2019.03.013

  • CrossRef
  • Google Scholar

• 5

  Atashpaz-GargariE.LucasC. (2008). “Imperialist Competitive Algorithm: An Algorithm for Optimization Inspired by Imperialistic Competition,” in IEEE Congress on Evolutionary Computation, Singapore, September 25-28, 2007, 4661–4667.

  • Google Scholar

• 6

  BarshandehS.HaghzadehM. (2020). A New Hybrid Chaotic Atom Search Optimization Based on Tree-Seed Algorithm and Levy Flight for Solving Optimization Problems. Eng. Comput.37, 3079–3122. 10.1007/s00366-020-00994-0

  • CrossRef
  • Google Scholar

• 7

  BeardwoodJ.HaltonJ. H.HammersleyJ. M. (1959). The Shortest Path through many Points. Math. Proc. Camb. Phil. Soc.55 (4), 299–327. 10.1017/s0305004100034095

  • CrossRef
  • Google Scholar

• 8

  BhattacharyyaT.ChatterjeeB.SinghP. K.YoonJ. H.GeemZ. W.SarkarR. (2020). Mayfly in harmony: a New Hybrid Meta-Heuristic Feature Selection Algorithm. IEEE Access8, 195929–195945. 10.1109/access.2020.3031718

  • CrossRef
  • Google Scholar

• 9

  BiJ.ZhouY.TangZ.LuoQ. (2021). Artificial Electric Field Algorithm with Inertia and Repulsion for Spherical Minimum Spanning Tree. Appl. Intelligence52, 195–214. 10.1007/s10489-021-02415-1

  • CrossRef
  • Google Scholar

• 10

  Bo JiangZ. (2009). Li. Research on Minimum Spanning Tree Based on Prim Algorithm. Comput. Eng. Des.30 (13), 3244–3247. 10.1115/MNHMT2009-18287

  • CrossRef
  • Google Scholar

• 11

  ChenL.XuC.SongH.JermsittiparsertK. (2021). Optimal Sizing and Sitting of Evcs in the Distribution System Using Metaheuristics: a Case Study. Energ. Rep.7 (1), 208–217. 10.1016/j.egyr.2020.12.032

  • CrossRef
  • Google Scholar

• 12

  ChenX.ZhouY.TangZ.LuoQ. (2017). A Hybrid Algorithm Combining Glowworm Swarm Optimization and Complete 2-opt Algorithm for Spherical Travelling Salesman Problems. Appl. Soft Comput.58, 104–114. 10.1016/j.asoc.2017.04.057

  • CrossRef
  • Google Scholar

• 13

  ChenY.YuJ.MeiY.WangY.SuX. (2016). Modified central Force Optimization (MCFO) Algorithm for 3D UAV Path Planning. Neurocomputing171, 878–888. 10.1016/j.neucom.2015.07.044

  • CrossRef
  • Google Scholar

• 14

  CrabbM. C. (2006). Counting Nilpotent Endomorphisms. Finite Fields Their Appl.12 (1), 151–154. 10.1016/j.ffa.2005.03.001

  • CrossRef
  • Google Scholar

• 15

  DerracJ.GarcíaS.MolinaD.HerreraF. (2011). A Practical Tutorial on the Use of Nonparametric Statistical Tests as a Methodology for Comparing Evolutionary and Swarm Intelligence Algorithms. Swarm Evol. Comput.1 (1), 3–18. 10.1016/j.swevo.2011.02.002

  • CrossRef
  • Google Scholar

• 16

  DhimanG.KumarV. (2019). Seagull Optimization Algorithm: Theory and its Applications for Large-Scale Industrial Engineering Problems. Knowledge-Based Syst.165, 169–196. 10.1016/j.knosys.2018.11.024

  • CrossRef
  • Google Scholar

• 17

  DiestelR. (2000). Graph Theory. Math. Gaz.173 (502), 67–128. 10.1007/978-3-662-53622-3

  • CrossRef
  • Google Scholar

• 18

  DinkarS. K.DeepK. (2018). An Efficient Opposition Based Lévy Flight Antlion Optimizer for Optimization Problems. J. Comput. Sci.29, 119–141. 10.1016/j.jocs.2018.10.002

  • CrossRef
  • Google Scholar

• 19

  EldemH.ÜlkerE. (2017). The Application of Ant colony Optimization in the Solution of 3D Traveling Salesman Problem on a Sphere. Eng. Sci. Technol. Int. J.20, 1242–1248. 10.1016/j.jestch.2017.08.005

  • CrossRef
  • Google Scholar

• 20

  FaramarziA.HeidarinejadM.StephensB.MirjaliliS. (2019). Equilibrium Optimizer: A Novel Optimization Algorithm. Knowledge-Based Syst.191, 105190. 10.1016/j.knosys.2019.105190

  • CrossRef
  • Google Scholar

• 21

  FaramarziA.HeidarinejadM.MirjaliliS.GandomiA. H. (2020). Marine Predators Algorithm: A Nature-Inspired Metaheuristic. Expert Syst. Appl.152, 113377. 10.1016/j.eswa.2020.113377

  • CrossRef
  • Google Scholar

• 22

  GibbonsJ. D.ChakrabortiS. (2011). Nonparametric Statistical Inference. Berlin: Springer, 977–979. 10.1007/978-3-642-04898-2_420

  • CrossRef
  • Google Scholar

• 23

  GoldbergD. E. (1989). Genetic Algorithm in Search, Optimization, and Machine Learning. Boston: Addison-Wesley Pub. Co.

  • Google Scholar

• 24

  GrahamR. L.HellP. (1985). On the History of the Minimum Spanning Tree Problem. IEEE Ann. Hist. Comput.7 (1), 43–57. 10.1109/mahc.1985.10011

  • CrossRef
  • Google Scholar

• 25

  GuoX.YanX.JermsittiparsertK. (2021). Using the Modified Mayfly Algorithm for Optimizing the Component Size and Operation Strategy of a High Temperature PEMFC-Powered CCHP. Energ. Rep.7 (11), 1234–1245. 10.1016/j.egyr.2021.02.042

  • CrossRef
  • Google Scholar

• 26

  HearnD. D.Pauline BakerM. (2004). Computer Graphics with Open GL. Beijing: Publishing House of Electronics Industry, 672.

  • Google Scholar

• 27

  HeidariA. A.MirjaliliS.FarisH.MafarjaI. M.ChenH. (2019). Harris Hawks Optimization: Algorithm and Applications. Future Generation Comp. Syst.97, 849–872. 10.1016/j.future.2019.02.028

  • CrossRef
  • Google Scholar

• 28

  HollandJ. H. (1992). Genetic Algorithms. Sci. Am.267 (1), 66–72. 10.1038/scientificamerican0792-66

  • CrossRef
  • Google Scholar

• 29

  HsinghuaC.PremkumarG.ChuC.-H. (2001). Genetic Algorithms for Communications Network Design - an Empirical Study of the Factors that Influence Performance. IEEE Trans. Evol. Computat.5 (3), 236–249. 10.1109/4235.930313

  • CrossRef
  • Google Scholar

• 30

  JainM.SinghV.RaniA. (2018). A Novel Nature-Inspired Algorithm for Optimization: Squirrel Search Algorithm. Swarm Evol. Comput.44, 148–175. 10.1016/j.swevo.2018.02.013

  • CrossRef
  • Google Scholar

• 31

  JosephK. B. (1956). On the Shortest Spanning Subtree of a Graph and the Traveling Salesman Problem. Am. Math. Soc.7 (1), 48–50.

  • Google Scholar

• 32

  KatagiriH.HayashidaT.NishizakiI.GuoQ. (2012). A Hybrid Algorithm Based on Tabu Search and Ant colony Optimization for K-Minimum Spanning Tree Problems. Expert Syst. Appl.39 (5), 5681–5686. 10.1016/j.eswa.2011.11.103

  • CrossRef
  • Google Scholar

• 33

  KennedyJ. (2003). “Bare Bones Particle Swarms,” in IEEE Swarm Intelligence Symposium, Indianapolis, IN, USA, April 26, 2003, 80–87.

  • Google Scholar

• 34

  KennedyJ.EberhartR. (1995). “Particle Swarm Optimization,” in Proc. of 1995 IEEE Int. Conf. Neural Networks, Perth, WA, Australia, November 27-December 01, 1995, 1942–1948.

  • Google Scholar

• 35

  KhisheM.MosaviM. R. (2020). Chimp Optimization Algorithm. Expert Syst. Appl.149, 113338. 10.1016/j.eswa.2020.113338

  • CrossRef
  • Google Scholar

• 36

  KirschsteinT.LiebscherS.BeckerC. (2013). Robust Estimation of Location and Scatter by Pruning the Minimum Spanning Tree. J. Multivariate Anal.120 (9), 173–184. 10.1016/j.jmva.2013.05.004

  • CrossRef
  • Google Scholar

• 37

  KirschsteinT.LiebscherS.PandolfoG.PorzioG. C.RagoziniG. (2019). On Finite-Sample Robustness of Directional Location Estimators. Comput. Stat. Data Anal.133, 53–75. 10.1016/j.csda.2018.08.028

  • CrossRef
  • Google Scholar

• 38

  LiJ.LiuZ.WangX. (2021). Public Charging Station Location Determination for Electric Ride-Hailing Vehicles Based on an Improved Genetic Algorithm. Sust. Cities Soc.74 (4), 103181. 10.1016/j.scs.2021.103181

  • CrossRef
  • Google Scholar

• 39

  LiS.ChenH.WangM.HeidariA. A.MirjaliliS. (2020). Slime Mould Algorithm: a New Method for Stochastic Optimization. Future Generation Comp. Syst.111, 300–323. 10.1016/j.future.2020.03.055

  • CrossRef
  • Google Scholar

• 40

  LiX.YinM. (2011). Hybrid Differential Evolution with Biogeography-Based Optimization for Design of a Reconfigurable Antenna Array with Discrete Phase Shifters. Int. J. Antennas Propagation2011 (1), 235–245. 10.1155/2011/685629

  • CrossRef
  • Google Scholar

• 41

  LiX.YinM. (2011). Design of a Reconfigurable Antenna Array with Discrete Phase Shifters Using Differential Evolution Algorithm. Pier B31, 29–43. 10.2528/pierb11032902

  • CrossRef
  • Google Scholar

• 42

  LiX.ZhangJ.YinM. (2014). Animal Migration Optimization: an Optimization Algorithm Inspired by Animal Migration Behavior. Neural Comput. Applic24, 1867–1877. 10.1007/s00521-013-1433-8

  • CrossRef
  • Google Scholar

• 43

  LiangX.KouD.WenL. (2020). An Improved Chicken Swarm Optimization Algorithm and its Application in Robot Path Planning. IEEE Access8, 49543–49550. 10.1109/access.2020.2974498

  • CrossRef
  • Google Scholar

• 44

  LiuY.ChaiY.LiuB.WangY. (2021). Bearing Fault Diagnosis Based on Energy Spectrum Statistics and Modified Mayfly Optimization Algorithm. Sensors21 (6), 2245. 10.3390/s21062245

  • CrossRef
  • Google Scholar

• 45

  LiuZ.JiangP.WangJ.ZhangL. (2021). Ensemble Forecasting System for Short-Term Wind Speed Forecasting Based on Optimal Sub-model Selection and Multi-Objective Version of Mayfly Optimization Algorithm. Expert Syst. Appl.177, 114974. 10.1016/j.eswa.2021.114974

  • CrossRef
  • Google Scholar

• 46

  LomnitzC. (1995). On the Distribution of Distances between Random Points on a Sphere. Bull. Seismological Soc. America85 (3), 951–953. 10.1785/bssa0850030951

  • CrossRef
  • Google Scholar

• 47

  MajumderA.MajumderA.BhaumikR. (2021). Teaching-Learning-Based Optimization Algorithm for Path Planning and Task Allocation in Multi-Robot Plant Inspection System. Arab J. Sci. Eng.46, 8999–9021. 10.1007/s13369-021-05710-8

  • CrossRef
  • Google Scholar

• 48

  MirjaliliS.MirjaliliS. M.LewisA. (2014). Grey Wolf Optimizer. Adv. Eng. Softw.69, 46–61. 10.1016/j.advengsoft.2013.12.007

  • CrossRef
  • Google Scholar

• 49

  NeumannF.WittC. (2007). Ant colony Optimization and the Minimum Spanning Tree Problem. Theor. Comp. Sci.411 (25), 2406–2413. 10.1016/j.tcs.2010.02.012

  • CrossRef
  • Google Scholar

• 50

  NezamivandS.AhmadB.NajafiF. (2018). PSOSCALF: A New Hybrid PSO Based on Sine Cosine Algorithm and Levy Flight for Solving Optimization Problems. Appl. Soft Comput.73, 697–726. 10.1016/j.asoc.2018.09.019

  • CrossRef
  • Google Scholar

• 51

  NingL.WangL. (2020). Bare-Bones Based Sine Cosine Algorithm for Global Optimization. J. Comput. Sci.47, 101219. 10.1016/j.jocs.2020.101219

  • CrossRef
  • Google Scholar

• 52

  OuyangX.ZhouY.LuoQ.ChenH. (2013). A Novel Discrete Cuckoo Search Algorithm for Spherical Traveling Salesman Problem. Appl. Math. Inf. Sci.7 (2), 777–784. 10.12785/amis/070248

  • CrossRef
  • Google Scholar

• 53

  RenH.LiJ.ChenH.LiC. (2021). Adaptive Levy-Assisted Salp Swarm Algorithm: Analysis and Optimization Case Studies. Mathematics Comput. Simulation181, 380–409. 10.1016/j.matcom.2020.09.027

  • CrossRef
  • Google Scholar

• 54

  SaremiS.MirjaliliS.LewisA. (2017). Grasshopper Optimisation Algorithm: Theory and Application. Adv. Eng. Softw.105, 30–47. 10.1016/j.advengsoft.2017.01.004

  • CrossRef
  • Google Scholar

• 55

  SinghA. (2009). An Artificial Bee colony Algorithm for the Leaf-Constrained Minimum Spanning Tree Problem. Appl. Soft Comput.9 (2), 625–631. 10.1016/j.asoc.2008.09.001

  • CrossRef
  • Google Scholar

• 56

  SinghP.MittalN. (2021). An Efficient Localization Approach to Locate Sensor Nodes in 3D Wireless Sensor Networks Using Adaptive Flower Pollination Algorithm. Wireless Netw.27 (3), 1999–2014. 10.1007/s11276-021-02557-7

  • CrossRef
  • Google Scholar

• 57

  SongX.-F.ZhangY.GongD.-W.SunX.-Y. (2020). Feature Selection Using Bare-Bones Particle Swarm Optimization with Mutual Information. Pattern Recognition112, 107804. 10.1016/j.patcog.2020.107804

  • CrossRef
  • Google Scholar

• 58

  StornR.PriceK. (1997). Differential Evolution – A Simple and Efficient Heuristic for Global Optimization over Continuous Spaces. J. Glob. Optimization11 (4), 341–359. 10.1023/a:1008202821328

  • CrossRef
  • Google Scholar

• 59

  UğurA.KorukoğluS.AliÇ.CinsdikiciM. (2009). Genetic Algorithm Based Solution for Tsp on a Sphere. Math. Comput. Appl.14 (3), 219–228. 10.3390/mca14030219

  • CrossRef
  • Google Scholar

• 60

  WangD.MengL.ZhaoW. (2016). Improved Bare Bones Particle Swarm Optimization with Adaptive Search center. Chin. J. Comput.39 (12), 2652–2667. 10.11897/SP.J.1016.2016.02652

  • CrossRef
  • Google Scholar

• 61

  Wikipedia (2012). Great circle. Available at: http://en.wikipedia.org/wiki/Great_circle.2012.

  • Google Scholar

• 62

  YanZ.ZhangJ.TangJ. (2021). Path Planning for Autonomous Underwater Vehicle Based on an Enhanced Water Wave Optimization Algorithm. Maths. Comput. Simulation181, 192–241. 10.1016/j.matcom.2020.09.019

  • CrossRef
  • Google Scholar

• 63

  YangX.-S. (2009). “Firefly Algorithms for Multimodal Optimization,” in Stochastic Algorithms: Foundations and Applications. SAGA 2009. Lecture Notes in Computer Science (Berlin: Springer), 5792, 169–178. 10.1007/978-3-642-04944-6_14

  • CrossRef
  • Google Scholar

• 64

  YapiciH.CetinkayaN. (2019). A New Meta-Heuristic Optimizer: Pathfinder Algorithm. Appl. Soft Comput.78, 545–568. 10.1016/j.asoc.2019.03.012

  • CrossRef
  • Google Scholar

• 65

  YueY.CaoLi.HuJ.CaiS.HangB.WuH. (2019). A Novel Hybrid Location Algorithm Based on Chaotic Particle Swarm Optimization for Mobile Position Estimation. IEEE Access7 (99), 58541–58552. 10.1109/access.2019.2914924

  • CrossRef
  • Google Scholar

• 66

  ZervoudakisK.TsafarakisS. (2020). A Mayfly Optimization Algorithm. Comput. Ind. Eng.145, 106559. 10.1016/j.cie.2020.106559

  • CrossRef
  • Google Scholar

• 67

  ZhangY.GongD.-w.SunX.-y.GengN. (2014). Adaptive Bare-Bones Particle Swarm Optimization Algorithm and its Convergence Analysis. Soft Comput.18 (7), 1337–1352. 10.1007/s00500-013-1147-y

  • CrossRef
  • Google Scholar

• 68

  ZhouG.GenM.WuT. (1996). “A New Approach to the Degree-Constrained Minimum Spanning Tree Problem Using Genetic Algorithm,” in 1996 IEEE International Conference on Systems, Man and Cybernetics. Information Intelligence and Systems (Cat. No.96CH35929), Beijing, China, October 14-17, 1996, 2683–2688.

  • Google Scholar

• 69

  ZhouY.WangR.ZhaoC.LuoQ.MetwallyM. A. (2019). Discrete Greedy Flower Pollination Algorithm for Spherical Traveling Salesman Problem. Neural Comput. Applic31 (7), 2155–2170. 10.1007/s00521-017-3176-4

  • CrossRef
  • Google Scholar