Abstract

Linguistic q-rung orthopair fuzzy numbers (Lq-ROFNs) are composed of a set of q-rung orthopair fuzzy numbers (q-ROFNs) with membership and nonmembership degrees of linguistic variables, which can be regarded as an extension of linguistic intuitionistic fuzzy numbers (LIFNs) and linguistic Pythagorean fuzzy numbers (LPFNs). Compared with LIFNs and LPFNs, Lq-ROFNs have a wider range of applicability because the value of q can have a wider range, thus allowing to handle more multiattribute group decision making (MAGDM) problems. Based on Lq-ROFS, this paper first proposes the Chebyshev distance metric, then develops the Chebyshev distance entropy for deriving objective solution of decision makers’ (DMs’) weight vector and attribute weight vector by combining the metric and entropy measure. The TODIM method has been widely valued by the society and has achieved good results in many MAGDM problems. Subsequently, the TODIM decision method in Lq-ROFS is presented and combined with the Chebyshev distance entropy model as a new decision method (CDE-TODIM) to solve the MAGDM problems. This method process is objective and direct and takes into account the decision maker’s preferences, and the decision results are more persuasive. Finally, the effectiveness and rationality of this MAGDM method are illustrated by a case as well as comparisons.

1. Introduction

In today’s human social activities, decision-making is becoming increasingly vital, so how to seek an optimal solution is the key to the problem [14]. When choosing the optimal solution, the choice is analyzed by collecting the evaluation information of a group of DMs, in which the decision making behaviors involving multiple criteria become multi-criteria group decision-making (MCGDM) [58] and those involving multiple attributes become MAGDM, and the problem studied in this paper focuses on MAGDM [911]. In the beginning, to express their subjective decisions, DMs often describe them by real numbers [12]. The size of real numbers represents the intensity of subjective will, but this means that the evaluation information is not so crisp. Therefore, to solve the situation mentioned above in the MAGDM problem, Zadeh [13] first proposed the concept of the fuzzy set (FS), which only considers the membership function. Once this method was proposed, it aroused widespread discussion in society [14, 15]. However, considering the membership function alone means that only membership degree (MD) information is included in the evaluation information so it will cause the loss of nonmembership degree (NMD) information, and it is not easy to guarantee the accuracy of the final decision result. Aiming at the fuzzy set (FS) proposed by Zadeh, Atanassov made further improvements to its deficiencies. Based on this, Atanassov [16] proposed the intuitionistic fuzzy set (IFS) in 1986; that is, the degree of membership and the degree of nonmembership are included in the evaluation information at the same time. The IFS requires two degrees cannot exceed one, that is, .

However, as the MAGDM problems are becoming more and more complex, this constraint is becoming too strong for decision-makers (DMs). For example, under the evaluation criteria, a DM’s preference for a certain alternative is recorded as 0.8 and NMD is recorded as 0.5, obviously, 0.8 + 0.5 > 1; so it cannot meet the conditions of IFS proposed by Atanassov. Thus, based on the predecessors, Yager and Abbasov [17] introduced Pythagorean fuzzy sets (PFSs), and the corresponding constraint conditions are relaxed compared with IFS and only need to meet . IFS and PFS have solved many MAGDM problems, respectively, have been recognized in the society and have also made some achievements in the academic field. In recent years, Yager [18] has further expanded fuzzy sets and proposed the concept of q-rung orthopair fuzzy sets (q-ROFSs); that is, MD and NMD only need to satisfy the following formula: . In fact, it is easy to see that q-ROFSs are much looser than PFSs, also can be better applied to more uncertain MAGDM problems than IFSs and PFSs. And, people found that if q = 1 in q-ROFS, q-ROFS becomes IFS, if q= 2, it becomes PFS. From here, it can be understood that q-ROFS is a kind of promotion of IFS or PFS [1922].

For different research objects, such as information aggregation and decision problems under uncertainty, Ullah [23] propose picture fuzzy set (PFS’), so as to study the MAGDM problem in this set, and Mahmood et al. [24] make further extensions based on PFS’ by proposing spherical fuzzy set (SFS) and T spherical fuzzy set (T-SFS) concept, which makes the decision process more explicit.

