#### Abstract

In a supply chain system, the prices with which the suppliers supply its local commodity to the retailers should satisfy the requirements of the retailers and the consumers. The supply and demand scheme satisfying these requirements is reduced into fuzzy relation inequalities (FRIs) with min-product composition. Due to the difference between the min-product composition and the classical max-t-norm one, we first study the resolution of such min-product FRI system. For optimization management in the supply chain system, we further investigate a maximin programming problem subject to the min-product FRIs. An algorithm is proposed to obtain the optimal solution based on the quasi-maximal matrix and corresponding index set. To illustrate the efficiency of our proposed algorithm, we provide a simple numerical example. The obtained optimal solution reflects an optimal pricing scheme, which maximizes the minimum prices of the commodity from the suppliers.

#### 1. Introduction

The classical fuzzy relational equation (FRE) system could be formulated asin which All the elements in , , and belong to . The composition in the equation is the classical* max-min* (). For the convenience of presentation, we always denote the index sets and from now on. Moreover, the components in the right hand vector are assumed to be in this paper.

Concept of FRE was first introduced by E. Sanchez [1–3]. Since E. Sanchez investigated the max-min FREs with application in medical diagnosis, some fuzzy mathematical scholars became to focus on FRE and its resolution. Efficient method for finding out all its solutions was one of the most interesting and important research topics [4–10]. The operations in a classical linear equations system are addition (+) and multiple (.). However, the operations in system (1) are () and (). Hence resolution of system (1) is much different from that of the classical linear equations. As known to everyone, solution set of the classical () is usually a convex set. However the solution set of system (1) (when it is consistent) is nonconvex in most cases. We call system (1) consistent (or solvable), iff , whererepresents its complete solution set.

It is well known that the consistency of (1) could be checked by its potential maximum (or greatest) solution. The necessary and sufficient condition, with which system (1) is consistent, is there exists a (unique) maximum solution. When the system is consistent, then its complete solution set could be formed and generated by the unique maximum solution and all the minimal ones. In general situation, a consistent system of max-min FREs has a finite number of minimal solutions. Although the number of minimal solutions is finite, it is difficult to find out all of them. The set of all minimal solutions could be found theoretically. But it is still a challenging problem until now. Due to its important role in solving (1), many scholars tried their bests to develop novel and efficient approaches to obtain the minimal solution set [11–17].

In the FRI system (1), the max-min composition could be extended to other ones, such as max-product and addition-min [18–22]. Data transmission on the Peer-to-Peer (P2P) network system was investigated by Yang et al. FRIs with addition-min composition were introduced to describe the transmission mechanism in such system [18, 20, 21]. The author established some optimization management models to the P2P network system. Specific methods were proposed to search one of the optimal solution(s), which represent(s) an optimal quantity of flow.

In a fuzzy relational equations system, the composition operator is indispensable. As mentioned above, the most common and typical one is* max-min* (). However, for wide application, the composition was generalized to max--norm one. Here is a continuous triangular norm. However we have to point out, the minimum operator, as a specific kind -norm, is irreplaceable and most frequently adopted in many application fields. Moreover, resolution method for a max--norm FRE or FRI system is similar to that for system (1). The structures of their complete solution sets are also the same.

In many practical application fields, mathematical programming with FREs or FRIs constraint was established and investigated, for describing the corresponding optimization model. Resolution of such optimization problems is usually related to the properties and structure of the feasible domain (i.e., the solution set of a FREs or FRIs system), and also related to the characteristic of the objective function.

