Abstract
Evolutionary algorithms have been widely used to solve large and complex optimisation problems. Cultural algorithms (CAs) are evolutionary algorithms that have been used to solve both single and, to a less extent, multiobjective optimisation problems. In order to solve these optimisation problems, CAs make use of different strategies such as normative knowledge, historical knowledge, circumstantial knowledge, and among others. In this paper we present a comparison among CAs that make use of different evolutionary strategies; the first one implements a historical knowledge, the second one considers a circumstantial knowledge, and the third one implements a normative knowledge. These CAs are applied on a biobjective uncapacitated facility location problem (BOUFLP), the biobjective version of the wellknown uncapacitated facility location problem. To the best of our knowledge, only few articles have applied evolutionary multiobjective algorithms on the BOUFLP and none of those has focused on the impact of the evolutionary strategy on the algorithm performance. Our biobjective cultural algorithm, called BOCA, obtains important improvements when compared to other wellknown evolutionary biobjective optimisation algorithms such as PAES and NSGAII. The conflicting objective functions considered in this study are cost minimisation and coverage maximisation. Solutions obtained by each algorithm are compared using a hypervolume S metric.
1. Introduction
Evolutionary algorithms (EAs) are an effective alternative to (approximately) solve several large and complex optimisation problems, as they are able to find good solutions for a wide range of problems in acceptable computational time. Although less studied, EAs for multiobjective optimisation (MO) problems, called evolutionary multiobjective optimisation (EMO), have demonstrated to be very effective. In fact, during the last two decades, several authors have focused their efforts on the development of several EMO algorithms to solve a wide range of MO problems. For instance, Maravall and de Lope [1] use a genetic algorithm (GA) to solve the multiobjective dynamic optimisation for automatic parking system. In [2] the authors propose an improvement to the wellknown NSGA algorithm (called NSGAII) based on an elitist approach. In [3] the author presents an EMO algorithm applied on a specific variation of the wellstudied capacitated vehicle routing problem (CVRP), where the author includes in the EMO algorithm an explicit collective memory method, namely, the extended virtual loser (EVL). Other wellknown EMO algorithms developed during the last two decades are PAES [4] and MO particle swarm optimisation [5]. More recently, hybrid techniques have been also applied to a large number of optimisation problems (see [6]). A comprehensive literature review related to EMO algorithms can be found in [7].
EMO algorithms have some problems that must be taken into account, though. For instance, they tend to fall into premature convergence with low evolution efficiency [8]. This is because of implicit information embodied in the evolution process and domain knowledge corresponding to optimisation problems which are not fully included in the solution approach [9]. To overcome these problems, one can make use of implicit evolution information. Reynolds [10] proposes an EA, called cultural algorithm (CA), which is inspired from human culture evolution process and that makes use of the implicit evolution information generated at each iteration. The CAs have a dual evolution structure which consists of two spaces: population and belief space. On the one hand, the population space works as in any other EA, that is, using evolutionary operators such as mutation and crossover. On the other hand, in the belief space, implicit knowledge is extracted from selected individuals in the population and stored in a different way. Then, they are used to guide the whole evolution process in the population space such that they can induce the population to escape from local optimal solutions. It has been proved that CAs can effectively improve the evolution performance [9]. Although less studied, CAs have been also used to solve MO problems. Coello et al. [11] and Coello et al. [12], two remarkable surveys on CAs, only mention Coello and Landa [13] work as example of a CA application solving MO problems. More recently, Zhang et al. [14] present a CA which is enhanced by using a particle swarm optimisation algorithm. This enhanced CA is applied to fuel distribution MO problem. Srinivasan and Ramakrishnan [15] present a MO cultural algorithm that is applied on data mining domain. In [16] the authors applied a CA to a biobjective portfolio selection problem using normative knowledge in the belief space. In [17] authors present a formal framework to implement MO cultural algorithms. To the best of our knowledge, [18] is the only article that uses CAs to solve the biobjective uncapacitated facility location problem (BOUFLP). Furthermore, we did not find any article which compares the performance of CAs using different evolutionary strategies at the belief space level. Thus, in this paper we present an extension of the biobjective cultural algorithm (BOCA) developed in [18]. We use two different strategies at the belief space level and compare the performance of our new algorithms with the performance of the previous one. We also compare its results with other wellknown EMO algorithms such as PAES and NSGAII. Obtained solutions were compared using the hypervolume S metric proposed in [19].
The remaining of this paper is organised as follows. Section 2 shows an overview on MO focused on EMO algorithms. Section 2.1 presents briefly the BOUFLP and some of its distinctive features. In Section 3, we describe our implementation for the BOCA algorithm. We describe the main differences between our implementation and the one in [18]. Section 3.2 presents the BOCA algorithm applied to a set of wellknown instances from the literature. Finally, Section 4 presents the conclusions of this work.
2. Overview
In this section we show an overview of topics related to this paper. In Section 2.1, MO concepts are presented, emphasizing EMO algorithms and its state of art. More specifically, we focus on the development of EMO algorithms for MO combinatorial optimisation (MOCO) problems. In Section 2.2 we present the BOUFLP formulation based on a costcoverage approach. Finally, in Section 2.3 we present an overview of CAs and its multiobjective extension. Details of our CA implementation are also presented at the end of this section.
2.1. (Evolutionary) Multiobjective Optimisation
In this section we briefly introduce the main principles of MO problems and, particularly, MOCO problems. For a comprehensive review of this topic see [20, 21]. In this paper we will make use of the following notation for the comparison of vectors (solutions). Let and . We say that if . Similarly, we will say that if but . Finally, we say that if . A solution , with being equal to the number of decision variables, is called an efficient solution and its image , with , a nondominated point of the MO problem if there is no , with , such that . In [22] the author describes several excellence relations. These relations establish strict partial orders in the set of all nondominated points related to different aspects of their quality. Previously, in [23, 24] the authors consider several outperformance relations to address the closeness of the set of nondominated points found by an algorithm to the actual set of nondominated points, called Pareto Frontier (PF). In [25] a comprehensive explanation of the desirable features of an approximation to the PF is presented. In this paper, we choose the metric which is properly explained in [24]. The metric calculates the hypervolume of a multidimensional region [19] and allows the integration of aspects that are individually measured by other metrics. An advantage of the metric is that each algorithm can be assessed independently of the other algorithms involved in the study. However, the values of two sets and , namely, and , respectively, cannot be used to derive whether either set entirely dominates the other. Figure 1(a) shows a situation where and set completely dominates set . Figure 1(b) shows a situation where but neither dominates nor dominates .
(a) Nondominated set completely dominates nondominated set 
(b) Neither dominates nor dominates 
In this paper we use EAs to solve the BOUFLP. An EA is a stochastic search procedure inspired by the evolution process in nature. In this process, individuals evolve and the fitter ones have a better chance of reproduction and survival. The reproduction mechanisms favour the characteristics of the stronger parents and hopefully produce better children guaranteeing the presence of those characteristics in future generations [26]. EAs have been successfully applied to a large number of both single and multiobjective optimisation problems. Comprehensive reviews of EMO algorithms are presented in [11, 24] and more recently in [12]. A more general review of hybrid heuristics solving MOCO problems, where EMO algorithms are also included, is presented in [27].
2.2. Biobjective Uncapacitated Facility Location Problem
Facility location problem (FLP) is one of the most important problems for companies with the aim of distributing products to their customers. The problem consists of selecting sites to install plants, warehouses, and distribution centres, allocating customers to serving facilities and interconnecting facilities by flow assignment decisions. Comprehensive reviews and analysis of the FLP are presented in [28–30].
In this paper we consider a twolevel supply chain, where a single plant serves a set of warehouses, which serve a set of end customers or retailers. Figure 2 shows this configuration.
Two main (conflicting) objectives can be identified in the FLP:(i)minimise the total cost associated with the facility installation and customer allocation and(ii)maximise the customers rate coverage.
Several works on singleobjective optimisation have been carried out considering these two objectives separately. On the one hand, uncapacitated FLP (UFLP) is one of most studied FLPs in the literature. In the UFLP the main goal is to minimise the locationallocation cost of the network. On the other hand, median FLPs are one of the most common FLPs among those that are focused on coverage maximisation. Most important FLP models are well described and formalised in [29]. MO FLPs have been also well studied in the literature during the last two decades. A survey on this topic can be found in [31].
As we mentioned before, in this paper we solve the BOUFLP. BOUFLP has been modelled with minisum and maxisum objectives (cost and coverage). The following model formulation is based on [32]. Let be the set of potential facilities and the set of customers. Let be the fixed cost of opening facility and the demand of customer . Let be the cost of assigning the customer to facility and the distance between facility and customer . Let be the maximal covering distance; that is, customers within this distance to an open facility are considered well served. Let be the set of facilities that could serve customer within the maximal covering distance . Let be 1 if facility is open and 0, otherwise. Let be 1 if the whole demand of customer is served by facility and 0, otherwise. Consider with and . Objective functions and are as follows:
Equation (6) represents total operating cost; the first term corresponds to location cost, that is, the sum of the fixed costs of all the open facilities, and the second term represents the allocation cost, that is, the cost of attending customer demand by an open facility. Equation (7) measures coverage as the sum of the demand of customers attended by open facilities within the maximal covering distance. Equations (2) and (3) ensure that each customer is attended by only one facility. Equation (3) also forces customer to be assigned to an open facility. Finally, equations (4) and (5) set decision variables as binary.
2.3. Biobjective Cultural Algorithm
The experience and beliefs accepted by a community in a social system are the main motivations for the creation of the CAs. Originally proposed by Reynolds [10], CAs model the evolution of cultural systems based on the principles of human social evolution. In this case, evolution is seen as an optimisation process [10]. The CAs guide the evolution of the population based on the knowledge. Knowledge acquired during previous iterations is provided to future generations, allowing accelerating the convergence of the algorithm to good solutions [33]. Domain knowledge is modelled separately from the population, because there is certain independence between them which allows us to work and model them separately, in order to enhance the overall algorithm performance. Figure 3 shows this interaction.
CAs are mainly characterised by presenting two inheritance systems: one at population level, called population space, and the other at knowledge level, called belief space. This key feature is designed to increase the learning rates and convergence of the algorithm and thus to do a more responsive system for a number of problems [34]. Moreover, it allows us to identify two significant levels of knowledge: a microevolutionary level (represented by the population space) and macroevolutionary level (represented by the space of beliefs) [35].
CAs have the following components: population space (set of individuals who have independent features) [35]; belief space (stored individuals acquired in previous generations) [34]; computer protocol, connecting the two spaces and defining the rules on the type of information to be exchanged between them by using the acceptance and influence functions; and finally knowledge sources which are described in terms of their ability to coordinate the distribution of individuals depending on the nature of a problem instance [35]. These knowledge sources can be of the following types: circumstantial, normative, domain, topographic, and historical.
The most distinctive feature of CAs is the use of the belief space which through an influence function affects future generations. For this reason, in this paper we focus on the effect on the algorithm performance of changes in such an influence function. To do this, we have considered results obtained previously in [18], where the authors used an influence function based on historical knowledge, and we compare those results with our BOCA implementation which considers two influence functions: the first one based on circumstantial knowledge and the second one based on normative knowledge. Algorithm 1 shows the general procedure of our BOCA algorithm.