In the methods like IFSs, PFSs, and q-ROFSs, MD and NMD given by DMs are expressed in the form of real numbers, but for DMs themselves, this group tends to use language to describe and express vague evaluation information. This is actually more convenient [25]. For example, when people select outstanding students to issue scholarships, the linguistic evaluation system is utilized to evaluate students’ ability such as “very good,” “ordinary,” and “poor”. In the early days, Zadeh [26] established the linguistic term set (LTS) and proposed some linguistic calculation models to model the evaluation information, thereby improving the accuracy [2729]. However, these linguistic calculation models only contain MD, without NMD. Based on this, Zhang [30] combined the established IFS and LTS to create a new concept named the linguistic intuitionistic fuzzy set (LIFS), which elements are called the linguistic intuitionistic fuzzy numbers (LIFNs). Each element in LIFSs, named LIFN, is composed of a pair of linguistic terms, including MD and NMD. And, it needs to meet the condition: the sum of the subscripts of MD and NMD must less than the cardinal number of LTS. The development of linguistic fuzzy sets is similar to the fuzzy sets. Garg [31] proposed a novel linguistic fuzzy set called linguistic Pythagorean fuzzy set (LPFS), and the elements in the set are called LPFNs. The analogy LIFN describes MD and NMD through a pair of linguistic terms, and LPFN is also the same. But the difference is that LPFN only needs to satisfy that the sum of MD and NMD subscripts less than the square of the cardinality number of the LTS. For easily understanding, we assumed that continuous LTS (CLTS) as . Thus, for LIFN, its related conditions are . And, as comparison, condition that LPFNs needs to meet . From here, people have concluded that all LIFNs will certainly meet the conditions of LPFNs, but the converse is not necessarily true. The conclusion comes naturally, that LPFNs cover a larger range than LIFNs.

In order to further expand the coverage of linguistic evaluation information for decision makers (DMs), Liu and Liu [32] proposed the linguistic q-rung orthopair fuzzy numbers (Lq-ROFNs). And condition is , . Over time, more and more MAGDM problems have been built on Lq-ROFS, and thus more and more authors have proposed decision methods for Lq-ROFNs [3349]. If the evaluation of DM is inscribed by selecting real numbers directly from the number field, it is not precise enough and does not aggregate all the information provided by DM, while Lq-ROFNs can precisely describe the linguistic information of DM, so this is the main reason why we choose to solve the MAGDM problem in Lq-ROFS in this paper.

It has been shown in the previous section that people collect evaluation information by using Lq-ROFNs, the information will be more accurate and the decision results obtained will be more convincing and valid. The TODIM method has not been used in Lq-ROFS yet, therefore, the decision making method proposed in this paper is mainly implemented in Lq-ROFS. The following section will illustrate the superiority of choosing the TODIM method.

Of course, in addition to the TODIM method, there is no shortage of methods to solve the MAGDM problem, such as VIKOR [50, 51], TOPSIS [52], ELECTRE [53, 54], and DEMATEL [55] (some of them have not yet been extended to Lq-ROFS, which will also be considered for future work). These methods have their suitable applications in different fields and different directions. However, it is a pity that all the above-given methods have a common shortcoming. That is, these methods do not consider the orientation of the DMs’ subjective to strengthen the decision-making results. For the time being, it is called the bounded rationality of the DMs. Gomes and Lima [56] considered the bounded rationality of DMs and proposed the TODIM method. This is the first time that the bounded rationality of DMs has been considered in the MAGDM problem, which is an unprecedented breakthrough [57]. TODIM method constructs the dominance matrix of pairwise comparison of schemes according to the different preferences of DMs for benefits and costs and then aggregates the individual dominance to form the dominance of each scheme. The ranking of schemes can be obtained by the size of comparative dominance. At present, this method has been widely valued by the society, and has achieved good results in many MAGDM problems. For example, Zhao et al. [58] applied the TODIM method in scientific and technological assessment in risk management, He et al. [59] proposed a shadowed set-based TODIM method, and applied it to large-scale group decision making, Wu et al. [60] developed a linguistic distribution behavioral MCGDM model integrating extended generalized TODIM method and quantum decision theory. Decision making using the TODIM method involves solving for different attribute weights. If the weights are unknown, the subjective assignment of weights tends to create a preference for a particular scheme, thus making it impossible to objectively evaluate the decision scheme, so an objective method of solving for the weights is considered. In the paper, Liu and Liu [32] propose the Hamming distance and combine it with the entropy method to solve the attribute weights. The superiority of Chebyshev distance will be discussed and the distance will be extended to distance entropy to solve for expert weights and attribute weights.

