Decision support model with Pythagorean fuzzy preference relations and its application in financial early warnings

Pythagorean fuzzy sets (PFSs) retain the advantages of intuitionistic fuzzy sets (IFSs), while PFSs portray 1.57 times more information than IFSs. In addition, Pythagorean fuzzy preference relations (PFPRs), as a generalization of intuitionistic fuzzy preference relations (IFPRs), are more flexible and applicable. The objective of this paper is to propose a novel decision support model for solving group decision-making problems in a Pythagorean fuzzy environment. First, we define the concepts of ordered consistency and multiplicative consistency for PFPRs. Then, aiming at the group decision-making problem of multiple PFPRs, a consistency improving model is constructed to improve the consistency of group preference relations. Later, a consensus reaching model is developed to reach the degree of group consensus. Furthermore, a decision support model with PFPRs is established to derive the normalized weights and output the final result. Holding these features, this paper builds a decision support model with PFPRs based on multiplicative consistency and consensus. Finally, the described method is validated by an example of financial risk management, and it is concluded that the solvency of a company is an important indicator that affects the financial early warning system.


Introduction
Decision-making implies that there is a range of choices. In this case, decision-makers (DMs) rank these choices from best to worst, and choose a goal that meets the expectations of DMs. In group decision-making (GDM) problems, DMs usually need to give a precise preference for a group of options using pairwise comparisons to express their preferred information and construct a preference relation judgment matrix [1][2][3]. Obviously, there are many methods to support decision-making [4][5][6]. Nevertheless, in many GDM problems, DMs may be challenging to depict their preference with a precise number. One reason is that DMs may not have an accurate or sufficient knowledge level of the problem, and the other is that DMs cannot clearly distinguish the extent to which one is superior to others [7,8]. In these cases, DMs may prefer to compare the imprecise judgment information in a matrix in pairs. To describe this vagueness and indeterminacy, researchers have introduced fuzzy sets (FSs) [9] theory to different preference relations (PRs). Then, fuzzy preference relations (FPRs) [10] and interval fuzzy preference relations (IFPRs) [11] emerged. Compared with traditional PRs, FPRs can provide reasonable feedback on the preferences given by DMs via a membership function [12][13][14]. With the growing complexity of the decision-making environment, there will be recognitions of affirmation, hesitation, or negation in decision-making. Therefore, IFPRs can feedback assessment intuitively with membership, nonmembership, and hesitation [15][16][17][18][19][20].
However, the domain inscribed by the IFSs [21,22] can be visualized as a triangular domain using a two-dimensional coordinate axis. Yager et al. [23,24] proposed PFSs, which are a further extension of IFSs. PFSs not only possess a membership function and non-membership function, but also extend the area composed of two functions to 1/4 circles. While maintaining the characteristics of IFSs, PFSs have a wider range. Therefore, PFSs are more prominent in handling decision-making problems. Currently, the research on PFSs mainly involves aggregation operators [25][26][27], similarity measures and distance measures [28][29][30], and extensions of classical methods [31,32]. To give a clear perspective on the different concepts, tools, and trends associated with their extension, Peng and Selvachandran [33] provided an overview of PFSs.
Given the development of PFSs, PFSs were introduced to PRs, the concept of Pythagorean fuzzy reference relations (PFPRs) [34,35] was proposed, and the GDM problem under PFPRs was studied.
The core of the decision-making problem based on PRs includes how to define the consistency of PRs and how to obtain the weight vector of PRs. Chu et al. [18] constructed a GDM framework with IFPRs based on additive consistency and the group consensus. Jin et al. [19] designed two methods to obtain the normalized priority weights from IFPRs. To represent preferences under uncertain weights, Lin and Wang [20] minimized extreme deviation from the relevant consistent IFPRs. Zhang [36] proposed a new approach index of Pythagorean fuzzy numbers (PFNs) and introduced a PFN ranking method based on the approach index. Additionally, the concept of the interval-valued Pythagorean fuzzy set (IVPFS) parallel to the interval intuitionistic fuzzy set was proposed [36]. In view of the imprecise modeling ability of PFS for human perception in multi-criteria decision-making, Zhang et al. [36] introduced the classic preference ranking organization method of enrichment evaluations (PROMETHEE) in the Pythagorean fuzzy environment. In Zhang et al.'s method [37], the weights and preferences are considered synthetically, which increase the range of the decision-maker's choice. Based on a new multiplicative consistency definition of interval fuzzy preference relations, Zhang [38] constructs a goal programming model.
However, when using additive consistency to estimate incomplete elements in FPRs, it is easy to exclude partial incomplete elements from [0, 1]. Thus, the final decision result may be unreliable. In some cases, additive consistency is not suitable for fuzzy preference relations. In contrast, the multiplicative consistency in FPRs can effectively compensate for this deficiency. Furthermore, Jin et al. [19] discussed multiplicatively consistent IFPRs and constructed a group decision model for solving practical problems. However, we know that the area of an intuitionistic fuzzy set (IFS) can be regarded as a triangle while the area of a Pythagorean fuzzy set (PFS) is a quarter circle. In short, the spatial membership of a PFS is greater than that of an IFS, which is equivalent to an intuitionistic fuzzy number (IFN) being a PFN. Otherwise, it does not hold [22,23]. Therefore, PFPRs can solve the problem that IFPRs are incompetent. Yang et al. [35] defined additive consistent PFPRs and Pythagorean fuzzy weighted quadratic (PFWQ) operators to construct group decision methods. In the GDM process, due to the complexity of the decision-making problem and the limitation of the evaluator's knowledge or background, not all PRs provided by each evaluator can meet the condition of additive consistency. Therefore, an optimization model with the minimum deviation between the given PFPR and the constructed multiplicatively consistent PFPR is established, and the objective optimization model is solved to obtain the fuzzy weight vector of the preference relation. It is worth affirming that Mandal and Ranadive [34] introduced PFPRs. However, the qualified consensus relation cannot be obtained through the given group collective PR [34]. He et al. [39] gave the multiplicative PFPR without considering the consensus among multiple DMs. Hence, it is only suitable for individual decision-making. In summary, how to generate a decision support model based on acceptably multiplicatively consistent PFPR and consensus is of great significance for GDM problems. Although it is worthwhile to discuss the multiplicative consistency under PFPRs, there are few related achievements. It is necessary to discuss the decision support model from the perspective of multiplicative consistency and consensus.
Taking the essential advantages of the Pythagorean fuzzy environment and multiplicative consistency to describe preference values, the main contributions of this paper are as follows.
1. We define the order consistency and multiplicative consistency with PFPRs. 2. We present three algorithms for ordered consistency, multiplicative consistency, and group consensus with PFPRs, respectively. 3. We demonstrate the proposed decision support model with an application in financial early warning to indicate its practicability and effectiveness.
The remainder of this paper is organized as follows. In "Preliminaries", some basic concepts, including FPRs, PFSs, and PFPRs, are briefly reviewed. "Method for deriving priority weight vectors with PFPRs" defines the multiplicative consistency and properties of PFPRs. "Multiplicative consistency improving algorithm" introduces consensus reaching with multiplicatively consistent PFPRs. Then, "Algorithm for group decision-making with PFPRs" develops a model for group decision-making with PFPRs. "Numerical example and comparison analysis" presents one illustrative example of financial early warning, and comparative analyses are implemented simultaneously on the proposed approach and other existing methods. Finally, conclusions are summarized in "Conclusions".