Wang [23] was the pioneer who investigated the fuzzy relation programming problem. In such optimization problem studied in [23], the target is the latticized linear function, and the constraint is formed by a group of FRIs composed by (). The constraint system was firstly discussed. To characterize the feasible domain, the authors hoped to find out all the minimal solutions to the constraint. Conservative path approach was the effective method that proposed in [23]. Characteristic matrix was defined for the constraint system, based on which all the conservative paths could be found. It was further proved that a conservative path always corresponds to a minimal solution. This sets up the relationship between the minimal solution and the conservative path (formed by some indices). Due to the monotonicity of the objective function, the optimal solution of the target problem should be one of the minimal solutions. Hence, after finding out all the minimal solutions by the conservative path approach, the optimal solution could be selected by simple comparison calculation. S.-C. Fang is another settler on fuzzy relation optimization problem. Together with G. Li [24] and J. Loetamonphong [25], he proposed and studied the minimization problem with FRE constraints. Different from the optimization problem presented in [23], Fang and Li focused on the classical linear function. They tried to minimize a linear function, i.e., , with the constraint of a group of max-min FREs. Corresponding problem was written asIn problem (4), the coefficients were classified into two kinds, negative and nonnegative. Correspondingly, each kind of coefficients formed a new subproblem. Two subproblems were generated according to these objective coefficients. One of the subproblems was constructed by the negative coefficients. This subproblem could be simply solved by the maximum solution of the feasible domain. But resolution of the other subproblem was much harder. It might take much more computation cost. The optimal solution of such nonnegative subproblem could be selected from the minimal solutions. However, for decreasing the computation cost, the authors [24] avoided to compute all the minimal solutions. Novel way was developed for searching the optimal solution. They defined some index sets, corresponding the so-called quasi-minimal solutions. Complete feasible domain was characterized by the index sets. The nonnegative subproblem was further converted into 0-1 integer programming. Hence, the optimal solution could be found step by step according the the typical branch-and-bound approach. In fact, the idea that dividing the main problem into a negative-coefficients subproblem and a nonnegative-coefficients one, was further deeply developed and widely used afterwards, in other relevant fuzzy relation optimization problems. It turns out to be the most important method for such kind of problems [26–31].

Optimizing a linear function, i.e., , which is restricted to the a FREs or FRIs system, becomes attractive and interesting research topic. Since the feasible domain is usually nonconvex, classical optimization method seems to be useless. To deal with such problem, some researchers improved the existing methods, or proposed some novel resolution methods [27–33]. In recent years, some researchers have turned their attentions to the fuzzy relation nonlinear optimization problems [34–41], especially to the fuzzy relation geometric programming problems [42–49].

Bingyuan Cao [42, 44, 45] was the one who proposed the so-called “fuzzy relation geometric programming” for the first time. Some special forms of the geometric function were first studied. In [42], Yang and Cao gave the mathematical formulae of the monomial geometric programming problem. In such optimization problem, the objective function was nonlinear. Corresponding mathematical model was described as follows:In the above problem, the coefficient and all the exponents are constant numbers. Notice that the feasible domain could be written as a union of finite close intervals, whose left endpoints are indeed the minimal solutions. There are finite number of closed intervals. The optimal solution of one of the subproblems could be chosen from the minimal solution set. One just need to compare all the function values of the minimal solutions. Similar problem was also studied by Shivanian and Khorram [46]. They further improved the resolution method in [42], for the following fuzzy relational optimization: In order to accelerate resolution procedures and decrease the computation, they deleted the components of , who were redundant and useless for searching the optimal solution.

With similar constraint, the mathematical modelsandwere investigated by Y.-K. Wu [43] and X.-G. Zhou [50], respectively. In problems (7) and (8), the coefficients are all nonnegative, i.e., and , . The authors did not consider the negative coefficients. To solve these problems, it was unnecessary to compute all the (potential) minimal solutions. Hence the computation was decreased in some degree. Polynomial algorithms were proposed for such problems [43, 50].

Recently, fuzzy relation inequality with min-product composition was defined and investigated [51–53]. It was first introduced to describe the pricing in the supply and demand scheme [53]. We view this supply and demand scheme as a simple situation of a supply chain system. In [53], the authors only considered the requirement of the retailers. In such case, the requirement of the suppliers might not be fulfilled and a group of min-product fuzzy relation inequalities was constructed. In order to satisfy both requirement of the retailers and the suppliers, Yang et al. [51, 52] further established and studied the min-product system with one more constraint . In [51, 52], checking method of consistency of a min-product system and structure of its solution set (when consistent) were investigated. Based on the specific structure of the complete solution set of min-product fuzzy relation inequalities system, Yang et al. [52] further studied a lexicographic optimization problem. Detailed resolution algorithm was proposed for obtaining the unique optimal solution of such problem. The optimal solution is indeed an optimal pricing scheme maximizing the profits of the suppliers in a fixed lexicographic order relation. Or in other words, the profits of the suppliers are maximized with fixed priority grade. Hence, the suppliers are treated distinguishingly. However in some cases, the suppliers should not be treated with fixed priority grade. They should be treated equally. Based on such equalitarianism consideration, we propose a new type of optimal model, i.e., the maximin optimization problem subject to min-product fuzzy relation inequalities, and investigate its resolution method in this work.