In this paper, we propose a new distance measure based on the Hamming distance in Lq-ROFS, which is the Chebyshev distance. It has similar properties to the Hamming distance. However, it is easier to define and considers fewer factors than the Hamming distance. For example, it considers only one between membership, nonmembership, and uncertainty, and it is more able to bring out the subjective tendency of DMs, which is objective mixed with subjective. And, to some extent, it is more effective than the metrics established by considering only subjective or objective. Furthermore, in most papers proposing methods for solving the MAGDM problem in Lq-ROFS, the authors only explain the steps of the method but do not indicate how to solve for the DMs’ weights and attribute weights. For example, Lin et al. [61] propose a new decision operator when the DM and attribute weights are known. Therefore, in this paper, we propose the Chebyshev distance entropy model for solving the DMs’ weights and attribute weights by using the definition of Chebyshev distance measure and entropy, which can explain the accuracy of the decision results and the effectiveness of the method more objectively.

The benefits and weaknesses of existing methods are listed in Table 1.

From the comparison in Table 1, it can be concluded that the innovation of this paper mainly focuses on extending the TODIM method to Lq-ROFS, defining the concept of Chebyshev distance entropy, proposing an objective method for solving expert weights and attribute weights, and finally combining the two to solve the MAGDM problem. And, we can also know that the traditional model for solving the MCDM problem has its own shortcomings, MABAC, RAFSI, and other methods are applicable to different environments, while for DMs’ linguistic evaluation to collect information, Lq-ROFS is more accurate.

Therefore, the advantages of the CDE-TODIM method are as follows:The first point is that the Chebyshev distance only considers the largest one of MD, NMD, and uncertainty, which can better express the preference of DM compared to the Hamming distance which directly takes the average of the three.The second point is to expand the Chebyshev distance to Chebyshev distance entropy, so as to solve the weight of DM and the attribute weight. The overall process is objective and direct, and the results are more convincing than those obtained by direct or subjective assignment.The third point is the extension of the TODIM method to Lq-ROFS, which is an unprecedented work that retains the advantageous nature of the TODIM method itself, i.e., considering the bounded rationality of DM and uses Lq-ROFNs for decision making, which makes the linguistic expressions of DMs transformed into real numbers and thus solves the MAGDM problem objectively.

For ease of reading, Table 2 will explain the abbreviations in the paper.

The rest of this paper is arranged as follows. In Section 2, some preliminary knowledge is introduced, including the concepts of q-ROFS, LTS, Lq-ROFS, and LqROFWA; in Section 3, how to use the Chebyshev distance entropy measure model to solve for DMs’ and attribute weight vector is proposed; in Section 4, we combine the Chebyshev distance entropy model with the TODIM method for solving the MAGDM problem in Lq-ROFS; in Section 5, a case is given to illustrate the effectiveness of the proposed method. Then, the paper solves the case using different methods and the results obtained are compared and analyzed. At the end of this paper, the practicability of CDE-TODIM method will be explained.

2. Preliminaries

The following concepts will be briefly introduced: q-ROFS, LTS, Lq-ROFS, and LqROFWA.

2.1. Q-Rung Orthopair Fuzzy Set (Q-ROFS)

Definition 1. (Liu and Wang [20]). Assuming that the given q-ROFS T is in finite universe , it can be expressed aswhere , respectively, represent the membership degree and the nonmembership degree of element belonging to the q-ROFS T. The corresponding constraints are as follows: for each , it holds the following results:Due to , the indeterminacy degree of q-ROFS can be expressed as , . Thereby, q-ROFN corresponds to , and it can be recorded as .

2.2. Linguistic Term Set (LTS)

Definition 2. (Herrera et al. [72]). LTS is a finite-ordered discrete with odd cardinality set, which is recorded as . The elements in A, like , can be defined according to different semantic environments, but in any case, the following conditions need to be met:(1)Ordered: if , it means , and is worse than (2)Negative operator: neg () = (3)Min operator: min () = , if i < j, k = i, else k = j(4)Max operator: max () = , if i < j, k = j, else k = iDue to the discrete language set cannot well include all experts decision information, Xu [73] extends it to a continuous set, which is recorded as CLTS (l is a positive integer), and its conditions are as same as LTS.

2.3. Linguistic Q-Rung Orthopair Fuzzy Set (Lq-ROFS)

Definition 3. (Liu and Liu [32]). A Lq-ROFN in K is recorded aswhere stand for linguistic MD and linguistic NMD, respectively. The condition holds for any . For this q-ROFS, q-ROFN is , and the linguistic indeterminacy degree of is .
The complete Lq-ROFNs on the basis of just recorded as .