Preliminaries
To illustrate this further, we introduce some basic concepts.

Fuzzy preference relations
Definition 1 [10] Assume that X = x 1 , x 2 , … , x n is a finite set of alternatives. A fuzzy preference relation R on X is defined by a reciprocal matrix R = r ij n×n with the circumstance that where r ij indicates the degree of preference of alternative x i relative to x j .
If r ij = 0.5 , there is no difference between alternative x i and alternative x j . If r ij > 0.5 , alternative x i is superior to alternative x j . If r ij = 1 , alternative x i is completely superior to alternative x j .

Intuitionistic fuzzy preference relations
Definition 2 [15] Assume that X = x 1 , x 2 , … , x n is a finite set. The intuitionistic fuzzy preference relation is denoted as where I (x i ) and I (x i ) represent the degree of membership and degree of non-membership of the elements x i ∈ X to the set I, respectively. Let = ( , Definition 3 [15] Let i = i , i , i = 1, 2 be double IFNs that satisfy the following sequential relations: Here, the score function Δ( i ) = i − i , and the accuracy function

Pythagorean fuzzy set
Definition 4 [23,24] Assume that X = x 1 , x 2 , … , x n is a finite set. The Pythagorean fuzzy set is denoted as where A (x i ) and A (x i ) represent the degree of membership and degree of non-membership of the elements x i ∈ X to the set A, respectively. Let Definition 5 [40] Let i = i , i , i = 1, 2 be PFNs. Then, we have the following: Here, the score function s( i ) = 2 i − 2 i , and the accuracy function

Pythagorean fuzzy preference relations
Definition 6 [34,35] Let p ij (i, j ∈ N) be the PFNs and P be the Pythagorean fuzzy preference relation defined on the alternative set X, which is characterized by Generally, ij represents the degree of preference of alternative X i relative to X j , and ij represents the degree of preference of alternative X j relative to X i . The indeterminate d e g r e e ij =

Method for deriving priority weight vectors with PFPRs
Multiplicative consistency and properties of PFPRs for all js, is called an order consistent PFPR.
Example 1 Let X = x 1 , x 2 , … , x n be a set of alternatives. A PFPR P = (p ij ) 4×4 , i, j ∈ N is as follows: From the above-mentioned PFPR P, it is apparent that p 1j ≥ p 2j ≥ p 3j ≥ p 4j , j ∈ {1, 2, 3, 4} , which denotes that alternative x 4 is inferior to x 3 , alternative x 3 is inferior to x 2 , and alternative x 2 is inferior to x 1 . Therefore, according to the ordered consistency of PFPRs, the rank of the alternative sets is Definition 8 Suppose that the PFPR P = (p ij ) n×n , i, j ∈ N is regarded as multiplicatively consistent, which satisfies the multiplicative transitivity For ij = ji , ij = ji , i, j = 1, 2, … , n in Definition 4, Eq. (6) is equivalent to the following: Based on Definition 3.2, we introduce the following theorem.

Proof
1. We first prove the necessity of Theorem 1.
Obviously, by Definition 3.2, we know that ij , i, j, k = 1, 2, … , n can be obtained. Furthermore, we have the equation H e n c e , f ro m E q . ( 9 ) , i t fo l l ows t h a t (p ij ) = (p ik ) ⋅ (p kj ), i, j, k = 1, 2, … , n. 2. Here, we continue to prove the sufficiency of Theorem 1.
S i n c e ij = ji , ij = ji , i, j = 1, 2, … , n a n d ⟩ n×n is regarded to be standardized, which meets the following conditions: Based on the above-mentioned normalized Pythagorean fuzzy weight vector w = w 1 , w 2 , … , w n T , a multiplicatively consistent PFPR G p = g ij n×n is constructed.
Theorem 2 A PFPR G p = g ij n×n is multiplicative consistent PFPR, and it is expressed as follows: 1. First, we prove that G p = g ij n×n is a PFPR.
Simultaneously, it can be deduced that the following equation holds: Therefore, G p = g ij n×n is a PFPR. 2. Second, we prove that PFPR G p is a multiplicatively consistent.
As for all i, j, k = 1, 2, … , n , the following equation holds: Therefore, PFPR G p is a multiplicatively consistent PFPR in light of Theorem 1.
In summary, this completes the proof of Theorem 2. ◻ Corollary 1 Assume that P = p ij n×n = ⟨ ij , ij ⟩ n×n is a PFPR. If there is a standardized Pythagorean fuzzy weight vector w = w 1 , w 2 , … , w n T with the following condition:

Models for constructing a multiplicative consistent PFPR
In this complex world, when handling decision-making problems with certain complexity and uncertainty, it is unimaginably difficult for decision-makers to establish PFPRs with multiplicative consistency. Therefore, Eq. (3.10) cannot be established, namely, at least one of the following cases exists: , one can obtain Hence, one can see that Next, we introduce two non-negative deviation variables d − ij and d + ij into the equation, which are expressed as follows

3
Obviously, the smaller the two non-negative deviation variables d − ij , d + ij are, the better the multiplicative consistency. Therefore, the optimization model can be constructed to generate the standardized Pythagorean fuzzy weight vector, which minimizes the sum of two non-negative deviation variables we introduce two nonnegative deviation variables e − ij and e + ij into the equation, which is expressed as follows: where e − ij ⋅ e + ij = 0. i, j = 1, 2, … , n. Similarly, the smaller the two non-negative deviation variables e − ij and e + ij are, the better the multiplicative consistency. Therefore, we can establish an optimization model and minimize the sum of the two non-negative deviation variables to obtain the normalized Pythagorean fuzzy weight vector From Definition 4, we know ij = ji , ij = ji . Then, Consequently, M-4.3 can be transformed into the model of M-4.4 by considering the upper triangular elements of PFPRs, which are expressed as follows: , and i, j = 1, 2, … , n.

Assuming that PFPR
In addition, according to Eq. (4.9), one can obtain 2. Support that ̃= 0 and i, j ∈ N , we can derive Through the above proof, we can see that there is a normalized Pythagorean fuzzy weight vector w = w 1 , w 2 , … , w n T , and the elements in the PFPR can be expressed by Eq. (14). According to Corollary 1, the PFPR P = p ij n×n = ⟨ ij , ij ⟩ n×n is multiplicatively consistent.
Therefore, the proof is completed. ◻ One can obtain that the optimal Pythagorean fuzzy weight vectors w = w 1 ,w 2 , … ,w n ⟩ n×n is multiplicatively consistent. Otherwise, we can improve the consistency of PFPRs using the following methods: n×n is a multiplicatively consistent PFPR, in which the elements satisfy Eq. (4.10).
Proof Based on the content mentioned above, we can know Furthermore, an algorithm that obtains the multiplicatively consistent PFPR by aggregating the optimal weight values is proposed, as shown below. (A-3.1)

S t e p 1 : C o n s t r u c t t h e P F P R m a t r i x
⟩ n×n from decision-making information given by the DM.
Step 2: Determine whether the PFPR matrix P is ordered consistent with Definition 7. If yes, the ranking result will be obtained directly. Otherwise, please proceed to the next step 3.
Step 3: The optimal standardized Pythagorean fuzzy weight vectors and deviation variables are obtained by solving the optimal model in connection with the decision-making problem.
Step 4: Calculate the score functions s(w i ) and accuracy functions h(w i ) by the following equations: Step 5: Rank the standardized Pythagorean fuzzy weight vectors w i according to the comparison method mentioned above in Definition 8.
Step 6: Determine the optimum alternative in terms of ranking results.
Step 7: Form a PFPR P = p ij n×n = ⟨̂i j ,̂i j ⟩ n×n with complete multiplicative consistency based on the improved equations.

Multiplicative consistency improving algorithm
In general, it is impractical to derive the conditions for complete multiplicatively consistent PFPRs. Therefore, the consistency index (CI) will be introduced to measure the individual consistency of PFPRs.
n×n is a multiplicatively consistent PFPR. The consistency index (CI) is denoted by Obviously, the smaller CI(A) is, the higher the degree of consistency of A. A is fully multiplicatively consistent when the condition CI(A) = 0 is satisfied.

D e f i n i t i o n 1 1 S u p p o s e t h a t P F P R
n×n and CI is the given consistency index threshold. A is called acceptably multiplicatively consistent PFPR, which is expressed as CI(A) ≤ CI.
To better improve the consistency of the PFPRs given by experts, the consistency improving algorithm (A-3.2) is constructed below. Through the effective improvement of the algorithm (A-3.2), the adjusted PFPRs can reach the preset consistency index threshold, increasing the reliability of the decision result.
Step 1: Let A l = (a l ij ) n×n = A = (a ij ) n×n , where l = 0.
Step 2: The optimal standardized Pythagorean fuzzy weight vectors ŵ l = ŵ l 1 ,ŵ l 2 , … ,ŵ l n T and deviation variables d l− ij ,d l+ ij ,ê l− ij ,ê l+ ij are obtained by solving the optimal model in connection with the decision-making problem.
Step 3: Construct a multiplicatively consistent PFPR matrix If

CI(A)
≤ CI , which indicates that A l = (a l ij ) n×n is the acceptable multiplicative consistency, continue to step 6. Otherwise, go to step 5. . Then, return to step 2.

S t e p 5 : L e t t h e i t e r a t i v e P F P R
Step 6: Let the acceptable multiplicatively consistency A = A l , whose consistency index is CI(A) and w * =ŵ l .
Step 7: The weights are sorted by calculating the score function s � w * i and accuracy function h � w * i , and then, the best alternative x * is selected.
Step 8: End. Proof From the content mentioned above, we can know Let PFPR Â (l) be multiplicatively consistent. With each iteration, we know Here, d l− ij ,d l+ ij ,ê l− ij ,ê l+ ij are the deviation variables and the optimal standardized Pythagorean fuzzy weight vectors
In conclusion, the proof of Theorem 5 is completed. ◻

Consensus reaching process with PFPRs
The model construction in the previous section only focuses on an individual PFPR. However, in actual decision-making, DMs usually give multiple Pythagorean fuzzy preference relations. To address this, we will conduct further research. It should be noted that the smaller the value of ICD(A k ) is, the higher the consensus of A k with other PFPRs, which indicates the higher the consensus between a decision-maker and other decision-makers.

Definition 13 Let
be a set of PFPR on the set X = x 1 , called the degree of group consensus of PFPRs, which is expressed as Similarly, the smaller GCD(A k ) is, the higher the group consensus of all PFPRs A 1 , A 2 , … , A m . That is, the higher the consensus of all decision-makers.

Definition 14 Let
be a set of PFPR given by DMs e k (k ∈ m) on the set X = x 1 , x 2 , … , x n . A c is denoted as a comprehensive Pythagorean fuzzy preference relation, which is expressed as where = 1 , 2 , … , m T is the corresponding weight vector of e k satisfying k ≥ 0, k ∈ M , and ∑ m k=1 k = 1.
The following lemma is introduced to facilitate the later discussion.

3
To summarize, A c is a PFPR, which completes the proof. ◻ However, in actual group decision-making, the conclusion of Theorem 6 has difficultly holding because of the different preferences of decision-makers. Therefore, we give further discussions on it.
Suppose that Corollary 1 is not valid, namely or Next, we introduce some non-negative deviation variables d k− ij , d k+ ij , e k− ij , e k+ ij into the equation, which are expressed as follows: Obviously, the smaller the non-negative deviation variables d k− ij , d k+ ij , e k− ij , e k+ ij , i, j ∈ N, k ∈ M are, the better the multiplicative consistency. Therefore, the optimization model can be constructed to generate the standardized Pythagorean fuzzy weight vector, which minimizes the sum of two non-negative deviation variables

3
Therefore, one can be obtained that the optimal Pythagor e a n f u z z y w e i g h t v e c t o r s

Theorem 7 Let
be a PFPR defined on alternative set X = x 1 , x 2 , … , x n , a n d After deforming Eq. 4.8, we can obtain Considering the above equation and CI(A c ) ≤ CI , the following process is derived: Overall, Theorem 7 holds, which completes the proof. ◻ Inference 1 Let A k = (a k ij ) n×n = ((a k ij, , a k ij, )) n×n (k ∈ m) be a PFPR defined on alternative set X = x 1 , x 2 , … , x n , and where CI is the given consistency index threshold. ) be a PFPR defined on alternative set X = x 1 , x 2 , … , x n , and