To initialise the population, we use a semirandom function. In its first phase, this function defines in a stochastic way the set of facilities that will be opened (selected facilities). Then, we allocate each customer to a selected facility minimising the cost function while avoiding minimising the coverage function . This strategy provides better results than using completely random initial populations, and its computational time additional cost is marginal.
To obtain the next generation, two parents are used in a recombination process. To avoid local optimal values, we do not overuse the culture. Thus, a parent is selected from the population to obtain diversity and the other parent is selected from the belief space to influence the next generation. The belief space keeps a list of all the individuals which meet some criteria. These criteria depend on what knowledge the algorithm implements. In this paper the circumstantial knowledge selects the best individuals found so far for each objective function. Thus, one individual will give us information on the best value found for and the other will do the same for . The historical knowledge stores a list of individuals with the best fitness value found so far. The fitness value is calculated as the hypervolume that is covered by an individual. Finally, normative knowledge considers a list of individuals which are pairwise nondominated with respect to the other individuals of their generation.
Let be the number of available facilities and let be the number of customers of our BOUFLP. In this paper decisions variables and are represented by a binary length vector and matrix, respectively. Comparing two different solutions (individuals) needs an evaluation criterion. In this paper we use the same criterion explained in [18].
3. Computational Experiments
In this section we present the set of instances that are used in this study as well as results obtained by our BOCA implementation.
3.1. Instances Presentation
The instances that were used in this paper correspond to random instances using a problem generator that follows the methodology from UflLib [36]. Previous works in the literature have also used this problem generator to create their test instances [18, 26].
The BOCA algorithm has several parameters that need to be set. As in [18], the number of generations considered in this paper is equal to 100. Population size is set equal to 100, mutation probability in the population space is equal to , and probability of mutation in the belief space is . Both and values are different from the values used in [18]. These values are chosen as they all together yield to the best performance of the algorithm given some test instances. Thus, although resulting values are different from that used in [18], the method we use to set them is the same as that used in that work. This is important in order to fairly compare the different BOCA algorithms.
3.2. Results and Discussion
In this section we compare the results obtained by the previous BOCA algorithm (see [18] for further details) and our approach. Moreover, a comparison between results obtained by wellknown EMO algorithms such as NSGAII and PAES and our BOCA algorithm is also presented in this section.
Tables 1 and 2 show the results obtained by the BOCA implementations using historical [18], circumstantial, and normative knowledge, respectively. In the same way, Tables 3 and 4 present the results obtained by the wellknown NSGAII and PAES algorithms. For each algorithm value (%), time (in seconds), and the number of efficient solutions have been included in these tables. As we mentioned before, we want to produce a set with a large number of efficient solutions , a value close to 100% (ideal), and a small . For the sake of easy reading, we have split the set of instances into two subsets (instances type and ).
We then compare our BOCA implementations with the one presented in [18]. Tables 5 and 6 show a comparison between those algorithms. As we can see, when compared in terms of its value (the bigger, the better), BOCA algorithm using historical knowledge () performs consistently better than the ones using circumstantial () and normative () knowledge. In fact obtains a value that is, in average, 5.8% bigger than the one obtained by and 6.5% bigger than the value obtained by . When compared in terms of the CPU time needed to reach the number of iterations (generations), is, in average, faster than both and algorithms. We can note that for instances, times required by and are, in average, quite similar (only 1.6% of difference). Finally, when we look at the number of efficient solutions found by each algorithm ( column), we can see that, again, outperforms both and algorithms. In this case, the average number of efficient solutions found by the algorithm is about 20% bigger than the one obtained by the other two approaches.
Results above are consistent with the good performance obtained by the approach in [18]. Moreover, results show that performance of the BOCA algorithm depends largely on the selected knowledge and it can make the difference in terms of value, time, and number of efficient solutions found by the algorithm. This is an important finding as it points out the relevance of the choice of a specific type of knowledge.
We now compare and algorithms to the wellknown NSGAII and PAES algorithms. Tables 7 and 8 show a comparison between our algorithm and the NSGAII and PAES algorithms. As we can see, although obtains, in average, a value 6.8% lower than the one obtained by the NSGAII algorithm, it is more than three times faster. Moreover, when is compared to PAES algorithm, the obtained values are, in average, equivalent while is around 30% faster than PAES. PAES obtains, in average, more efficient points than though (9.32%).
Finally, Tables 9 and 10 show a comparison between and NSGAII and PAES algorithms. performs quite similar to PAES algorithm with respect to both value and the number of obtained efficient solutions. However, is faster than PAES. Similar situation occurs when is compared to NSGAII algorithm. Although NSGAII obtains better values for both and , is much faster than NSGAII. This situation can be explained by the very fast performance that our algorithm obtains for the set of small instances. When we look further at the results, we can note that if we only consider both medium and large size instances, execution times obtained by both algorithms are quite similar to each other. This result confirms what is outlined in [18] in the sense of the good performance that the BOCA algorithm shows. Furthermore, our results confirm this good performance with respect to other wellknown EMO algorithms does not depend on which type of knowledge is considered. However, as we mentioned before, the choice of the knowledge used on the BOCA algorithm is an important issue and it has an impact on the algorithm performance.
4. Conclusions and Future Work
Evolutionary algorithms are a very good alternative to solve complex combinatorial optimisation problems. In this paper we have implemented a biobjective cultural algorithm to solve the wellknown BOUFLP. We have considered two different sources of knowledge, namely, circumstantial and normative, and compare them with a previously implemented historical knowledge. Furthermore, we compare our BOCA approaches with two wellknown EAs, namely, NSGAII and PAES.
Although BOCA approaches using both normative and circumstantial knowledge could not improve the results obtained by the BOCA algorithm with the historical knowledge, results pointed out that performance of the BOCA algorithm depends largely on the selected knowledge and it can make the difference in terms of value, time, and number of efficient solutions found by the algorithm. This is an important finding as it points out the relevance of the choice of a specific type of knowledge. Moreover, our results also confirm the good performance showed by the BOCA algorithm with respect to other wellknown EMO algorithms such as NSGAII and PAES algorithms. The BOCA algorithm is very competitive when compared to those EMO algorithms independently of the type of knowledge implemented.
As a future work we think that more investigation is needed in order to find patterns that allow us to get the right knowledge implemented depending on the problem features. As we mentioned before, the knowledge choice has an impact on the performance of the BOCA algorithm and therefore it must be studied in depth. Also, as future work, hybrid knowledge could be implemented in order to exploit the advantages of each kind of knowledge at the same time. Moreover, our BOCA algorithm can be used to solve other interesting MOPs arising in the logistic field, such as routing or scheduling problems.
Appendix
Result Tables
In this appendix section, obtained results are presented. Columns show the value obtained by algorithm as %. Algorithms are indexed as follows. The original BOCA algorithm is indexed by . BOCA algorithms using circumstantial and normative knowledge are indexed by and , respectively. Finally, the other EAs considered in this paper, namely, NSGAII and PAES, are indexed by and , respectively. Columns show the time obtained by each algorithm in seconds. Columns show the number of efficient solutions found by the corresponding algorithm. Finally operator shows a value that is equivalent to .
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.