Definition 4. (Liu and Liu [74]). Suppose , are two Lq-ROFNs, , p is a constant, and p > 0.

2.4. LqROFWA (Lq-ROFNs Weighted Average Operator)

Definition 5. (Lin et al. [61]). Let be n Lq-ROFNs, the LqROFWA has the following equation:where denotes the weight of that satisfies and .

3. Chebyshev Distance Entropy Measure Model

In this section, the definition of Chebyshev distance between Lq-ROFNs, Lq-ROFV (Linguistic Q-rung Orthopair Fuzzy Vectors), and Lq-ROFM (Linguistic Q-rung Orthopair Fuzzy Matrices) will be introduced, and then, the methods to solve for DMs’ weights and attribute weights by Chebyshev distance entropy model will be presented.

3.1. Chebyshev Distance

Definition 6. The Chebyshev distance between any two Lq-ROFNs , can be defined as follows:

Theorem 1. The Chebyshev distance between any two Lq-ROFNs , satisfies the following properties:(1)(Symmetry) (2)(Non-negativity) (3)(Triangle rule)

Proof. (Symmetry) :(Non-negativity).
Due toSo, .
And . Reasons are as follows:(Triangle rule) :Thus, .Comparing to the Hamming distance presented in paper [74],The Chebyshev distance also operates with the subscripts of Lq-ROFN and converts Lq-ROFN to real number in a similar way, which means that they have similar properties. However, compared to Hamming distance, which averages the sum of subscript differences, Chebyshev distance considers only the one with the largest difference among membership, nonmembership, and uncertainty, and is more able to take into account the subjective tendencies of decision makers.

Definition 7. The Chebyshev distance between any two Lq-ROFVsAnd can be defined as follows:

Theorem 2. still satisfies symmetry, non-negativity, and triangle rule.

Definition 8. The Chebyshev distance between any two Lq-ROFMs:And can be defined as follows:where .

Theorem 3. also satisfies symmetry, non-negativity, and triangle rule.

3.2. Chebyshev Distance Entropy for DMs’ Weights

Entropy measures the uncertainty of information based on probability theory, which reveals the law that the more dispersed the distribution of information, the greater the uncertainty. After being inspired by the definition of the entropy power method based on the linguistic intuition fuzzy number in one paper (2017), the definition of Chebyshev distance entropy will be subsequently proposed in the following steps.

Before presenting the method, first assume that evaluate each scheme based on every attribute , and each appraised value can be recorded as Lq-ROFN variable . Here, we also assumed that the DMs’ weight vector is and the attribute weight vector is .

In this way, we can obtain the Lq-ROFS decision matrix , with , where i= 1,…,n, and j= 1,…,m. Then, the specific steps of Chebyshev distance entropy to solve for DMs’ weight vector are as follows:Step 1: calculate the Chebyshev distance between the decision matrices given by each DM, then the DMs’ distance matrix can be obtained such thatStep 2: calculate the deviation of the from each of the other DMs:Step 3: calculate the entropy value of each to represent the DM’s information:Step 4: calculate the variance of each :Step 5: calculate the weights corresponding to each :

3.3. Chebyshev Distance Entropy for Attribute Weights

The following section will describe how to solve for attribute weights using the Chebyshev distance entropy model:Step 1: calculate the deviation of the scheme from each of the other schemes under attribute :Here, indicates the Chebyshev distance between Lq-ROFNs.Step 2: calculate the entropy value of each attribute to represent the decision information:Step 3: calculate the variance of attribute :Step 4: calculate the weights corresponding to each attribute :

As a way to obtain attribute weights objectively, the weights solved are more representative and can make the decision results more accurate and effective, so it will be a powerful tool for solving MAGDM problems.

The following will introduce how distance entropy will be combined with the TODIM method.

4. Lq-ROFS CDE-TODIM

In this section, based on the intuitionistic linguistic TODIM method, we extend this method to Lq-ROFS. The followings are the specific steps of this method.