Inference 2 Let
can be obtained, which shows that if all PFPRs A k are multiplicative consistent, then the comprehensive PFPR A c is also multiplicatively consistent.

Consensus reaching algorithm
To improve the consensus of all decision-makers in the decision-making problem, this section designs an algorithm to improve the consensus level, which is realized by filtering the preference relation with the worst consensus and adjusting it to reach the consensus index threshold. (A-4.1) Input: A set of PFPRs A 1 , A 2 , … , A m , a group consensus index threshold GCD , and parameter (0 < < 1).
Output: Adjusted PFPR A k * (k ∈ M) , the degree of group consensus GCD * , and iterations l.
Step 2: The degrees of individual consensus ICD (l) and degrees of group consensus GCD (l) of PFPRs A k(l) (k ∈ M) are calculated by Eq. 4.1 and Eq. 4.2, respectively.
Step 3: Check for a consensus for PFPRs. If GCD (l) ≤ GCD , we continue to step 5. Otherwise, go to step 4.
Step 4: Focus on the PFPR with the worst consensus and adjust it. Assume that ICD A (l) = max k ICD A k(l) , and then , a k(l+1) ij, n×n , which satisfies Then, let l = l + 1 , return to step 2.
Step 5: Let the iterative PFPR A k * = A k(l) (k ∈ M) and GCD * = GCD (l) . The output adjusted PFPR A k * (k ∈ M) , the degree of group consensus GCD * , and iterations l.