The structure of the rest part of our work is as follows. In Section 2 we introduce the min-product FRIs system. Relevant properties and results are presented. Matrix-based method is proposed to find out all the solutions of such system in Section 3. In Section 4, we study the maximin programming problem subject to min-product FRIs system. Resolution algorithm is developed to search its optimal solution, with illustrating example. Sections 5 and 6 are simple discussion and conclusion respectively.

#### 2. Preliminaries

##### 2.1. Pricing Relation and the Corresponding Mathematical Model in a Supply Chain

This subsection provides simple application background for our studied min-product FRIs, which is adopted from [51, 52].

As shown in Figure 1. We assume that the supply chain is composed by suppliers and retailers. The suppliers were denoted by , while the retailers were denoted by . Suppose the suppliers will supply a sort of commodities to the retailers. The trade is free. We consider the subjective factors whether the suppliers are willing to supply its local commodities. What influence the transaction is the price of the commodities. For convenience to express, we define two index sets, i.e., Suppose the selling price of the commodities from is . Considering the transportation and other costs, when the commodities are received by , the price should be bigger than . We denote this received price by , where is a real number bigger than 1. There is no doubt that the retailer will choose the suppler with the cheapest received price to stock the commodities. Assume that the price requirement of the commodity of retailer is no bigger than , , then we have

On the other hand, if the prime (cost) prices of the commodities for the suppliers are , respectively, then to ensure the supplier is profitable, it holds thatCombining these two aspects, the prices should satisfy the following system:

In system (12), since the requirements of all the retailers are satisfies, this indicates all the retailers are supplied with the commodities by at least one supplier. However, this is not able to make all supplier to supply its local commodities to at least one retailer. That is to say, although some price scheme satisfies system (12), there might exists some supplier who has not been chosen by the retailer. In order to make all suppliers able to supply their local commodities, a feasible price scheme should also fulfill at least one of the inequalities below,Satisfying at least one of the inequalities in system (13) is in fact equivalent to

Combining both inequalities (12) and (14), we get

After standardization of all the variables and parameters, system (15) becomes a system of min-product FRIs, with , . Moreover, in this paper we always denote , . Then the matrix form of (15) turns out to beIn system (16), the entries of matrix and vectors are respectively.

##### 2.2. Consistency Checking for System (16)

In this subsection we provide some basic concepts and existing results of (16). The min-product FRI system has been studied in [52]. The authors investigated some basic properties of the system, including consistency checking and structure of its solution set. Besides, lexicographic maximum solution, as a specific solution, was defined and introduced to maximize the profits of the suppliers under some fixed priority grade [52].

In the rest we always denote for convenience. Moreover, the solution set of (16) is denoted byin which , . System (16) is called* consistent* unless . Otherwise, it is said to be* inconsistent*.

*Definition 1 (see [51, 52]). *The minimum (or minimal) element in is called minimum (or minimal) solution of system (16). Analogously, the maximum solution and maximal solution could be defined in the same way.

The minimum solution and maximal solution are important for generating all the solutions of (16). Before computing all the solutions of (16), we present method for checking the system’s consistency of (16).

Theorem 2 (see [51, 52]). *System (16) is consistent iff .*

It is clear that, if system (16) is consistent, then the vector should be its minimum solution.

Define the matrix , where

Theorem 3 (see [51, 52]). *System (16) is consistent if and only if there exists a nonzero element in every row and every column in the matrix .*

Both Theorems 2 and 3 could be applied to check the consistency of (16).

When system (16) is consistent, i.e., its solution set is nonempty, then the special structure of this solution set is shown in the following Theorem 4.

Theorem 4 (see [51, 52]). *For the consistent system (16), suppose represents its solution set. Then we have Here and is the collection of all its maximal solutions.*

Since is self-evident, solving is equal to computing . In next section we will provide matrix-based method for solving , i.e., the solution set of (16).

Illustrative example for the above-provided consistency-checking theorems could be found in References [51, 52].

#### 3. Matrix-Based Resolution for System (16)

In this section we provide a novel method for obtaining the complete solution set of system (16). The solution method is based on the below-defined concept of quasi-maximal Matrix. Thus we call it matrix-based method.

Theorem 5. *Take arbitrary . Suppose . Then if and only if for any , there exists a corresponding index such that .*

*Proof. *() Since , it holds for any that Notice that is a finite set. For any there exists such that Thus we have .

() If for any , then Combining , it is clear that .

Based on the above-provided matrix , we define concept of quasi-maximal matrix as follows.