We know that attributes can be divided into two types according to their properties, namely, benefit-type attributes and cost-type attributes. Therefore, the first step should be normalization, which is the following:Step 1: the decision matrix given by each DM is normalized to give the matrix .When the elements in the decision matrix K are normalized, there are following conclusions:Step 2: assuming that the DMs’ weights are derived by the Chebyshev distance entropy model, or are given in the case. Then, the DMs’ decision matrices can be assembled by LqROFWA operator; in this way, they will turn into a new decision matrix containing all the wishes of each DM.Step 3: assuming that attribute weights are derived by the Chebyshev distance entropy measure model or given in the case, which are noted as , and each , . Then, calculate the relative weight of each attribute relative to the highest weight .Step 4: calculate the overall dominance of scheme over .wherewhere represents the Chebyshev distance between two Lq-ROFNs, and represents profit or 0, while indicates costs. θ represents the attenuation factor of the losses. When θ> 1, the DMs are risk averters, the higher value of θ, the higher the degree of avoidance; when 0 <θ< 1, DMs are risk appetite. From paper [75], it suggests DM give the value of θ between [1, 2.5]. And θ= 1 or θ= 2.5 are common values.Then, prospect value function of scheme is shown as follows:Step 5: normalize :Step 6: the scheme is sorted according to the numerical size of . The larger , the better scheme is.

In Figure 1, we can see the specific method flow, which is presented in a more macroway.

5. Case Study

In this session, a case of decision making utilizing the CDE-TODIM method under Lq-ROFS will be given, and the results will be profiled and illustrated and then compared with other decision making methods on this basis. As a result, the validity of the CDE-TODIM method will be analyzed from the above-given solution conclusions.

5.1. Case Study Based on the CDE-TODIM Model

The epidemic caused by COVID-19 is a great concern for humans worldwide at this time, and while hoping for a speedy improvement in the epidemic, business opportunities have been identified in home isolation. This case will apply the CDE-TODIM method to decide the best investment option for investors and to illustrate the final results. An investor currently has a sum of money on hand and while isolating his home, he has identified several business opportunities and intends to invest in one of these options, which are as follows: is a pharmaceutical company, is a mask manufacturing plant, is an online medical platform, and is an online educational platform. In order to find the optimal investment solution, the investor asks three decision makers to make decisions and takes into account all the wishes of these three DMs. The DMs considered the following four main evaluation criteria in selecting the schemes: profitability , market prospects , risk , and cost . The first two attributes are benefit-type and last two are cost-type; therefore, these can be divided into two classes and . On the basis of the CLTS , the DMs give the decision result (Lq-ROFR) for each scheme under each attribute .

S = { = extremely poor,  = very poor,  = poor,  = slightly poor,  = general,  = slightly good,  = good,  = very good, and  = extremely good}

The decision making steps using the CDE-TODIM method are as follows:Step 1. The Lq-ROFR given by the DMs is in Tables 35.Step 2. The results of normalizing the Lq-ROFR in Tables 68.Step 3. The Chebyshev distance entropy model is used to find the DMs’ weights as Table 9.Step 4. LqROFWA method of DM’s decision matrices for aggregation (q = 3).Step 5. Applying Chebyshev distance entropy measure method to solve for attribute weights. And calculated attribute weights are shown in Table 10.Step 6. Calculate the relative weight of each attribute relative to the highest weight .Step 7. Calculate the total dominance of scheme over scheme (As a result of the conservative decision analysis, first take ).Step 8. Calculate .Step 9. Normalize , then we can obtain .Step 10. Sort schemes: . Therefore, the optimal scheme is .Step 11. End.

Tables 915 show the intermediate results obtained from the model calculations, which more visually reflect the feasibility of the method.

5.1.1. Sensitivity Analysis

The scores of different schemes when q is taken from 1 to 20 are displayed in the figure.

When is taken from 1 to 10 ( = 3), the figure of the change in scheme scores.

It can be summarized from Figure 2 that when is fixed at 2.5, which means that DMs are risk averse, then for q from 1 to 20 (different Lq-ROFS environments), it can be found that the optimal solution is when , followed by , followed by , and the worst is , but when , the optimal solution becomes , and when q= 16, 17, 18, the optimal solution becomes again, and when q > 18, the optimal solution is . The above mentioned discussion shows that the optimal solution fluctuates between and and is influenced by the value of q.

From Figure 3, when q is taken as 3 and takes different values, it can be seen that when , the scheme ranking is , when , the scheme ranking is , which indicates that the value of also has an effect on the scheme ranking, and the optimal scheme also fluctuates between and .

Therefore, it can be seen from the figure that the scheme ordering changes significantly when the q value changes, while the scheme does not change much when the value changes, which also indicates that the choice of which Lq-ROFS has a greater impact on the scheme.