Algorithm for group decision-making with PFPRs
This section will integrate the contents of the third and fourth sections to build a decision support model based on consistency and consensus, which can ensure that all PFPRs reach the preset consensus level under the premise of satisfying the acceptable multiplicative consistency. In addition, the model is used to calculate the standardized weight vector  Output: The best alternative x ▵ .

Stage 1: Check order consistency
If the PFPRs satisfy the constraints in Definition 7, then they are order consistent and the ranking results can be directly output. Otherwise, go to the next stage.

Stage 2: Consistency control process
The improved consistency algorithm (A-3.2) is used to adjust and obtain the PFPRs A 1 , A 2 , … , A m with acceptable multiplicative consistency.

Stage 3: Consensus building process
The adjusted PFPRs A 1 * , A 2 * , … , A m * of acceptable multiplicative consistency are calculated by the consensus improving algorithm (A-4.1), so that the degree of group consensus reaches the consensus index threshold.

Stage 4: Alternatives selection process
Equation (30) is used to aggregate the PFPRs A 1▵ , A 2▵ , … , A m▵ into comprehensive PFPR A c▵ . According to the ranking of the standardized weight vector calculated by the optimization model (M-4.2), the alternatives are sorted to select the best one x ▵ (Fig. 1)

Numerical example and comparison analysis
In this section, we use a numerical example to illustrate the proposed GDM method and compare the experimental results of the proposed GDM method with those of the GDM methods.

