Solving fuzzy multi-objective shortest path problem based on data envelopment analysis approach

The shortest path problem (SPP) is a special network structured linear programming problem that appears in a wide range of applications. Classical SPPs consider only one objective in the networks while some or all of the multiple, conflicting and incommensurate objectives such as optimization of cost, profit, time, distance, risk, and quality of service may arise together in real-world applications. These types of SPPs are known as the multi-objective shortest path problem (MOSPP) and can be solved with the existing various approaches. This paper develops a Data Envelopment Analysis (DEA)-based approach to solve the MOSPP with fuzzy parameters (FMOSPP) to account for real situations where input–output data include uncertainty of triangular membership form. This approach to make a connection between the MOSPP and DEA is more flexible to deal with real practical applications. To this end, each arc in a FMOSPP is considered as a decision-making unit with multiple fuzzy inputs and outputs. Then two fuzzy efficiency scores are obtained corresponding to each arc. These fuzzy efficiency scores are combined to define a unique fuzzy relative efficiency. Hence, the FMOSPP is converted into a single objective Fuzzy Shortest Path Problem (FSPP) that can be solved using existing FSPP algorithms.


Introduction
Network optimization models can be seen in a wide range of application such as transportation systems, communication systems, pipeline distribution systems, fluid flow systems and neural decision systems [1]. In many real-life applications, it is required to obtain a best or shortest paths from one node to another. This kind of problem as one of the important types of network models is called shortest path problem (SPP). However, the best path cannot be mostly measured by a single objective alone, but by some or all of them together. In such cases, we will encounter a multi-objective shortest B Ali Ebrahimnejad aemarzoun@gmail.com; a.ebrahimnejad@qaemiau.ac.ir 1 path problem (MOSPP). In these problems, multiple, conflicting and incommensurate objectives such as optimization of cost, profit, time, distance, risk, quality of service etc. may arise. There are various approaches to solve MOSPPs. For example, Martins and Santos [2] analyzed the labeling algorithm for the MOSPP by the generalization of the classical shortest path problem. Sastry et al. [1] proposed three algorithms for MOSPP and an algorithm to detect negative cycle in a network. Using their proposed algorithms, MOSPP in a cyclic and acyclic network having weights either positive or negative or both can be solved. Sauvanet and Neron [3] focused on the determination of a well-balanced trade-off path between the objectives. They recalled the original Best Compromise A* method and then presented two improvements to speed up the search of the best compromise solution. Pulido et al. [4] introduced an exact label-setting algorithm that returns the subset of Pareto optimal paths that satisfy a set of lexicographic goals, or the subset that minimizes deviation from goals if these cannot be fully satisfied. Shi et al. [5] proposed an exact method for finding all the Paretooptimal paths for a multi-criteria constrained shortest path problem. Cintrano et al. [6] proposed a robust model for the bi-objective shortest path problem applied in a smart mobility context. Zajac and Huber [7] proposed a classification scheme in which every identified objective (in a given multi-objective routing problem) is sorted into a category and related to problem elements. Ajeil et al. [8] solved a path planning problem for an autonomous mobile robot in static and dynamic environments. Dehghani et al. [9] proposed a multi-objective optimization model for reliable communication flow networks, including maximizing the network reliability, minimizing total cost, and maximizing network flow, simultaneously.
Among the aforementioned methods, there are a number that are based on data envelopment analysis (DEA). DEA is a widely used technique to evaluate the performance of a set of homogeneous decision making units (DMUs). Nowadays, extensive researches have been performed in DEA and its scope has been broadened. "One reason is that DEA has opened up possibilities for use in cases which have been resistant to other approaches because of the complex (often unknown) nature of the relations between the multiple inputs and multiple outputs involved in many of these activities" [10]. The DEA-based approaches to make a connection between the MOSPP and efficiency are more flexible to deal with the real-world applications. Among these approaches, the work of Masoumi et al. [11] can be mentioned that considered each possible path from the source node to the destination one as a DMU. They proved that each nondominated path in the MOSPP is equivalent to the strongly efficient DMU under the Variable Returns to Scale (VRS) FDH model, and in this way, they tried to determine the set of all non-dominated paths on the network. Amirteimoori et al. [12] obtained a unique efficiency measure corresponding to each arc using DEA models. Then he defined the overall efficiency for each possible path and determined an efficient path, a path with maximum overall efficiency, as an optimal solution to the MOSPP.
In this paper, we develop a DEA-based approach to solve the MOSPP with fuzzy parameters. As you know, for various reasons such as shortage of information, insufficient data, lack of evidence, and so forth, the decision maker usually does not have exact information about the parameters in the real-life applications of MOSPP. In fact, in classical MOSPP it is assumed that the costs and time of paths are exactly known. In real-life shortest path cases, decision-makers may face many uncertainties on the cost and time of paths because of aesthetic quality of the scenery and changing weather, social, or economic conditions. In other words, the decisionmakers cannot exactly know the parameters of a MOSPP in reality. Additionally, these parameters are not stable since this imprecision may follow from the lack of exact information or data, uncertainty in judgement and high information cost. This imprecision embedded into the parameters may not be of a probabilistic type, so this type of uncertainty can be handled via fuzzy parameters. The MOSPP with at least one fuzzy parameter is called Fuzzy MOSPP (FMOSPP). There are few methods to solve the FMOSPPs. For example, Rani and Reddy [13] discussed the multi-objective fuzzy shortest path selection, where the arc lengths are expressed as trapezoidal fuzzy numbers. Zero et al. [14] proposed two algorithms to the Bi-objective Shortest Path (BSP) problem considering a fuzzy objective. Hasuike [15] proposed a model for a robust shortest path problem with random arc costs. The author formulated the proposed model as a bi-criteria optimization problem to minimize the total cost of the path from a source node to a destination node as well as to maximize the range of the confidence interval satisfying the condition that the total cost of any path in the region is less than a target value in terms of robustness. Patle et al. [16] presented the efficient obstacle avoidance mechanism in a complex environment based on the fuzzy optimized decision function.
Abbaszdeh Sori et al. [17] proposed an approach using fuzzy inference system to find the most desirable path of a fuzzy constrained SSP based on three factors of cost, time, and risk which the length of path does not exceed a predetermined value. Majumder et al. [18] presented an uncertain MOSPP for a weighted connected directed graph, where every arc weight is associated with uncertain cost and time. They have formulated the expected value model and chanceconstrained model of the proposed uncertain MOSPP and employed two multi-objective genetic algorithms for solving these models. Nonetheless, to the best of our knowledge, there is very limited literature available based on fuzzy DEA to solve the FMOSPPs.
Note that each arc in a FMOSPP corresponds to multiple fuzzy attributes. In this paper, using Fuzzy DEA (FDEA) models, each batch of multiple fuzzy attributes associated with the arcs is converted into a unique fuzzy attribute. To do this, each arc in a FMOSPP is considered as a DMU with multiple fuzzy inputs and outputs. Indeed, the objective functions that need to be maximized are used to define outputs and those ones that need to be minimized are used to define inputs. Then two fuzzy efficiency scores are obtained corresponding to each arc. These fuzzy efficiency scores are combined to define a unique fuzzy attribute. By considering the unique fuzzy relative efficiency as the unique attribute of the arcs, the FMOSPP under consideration is converted into a single objective fuzzy shortest path problem (FSPP) that can be solved using the existing algorithms for solving fuzzy FSPPs such as the generalized Dijkstra algorithm proposed by Deng et al. [19]. The artificial bee colony (ABC) algorithm proposed by Ebrahimnejad et al. [20] can also be used to solve the fuzzy SP (FSP) problems with fuzzy arc weights. The proposed algorithm by Tajdin et al. [21] can compute a shortest path in a network having various types of fuzzy arc lengths. Moreover, the fuzzy-based modified particle swarm optimization algorithm [22] and elite artificial bees' colony algorithm [23] can be used for solving FSP problem. An overview of efficient algorithms for solving SP problems with various types of input data in junction with fuzzy, intuitionistic, vague, interval fuzzy, interval-valued intuitionistic fuzzy and neutrosophic sets can be found in Broumi et al. [24].
It should be mentioned that there are various FDEA methods to measure the efficiency of DMUs with fuzzy inputs and outputs. For example, the tolerance approach developed by Sengupta [25] is one of the first fuzzy DEA models. Chen [26] modified the α-level approach and proposed an alternative fuzzy DEA to handle both the crisp and fuzzy data. Saati and Memariani [27] developed a fuzzy slack-based measure (SBM) based on the α-level approach. Hatami-Marbini et al. [28] used the provided ranking method in Asady and Zendehnam [29] to propose a fuzzy DEA model for assessing the efficiency scores in the fuzzy environment. Sheth and Triantis [30] introduced a fuzzy goal DEA framework to measure and evaluate the goals of efficiency and effectiveness in a fuzzy environment. Azar et al. [31] developed a new model for calculating the common set of weights (CSW) in fuzzy environment. Ebrahimi et al. [32] proposed a pair of DEA models to find the lower and upper bound efficiency scores in the presence of interval and weak ordinal data. Heydari et al. [33] presented a fully fuzzy network DEA-Range Adjusted Measure model for evaluating airlines and extends the network DEA-RAM model in the fully fuzzy framework. Arana-Jimenez et al. [34] proposed slacks-based additive inefficiency measure to deal with the problem of efficiency assessment when the input and output data are given as fuzzy sets. Another method that underpins this paper is the work of Saati et al. [35]. They presented a fuzzy version of CCR model [36] with the asymmetrical triangular fuzzy number and suggested a procedure for its solution. The basic idea in this procedure is to transform the fuzzy CCR model into a crisp linear programming problem by applying an alternative α-cut approach. In this way, the procedure converts the problem to an interval programming in which, instead of comparing the equality (or inequality) of two intervals, a variable is defined in the interval, not only satisfies the set of constraints but also maximizes the efficiency value.
On such motivation basis, the main contributions of this study are summarized as follows: (1) To the best of our knowledge, this study is the first attempt for solving FMOSPP using FDEA approach, (2) The use of the FDEA approach allows the simultaneous use of maximization and minimization functions in modeling the FMOSPP, (3) The proposed FDEA approach gives an efficient plan in converting FMO-SPP into a crisp single objective SPP, (4) In contrast to the existing multi objective optimization solution approaches such as goal programming, the proposed solution technique keeps the structure of SPP for ease of solution and implementation, (5) The FDEA approach proposed in this paper provided an efficient plan for the FMOSPP without explicit the knowledge of the decision maker's fuzzy utility function.
The rest of the paper is organized as follows: Section 2 presents the necessary mathematical background including some basic definitions about the fuzzy numbers along with the preliminaries about the FDEA. A brief introduction to the FMOSPP is presented in Sect. 3. In Sect. 4, the concept of the fuzzy efficient path is introduced and a FDEA based approach is proposed to solve the FMOSPP. The proposed approach is illustrated with two numerical examples in Sect. 5. In Sect. 6, the proposed approach is compared with the extended version of goal programming (GP) approach. Finally, conclusions are presented in Sect. 7.