5.2. Comparison of Effectiveness of Multiple Approaches

To illustrate the effectiveness and feasibility of the CDE-TODIM method, the results were compared with the use of other methods to solve this case, and the comparison is presented in Table 16 to obtain an evaluation of the CDE-TODIM method.

The results of other methods for solving this example show that the difference between the decision results obtained by the CDE-TODIM method and other methods is not significant, which shows that the CDE-TODIM method works.

For this case, the optimal solution of decision making using the CDE-TODIM method is better in investment in mask manufacturing plant and online education platform when the value of q is less than 4, and the education platform solution is better when it is greater than 4. The optimal solution found by the rest of the methods is investment in mask manufacturing plant, which indicates that the decision results of different solutions tend to be consistent and the results are more credible. According to the data, more people in real life choose to invest in the mask manufacturing plant and a small number of people choose to invest in the education platform. This means that the best way to measure the decision results is to find similar cases from real life and compare them with the results, and if the deviation is large, further careful consideration is needed, and if the deviation is small, the optimal result of the decision can be chosen.

6. CDE-TODIM Model Effectiveness Evaluation

After a comparative analysis with different methods, other advantages of the CDE-TODIM method, in addition to the previously proposed advantages, are as follows:(1)Considering the solutions of all DMs and using Chebyshev distance entropy method to solve the DMs’ weights, and then assembling the decision matrix with the help of LqROFWA operator, which makes the results more objective and comprehensive.(2)By applying the extended TODIM method under Lq-ROFS, not only the original advantages of the method are retained, but also the combination of language sets makes the evaluation system more realistic and accurate.(3)Improved Hamming distance to Chebyshev distance for simplicity and ease of use

In general, the CDE-TODIM method is comprehensive and objective in its consideration of aspects that lead to better decision making results, and the method is simple and easy to follow with fewer steps.

7. Conclusions

In the context of the new crown epidemic and during home quarantine, the idea of how to find business opportunities, earn wealth, and subsidize family is on the minds of many people, but the exact aspect to start with is a difficult one. The method provided in this paper is to solve the realistic problem of difficult decisions and select the most suitable solution to achieve the goal.

In this paper, we have developed Chebyshev distance under Lq-ROFS and proposed Chebyshev distance entropy to objectively determine the DMs’ weights and attribute weights based on the above distance and entropy. Then, we have combined the proposed weight method and the TODIM method to solve the MAGDM with Lq-ROFS problem whose decision matrix elements are Lq-ROFN. Moreover, a case study and a comparative analysis have been provided to verify the feasibility of the the CDE-TODIM MAGDM method.

The contribution of this paper is proposing the CDE-TODIM method for solving the investment option selection problem in the context of epidemic.

In the future, CDE-TODIM method will be extended to PFS’ for research, and for Lq-ROFS will consider interval language, probabilistic language, etc. In addition to the investment case covered in this paper, extension applications in other fields such as health care assessment will also be considered.

Data Availability

The data used to support the study are included within the article.

Ethical Approval

This article does not contain any studies with human participants or animals performed by any of the authors.

All authors agreed with the content of the manuscript and the accepted submission and agreed to be accountable for all aspects of the work.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Yangcan Wang conceptualized the study, curated the data, formally analysed the study, carried out the funding acquisition, proposed the methodology, and wrote the original draft; Xin Tang, Yutong Pan, and Mingqian Cai supervised the study; and Yangcan Wang, Xin Tang, Yutong Pan, and Mingqian Cai wrote the review and edited.

Acknowledgments

The work was supported by the National Natural Science Foundation of China (Grant nos. 72171002, 71771001, 71701001, 71871001, 71901001, 71901088, 72071001, and 72001001), Natural Science Foundation for Distinguished Young Scholars of Anhui Province (Grant no. 1908085J03), Research Funding Project of Academic and Technical Leaders and Reverse Candidates in Anhui Province (Grant no. 2018H179), Top Talent Academic Foundation for University Discipline of Anhui Province (Grant no. gxbjZD2020056), Anhui Provincial Natural Science Foundation (Grant no. 1808085QG211), College Excellent Youth Talent Support Program (Grant no. gxyq20119236), Key Research Project of Humanities and Social Sciences in Colleges and Universities of Anhui Province (Grant no. SK2019A0013), Statistics and Science Research Foundation of China (Grant no. 2017LZ11), and College Student Innovation and Entrepreneurship Training Program (Grant nos. 202110357010 and S202110357021).