Numerical example
In recent years, the rapid development of the economy has brought opportunities to various industries, but has also brought many challenges. For example, increasingly more enterprises are facing financial difficulties, which makes it important to establish an effective financial distress early warning model to prevent financial crises and risk management.
To achieve a better foothold and long-term development, a listed company establishes a perfect financial early warning system. Through comprehensive analysis, it is found that the main indicators that affect the financial situation of the enterprise are the following: solvency x 1 , operations x 2 , profitability x 3 , and development x 4 . To explore the importance of these four indicators to the financial early warning system, the listed company now invites three experts in related fields to evaluate and analyze the above four impact indicators systematically. The weight vector of the experts is = (0.3, 0.4, 0.3) T . Besides, to simplify the decision-making process, preference relations are useful tools that only require the DMs to judge two alternatives at a time. Furthermore, PFVs are flexible in expressing the DMs' hesitancy, where there is more than one value for judgments. Considering these two aspects, PFPRs are used for this GDM problem.
Therefore, experts use PFNs to evaluate the above four impact indicators according to professional knowledge and financial experience. The final PFPRs are as follows: Some preset parameters are given as follows: consistency index threshold level CI = 0.1 , group consistency index threshold level GCD = 0.1 , consistency adjustment parameter = 0.2 , and consistency adjustment parameter = 0.3 . The detailed stages on the latter are shown in "Appendix".
As international competition intensifies and the market environment changes, enterprises will inevitably suffer from various financial crises. Establishing a financial early warning system is conducive to preventing and overcoming financial crises. To guard against financial risks, the listed company prefers to establish an effective financial early warning system based on several indicators affecting its financial situation.
From the stages in "Appendix", we obtain s x 1 = −0.0747, s x 2 = −0.4753, s x 3 = −0.6741, and s x 4 = −0.7762 via Definition 5. Then, the normalized Pythagorean fuzzy weight vector is ranked w 1 >w 2 >w 3 >w 4 . Therefore, the rank results of importance of the influencing factors are solvency > operations > profitability > development (x 1 > x 2 > x 3 > x 4 ) , indicating that the most important interfering factor is solvency x 1 . Moreover, the comprehensive PFPR A c▵ also shows that the preferred factor x 1 is higher than the other three factors by the three experts x j (j = 2, 3, 4) , such a s a c▵ 12 = (0.7299, 0.5672),a c▵ 13 = (0.7784, 0.5129) a n d a c▵ 14 = (0.7762, 0.4833) . Furthermore, most of the raw information in the PFPRs given by the three experts is robust enough to reflect that the preference value of x 1 is superior to those of the others. E.g., a 2 12 = (0.75, 0.55) represents the degree of preference of solvency x 1 relative to operations x 2 . Obviously, this is consistent with the results obtained in this paper, which side by side proves that the method of this paper is robust and reasonable.
Therefore, under the analysis and judgment of the three experts in the related field, it is concluded that solvency, operations, profitability, and development are all important indicators reflecting financial status; however, to avoid risks effectively, listed companies should give primary consideration to their own solvency when establishing a financial early warning system.
Generally, solvency is the ability to pay debts as they fall due. If there is a problem with solvency, it will definitely pose a threat to the entire financial situation of an enterprise, and may even cause the enterprise to experience a financial crisis. Therefore, in the fierce market competition, the first goal of an enterprise to survive is to cover its expenses with revenues and pay off its debts on time.