*Definition 6 (Quasi-maximal matrix). *A matrix is called a quasi-maximal matrix, if it fulfills the following conditions:(i);(ii)for any , , where .It is shown in Definition 6 that the element in a quasi-maximal matrix is either or . Moreover, in each row in the quasi-maximal matrix , there exists and only exists one nonzero (or positive) element. By Theorem 3 and Definition 6 we quickly get the following Corollary 7.

Corollary 7. *If system (16) is consistent, then there exists at least a quasi-maximal matrix.*

For system (16), we denote the set of all quasi-maximal matrices by . Following Corollary 7, when (16) is consistent, it holds that .

For a given quasi-maximal matrix , defineThen we are able to construct a corresponding vector, denoted by , as follows:

Theorem 8. *If is a quasi-maximal matrix of system (16), then the vector defined by (25) is a solution of (16), i.e., .*

*Proof. *(i) . Take arbitrary . If , then . If , then for any , . That is . Thus . On the other hand, it is clear that . Hence and then(ii) . Take arbitrary . If , then . If , then for any , similarly we get Hence(iii) Since is a quasi-maximal matrix, it follows from Definition 6 that and , for arbitrary . It is reasonable to assume that , . Then it holds that . Hence and .

Notice that . We have According to Theorem 5, the above-proved (i), (ii) and (iii) contribute to .

Theorem 9. *Let be an arbitrary solution of system (16). There exists a quasi-maximal matrix , such that .*

*Proof. *Since , is evident. According to Theorem 5, for each , there exists a such thatUntil now we have found indices . Inequality (30) indicates . By (18), we haveDefine , in whichAccording to Definition 6, it is easy to check that , i.e., is a quasi-maximal matrix of system (16). To complete the proof, we just need to verify that . Take arbitrary . Next we check in two cases.

(i) If , then . By (24), for any , it holds that . By (32), it holds thatandIt is obvious that , . Combining (30), we haveFor , it holds thatHence .

(ii) If , then . It is obvious that .

Combining the above Theorems 8 and 9, the following Theorem 10 is self-evident.

Theorem 10. *Suppose (16) is consistent system, with minimum solution . Besides, the set of all its quasi-maximal matrices is denoted by . Then we have*

As shown in Theorem 8, for arbitrary , the vector defined by (25) is indeed a solution of (16). We call the quasi-maximal solution corresponding to the . Hence the set of all quasi-maximal solutions isBased on the definition of quasi-maximal matrix, it is easily found that is a finite set. Consequently, is also a finite set. Moreover, the relation between and the maximal solution set is as described below.

Proposition 11. *In system (16) it holds that .*

To get the maximal solution set of system (16), denoted by previously, we may delete all the quasi-maximal solutions which are not maximal from the set .

In fact, Theorem 10 indicates a resolution approach for obtaining the complete solution set of system (16). Based on the matrix , we can compute all the quasi-maximal matrices. Then we further get the set of all quasi-maximal matrices, i.e., , and the set of all quasi-maximal solution, i.e., . At last we are able to find the solution set by Theorem 10 or by Theorem 4.

#### 4. Maximin Programming with Min-Product FRIs Constraint

In this section we study the maximin programming problem subject to system (16), i.e.,In problem (39), the constraint is system (16).

As pointed out previously, the the supplier supply its local commodity with price . The bigger value take, the higher profit gain. To enhance the profit of the suppliers, we should maximize the prices . However, in most common cases, one is not able to maximize all the valuables simultaneously. In this paper, based on the fairness consideration, we try to ensure the worst profit of all the suppliers. To reach such goal, we should maximize the minimum price, i.e., . As a consequence, we construct the corresponding mathematical model as (39).

Next we investigate some properties of problem (39) and try to develop an algorithm to find out one of its optimal(s).

Theorem 12. *If system (16) is consistent, then problem (39) has at least one optimal solution. Furthermore, there exists a maximal solution , such that is an optimal solution of (39).*

*Proof. *If system (16) is consistent, then is a finite set and . Obviously is also a finite set. We may assume that Let , and Obviously there exists such that i.e., Next we verify that is an optimal solution of problem (39). Notice that is obviously a feasible solution of (39). It corresponding function value is . To finish the proof, we should further check that holds for any .

Take arbitrary . According to Theorem 4, there exists a maximal solution such that . Since , we have The proof is completed.

Proposition 13. *If is an optimal solution of problem (39), then any vector fulfilling is also an optimal solution of (39).*

Theorem 12 allows us to gain an optimal solution of problem (39) by comparing the function value of all the maximal solutions of system (16). However, it is difficult to obtain all the maximal solutions. Hence we usually do not recommend this resolution approach. Next we will provide another effective resolution algorithm for problem (39).

Based on the matrix defined by (18), we further define the index sets as follows. For each , letandSuppose denotes the Cartesian product of .

For a given , define matrix , in which

Proposition 14. *If system (16) is consistent, then holds for all .*

*Proof. *If system (16) is consistent, it follows from Theorem 3 that each row in the matrix has at least one nonzero element. The nonzero element should be positive. That is, for any , there exists some . Hence .

Proposition 15. *If system (16) is consistent, then for any , the matrix defined by (48) is a quasi-maximal matrix.*

*Proof. *The proof is trivial according to expression (48) and Proposition 14.

Theorem 16. *Suppose system (16) is consistent, and is defined by (48) based on . Then is an optimal solution of problem (39).*

*Proof. *Denote .

(i) According to Proposition 15 and Theorem 8, is a feasible solution of problem (39).

(ii) Take arbitrary . For any , it follows from Theorem 5 that there exists such that Thus . On the other hand, it is clear that . Then it holds that Hence Due to the arbitrariness of , we have . On the other hand, indicates . It is obvious that . So we get (iii) Take arbitrary .*Case 1*. If , then .*Case 2*. If , thenandInequality (54) indicates . Considering formulae (48), we further get Combining Cases 1 and 2, it holds for any that . So we have .

On the other hand, the above-proved point (ii) implies that since (see (i)). Hence it turns out to be .

Points (i), (ii) and (iii) contribute to the optimality of . That is to say, is an optimal solution of problem (39).

To find an optimal solution of problem (39) without computing all the maximal solutions of (16), we now propose a novel resolution algorithm as follows.

*Algorithm for Solving Problem (39)*

*Step 1. *Compute the matrix by (18).

*Step 2. *Check the feasibility of problem (39). If there exists at least a nonzero element in each row and also in each column in the matrix , then system (16) is consistent and problem (39) is feasible, continue to Step 3. Otherwise, and problem (39) has none optimal solution, stop.

*Step 3. *Compute by (46).

*Step 4. *Compute by (47).

*Step 5. *Construct the Cartesian product and take arbitrary index vector .

*Step 6. *Based on the index vector , construct the quasi-maximal matrix by (48).

*Step 7. *Compute the index sets by (24).

*Step 8. *Compute the quasi-maximal solution corresponding to by (25). Then is an optimal solution of problem (39) according to Theorem 16.

Next we show the computational complexity of our proposed resolution algorithm.

*(i) Computational Complexity* Number of variables; Number of inequalities.

Computing the matrix in Step 1 costs operations. In Step 2, checking whether the elements in each row and in each column of are all zeroes costs operations. Computing in Step 3 costs operations, while computing in Step 4 costs operations. In Steps 5 and 6, it costs operations for obtaining the quasi-maximal matrix . Computing the index sets in Step 7 costs operations. At last, Step 8 costs operations for generating the optimal solution of problem (39), based on the quasi-maximal matrix . As a consequence, the above-proposed algorithm costsoperations in total. Hence the computational complexity of the algorithm is . It also indicates our proposed algorithm could be achieved in polynomial time.

For illustrating the feasibility of the above-presented algorithm, we provide a numerical example as below.

*Example 1. *Letbe the matrix form of a system of min-product FRIs. In system (57), the matrix A and the vectors are as Check whether system (57) is consistent. If it is consistent, find an optimal solution of the maximin programming subject to system (57), i.e.,*Solution**Step 1-2*. In fact, the matrix is [51] Following Theorem 3, system (57) is consistent and problem (59) is feasible.*Step 3*. Compute by (46).*Step 4*. Compute by (47).*Step 5*. Construct the Cartesian product and take arbitrary index vector . We take the index vector .*Step 6*. Based on the index vector , we are able to compute the elements in by (48) as follows: The quasi-maximal matrix corresponding to is *Step 7*. Compute the index sets by (24). *Step 8*. Compute the quasi-maximal solution corresponding to by (25). Firstly, Since , we have In addition, Hence we find an optimal solution of problem (59) as

#### 5. Discussion

##### 5.1. Optimal Solution Set of Problem (39)

In Example 1, we have found an optimal solution of problem (59) as The corresponding optimal function value is In Step 5, we take the index vector as . Following Steps 6–8, we get the vector . It is easy to check that satisfies the constraint in problem (59) and moreover, . Therefore, is also an optimal solution of problem (59).

It has shown above that the optimal solution of problem (59) is not unique. Hence we further discuss how to obtain all the optimal solutions of our proposed maximin optimization problem in this subsection. The following Theorem 2 gives exactly description of the set of all optimal solutions to problem (39), according to one of the optimal solutions obtained by our proposed algorithm in Section 4.

Theorem 2. *Let be an optimal solution to problem (39), with corresponding optimal value . Then the set of all optimal solutions to problem (39) is the solution set of the following inequalities,*

*Proof. *We complete the proof from two aspects as follows:

(i) Take arbitrary satisfying system (73). It holds that and , . Hence is a solution of system (16), i.e., a feasible solution of problem (39).

Since is a feasible solution of problem (39), while is an optimal solution with objective function value , we haveOn the other hand, since satisfies system (73), we have , .Inequalities (74) and (75) contribute to .

Hence is an optimal solution of problem (39).

(ii) Suppose is an optimal solution of problem (39). Next we have to check that satisfies system (73). It is obvious that satisfies the constraint of (39), i.e., system (16). On the other hand, due to the optimality of , we have This indicates for all . Hence , , and is a solution of system (73).

It is clear that system (73) is indeed a system of min-product fuzzy relation inequalities. In this sense, the set of all optimal solutions to problem (39) (when not unique), is exactly the solution set of a min-product fuzzy relation inequalities system. According to Theorem 10, there exist a unique minimum optimal solution and a finite number of maximal optimal solutions.

##### 5.2. Comparison to the Existing Works

In this subsection we make some comparison between the results presented in this paper and those in the existing works [51–53].

In [53], only requirement of the retailers was considered. Both solutions and strong solutions were studied in [53]. The authors focused on the strong solutions. Structure and resolution method of the strong solution set are the major part of this work. A path-based method was proposed to find out all the strong solutions of a system of min-product fuzzy relation inequalities. No optimization problem was considered in this paper.

In [51], although both requirements of the retailers and the suppliers were considered, the authors only discussed the method of consistency checking of the min-product fuzzy relation inequalities system. In such system, a new constraint, i.e., , was added. The consistency of such system could be checked by its potential minimum solution , or by its distinguishing matrix . Besides, structure of the solution set to such system was simply presented. No optimization problem was considered in this paper.

In [52] optimization problem under lexicographic order was considered. The constraint of the problem was exactly the min-product system proposed in [51]. In such optimization problem, the target was to maximize the profits of the suppliers. Under a specific lexicographic order, the suppliers were treated with fixed priority grade. This is conflict with the equalitarianism. The authors proposed a cyclic algorithm to obtain the unique optimal solution. There were procedures in a cycle. In each procedure, a new vector was generated and its feasibility was checked. Moreover, procedures produced components of the optimal solutions.

Both of the studied problem and the resolution method are different from those in the above-mentioned existing works. In this paper, we mainly study the maximin problem subject to min-product fuzzy relation inequalities, considering both requirements of the retailers and the suppliers. Quasi-maximal matrix is defined for investigating the structure of the complete solution set of the proposed min-product system. Based on the special structure of the solution set, a novel approach based on the quasi-maximal matrix and its corresponding index sets is proposed to find an optimal solution of the maximin optimization problem. The motivation of our proposed maximin optimization problem has been stated in Section 1.

#### 6. Conclusion

FRIs with min-product composition could be applied to describe the price conditions in a supply chain. Consistency checking for the min-product system (16) has been discussed in [52]. In this work we further investigate the resolution method of (16). In fact, each solution of (16) represents a feasible pricing scheme of the suppliers in the supply chain. Consistency and solution-set characteristic of a system of min-product fuzzy relation inequalities are discussed in detail. The consistency could be checked in two ways. When consistent, the complete solution set of a system of min-product fuzzy relation inequalities is fully determined by its unique minimum solution and finite maximal (or quasi-maximal) solutions. Matrix-based method is proposed to compute all the solutions of (16). However, in practical application, it is no need to compute all the solutions in some cases. Optimal solution might be more meaningful, considering some specific optimization management objectives. In this paper, maximin programming is established for improving the minimal profit of the suppliers. In such maximin optimization problem, the objective is to maximize the maximin function, i.e., . We propose a novel algorithm to search an optimal solution, in case that the feasible domain is nonempty. The algorithm is developed step by step and illustrated by a detailed numerical example.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Authors’ Contributions

Hai-Tao Lin and Xiao-Bin Yang contributed equally to this work.

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China (61877014) and other funds provided by Guangdong Province and Hanshan Normal University (2016A030307037, LQ201702, and QD20171001).