Some basic definitions on fuzzy numbers
Fuzzy sets theory was first introduced by Zadeh [42] to deal with the uncertainty of human judgment (Wang et al. [43]). Each fuzzy set is an extension of a crisp set that allows partial membership between 0 and 1. A fuzzy set on the real line that satisfies some conditions is called a fuzzy number. Here, some basic definitions and basic fuzzy arithmetic operations on the fuzzy numbers are reviewed [44].

Definition 2.1 A fuzzy set A in R (real line) is defined to be a set of ordered pairs
is called the membership function for the fuzzy set.

Definition 2.3 A fuzzy set A on
R is called convex if for any x, y ∈ R and any λ ∈ [0, 1] the following inequality is satisfied: (1)

Definition 2.4
A fuzzy setÃ on the real line that satisfies the conditions of normality and convexity is called a fuzzy number.

Definition 2.5 A fuzzy numberÃ, denoted bỹ
, is called a triangular fuzzy number if its membership function is given by relation (2): Definition 2.6 A triangular fuzzy numberÃ = (a l , a m , a u ) is said to be a non-negative (positive) triangular fuzzy number if and only if a l ≥ 0 (a l > 0).

Definition 2.7 Two triangular fuzzy numbers
be two positive triangular fuzzy numbers and t is a non-negative scaler. Then the basic fuzzy arithmetic operations along with the scaler multiplication are defined as follows: Scaler multiplication: tÃ = (ta l , ta m , ta u ). (3)

Definition 2.9
The α-cutof fuzzy setÃ on R, denoted by [Ã] α is a crisp set that contains all the elements of R whose membership degree inÃ are greater than or equal to the value of α, that is,

Theorem 2.10
The α-cut of the triangular fuzzy numberÃ = (a l , a m , a u ) is the closed interval as relation (5) [37]:

Mathematical model of FDEA
The DEA is a useful non-parametric method that allows handling multiple outputs and inputs at the same time. Mathematically, DEA is a linear programming-based technique for evaluating the relative efficiency of a set of DMUs with multiple inputs and outputs. A major stated advantage of the DEA is that it does not require any prior assumptions on underlying functional relationships between inputs and outputs [45].
In classic DEA models, input and output data are assumed to be crisp values. However, in many real-world applications, the available data are imprecise or vague. This issue may result from unquantifiable, incomplete, missed information, etc. In recent years, much effort has been made on the theoretical development of Fuzzy DEA (FDEA) to deal with this impreciseness and ambiguity in DEA.
Suppose we wish to evaluate the efficiencies of n homogeneous DMUs. Each DMU j ( j = 1, . . . , n) produces s different fuzzy outputsỹ j = (ỹ 1 j , . . . ,ỹ s j ), using m different fuzzy inputsx j = (x 1 j , . . . ,x m j ). The multiplier form of the input-oriented CCR model [36] for evaluating the relative efficiency of DMU p is given by model (6): where '˜' indicates the fuzziness and u r s (r = 1, . . . , s) along with the v i s (i = 1, . . . , m) are the weights assigned to the outputs and inputs, respectively. In this paper, to transform the Fuzzy Multi-Objective Shortest Path Problem (FMOSPP) into a single objective Fuzzy Shortest Path Problem (FSPP), some FDEA models will need to be solved. There exist several methods to solve the FDEA model (6) (e.g., [27,28,30,31]). Nonetheless, we use the proposed approach by Saati et al. [35] which is a fuzzy version of the CCR model [36] with the asymmetrical triangular fuzzy numbers. The basis of this approach is to transform each FDEA model into an interval model. In fact, using this approach, the FDEA model (6) is converted into a possibilistic programming problem and then transformed into a crisp linear programming model based on α-cut. In this way, by solving the linear programming problem for a given α-cut, it will be possible to generate a reliable and robust solution for the possibilistic mathematical programming problem.

Fuzzy multi-objective shortest path problem (FMOSPP)
The aim of the MOSPP is to find the shortest path between a pair of source and destination node in a network with respect to the multiple predefined objectives. In other words, there are two or more objective functions that must be optimized when selecting a path. When at least one of the parameters of a MOSPP is represented in terms of fuzzy numbers, it is called FMOSPP. The fuzzy numbers can appear in different parameters simultaneously and also be of different types such as normal, triangular or trapezoidal fuzzy numbers along with crisp numbers. Throughout the paper, without loss of generality, we consider all the parameters of the FMOSPP as non-negative triangular fuzzy numbers.
For more details, consider a directed graph G = (V , A), where the vertices or nodes of the graph are contained in the Corresponding to each arc (i, j), there are p associated non-negative fuzzy attributesc k i j (k = 1, . . . , p). Thus, the FMOSPP has p fuzzy objectives that need to be minimized and maximized. The problem is to find a best feasible path from the source node 1 to the destination node n such that each node in this path must be considered just one time. Mathematical model of this FMOSPP can be written as model (7): where x i j = 1 or x i j = 0 indicates that the arc (i, j) is in the path or not, respectively. There are a few methods to solve the FMOSPP (7) (e.g. [13][14][15]). In Section 4, we propose a DEA-based approach to solve the FMOSPP (7).

Fuzzy efficient path
Consider the FMOSPP (7) with p fuzzy objectives corresponding to the undirected network G = (V , A) that needs to be minimized and maximized. To find an optimal solution to this problem, each arc of the network is considered as a DMU. The arc attributes corresponding to the fuzzy objectives that should be minimized are considered as fuzzy inputs and the attributes corresponding to the fuzzy objectives that should be maximized are considered as fuzzy outputs of the DMUs. The aim of the proposed approach is to convert the p fuzzy attributes associated with each arc into a unique fuzzy attribute. In this way, the given FMOSPP (7) is transformed into a single objective fuzzy SPP. Then an efficient path with the maximum fuzzy efficiency is determined as an optimal solution to this problem. To this end, based on the FDEA, two fuzzy efficiency scores are obtained corresponding to each arc as a criteria of its relative performance. After that, the mean of the fuzzy efficiency scores is used to derive a new fuzzy attribute for each arc.
For a detailed description of the proposed method, let is also a non-negative triangular fuzzy number for all r = 1, . . . , s. In other words, each arc (i, j) ∈ A has k fuzzy inputs and s fuzzy outputs such that p = k + s. Figure 1 shows the structure of each arc as a DMU. Now, to solve the FMOSPP (7), the following steps should be taken: Step 1. For each arc (i, j) ∈ A, consider all the existing arc (i, h) ∈ A. Then with the node i as a target and changing the h, subject to (i, h) ∈ A, obtain the fuzzy relative efficiency of the arc (i, j) ∈ A, i.e.Ẽ 1 i j , as the optimal objective value of the FDEA model (8): By solving model (8) for all i and j, obtain the fuzzy relative efficiency scoresẼ 1 i j s related to all of the existing arcs in A. How to solve model (8) will be discussed later.
Step 2. For each arc (i, j) ∈ A, consider all the existing arc (h, j) ∈ A. Then, with the node j as a target and changing the h, subject to (h, j) ∈ A, obtain the fuzzy relative efficiency of the arc (i, j) ∈ A, denoted byẼ 2 i j , as the optimal objective value of the FDEA model (9): Step 3. By having the efficiency scoresẼ 1 i j andẼ 2 i j , obtain the unique fuzzy efficiency scoreẼ * i j for the arc (i, j) from the relation (10) as follows: Step 4. Replace the existing multiple fuzzy attributes To solve the FDEA model (8) and (9), we follow the Saati et al.'s [35] approach. For this purpose, consider model (8) in its detailed form as model (12): where 1 l = 1 m = 1 u is considered equal to 1 and 0 l = 0 m = 0 u = 0, to ensure that the efficiency values are not greater than 1. Model (12) is a possibilistic linear programming. By introducing α-cuts of the objective function and constraints, this model can be written as model (13): Since all of the coefficients and RHSs are stated as intervals, model (13) involves treating an interval problem and, hence, cannot be treated by standard methods without further transformations. There are several methods to solve this model that most of them compare the intervals in both sides of the constraints with each other. However, according to Saati et al. [35], instead of comparing the intervals, the variables in the intervals are defined such that they satisfy the set of constraints, and, at the same time the objective function is maximized. To this end, consider the variablesx t ih andŷ r ih according to the relation: By substituting the variablesx t ih andŷ r ih , model (13) is transformed to model (15) as follows: y r ,l i j + (y r ,m i j − y r ,l i j )α ≤ŷ r i j ≤ y r ,u i j − (y r ,u i j − y r ,m i j )α, ∀h; (i, h) ∈ A, r = 1, . . . , s, Model (15) is a non-linear programming problem. To transform this model to a linear programming model, the variable substitutions can be used as the relation: According to relation (16), an optimization point of model (13) can be obtained from the parametric linear programming problem (17), while α ∈ (0, 1] is its parameter: Similarly, model (9) can be solved following the Saati et al.'s [35] approach. Therefore, by solving model (17) for all (i, j) ∈ A, the fuzzy efficiency score E * ,α i j can be obtained from the relation (18), corresponding to each particular α ∈ (0, 1]: Finally, an efficient path corresponding to the desired 'α' can be determined by solving model (19):   (17) at α = α 1 and α = α 2 , respectively. . . . , s, t = 1, . . . , k) be an optimal solution of model (17) at α = α 2 . Therefore, relations (20) to (22) hold: On the other hand, for α 1 ≤ α 2 , relations (23) and (24) clearly hold: Thus, relations (23) and (24) lead to the relation (25): Therefore, from relations (20), (22) and (25), it is concluded that u * r , v * t ,ȳ r * ih andx t * ih (∀h; (i, h) ∈ A, r = 1, . . . , s, t = 1, . . . , k) is a feasible solution for model (17) at α = α 1 ; then, it is clear that Note that Theorem 4.1 is also true for E 2,α 1 i j and E 2,α 2 i j . Thus, from Theorem 4.1 and relation (18), it is immediately concluded that for any α 1 , α 2 ∈ (0, 1] such that α 1 ≤ α 2 , the inequality E * ,α 2 i j ≤ E * ,α 1 i j holds.

Note 4.2 The objection function
In this way, model (19) is transformed into a classic single objective SPP with a minimization objective function that can be solved using any standard algorithm for solving the SPP such as Dijkstra's algorithm [48].

Note 4.3
Using the mentioned α-cut approach, instead of a unique fuzzy efficient path, the efficient path table corresponding to different α ∈ (0, 1] will be provided for the decision-maker.

Numerical examples
In this section, the proposed fuzzy DEA approach for solving fuzzy multi-objective shortest path problem is illustrated using two numerical examples.  by (a l , a m , a u ), are given in Table 1. In this case, the following fuzzy three-objective shortest path problem is formulated to achieve the fuzzy efficient path: For solving model (26), the fuzzy DEA approach proposed in this study is used. Thus, according to this approach, flyingcost and flying-delays are considered as inputs and flyingsafety is considered as output. Based on this approach, to achieve a fuzzy efficient path, the values of the efficiency scores E 1,α i j and E 2,α i j must first be obtained. To this end, the corresponding model (8) and (9) which is equivalent to model (28) according to the data of Table 1: s.t.x 1 13 +x 2 13 = 1, y 1 12 −x 1 12 −x 2 12 ≤ 0, y 1 13 −x 1 13 −x 2 13 ≤ 0, Similarly, for a given α ∈ (0, 1], the value of E 2,α 13 can be obtained as the optimum objection value of model (29): By obtaining the values of the E 1,α i j s and E 2,α i j s for all i and j, the values of the E * ,α i j s can be computed according to Relation (18). Tables 2, 3 and 4 represent the obtained results for the E 1,α i j s, E 2,α i j s and E * ,α i j s corresponding to each α ∈ {0.1, 0.2, 0.3, . . . , 0.9, 1}, respectively.

Comparison of results with fuzzy goal programming
In this section, the performance of the proposed fuzzy DEA approach for solving FMOSPP (7) is compared with the extended version of goal programming (GP) approach.
GP is a popular approach for solving the multi-objective optimization problem by reducing it to a single objective one. The idea of GP is to minimize the distance between objective functions and an aspiration level vector either determined by the decision-maker or equalsf = f * Assume thatñ t = n l t , n m t , n u t andp t = p l t , p m t , p u t are the under deviations (achievement) and over deviations of the objectivesf t = (i, j)∈Ac t i j x i j from their aspiration values . . , k + s) be those objective functions of FMOSPP (7) that should be minimized and maximized, respectively. Accordingly, the FMOSPP (7) is converted by GP into a minimization problem of the deviational variables which minimizes the sum of deviational variables:     Table 6 The values of the E * ,α i j s for different levels of α ∈ (0, 1] in Model (32) is a linear programming problem that can be solved using the simplex method. It is worthwhile mentioning that both fuzzy DEA approach proposed in this study and the extended fuzzy GP approach reduces the FMOSPP (7) into a single objective problem. However, there are several important reasons for using our proposed method compared with the GP technique: -The classical LP problem (32) applied for solving FMO-SPP (7) is not a shortest path structured problem, whereas problem (19) for solving FMOSPP (7) are classical shortest path problem. -The classical LP problem (32) applied for solving FMO-SPP (7) has 9(k +s)+m constraints (without considering the non-negative constraints) and 3(k + s) + n variables (without considering slack variables), whereas problem (19) has only m constraints and n variables. Hence, from a computation point of view the proposed fuzzy DEA approach is preferable to the fuzzy GP technique for solving the FMOSPP (7). -By use of the fuzzy DEA approach proposed in this study, the FMOSPP (7) is converted into a single objective fuzzy SPP keeping the structure of the SP problem, whereas using GP approach increases the number of constraints and variables of the problem by adding new constraints to the shortest path problem. Therefore, utilizing the proposed technique to is highly economical in comparison with the goal programming technique from a computational viewpoint, considering the number of constraints and variables. -In the case that the problem requires an integer optimal solution, the fuzzy DEA approach can obtain an optimal solution to the SPP by ignoring the integrality restrictions, while fuzzy GP approach is not able to find integer solutions without adding the integrality restrictions.
Anyway, by solving model (30) Fig. 3. Therefore, it can be said, the proposed method is at least preferable considering the safety of the shortest path. Furthermore, although both optimal solutions obtained using the GP technique and the proposed method are Pareto optimal solutions, the efficiency scores of shortest paths obtained from the proposed method are more than the one derived from the GP technique. Finally, model (32) is an integer linear programming problem while model (19) is a linear programming problem. Solving an integer linear programming problem is harder than a linear programming problem.
Similarly, by solving model (30) corresponding to Example 5.2, aspiration level vector is obtained asf * = f * 1 ,f * 2 ,f * respectively. Again, although both optimal solutions obtained using the fuzzy GP technique and the proposed method are Pareto optimal solutions, the efficiency score of shortest path obtained from the proposed method is more than the one derived from the GP technique.

Conclusions
In this paper, a DEA-based approach was proposed to solve the SPP with multiple fuzzy attributes according to each arc called "FMOSPP". The idea behind the proposed approach was to convert the multiple fuzzy attributes into a unique fuzzy attribute. To this end, each arc in the given FMOSPP was considered as a DMU. Then those objective functions that need to be maximized were used to define the outputs of DMU and those ones that need to be minimized were used to define the inputs of DMU. In the following, using two fuzzy efficiency scores, the unique fuzzy relative efficiency was defined to each arc. By considering the unique efficiency as the only attribute, the FMOSPP was converted into a single objective FSPP that can be solved to determine a path with maximum efficiency as an efficient path. It is worth mentioning that all of the multiple fuzzy attributes were considered as non-negative triangular fuzzy numbers. Accordingly, we used the proposed method by Saati et al. [35] to obtain the mentioned fuzzy efficiency scores based on the α-cut approach. With this method, instead of facing a FSPP, a classic SPP with crisp parameters would be obtained for any selected α ∈ (0, 1] that can be solved using the standard algorithms for solving the SPP such as Dijkstra's algorithm [48]. One of the main advantages of the proposed fuzzy DEA approach relies on its computational simplicity, allowing for its implementation within complex network structures. In addition, a simple technique to solve a fuzzy multi objective shortest path problem would be to optimize the fuzzy utility function of decision maker. However, the mathematical representation of the fuzzy utility function may not be determined easily in many problems. But, the fuzzy DEA approach proposed in this paper provided an efficient plan for the FMOSPP without explicit the knowledge of the decision maker's fuzzy utility function. In sum, the proposed fuzzy DEA approach applied in this study not only gives the efficient path without any knowledge of decision-maker, but also keeps the structure of the shortest path problem.
At the end, it should be noted that this study links the fuzzy MOSPP and fuzzy DEA together. As a future research it would be worthwhile to solve the stochastic multi-objective capacitated transportation problem [49] via stochastic DEA approaches. As another research, someone may be interested in establishing links between MOSPP and DEA in the case that the arc weights are represented in terms of intuitionistic fuzzy numbers [50][51][52] and single-valued neutrosophic num-bers [53]. Moreover, the development of the proposed method for solving multi-objective all-pairs shortest path problem in an interval-valued fuzzy network [54] and an interval-valued Pythagorean fuzzy network [55] is left to the next study. Finally, the generalization of the proposed method for determining the efficient solution of multi-objective constrained shortest pat problem using new forms of fuzzy inference systems such as type-1 fuzzy logic systems [56,57], interval type-2 fuzzy logic systems [58] and generalized type-2 fuzzy logic system [57] is an interesting topic for future research.