Comparison analysis
The following is the comparison and analysis of methods of [18-20, 34, 35, 39], which includes the following main aspects.
1. To highlight the excellent properties of multiplicative consistency, we will compare multiplicative consistency with the method of additive consistency. Therefore, the methods proposed by Chu et al. [18] and Yang et al. [35] are used to find the optimal solution to the above problem. 2. To determine whether the decision-making information of the Pythagorean fuzzy environment is more comprehensive than that of the intuitionistic fuzzy environment, we will compare the proposal with the IFPR methods proposed by Chu et al. [18], Jin et al. [19] and Lin and Wang [20]. 3. To accentuate the consensus in GDM, we will compare the proposal with the individual decision-making approach by He et al. [39] and the approach that consid-   [19], Mandal and Ranadive [34], and Lin and Wang [20].
Considering the relationship between PFPRs and IFPRs, we will compare the methods [18][19][20] under an intuitionistic fuzzy environment to highlight the advantages of PFPRs. First, according to the transformation between them, the following IFPRs are obtained: Considering that the method proposed by He et al. [39] solves an individual decision-making problem, the original PFPRs need to be converted to an overall decision-making matrix. The overall PFPR A ⊙ can be obtained by utilizing the method in [37]. We need to convert the above PFPRs A k = a k ij 4×4 (k = 1, 2, 3) of these four factors for the financial early warning system into the following overall PFPRs A ⊙ : Based on the example of financial crisis and risk management, the calculation results of the above five methods are shown in Table 1. To make the analysis not too redundant and to better highlight the results, the specific calculation procedure is omitted.
First, compared with the method proposed by He et al. [39], we find that the results obtained by the method in [39] are different from those obtained by our proposed method. The method in [39] yields a ranking result of x He 2 > x He 1 > x He 4 > x He 3 , while our method produces a ranking result of x 1 > x 2 > x 3 > x 4 . However, from the overall PFPR matrix A ⊙ , we have A ⊙ 12 = (0.5250, 0.4500) , which means the degree of certainty that solvency x 1 is preferred to operations x 2 is 0.525, and the degree of certainty that solvency x 1 is non-preferred to operations x 2 is 0.45. Obviously, solvency x 1 cannot be inferior to operation x 2 . Similarly, profitability x 3 should not be inferior to development x 4 by A ⊙ 34 = (0.4750, 0.4500) . Therefore, our method is more robust than that of [39].
Second, from the four indicator rankings obtained by our method with those of Chu et al. [18], Jin et al. [19], Lin and Wang [20], Yang et al. [35], and Mandal and Ranadive [34], we can see that the ranking results are consistent, i.e., solvency > operations > profitability > development , but with some differences in approaches.
1. Compared with the methods proposed by Chu et al. [18], Yang et al. [35] found the optimal solution to the above problem. We know that when additive consistency is used to describe the incomplete information in FPRs, an overflow range of 0 to 1 may be generated. Compared with additive consistency [18,35], multiplicative consistency does not have information spillover in FPRs, or PFPRs. Therefore, the advantages of our decision support model are obvious. In response to PFs that do not achieve acceptable consistency, Chu et al. [18] built non-preferred relations for adjustment. This method can consider the additive consistency of the nonpreferred relations, but the process is complicated and information is easily lost. Our decision support model is directly based on the decision-making information given by experts, making full use of the original decision-making information for decision-making, which can reduce the information loss in the decision-making process and increase the reliability of the final results. By comparison, our method is simpler and more convenient. 2. As above, the decision results by Chu et al. [18], Jin et al. [19] and Lin and Wang [20] are also consistent with our final result, which again proves the excellent performance of the decision support model proposed in this paper. Although Jin et al. [19] constructed an improved model of the multiplicative consistency of IFPRs, our integrated environment is based on PFPRs, which is more extensive and sweeping than IFPRs in describing decision information. Moreover, the group decision model based on multiple IFPRs is just a simple extension of the decision model based on individual IFPRs. The model does not consider the consistency and compatibility of IFPRs, but our group decision model is the opposite. In our decision support model, we not only discuss how to improve the consistency of PFPRs, but also consider how to enhance the compatibility between them. Therefore, the virtues of this method are self-evident. 3. We compare the proposed method with the individual decision-making approach by He et al. [39] and the approach that considers only consistency proposed by Jin et al. [19], Mandal and Ranadive [34], and Lin and Wang [20]. He et al. [39] focused on multiplicatively consistent PFPRs, which is similar to this paper. Their method only considered the individual decision problem, while our method discusses the GDM consensus problem. Although Jin et al. [19] considered GDM, GDM was only an aggregation of individual decisions and did not consider the problem of non-consensus opinions among experts in GDM. As for Mandal and Ranadive [34], a consensus relationship was constructed. However, there is a problem with the definition of the relevant equation, which makes the value on the diagonal of the constructed consensus relation not equal to one. Thus, our method is more practical.
Through the introduction of the above different methods of [18-20, 34, 35], it is found that the ranking results of the impact of factors are consistent (see Table 1), which shows that the decision support model constructed in this paper is reasonable and effective. Furthermore, the decision support model proposed in this paper integrates a variety of algorithm models such as a consistency improving algorithm and a consensus building algorithm. Therefore, our research methods consider decision information from multiple perspectives and solve decision-making problems at different levels. Besides, there is no information overflow in multiplicative consistency, which is more suitable for PFPRs than additive consistency. However, the degree of spatial membership of intuitionistic fuzzy environment is less than that of Pythagorean fuzzy environment, which is equivalent to the IFNs must be the PFNs. Otherwise, this does not hold. Therefore, the PFPRs can solve the problem that IFPRs are incompetent. In conclusion, it is valuable to build a decision support model with PFPRs based on multiplicative consistency.
To make a comparison, we count the Pearson correlations of rankings similarity via SPSS, as shown in Table 2.
According to the data of correlations of different methods with the same ranking, there are some correlations among these methods, and the difference is only the strength of the correlations. The correlation coefficients of our method and methods in [19] and [20] are less than 0.05, and are significantly correlated. In particular, the correlation coefficient between our method and the method in [34] is 0.006 < 1, indicating a strong correlation between them.

