Optimal Inconsistency Repairing of Pairwise Comparison Matrices Using Integrated Linear Programming and Eigenvector Methods

Satisfying consistency requirements of pairwise comparison matrix (PCM) is a critical step in decision making methodologies. An algorithm has been proposed to find a new modified consistent PCM in which it can replace the original inconsistent PCM in analytic hierarchy process (AHP) or in fuzzy AHP. This paper defines the modified consistent PCM by the original inconsistent PCM and an adjustable consistent PCM combined. The algorithm adopts a segment tree to gradually approach the greatest lower bound of the distance with the original PCM to obtain the middle value of an adjustable PCM. It also proposes a theorem to obtain the lower value and the upper value of an adjustable PCM based on two constraints. The experiments for crisp elements show that the proposed approach can preserve more of the original information than previous works of the same consistent value. The convergence rate of our algorithm is significantly faster than previous works with respect to different parameters. The experiments for fuzzy elements show that our method could obtain suitable modified fuzzy PCMs.


Introduction
Analytic Hierarchy Process (AHP) is developed by Saaty [1], which is a multicriterion decision-making methodology widely used in many real problems [2,3].The AHP methodology expresses the relative importance of criteria by pairwise comparisons and converts the values of pairwise comparisons to priorities.Fuzzy AHP methodology [4] is an advanced AHP methodology, which is used to tackle the uncertainty and inaccurate problems in multicriteria decision-making process.Fuzzy AHP derives the fuzzy priorities of criteria from pairwise comparisons matrix with triangular (or trapezoidal) fuzzy elements.To make sure the priorities of each criterion are accurate and sensible, consistency of pairwise comparison matrix (PCM) with crisp or fuzzy elements must be achieved.
Several works [5][6][7][8] focus on reducing inconsistent PCMs with crisp numbers.Karapetrovic and Rosenbloom [5] revised the single entry of a ratio's value till the consistency of relative matrix was at an acceptable level.Xu and Wei [6] preserved the initial ratios' value in the pairwise comparison matrix while obtaining satisfactory consistency requirements.Cao et al. [7] developed a heuristic approach, which can preserve more of the original information compared to Xu and Wei [6].However, for these three works, when the consistency requirement is increased, the computing times will be largely increased and the information on the original matrix cannot be well preserved.Anholcer's [8] work is used to minimize the distance between inconsistent PCMs and their corresponding consistent PCM.The aim of this work is to find out a new modified PCM which is consistent and has the closest resemblance to the optimal one.But the modified new matrix has a long distance with the original matrix according to the analysis of the parameters given in [6,7].
Some works used to solve the inconsistency of PCMs with fuzzy elements are given in [9][10][11][12].Xu and Wang [9] repaired incomplete and inconsistent fuzzy preference relations by finding out the unusual and false element until the consistency ratio was at a satisfactory level.Leung and Cao [10] proposed a new definition of fuzzy positive reciprocal matrix by setting deviation tolerances based on an idea of allowing inconsistent information.Morteza and Bafandeh [11] further discussed Leung and Cao's work and proposed a new method of fuzzy consistency tests by direct fuzzification of a QR (quick response) algorithm, which is one of the methods for the eigenvalues calculation of an arbitrary matrix.Wang and Chen [12] applied fuzzy linguistic preference relations to construct consistent PCMs by considering reducing the number of pairwise comparisons.However, these works do not have standard parameters to verify the reliability of their theory so far.Therefore, it is very important to prove the feasibility of a methodology that can reduce the inconsistency of the original matrix and preserve the original matrix's information as much as possible.
In this research work, we propose a modified consistent PCM as a combination of original inconsistent PCM and an adjustable consistent PCM.In order to achieve the modified PCM, this paper is structured in the following way: Section 2 and Section 3 give the basic concepts of the PCM with fuzzy and crisp elements, consistency indices for crisp and fuzzy elements, and parameters to judge the effectiveness of the modified matrix.Section 4 and Section 5 propose an algorithm to obtain the middle value, upper value, and lower value of the adjustable matrix.The main idea of our algorithm is to find the optimum priority vector by solving a linear programming problem and use the eigenvector method to obtain the adjustable matrix based on the optimum priority vector.In Section 6, the algorithm was applied to obtain the adjustable PCM by adopting the same illustrating PCM in [6,7]; the comparison results show that our algorithm can preserve more original information than Cao et al. [7] and Xu and Wei [6].In Section 7, the algorithm was used to find the sensible and the closest consistent modified PCM with fuzzy elements.The experiments show that the new modified matrix can satisfy consistent indices' requirements of NI [13] and CCI [14,15] and also can preserve more original information under parameters of  fuzzy and  fuzzy .Section 8 concludes with effective and efficient analysis to show that our algorithm can be performed easily.

State of the Art
where for each element ã = ( Where the operator ⊙ is one of the operation rules of triangular fuzzy elements, this operation can be calculated by the following equation: When positive reciprocal matrix Ã is crisp numbers , then the consistent condition is   ⋅   =   , {∀, ,  | 1 ≤ , ,  ≤ }. However, this definition is too strict, because it is unrealistic to reach perfect consistency of a PCM (crisp or fuzzy elements).Some works [10,13,[17][18][19][20][21][22][23] have developed consistency indices to accept a certain level of deviations.We will adopt consistency indices in Section 2.2 to determine whether the current PCM matrix is at an acceptable consistency level.

Consistency Indices.
Several consistency indices have been proposed for crisp numbers: for instance, geometric consistency index [17], singular value decomposition method [18], and harmonic consistency index [19].However, it was proven that all these consistency indices are linear or nonlinear transformations of Satty's CR [20].Therefore, we use Satty's CR in this paper for measuring consistency of crisp numbers.
The consistency index of CR [20] is defined as follows: where  max is the principle eigenvalue of .RI is random index, which can be gotten by searching a defined table.When the value of 0 < CR < 0.1, the consistency can be accepted.Several important works focus on the consistency of pairwise comparison matrix with fuzzy elements.The first one is Leung and Cao [10], who proposed a notion with consideration of a tolerance deviation.However, the notion is strongly related to Satty's CR and it has shortcomings to calculate consistency of pairwise comparison matrix with fuzzy elements [11].The second one is Ramík and Korviny's [13] work, which proposes a new consistency index NI to examine fuzzy elements based on the distance of the matrix to a special ratio matrix and compare the properties with CR.This work has been further studied and has been used by several important works [21][22][23].We have tested this work's performance.The results indicated that it can satisfy reasonable results with fuzzy elements, although it has some shortcomings [23].
The consistency index of NI [13] is defined as follows: where    is a normalization factor and the values of    ,    , and    can be obtained from    ,    , and    .Another successful index that can be extracted from [14,15] is that they extend GCI (geometric consistency index) [24] to CCI (centric consistency index) to deal with PCM with triangle fuzzy elements.The consistency index CCI is defined as follows: where Ã = (   ,    ,    ) are the elements for fuzzy PCM; ω = (ω  ) is a priority vector derived by logarithmic least squares.When CCI( Ã) = 0, Ã is considered fully consistent.Thresholds remain identical with index GCI as CCI being a fuzzy extension of GCI.The thresholds are provided as follows: CCI = 0.3147 for  = 3, CCI = 0.3526 for  = 4, and CCI = 0.370 for  > 4 based on Aguarón and Moreno-Jiménez [24].
In conclusion, we will adopt CR to examine the consistency of crisp elements and use NI and CCI to examine the consistency of fuzzy elements.

Parameters to Judge the Effectiveness of Modified Matrix.
Next, we are going to describe the necessary parameters which can be used to measure the effectiveness of modified PCM .Xu and Wei [6] have given two parameters in the following: and  are used as the parameters of modificatory effectiveness.The authors of [6,7] argue that a modified matrix that preserves the most information of the original one must satisfy the following condition: 0 <  < 2 and 0 <  < 1. Extending these two parameters to be suitable to judge modificatory effectiveness of fuzzy elements so the range is identical with  and ,  fuzzy and  fuzzy are Besides  and , two parameters proposed by Xu and Wei [6], we propose a third parameter should be added, which is Condition of Order Preservation (COP) (Table 3) [25].For example, suppose the original matrix  has alternatives ( 1 ,  2 ,  3 , and  4 ); it has that the relationship  1 dominates  2 and  3 dominates  4 , and the judgments indicate that the extent to which  1 dominates  2 is greater than the extent to which  3 dominates  4 ; then the priority vector  should satisfy ( 1 ) > ( 2 ) and ( 3 ) > ( 4 ) (preservation of order of preference) and ( 1 )/( 2 ) > ( 3 )/( 4 ) (preservation of order of intensity of preference).

Distance Analysis between Original Matrix and Modified
Matrix.To measure the distance from original PCM  (or Ã), several methods can be used [26][27][28][29][30].The PCM Ã can be considered as the combination of three matrixes of PCM .If the measuring distance methodologies for crisp  have been properly handled, then they can be suitable for fuzzy Ã as well.We will first discuss the distance methods for crisp numbers.The measurement has logarithmic least square method (LSM) [26,31], eigenvector method [27,28], and least squares method [29,30].The optimum eigenvector should be as close as possible to the original eigenvector (derived from the original matrix).The adjustable matrix could strongly resemble the original matrix once the optimum eigenvector has been determined.Therefore, eigenvector can be used to calculate adjustable PCM   ( B ).
According to [27,28,32], if matrix  is consistent, then we could find positive weights  = { 1 ,  2 , . . .,   }  which can satisfy such condition   =   /  , ∀,  ∈ (1, 2, . . ., ).Therefore, if matrix  is close to consistent, then it must have (  /  ) −   ≈ 0. On the basis of this idea, we want to obtain the minimum value of the fastest distance between   /  and   .The fastest distance between   /  and   is defined as follows: The issue of finding the optimal modified matrix can be portrayed as resolving the minimum value of function , which can be expressed in the following equation: This equation can reach the absolute minimum value (lowest point) when even the worst situation of proportion of   /  is closest to   .The method to obtain the optimum positive eigenvector  is to find out all possible constraints and obtain feasible solutions by selecting a suitable linear programming pattern.

Constraints Analysis for Achieving Modified Matrix.
The new modified PCM  (or B) should satisfy three conditions: (1) the consistent value should be at an acceptable level; (2) the farthest distance between new PCM  (or B) and orginal PCM  (or Ã) should be as small as possible; (3) the obtained new matrix should have a strong similarity with the original matrix.These three conditions can guarantee a new consistent PCM  (or B) which has the closest and maximal similarity with the orginal inconsistent PCM  (or Ã).To suit these conditions, we provide an adjustable PCM which can reduce the inconsistency of original PCM as much as possible.On the basis of the adjustable matrix, we propose the definition of the modified matrix.
Definition 2. The new modified PCM  is defined as a combination of the original PCM  and an adjustable matrix called   , which is derived from ; the modified PCM  is defined in For fuzzy elements, the new modified PCM B is defined as follows: The properties of the adjustable PCM  (or B) and the way to obtain them will be studied in Section 4.

Calculation Processes of Obtaining
Adjustable Matrix On the basis of this idea, we build an adjustable matrix B = [ b  ] × which can reduce the inconsistency of the original matrix.The new modified matrix is constructed by two parts.(1) One is the original matrix.The function of this matrix is to keep the original information and make sure the two matrixes are in an acceptable distance.(2) One is the adjustable matrix.The function of this matrix is to modify the inconsistency level to make sure the new modified matrix's consistency is based on consistency indexes.The mathematical expression objective is as follows: The mathematical expression objective can be developed as This equation is also suitable when    =    =    , which is for crisp elements.

Stage 1: Specify Formulas to Obtain the Middle
(The calculation steps of ( 19) are given in the appendix.)Assume the value of  is given, and then the constraint (18) can be rewritten as follows: Take the reciprocal of constraint (22); then the new constraints (23), (24), and (25) can be gotten: Combine constraint (21) and constraint (25) and remove one of the inequalities, and then the new inequality can be gotten: Analogously, a similar inequality can be gotten: Next, we add slack variable , , and objective function  =  + , to change the constraints ( 26)-( 27) to equality constraints Now we specify constraints, propose formulas to calculate    ,    , and add additional stopping parameter .Then constraints (28) correspond to the following linear programming problem: 0 ≤ , , , . (34) In constraint (29), we define   1 = 1 to normalize vector ω .If the stopping parameter  = 0, then it means constraint (29) has a solution, and the problem ( 29)-(34) has a feasible solution; the value of    is the optimal solution.If the stopping parameter is  = 1 which contradicts with constraint (32), then the problem (29)-(34) is inconsistent; then the value of    is not the optimal solution.(2) The equalities ( 29)-(34) can be solved by simplex algorithm [33].The main idea of this algorithm is to walk along edges of the polytope to find out extreme points with lower and lower objective values till the minimum value is reached or an unbounded edge is visited.If the extreme point is reached, then the problem ( 29)-(34) has feasible solutions.
(3) If  =   can make problem ( 29)-( 34) have a feasible solution, then it must have  ≥   that can also make the problem ( 29)-(34) have a feasible solution.In order to find the greatest lower bound of , we investigate it in Section 4.3.
Figure 1: The structure of the segment tree.
(5) In this section, we aim to find out the least absolute worst distance by setting a more precise priority weights range and adding slack variables, and an iterative way to obtain feasible solutions will be proposed in Section 4.3.

Stage 2: Find Feasible Solutions to Obtain the Middle
Value     of Adjustable Matrix B .Next, we focus on how to find out the greatest lower bound of .The problem can be described as storing intervals of [0,  max ], analyzing the corresponding  value, and finding the greatest lower bound of  that makes  = 0.This problem can be solved by segment tree.
Substep 2 (builds a segment tree by using interval [0, ]).For example, in Figure 1, [0, ] is the interval, and P1, P2, P3, P4 is the list of distinct interval endpoints.We separate intervals into two parts in every division and terminate this process till the value of the interval is less than the accuracy level ().
Then obtain the  value in problem ( 29 Then the value of     and     is defined as follows: Proof.The value of   and   should satisfy two conditions: one is minimum fuzziness of    , and the other one is maximally maintaining similarity of original matrix Ã.
By the first condition, we can get By the second condition, we need to consider the distance between    ,    , and    , maintain the relationship among the original matrix, and make the new matrix closest to the original matrix's pattern.It means to find out the smallest coefficient between    and    and the smallest coefficient between    and    .We can get the following formulas based on the second condition: Therefore, we can prove formula (35) exists and prove that it is correct by (38).To be more precise, for instance, in Figure 2, ( b )   2 is not satisfied with condition 1, and then it is not the minimum fuzziness of associated weights.If ( b )   1 is not satisfied with condition 2, then it does not have the same pattern with matrix Ã. ( b )   3 is the optimal solution.

An Algorithm to Obtain Modified PCM B
After analysis in Sections 4.2, 4.3, and 4.4, we propose an algorithm to conclude how to obtain the modified PCM B (see Algorithm 1).

Numerical Illustration and Comparison
with Crisp Numbers

Calculation of an Illustration by Using Proposed Algorithm.
We run the experiments by software Matlab (R2009a) on a personal computer with Intel Core 2.2 GHZ and 4 G RAM.
The value of CR is more than 0.1.Therefore, we will adopt Algorithm 1 to obtain the new consistency matrix.

Calculation Step 1. Input of the initial value of Ã and initial value of 𝜉.
(1) The original matrix Ã = (   ,    ,    ); here, crisp number is a specific case of fuzzy number in our model.Set    =    =    = .(2) The acceptable precise degree  = 0.1.(Here, we adopt 0.1 as an example.)
(  (40) Now, we go to discuss the results.When in the 7th iterative time, the gap between  min and  max is 0.0787, which is less than 0.1; the process stops at this point.In order to better understand, we use the segment tree to express the changing process in Figure 3.
In Figure 3, the current  value is colored with red and bold.The first  value is 10.0725, and  = 0, which means there is a solution for inequalities ( 29)-(34); the second  value is 5.0362, and  = 0; by combining the first result, it must have a solution when Z value is between 5.0362 and 10.0725.Therefore, the program goes to search the left subtree of the node; the third  value is 2.5181, and  = 1, which means there is not a solution for inequalities ( 29)-(34); the program needs to search the right subtree of this node, because the right child node is bigger than 2.5181.We continue the process till the gap between two child nodes is less than .
On the basis of 3, we construct the new consistency matrix   as in the following: ) .
(42) [6] defined the original matrix  (the elements are   ) that can be replaced by the new matrix  (the elements are   ), which is showed in the following equation:

Comparison with References. Xu and Wei
where  is the positive value which is less than, but approaching to, 1. Cao et al. [7] proposed an equation to obtain the new matrix , which is showed by the following: where ∘ is the symbol of Hadamard product.For example,  =  ∘  means   =   ×   , ∀,  = 1, . . ., , and   is the modified deviation matrix, which is showed in the following equation: where  is the deviation matrix and DI is a zero deviation matrix when [  ] = 1.The value of  is between 0 and 1.
and  has different meanings for two papers, but the two parameters should be as close to 1 as possible.In two papers, they mentioned that  =  = 0.98 is the most suitable value to get the optimal new matrix.We will discuss 's meaning and the relationship with  and  in Section 8. We compare with two references in two situations: one is the required critical ratio (CR) less than 0.1 (Table 2); the other one is the required critical ratio (CR) close to 0 (Table 4).
From Table 2, we reach the outcome that the CR value is lower than [6,7]; at the same time, the method could achieve lower values of  and  in short iterative times.It means that our method can preserve more original information and obtain a more consistent new matrix in short iterative times.When we rank the priority weight  derived from , the ranking results are the same except when  = 0.1 in two similar weights, which means the priority weight  which is derived from our method is acceptable.
On the basis of Table 2, we go to discover the difference of COP parameter in our method and two references.The preference value between every two alternatives is gotten from priority weight  of the column "priority weight" in Table 2.For example, the value of { 1 ,  2 } is obtained in  XU = (1.00000.3241 1.0329 0.1034 0.1843 0.2172 0.9804 1.9353), which is 1/0.3241= 3.0855.
The modified new matrix  XU ,  Cao ,  Our (0.05),  Our (0.1), and  Our (0.2) is in the following: The gap of elements in every two new matrixes ( XU ,  Cao ,  Our (0.05),  Our (0.1) and  Our (0.2)) is less than 1.5 < 2, which proves these five matrixes are considerably close with each other; but, in Table 2, the CR value of our method is less than the other two methods, which means  Our (0.05),  Our (0.1), and  Our (0.2) are more consistent than  XU and  Cao .Meanwhile, the value of  and  for  Our (0.05) is the lowest, which means  Our (0.05) is the best solution for original matrix .Next, we go to discuss the second situation when CR is close to 0. The analysis results are in Table 4.
The iterative times for [6,7] 820 and 1619 are quite large, which will cost the running time.The value of  and  in these two references is far away from the acceptable value, which means the priority weight of  XU and  Cao is not acceptable at all; but the iterative times of our method are less than 10 times, the value of  is less than 4, and the value of  is less than 2, which means our result is acceptable to some extent.

Numerical Illustration and Comparison with Fuzzy Numbers
To judge the effectiveness of Algorithm 1 for fuzzy PCMs, we analyze several parameters including the inconsistency index NI [13], consistency index CCI [14,15],  fuzzy ,  fuzzy , COP, iteration times, and running time.
For this matrix, the modified matrices B are quite similar when the accurate degrees are  = 0.1 and  = 0.01, and, therefore, we only give the modified matrix B when  = 0.1 and  = 0.001 under  = 0.5 and  = 0.6: In Table 5, the inconsistency index NI is less than 0.1 and consistency index CCI is less than 0.3147, which indicates the modified matrix can reach a consistent requirement;  fuzzy satisfies 0 <  fuzzy < 2 and  fuzzy satisfies 0 <  fuzzy < 1, which shows the obtained new matrix is in an acceptable distance with Ã; the value of COP, which is gotten from priority vector ω of B, can preserve order of preference and order of intensity preference, which presents B can maintain the pattern of Ã (similarity).Meanwhile, Algorithm 1 has high convergence speed based on the less iteration times and running time.The gap among these four matrixes is quite narrow; B ( = 0.1,  = 0.5) is the best choice when we select the smallest value of the consistency index (NI, CCI),  fuzzy , and  fuzzy .

Discussion
Tables 2 and 4 conclude the effectiveness of Algorithm 1 by comparing with Cao et al. [7] and Xu and Wei [6] based on PCMs with crisp elements.Algorithm 1 can retain more original information and achieve lower value of  and  when both  [6] and  [7] approach 1. Table 5 summarizes the effectiveness of Algorithm 1 for PCMs with fuzzy elements.The modified PCMs with fuzzy elements can maintain original information and reach acceptable consistency level as well.
In Algorithm 1,  has different meaning with  and .In Cao et al. [7] and Xu and Wei [6], the results gotten by using  and  equal to 0.98 are significantly better than those by using  and  less than 0.98.Yet, it is not always the case for .The true meaning of  is the proportion of the original matrix to adjustable matrix.The approximate prefect range of  is [0.25, 0.8] based on the experiments' results, and there is no clue which value is the best one for  as the inconsistency level of PCMs is different.As a matter of fact, changing 's value will not add running time, because this operation runs in constant time.As a future research, it would be interesting to figure out the exact prefect range of .In this section, we focus on efficiency analysis of Algorithm 1 as effective analysis is done in Sections 6 and 7.   the average running time and iteration times for crisp and fuzzy elements.The elements of crisp matrix  and fuzzy Ã are gotten from randomly generated numbers by programming.The size () of matrix  and Ã varies from 3 to 10 (in real life problems, the size is usually no more than 10).We adopt twenty matrixes for each size.The assumed accuracy level  is the following three numbers: 0.1, 0.01, and 0.001.The value of  in the final step ( =    (   ) 1− ) will not affect the running time.The average running time of crisp elements is shown in Figure 4(a).The average running time of fuzzy elements is shown in Figure 4(b).In all cases, the running times are less than half a second, which is acceptable in real experiments.As expected, the higher the accuracy degree, the longer the running time in the same matrix size.The figure shows the growth rate of fuzzy matrixes' running time is more irregular than crisp matrixes.

Efficiency Analysis: Convergence Rate.
The speed of convergence is one important factor of the efficiency for an iterative method.In Algorithm 1, the only parameter that can influence iteration times is the accuracy degree (). and  are the parameters which can influence the number of iterations for [6,7], respectively.According to the numerical illustration in Section 6, we study the convergence rate with respect to the values of , , and , by using Algorithm 1 and [6,7], respectively.Meanwhile, we examine the convergence rate for fuzzy elements based on the data in Section 7. Figure 5 shows the results.
The convergence rate of Cao et al. [7] is slower than Xu and Wei [6] when  = .For these two methods, the lower the parameters value, the higher the rate of convergence.Algorithm 1 can approach its limits faster than Cao et al. [7] and Xu and Wei [6].The maximum iteration time for Algorithm 1 is less than twenty when  = 0.0001, while the maximum iteration times for Cao et al. [7] and Xu and Wei [6]  are 1619 and 820, respectively.For fuzzy elements, the convergence rate of Algorithm 1 has low growth rate as well.

Conclusions
In this work, an algorithm has been proposed to derive a consistent PCM with crisp or fuzzy elements from an inconsistent one.The presented approach used the same numerical example used by [6,7,13].The experiments reveal that the proposed approach could retain more original information and the convergence rate is faster than Cao et al. [7] and Xu and Wei [6] in different values of CR.Two effective criteria,  and , show that the distance between modified PCM and original PCM is acceptable for both crisp and fuzzy elements.A new effectiveness criterion COP reflects that the modified PCM resembles the original PCM (crisp and fuzzy elements).In conclusion, this approach could enhance the quality of vague and inaccuracy data for decision makers and also could better handle the inconsistency problem of AHP and fuzzy AHP.
In the future work, we will apply this approach to different applications.Meanwhile, efforts should be made to explore the prefect range of  and study the issue of selecting the situation where we can supply a consistency PCM to an inconsistency one based on real-life meaning.

Conclusion
Cao et al. [7] are closer to the original matrix , but the gap among Cao et al. [7], Xu and Wei [6], and our method is very narrow (less than 1), the three methods almost in the same extends to close to original .Thus, our method is acceptable in this parameter.
Given a vector of fuzzy weights ω = (ω 1 , ω2 , . . ., ω ), which is the corresponding eigenvector matrix ω closest to original matrix Ã, the value of eigenvector matrix ω should satisfy the following formula:

Figure 3 :
Figure 3: An example of iterative times in segment tree structure.

Figure 4 :
Figure 4: Average running time with fuzzy and crisp elements (in seconds).
,    ,    ),    is the lower value,    is the middle value, and    is the upper value.In the particular situation, if  = , when the following condition is satisfied, then Ã is a reciprocal matrix: The aim of this paper is to find out a consistent matrix B = [ b ] × which can preserve the most information of the original matrix Ã = [ã  ] × and be the closest to Ã.One method to achieve this aim is that setting one element of B is b = (   ,    ,    ) and one element of Ã is ã = (   ,    ,    ), by making every element of B closely resemble the elements of Ã, which is    ≐    ,    ≐    , and    ≐    .Therefore, we go to find    ,    ,    separately.After finding out    ,    ,    , combine them together to get the new matrix B.
)-(34) by setting the current  value.If  = 0, then the next  value is equal to the lower bound of the current node.adapted to counting queries, and it helps us to query the number of segments that contain a given point.combination of the adjustable matrix and the original matrix based on Definition 2; then the adjustable matrix should have the minimum fuzziness and maximum preservation of the original matrix's pattern.If B is minimum fuzziness, then fuzziness of B will mostly come from Ã, and Bwill be more similar with Ã.In fact, the minimum fuzziness of B could reduce uncertainty factors of B. If B could maximally preserve the pattern of Ã, then the combination of B and Ã could reach the most potential of similarity with Ã.We propose Theorem 3 to obtain the value of    = 1, 2, . . ., .Set   ,   as arbitrary positive constants.Define the value of    ,    in the following formulas: Substep 3. Select the greatest lower value of  when  = 0; obtain the feasible solution of    ,    to achieve     .Note 2. The segment tree is special for storing intervals.The built time is ( log ) for  intervals, and it uses ( log ) storage.The reason we adapted to the segment tree is because the segments can be stored in any arbitrary manner, it can easily be

Table 1 :
The values of parameters for each iterative time.
) First time  = 10.0725;substitute  into Algorithm 1; The value of  can determine how to change the current  value.If  = 0, it means there is a solution for problem (29)-(34); then save the value of vector  and set  max = ; else  = 1; it means there is not a solution for problem (29)-(34); then set  min = ,  =  max .
: Algorithm Running Time.To demonstrate the overall efficiency of Algorithm 1, we measure

Table 3 :
COP (condition of order preservation) parameter value.

Table 5 :
Effective analysis for Algorithm 1 based on nine parameters.