Conclusions
Considering that different fuzzy environments, consistency improvement, consensus reaching, and the method of ranking weights are the key research topics, this paper focuses on the preference relations in a Pythagorean fuzzy environment.
This paper provides several significant contributions to group decision-making with PFPRs as follows: 1. A novel concept of multiplicatively consistent PFPRs is introduced. The advantage of the proposed multiplicative consistency is that it can overcome the malpractice of additive consistency with PFPRs.

A consistency measure index of PFPRs is introduced.
Moreover, when the CI does not reach the expected threshold, a consistency improving model is designed to achieve the expected threshold. 3. The ICD and GCD are developed to measure the consensus level. If the GCD is not within the expected threshold, the consensus reaching model is modified to reach the expected threshold. 4. For the individual decision-making problem, order consistency is defined to evaluate it efficiently. Furthermore, a decision support model, which is based on a consistency improving model and consensus reaching model, is constructed for the GDM problem.
However, we do not conduct deep discussions on how the iteration parameters will affect the iterative speed and final decision results, which will be considered in the future. In addition, the weight of the experts in this paper is fixed, without considering the complexity of practical problems. The acquisition of incomplete weight information will be included in the following research. 1 3 whose consistent indexes are CI A 1� = 0.0911 < CI and CI A 3� = 0.1011 > CI , respectively. Therefore, A 3′ is adjusted to meet the condition that its consistency index is less than 0.1. Then, let A 3� = A 3�� . One can obtain with CI A 3�� = 0.0786 < CI . Through this previous process of improving the consistency, all PFPRs are acceptably multiplicatively consistent. Then, let the PFPRs A 1 * = A 1� , A 2 * = A 2 and A 3 * = A 3�� , and proceed to the next stage 3.

Stage 3: Consensus building process
Based on the adjusted PFPRs A 1 * , A 2 * , … , A m * with acceptable multiplicative consistency, we will check their consensus.
First, let l = 0, and A k * (l) = A k * = a k * ij 4×4 , k = 1, 2, 3 . The degree of individual consensus ICD and the degree of group consensus GCD with PFPRs are calculated as follows: