1 Introduction

Uncertainty is a crucial concept for decision-making problems. It is not easy to make precise decisions in life since each information contains vagueness, uncertainty, imprecision. Fuzzy Set(FS) Theory, Zadeh’s [49] pioneering work, proposed a membership function to solve problems such as vagueness, uncertainty, imprecision, and this function took value in the range of [0,1]. FS Theory had solved many problems in practice, but there was no membership function in real life, which only includes acceptances. Rejection is as important as acceptance in real life. Atanassov [3] clarified this problem and posed the Intuitionistic Fuzzy Set(IFS) Theory using the membership function as well as the non-membership function. In IFS, the sum of membership and non-membership grades is 1. This condition is also a limitation for solutions of vagueness, uncertainty, imprecision. Yager [46] has presented a solution to this situation and suggested Pythagorean Fuzzy Sets(PFS). PFS is more comprehensive than IFS because it uses the condition that the sum of the squares of membership and non-membership grades is equal to or less than 1. PFS is also a particular case of the Neutrosophic Set initiated by Smarandache [43]. There are many studies in the literature on FS, IFS, and PFS theories [1, 2, 13,14,15,16,17,18,19,20, 23,24,25,26,27, 29, 30, 38,39,40,41, 43,44,45, 47, 48, 51]. Despite all the possible solutions, these theories have limitations. How to set the membership function in each particular object and the deficiencies in considering the parametrization tool can be given as examples of these limitations. These limitations handicap decision-makers from making a correct decision during the analysis.

A new method, called Soft Set, was proposed by Molodtsov [33], in which the preferences for each alternative were given in distinct parameters, and thus a solution was found for the limitations expressed. Immediately after the occurrence of SS theory, Fuzzy Soft Sets [29] and Intuitionistic Fuzzy Soft Set(IFSS) [30] were defined and their various properties were studied [31, 32]. Pythagorean Fuzzy Soft Set(PFSS) is defined by Peng et al. [38]. PFSS is a natural generalization of IFSS and is a parameterized family of PFSs. In [4, 5, 22, 34, 42], the main features of PFSS were examined and applied to various areas such as medical diagnosis, selection of a team of workers for business, stock exchange investment problem. The benefit of these extended theories is that they are capable of simplifying the characterization of real-life cases with the help of their parameterized feature.

It is possible that people may hesitate during decision-making. In order to avoid human hesitations from adversely affecting the decision-making process, hesitation value is also taken into account in PFS, just like IFS. Thus, experts may have hesitations about membership grades. If the expert participating in the decision process is only one person, this expert’s error or bias will affect the process negatively. However, hesitation is subjective and the expert’s hesitation can be directed by her/his own perceptions. In this case, enriching the decision process, making the evaluation with alternative decisions more meaningful, combining the subjective evaluations of more than one expert instead of the subjective evaluation of a single expert will provide a healthier decision-making process. With this in mind, Agarwal et al. [1] defined Generalized Intuitionistic Fuzzy Soft Sets(GIFSS). Feng et al. [14] identified some problems and difficulties in the definition of GIFSS and operations related to GIFSS in the manuscript of Agarwal et al [1]. Kirisci [27] defined GPFSS, considering the fixes in [14].

GPFSS ensures the frame for evaluating the reliability of the info in the PFSS so as to compensate for any distortion in the info given. The most important benefit of incorporation of the generalized parameter into the analysis is to decrease the likelihood of errors induced by the imprecise info by taking the chairperson’s view on the same. For example, a patient may give wrong information to a physician about her/his symptoms. If the physician does not notice this wrong information, errors in diagnosis and treatment will occur. In this case, an experienced physician can measure the reliability of the information given by the patient with a generalization parameter. So, there is a requirement for a generalization parameter, demonstrating an expert’s level of confidence in the reliability of the info, respectably making the approach quite close to real-world cases. This assists in extracting the singular bias from the input data and gets more credibility to the final decision. GPFSS has a generalization parameter to states the uncertainties.

Usually, the most present multiple criteria decision analysis (MCDA) methods goal to discover a unique solution that maximum achieves the general criteria as much as possible. The VIKOR method is skilled in specifying a set of compromise solutions(CS) in the existence of contradictory criteria according to the choices of the decision-makers. The VIKOR method developed by Opricovic and Tzeng [35, 36] focuses on ranking and selecting from a set of alternatives, and determines compromise solutions for a problem with conflicting criteria, which can help the decision makers to reach a final decision. This abbreviation comes from the Serbian name of the method: Višekriterijumska Optimizacija I Kompromisno Resenje(multicriteria optimization and compromise solution). Here, the CS is a feasible solution that is the closest to the ideal solution, and a compromise means an agreement established by mutual concessions [21, 35, 37].

VIKOR focalizes on the specification of a set of CSs that are built by reciprocal compensations [6, 7, 28, 35, 37] to efficiently help MCDA and improve the quality of DM [50]. In recent years, VIKOR method studies have acquired popularity in varied MCDA areas because of their large ability to compromise rank performance by comparing the measure of closeness to the ideal solution [21]. Because of these properties, the VIKOR method has become popular and has been used in many areas [21]. The VIKOR method usually utilizes the positive- and negative-ideal solutions to serve as points of reference, and the specification of these ideal solutions is immediately produced from the evaluative ratings of all alternatives according to various criteria [35, 37]. Chen [12] defined the VIKOR method based on distance index for Pythagoras fuzzy clusters and applied it to various problems.

PF sets include more uncertainties and are generally able to accommodate higher degrees of uncertainty. To efficiently administrate the complex uncertainty in human cognitive and DM activities, it is necessary to expand the valid VIKOR methods that previously focalized on certain settings and fuzzy environments to the PF context. Although the VIKOR method is a very qualified approach in terms of DM processes, it has been used very little in medical DM, health services, and child development.

The aim of this study is to develop a new approach to solve MADM problems by integrating the VIKOR method with GPFSS. It is known that the importance of the VIKOR method is to determine the CS and to improve the standard of the DM process. The VIKOR method is balancing the majority’s maximum group utility and the opponent’s minimum individual regret. In the VIKOR method, the best and worst solutions serve as the point of reference and the distance between the best solution and evaluative ratings does not have an upper bound. The new idea of remoteness index is considerable because it ensures the upper bounds by dividing the distance between the best solution and evaluative ratings by the distance between the best solution and the worst solution. The novel ranking indexes based on the new idea of remoteness indexes are debated to ensure more efficient compromise rankings. Thus, a remoteness index-based VIKOR method for GPFSS is improved.

The negatives caused by early childhood children staying at home during quarantine times are measured with new algorithms. All areas of development are in interaction with each other, and the problem that may occur in one area has reflections on other areas. In this context, problems in the cognitive field may cause problems in other areas, as well as problems in other areas that may reflect on cognitive development. Therefore, it is necessary to know the development characteristics of each field well and to identify the problems that may arise in a timely and accurate manner. The pre-school period has a very important place in human life. In this period, the child’s personality development and emotional, social, cognitive, language, and academic development are largely completed. The child completes her/his cognitive and behavioral development in a healthy way with the contribution of her family and school environment. That’s why new algorithms have been applied to such an important subject.

The originality of this study can be expressed as follows: New distance measures such as Euclidean, Hamming, and Generalized for GPFFS are defined. The properties of these distance measures were examined, and an example was given. The displaced and fixed ideals are defined for PFVs. In addition, the displaced and fixed remoteness indexes are given for PFVs. Two types of weights called precise weights and PF weights are defined and discuss in their methods of generation. The precise importance weights are obtained by the expectation score function. The PF importance weights represented the importance and unimportance degrees of the criteria. The displaced and fixed ideals, the displaced and fixed remoteness indexes are identified as ranking indexes (four groups). Based on these definitions, four algorithms expressing VIKOR methods have been proposed. The problem of the cognitive and behavioral development of early childhood children in quarantine is examined with the suggested remoteness-based VIKOR method.

2 Preliminaries

In this section, basic definitions, theorems, propositions and properties that will be used in the whole study will be given. Throughout the study, the initial universe, parameters sets will denote \(\Upsilon\) , \(\Sigma\), respectively.

2.1 Pythagorean fuzzy environment

The theory of IFS was developed by Atanassov [3] and is a natural generalization of the theory of FS. IFS is defined by MS and NMS degrees and, for this reason, can indicate the fuzzy character of data in more elaboration comprehensively. The significant qualification of IFS is that it appoints to each element MS and NMS degrees (\(MS + NMS\le 1\)). However, in some DM processes, the sum of the MS and NMS values obtained may be greater than 1. In this case, we can take the sum of the squares of these MS and NMS values obtained, which will be less than or equal to 1.

As an original idea, PFSs were created by Yager [45]. PFS is a very useful tool for uncertainty. PFS offers good results especially in solutions where IFS is insufficient. The differences between PFSs and IFSs can be mentioned as follows: For IFS, \(MS+NMS\le 1\), \(0\le MS+NMS \le 1\), \(\pi =1-(MS+NMS)\) and \(\pi +MS+NMS=1\).

For PFS, \(MS+NMS\le 1\) or \(MS+NMS\ge 1\), \(0\le MS^{2}+NMS^{2} \le 1\), \(\pi =\sqrt{1-(MS^{2}+NMS^{2})}\) and \(\pi ^{2}+MS^{2}+NMS^{2}=1\).

For, the initial universe set \(\Upsilon\), the function \(u_{A}(x):\Upsilon \rightarrow [0,1]\) is called FS on \(\Upsilon\). The FS can be indicated by

$$\begin{aligned} A=\left\{ (a_{i},u_{A}(a_{i})):u_{A}(a_{i})\in [0,1]; \forall a_{i}\in \Upsilon \right\} . \end{aligned}$$
(2.1)

Choose the MS function \(u_{B}:\Upsilon \rightarrow [0,1]\) and the NMS function \(v_{B}:\Upsilon \rightarrow [0,1]\). Let’s assume that the condition \(0\le u_{B}(a)+v_{B}(a)\le 1\) for any \(a\in \Upsilon\) is satisfied. Then, the set \(B=\{( a, u_{B}(a), v_{B}(a)): a\in \Upsilon \}\) is said to be an IFS B on \(\Upsilon\) . In this case, the condition \(\pi _{B}=1-u_{B}(x)-v_{B}(x)\) holds [3].

Suppose that condition \(0\le [u_{C}(x)]^{2}+[v_{C}(x)]^{2}\le 1\) satisfies for \(u_{C}:\Upsilon \rightarrow [0,1]\) and \(v_{C}:\Upsilon \rightarrow [0,1]\). Then, an PFS C in \(\Upsilon\) is defined by \(C=\{( x, u_{C}(x), v_{C}(x)): x\in \Upsilon \}\). In this case, the condition \(\pi _{C}=\sqrt{1-[u_{C}(x)]^{2}-[v_{C}(x)]^{2}}\) holds [45, 46, 48].

Let C be a PFS over \(\Sigma\). In this definition, (FC) can be represented as

$$\begin{aligned} (F, C)= & {} \bigg \{\big (a, u_{F}(a), v_{F}(a) \big ): a\nonumber \\\in & {} {\mathscr {P}}, u_{F}(a)\in [0,1], v_{F}(a)\in [0,1]\bigg \}. \end{aligned}$$
(2.2)

The set of all PFS on \(\Upsilon\) will be denoted by \(\Omega (\Upsilon )\).

Let \(\Upsilon\) ,\(\Sigma\) be an initial universe and a parameters sets, respectively. For \(D\subseteq \Sigma\), consider a set-valued mapping \(F:D\rightarrow \omega (\Upsilon )\), where the power set of \(\Upsilon\) is showed by \(\omega (\Upsilon )\). Therefore, a pair (FD) is called a SS on \(\Upsilon\).

For \(E\subseteq \Sigma\), choose \(F:E\rightarrow \omega (\Upsilon )\), where the set of all PFSs over \(\Upsilon\) is indicated by \(\omega (\Upsilon )\). Then, a pair (FE) is called a PFSS on \(\omega (\Upsilon )\) [38].

Choose the Pythagorean fuzzy numbers (PFNs) \(R=(u_{R}, v_{R})\) [51]. For three PFNs \(\theta = (u,v), \theta _{1}=( u_{1}, v_{1}), \theta _{2}=( u_{2}, v_{2})\), some fundamental operations are as follows [45, 48]: For \(\alpha >0\), \({\bar{\theta }}=( v, u )\), \(\theta _{1} \oplus \theta _{2} = ( \sqrt{u_{1}^{2}+u_{2}^{2}-u_{1}^{2}u_{2}^{2}} , v_{1}v_{2})\), \(\theta _{1} \otimes \theta _{2} = ( u_{1}u_{2}, \sqrt{v_{1}^{2}+v_{2}^{}2-v_{1}^{2}v_{2}^{2}} )\), \(\theta _{1} \wedge \theta _{2} = ( \min \{u_{1}, u_{2}\} , \max \{v_{1}, v_{2}\})\), \(\theta _{1} \vee \theta _{2} = ( \max \{u_{1}, u_{2}\} , \min \{v_{1}, v_{2}\})\), \(\alpha .\theta = ( \sqrt{1-(1-u^{2})^{\alpha }},v^{\alpha } )\), \(\theta ^{\alpha } = ( u^{\alpha }, \sqrt{1-(1-v^{2})^{\alpha }})\). The PFS defined by MS and NMS which satisfies the condition \(MS^{2} + NMS^{2} \le 1\) is initiated by Yager [45]. For example, choose \(p=(0.87, 0.42)\). Since \((0,87) + (0.42) > 1\), these values are not suitable for IFS. However, since \((0,87)^{2}+(0.42)^{2}<1\), these values can be used with PFS.

A new PFN formula consisting of the strength of commitment \(r_{p}\) and the direction of commitment \(d_{p}\) was also given by Yager [45]. The new PFN is denoted by \(P=(r_{p}, d_{p})\), where \(r_{p}\in [0,1]\). The stronger commitment is defined by the larger the value of \(r_{p}\), in this case the lower the uncertainty of the commitment. \(r_{p}\) and \(d_{p}\) are characterized by \(u_{p}=r_{p}\cos (\theta _{p})\), \(v_{p}=r_{p}\sin (\theta _{p})\), i.e., these values are related to MS and NMS grades. Here \(\theta _{p}\) is expressed in radian and calculate \(\theta _{p}=(1-d_{p})(\pi /2)\). If \(P=(r_{p}, d_{p})\) is defined in terms of polar coordinates, it is denoted as \((r_{p}, \theta _{p})\). In this case, \(d_{p}\) has a formula in the form of \(d_{p}=1-(2.\theta _{p}/\pi )\). It can then be said that PFNs consists of \(u_{p}, v_{p}, \pi _{p}, r_{p}, d_{p}\) and \(\theta _{p}\) parameters. If \(u_{p}, v_{p} \in [0,1]\), then it is clear that \(u_{p}^{2}\le u_{p}\), \(v_{p}^{2}\le v_{p}\).

2.2 Generalized Pythagorean fuzzy soft sets

Consider the \(L=\{(a,b): a,b \in [0,1], a^{2}+b^{2}<1\}\) with the partial order \(\le _{L}\). The complete lattice \((L, \le _{L})\) is defined by \((a,b)\le _{L} (i,j) \quad \Leftrightarrow \quad a\le i \quad and \quad b \ge j\), for all \((a,b), (i,j)\in L\).

The \((a,b)\in {\mathscr {L}}\) given in this way is called Pythagorean fuzzy value (PFN)(or Pythagorean fuzzy number(PFN)) [16].

Definition 2.1

[25] Take the PFS C. Let the PFS \(f:C \rightarrow L\). If the set

$$\begin{aligned} (F,C,f)= & {} \bigg \{\big (a, u_{C}(a), v_{C}(a), f_{C}(a) \big ): a\nonumber \\\in & {} \sum , f_{C}(a)\in L, u_{C}(a)\in [0,1], v_{C}(a)\in [0,1]\bigg \} \end{aligned}$$
(2.3)

is a PFSS, the (FCf) is called Pythagorean fuzzy parameterized Pythagorean fuzzy soft set (generalized Pythagorean fuzzy soft set-GPFSS), where (FC) is a fundamental PFSS (FPFSS). Furthermore, \(f_{C}\) is said to be the parametric PFS(PPFS) and is indicated by \(u_{f}, v_{f}\).

Example 2.2

Choose the patients set \(H=\{h1, h2, h3\}\) and the symptoms set \(B=\{b_{1}, b_{2}, b_{3}, b_{4}\}\) and \(C=\{b_{1}, b_{3}, b_{4}\}\subset B\). The values \(f_{C}(b_{1})=(0.77, 0.35)\), \(f_{C}(b_{3})=(0.65, 0.42)\) and \(f_{C}(b_{4})=(0.88, 0.24)\) express the diagnostic view of a physician. The (FCf) is shown in Table 1.

Table 1 The GPFSS (FCf)

For all \(R=(u_{R}, v_{R})\in L\) , the score function for PFNs ie described

$$\begin{aligned} SC_{R}=u_{R}^{2}-v_{R}^{2} \end{aligned}$$
(2.4)

where \(SC:L\rightarrow [-1,1]\) [1, 51].

For any two PFNs RT,

$$\begin{aligned} \left\{ \begin{array}{ccl} R \prec T, \quad \textit{if} \quad &{} \mathscr {SC}(R)< \mathscr {SC}(T), \\ R \succ T, \quad \textit{if} \quad &{} \mathscr {SC}(R)> \mathscr {SC}(T), \\ R \sim T, \quad \textit{if} \quad &{} \mathscr {SC}(R) = \mathscr {SC}(T); \end{array} \right. \end{aligned}$$
(2.5)

As can be seen from (2.4), the larger the score \(s_{R}\), the greater the PFN R [16]. It should be noted that \(\mathscr {SC}\) cannot differentiate some evidently distinct PFNs that have the same score. To explain this situation, we can give examples using (2.5): Take two PFNs RT as \(R=(0.481, 0.402)\) and \(T=(0.527, 0.456)\). Then, \(SC_{R}=0.0697\) and \(SC_{T}=0.0697\). Again, for \(R=(0.123, 0.123)\) and \(T=(0.456, 0.456)\) we can write \(R \sim T\). Therefore, if only the scoring function is used for comparison, it is not possible to make a comparison between these numbers. To overcome this problem, we can define a new function [1], as follows:

The mapping \(\mathscr {AC}:{\mathscr {L}}\rightarrow [0,1]\) is called accuracy function, if

$$\begin{aligned} \mathscr {AC}_{R}=u_{R}^{2}+v_{R}^{2} \end{aligned}$$
(2.6)

for all \(R=(u_{R}, v_{R})\in {\mathscr {L}}\) [39].

Using the (2.4) and (2.6), for comparing PFVs, the following method is presented by Agarwal et al. [1]. For any two PFNs RT, if \(\mathscr {SC}(R) = \mathscr {SC}(T)\), then

$$\begin{aligned} \left\{ \begin{array}{ccl} R>T, \textit{if} \quad &{} \mathscr {AC}(R)> \mathscr {AC}(T), \\ R<T, \textit{if} \quad &{} \mathscr {AC}(R)< \mathscr {AC}(T) \\ R \sim T, \quad \textit{if} \quad &{} \mathscr {AC}(R)= \mathscr {AC}(T); \end{array} \right. \end{aligned}$$
(2.7)

For a binary relation \(\le _{(\mathscr {SC}, \mathscr {AC})}\in {\mathscr {L}}\) and \(R,T \in {\mathscr {L}}\), it can be written as

$$\begin{aligned} R\le & {} _{(\mathscr {SC}, \mathscr {AC})} T \quad \Leftrightarrow \quad (\mathscr {SC}_{R}< \mathscr {SC}_{T})\vee (\mathscr {SC}_{R}\nonumber \\= & {} \mathscr {SC}_{T} \wedge \mathscr {AC}_{R} \le \mathscr {AC}_{T}). \end{aligned}$$
(2.8)

From [14], we can give the new definition:

Definition 2.3

[27] The mapping \(\mathscr {ES}: {\mathscr {L}}\rightarrow [0,1]\) is called expectation score function such that for all \(R=(m_{R}, n_{R})\in {\mathscr {L}}\)

$$\begin{aligned} \mathscr {ES}_{R}=\frac{u_{R}^{2}-v_{R}^{2}+1}{2}. \end{aligned}$$
(2.9)

Proposition 2.4

[27] Let \(\mathscr {ES}: {\mathscr {L}}\rightarrow [0,1]\) and \(R=(u_{R}, v_{R}) \in {\mathscr {L}}\). Then, we have

  1. i.

    \(\mathscr {ES}(0,1)=0\) and \(\mathscr {ES}(1,0)=1\),

  2. ii.

    \(\mathscr {ES}(u_{R}, v_{R})\) is increasing with respect to \(u_{R}\),

  3. iii.

    \(\mathscr {ES}(u_{R}, v_{R})\) is decreasing with respect to \(v_{R}\).

Definition 2.5

[27] For two PFNs \(R,T \in {\mathscr {L}}\) and the relation \(\le _{(u,\mathscr {ES})}\) on \({\mathscr {L}}\), we have

$$\begin{aligned} R\le & {} _{(u,\mathscr {ES})} T \quad \Leftrightarrow \quad (u_{R}< u_{T})\vee (u_{R}\nonumber \\= & {} u_{T} \wedge \mathscr {ES}_{R} \le \mathscr {ES}_{T}). \end{aligned}$$
(2.10)

Theorem 2.6

[27] The relation \(\le _{(u,\mathscr {ES})}\) is a partial order on \({\mathscr {L}}\).

By replacing the approval rates in Definition 2.5, the other relation \(\le _{(\mathscr {ES},u)}\) is written as follows:

Definition 2.7

[27] For two PFNs \(R,T \in {\mathscr {L}}\) and the relation \(\le _{(\mathscr {ES},u)}\) on \({\mathscr {L}}\),

$$\begin{aligned} R \le _{(\mathscr {ES},u)} T \quad \Leftrightarrow \quad (\mathscr {ES}_{R}< \mathscr {ES}_{T})\vee (\mathscr {ES}_{R}= \mathscr {ES}_{T} \wedge u_{R} \le v_{T}). \end{aligned}$$
(2.11)

The relationship with each other of \(\le _{(\mathscr {SF},\mathscr {AC})}\) and \(\le _{(\mathscr {ES},u)}\) can be given as follows:

Proposition 2.8

[27] Let R and T be PFNs in \({\mathscr {L}}\). Then, \(R \le _{(\mathscr {SC},\mathscr {AC})} T\) \(\Leftrightarrow\) \(R \le _{(\mathscr {ES},u)} T\).

If we consider the debates of [14, 44], then the following proposition is given:

Proposition 2.9

[27] Take two PFVs \(R=(u_{R}, v_{R})\) and \(T=(u_{T}, v_{T})\) in \({\mathscr {L}}\). Then, the following conditions are equivalent:

  1. i.

    \(\mathscr {SC}_{R}=\mathscr {SC}_{T} \wedge \mathscr {AC}_{R} \le \mathscr {AC}_{T}\),

  2. ii.

    \(\mathscr {ES}_{R}=\mathscr {ES}_{T} \wedge u_{R} \le u_{T}\),

  3. iii.

    \(\mathscr {ES}_{R}=\mathscr {ES}_{T} \wedge v_{R} \le v_{T}\),

  4. iv.

    \(\mathscr {SC}_{R}=\mathscr {SC}_{T} \wedge v_{R} \le u_{T}\),

  5. v.

    \(\mathscr {SC}_{R}=\mathscr {SC}_{T} \wedge v_{R} \le v_{T}\).

Definition 2.10

[51] Let \(R=(u_{R}, v_{R}), T=(u_{T}, v_{T})\in {\mathscr {L}}\) be two PFVs. Then, for \(\alpha >0\), we have the following operations:

$$\begin{aligned} R \oplus T= & {} \left( \sqrt{u_{R}^{2}+u_{T}^{2}-u_{M}^{2}u_{T}^{2}}, v_{R}v_{T}\right) , \end{aligned}$$
(2.12)
$$\begin{aligned} R \otimes T= & {} \left( u_{R}u_{T}, \sqrt{v_{R}^{2}+v_{T}^{2}-v_{R}^{2}v_{T}^{2}} \right) , \end{aligned}$$
(2.13)
$$\begin{aligned} \alpha R= & {} \left( \sqrt{1-(1-u_{R}^{2})^{\alpha }}, (v_{R})^{\alpha }\right) , \end{aligned}$$
(2.14)
$$\begin{aligned} R^{\alpha }= & {} \left( (u_{R})^{\alpha }, \sqrt{1-(1-v_{R}^{2})^{\alpha }}, \right) . \end{aligned}$$
(2.15)

Theorem 2.11

[51] For \(R=(u_{R}, v_{R}), T=(u_{T}, v_{T})\in {\mathscr {L}}\) and \(\alpha , \alpha _{1}, \alpha _{2}>0\),

  1. i.

    \(R \oplus T = T \oplus R\),

  2. ii.

    \(R \otimes T= T \otimes R\),

  3. iii.

    \(\alpha (R+T)=\alpha R \oplus \alpha T\),

  4. iv.

    \(\alpha _{1}R \oplus \alpha _{2}R = (\alpha _{1}+ \alpha _{2})R\),

  5. v.

    \((R \otimes T)^{\alpha }=R^{\alpha }\otimes T^{\alpha }\),

  6. vi.

    \(R^{\alpha _{1}} \oplus R^{\alpha _{2}}=R^{(\alpha _{1}+\alpha _{2})}\).

Proposition 2.12

[27] For \(R=(u_{R}, v_{R}), T=(u_{T}, v_{T}), V=(u_{V}, v_{V})\in {\mathscr {L}}\), then we have,

$$\begin{aligned} R \le _{(u,\mathscr {ES})} V \Rightarrow T \oplus R \le _{(u,\mathscr {ES})} T \oplus V. \end{aligned}$$
(2.16)

Corollary 2.13

[27] Let For \(R=(u_{R}, v_{R}), T=(u_{T}, v_{T})\in {\mathscr {L}}\) and \(\alpha , \alpha _{1}, \alpha _{2}\in {\mathbb {R}}^{+}\). Then,

  1. i.

    If \(R \le _{(u,\mathscr {ES})} T\),    then    \(\alpha R \le _{(u,\mathscr {ES})} \alpha T\),

  2. ii.

    If \(\alpha _{1} \le \alpha _{2}\),    then    \(\alpha _{1} R \le _{(u,\mathscr {ES})} \alpha _{2} R\).

3 Distance measures for GPFSSs

The Hamming, Euclidean and generalized distance measures for GPFSSs are given in this section. Some properties of distance measures are discussed. The new distance measures offered in this section will be then applied to build useful concepts of remoteness indices, consisting of displaced positive- and negative-ideal remoteness indices as well as fixed positive- and negative-ideal remoteness indices.

Let \(\Upsilon =\{x_{1}, x_{2}, \cdots , x_{m}\}\), \(\sum =\{e_{1}, e_{2}, \cdots , e_{n}\}\) be universe set and parameter set, respectively. Choose the two GPFSSs \(\Omega _{1}, \Omega _{2}\). For convenience, take \(\triangle X(x_{i})=X(\Omega _{1})_{ej}(x_{i})-X(\Omega _{2})_{ej}(x_{i})\).

Definition 3.1

For two GPFSSs \(\Omega _{1}=(F, A, f)\), \(\Omega _{2}=(G, B, g)\), the normalized Hamming distance, normalized Euclidean distance and generalized distance measures between \(\Omega _{1}\) and \(\Omega _{2}\) are defined, respectively, as follows:

$$\begin{aligned} U_{H}(\Omega _{1}, \Omega _{2})= & {} \frac{1}{5mn}\begin{Bmatrix} \sum\limits _{i=1}^{m}\sum\limits _{j=1}^{n}\bigg [|\triangle u(x_{i})|+|\triangle v(x_{i})|+|\triangle \pi (x_{i})|+|\triangle r(x_{i})|+|\triangle d(x_{i})|\\ +|u(f)_{e_{j}}(x_{i})-u(g)_{e_{j}}(x_{i})|+|v(f)_{e_{j}}(x_{i})-v(g)_{e_{j}}(x_{i})|+|\pi (f)_{e_{j}}(x_{i})-\pi (g)_{e_{j}}(x_{i})|\\ +|r(f)_{e_{j}}(x_{i})-r(g)_{e_{j}}(x_{i})|+|d(f)_{e_{j}}(x_{i})-d(g)_{e_{j}}(x_{i})|\bigg ] \end{Bmatrix} \end{aligned}$$
(3.1)
$$\begin{aligned} U_{E}(\Omega _{1}, \Omega _{2})= & {} \begin{Bmatrix} \frac{1}{5mn}\sum\limits _{i=1}^{m}\sum\limits _{j=1}^{n}[(\triangle u(x_{i}))^{2}+(\triangle v(x_{i}))^{2}+(\triangle \pi (x_{i}))^{2}+(\triangle r(x_{i}))^{2} \\+(\triangle d(x_{i}))^{2} +(u(f)_{e_{j}}(x_{i})-u(g)_{e_{j}}(x_{i}))^{2}+(v(f)_{e_{j}}(x_{i})-v(g)_{e_{j}}(x_{i}))^{2}\\ +(\pi (f)_{e_{j}}(x_{i})-\pi (g)_{e_{j}}(x_{i}))^{2}+(r(f)_{e_{j}}(x_{i})-r(g)_{e_{j}}(x_{i}))^{2}\\ +(d(f)_{e_{j}}(x_{i})-d(g)_{e_{j}}(x_{i}))^{2}] \end{Bmatrix}^{1/2}. \end{aligned}$$
(3.2)
$$\begin{aligned} U_{G}(\Omega _{1}, \Omega _{2})= & {} \left\{ \begin{array}{l} \frac{1}{5mn}\sum\limits _{i=1}^{m}\sum\limits _{j=1}^{n}(\triangle u(x_{i}))^{\lambda }+(\triangle v(x_{i}))^{\lambda }+(\triangle \pi (x_{i}))^{\lambda }\\ +(\triangle d(x_{i}))^{\lambda }+(\triangle r(x_{i}))^{\lambda }+(u(f)_{e_{j}}(x_{i})-u(g)_{e_{j}}(x_{i}))^{\lambda }\\ +(v(f)_{e_{j}}(x_{i})-v(g)_{e_{j}}(x_{i}))^{\lambda }+(\pi (f)_{e_{j}}(x_{i})-\pi (g)_{e_{j}}(x_{i}))^{\lambda }\\ +(r(f)_{e_{j}}(x_{i})-r(g)_{e_{j}}(x_{i}))^{\lambda }+(d(f)_{e_{j}}(x_{i})-d(g)_{e_{j}}(x_{i}))^{\lambda } \end{array} \right\} ^{1/\lambda } \end{aligned}$$
(3.3)

where \(\lambda \ge 1\).

Definition 3.2

The weighted normalized Hamming distance, weighted normalized Euclidean distance and weighted generalized distance measures between \(\Omega _{1}\) and \(\Omega _{2}\) are defined, respectively, as follows:

$$\begin{aligned} U_{WH}(\Omega _{1}, \Omega _{2})= & {} \frac{1}{5mn}\left\{ \begin{array}{ccl} \sum\limits _{i=1}^{m}\sum\limits _{j=1}^{n}\omega _{j}|\triangle u(x_{i})|+|\triangle v(x_{i})|+|\triangle \pi (x_{i})|+|\triangle r(x_{i})|\\ +|\triangle d(x_{i})|+|u(f)_{e_{j}}(x_{i})-u(g)_{e_{j}}(x_{i})|\\ +|v(f)_{e_{j}}(x_{i})-v(g)_{e_{j}}(x_{i})|+|\pi (f)_{e_{j}}(x_{i})-\pi (g)_{e_{j}}(x_{i})|\\ +|r(f)_{e_{j}}(x_{i})-r(g)_{e_{j}}(x_{i})|+|d(f)_{e_{j}}(x_{i})-d(g)_{e_{j}}(x_{i})| \end{array} \right\} \end{aligned}$$
(3.4)
$$\begin{aligned} U_{WE}(\Omega _{1}, \Omega _{2})= & {} \left\{ \begin{array}{ccl} \frac{1}{5mn}\sum\limits _{i=1}^{m}\sum\limits _{j=1}^{n}\omega _{j}(\triangle u(x_{i}))^{2}+(\triangle v(x_{i}))^{2}+(\triangle \pi (x_{i}))^{2}+(\triangle r(x_{i}))^{2}\\ +(\triangle d(x_{i}))^{2}+(u(f)_{e_{j}}(x_{i})-u(g)_{e_{j}}(x_{i}))^{2}\\ +(v(f)_{e_{j}}(x_{i})-v(g)_{e_{j}}(x_{i}))^{2}+(\pi (f)_{e_{j}}(x_{i})-\pi (g)_{e_{j}}(x_{i}))^{2}\\ +(r(f)_{e_{j}}(x_{i})-r(g)_{e_{j}}(x_{i}))^{2}+(d(f)_{e_{j}}(x_{i})-d(g)_{e_{j}}(x_{i}))^{2} \end{array} \right\} ^{1/2} \end{aligned}$$
(3.5)
$$\begin{aligned} U_{WG}(\Omega _{1}, \Omega _{2})= & {} \left\{ \begin{array}{ccl} \frac{1}{5mn}\sum\limits _{i=1}^{m}\sum\limits _{j=1}^{n}\omega _{j}(\triangle u(x_{i}))^{\lambda }+(\triangle v(x_{i}))^{\lambda }+(\triangle \pi (x_{i}))^{\lambda }+(\triangle r(x_{i}))^{\lambda }\\ +(\triangle d(x_{i}))^{\lambda }+(u(f)_{e_{j}}(x_{i})-u(g)_{e_{j}}(x_{i}))^{\lambda }\\ +(v(f)_{e_{j}}(x_{i})-v(g)_{e_{j}}(x_{i}))^{\lambda }+(\pi (f)_{e_{j}}(x_{i})-\pi (g)_{e_{j}}(x_{i}))^{\lambda }\\ +(r(f)_{e_{j}}(x_{i})-r(g)_{e_{j}}(x_{i}))^{\lambda }+(d(f)_{e_{j}}(x_{i})-d(g)_{e_{j}}(x_{i}))^{\lambda } \end{array} \right\} ^{1/\lambda } \end{aligned}$$
(3.6)

where \(\lambda \ge 1\), \(0\le \omega _{i} \le 1 (i=1,2,\cdots , n)\) and \(\sum _{i=1}^{n}\omega _{i}=1\).

Theorem 3.3

Let \(({\mathscr {L}}, \le _{{\mathscr {L}}})\) be a partially ordered set. For the GPFSSs \(\Omega _{1}, \Omega _{2}, \Omega _{3}\) belonging to L with the partial order \(\le _{L}\). A generalized distance measure of \(U_{G}\) satisfies the following conditions:

\((U_{1})\):

\(U_{G}(\Omega _{1}, \Omega _{1})=0\) (reflexivity)

\((U_{2})\):

\(U_{G}(\Omega _{1}, \Omega _{2})=0\) if \(\Omega _{1}=\Omega _{2}\)

\((U_{3})\):

\(U_{G}(\Omega _{1}, \Omega _{2})= U_{G}(\Omega _{2}, \Omega _{1})\) (symmetry)

\((U_{4})\):

\(0 \le U_{G}(\Omega _{1}, \Omega _{1}) \le 1\) (boundedness)

\((U_{5})\):

\(U_{G}(\Omega _{1}, \Omega _{2})\le U_{G}(\Omega _{1}, \Omega _{3})\) and \(U_{G}(\Omega _{2}, \Omega _{3})\le U_{G}(\Omega _{1}, \Omega _{3})\) if \(p_{1} \le _{L} p_{2} \le _{L} p_{3}\).

The proof of this theorem is obvious.

Note 1: If \(\lambda =1\) and \(\lambda =2\), \(U_{H}=U_{G}\) and \(U_{E}=U_{G}\), respectively. Then \(U_{H}\) satisfies the conditions of Theorem 3.3 with separability and triangle inequality. Therefore, the Hamming distance measure \(U_{H}\) is a metric. Further, \(U_{E}\) satisfies the conditions of Theorem 3.3 with separability. Fulfilling these conditions indicates that the \(U_{E}\) is a semi-metric. The conditions separability and triangle inequality conditions, respectively, are as follows:

(\(A_{1}\)):

\(U_{G}(\Omega _{1}, \Omega _{2})=0 \Leftrightarrow \Omega _{1}=\Omega _{2}\), for \(\lambda =1\) and \(\lambda =2\).

(\(A_{1})\):

\(U_{H}(\Omega _{1}, \Omega _{3}) \le U_{H}(\Omega _{1}, \Omega _{2})+ U_{H}(\Omega _{2}, \Omega _{3})\)

Theorem 3.4

For the GPFSSs \(\Omega _{1}, \Omega _{2}\), \(U(\Omega _{1}, \Omega _{2})=U(\Omega _{1}^{t}, \Omega _{2}^{t})\).

Proof

We will only prove for normalized Hamming distance. It can be proved similarly for normalized Euclidean distance and generalized distance measures. The distance between two GPFSSs \(\Omega _{1}, \Omega _{2}\) can be written as follows:

$$\begin{aligned} U_{H}(\Omega _{1}, \Omega _{2})= \frac{1}{5mn}\left\{ \begin{array}{ccl} \sum\limits_{i=1}^{m}\sum\limits_{j=1}^{n}[|u(\Omega _{1})_{ej}(x_{i})-u(\Omega _{2})_{ej}(x_{i})|+|v(\Omega _{1})_{ej}(x_{i})-v(\Omega _{2})_{ej}(x_{i})|\\ +|\pi (\Omega _{1})_{ej}(x_{i})-\pi (\Omega _{2})_{ej}(x_{i})|+|r(\Omega _{1})_{ej}(x_{i})-r(\Omega _{2})_{ej}(x_{i})|\\ +|d(\Omega _{1})_{ej}(x_{i})-d(\Omega _{2})_{ej}(x_{i})|+|u(f)_{e_{j}}(x_{i})-u(g)_{e_{j}}(x_{i})|\\ +|v(f)_{e_{j}}(x_{i})-v(g)_{e_{j}}(x_{i})|+|\pi (f)_{e_{j}}(x_{i})-\pi (g)_{e_{j}}(x_{i})|\\ +|r(f)_{e_{j}}(x_{i})-r(g)_{e_{j}}(x_{i})|+|d(f)_{e_{j}}(x_{i})-d(g)_{e_{j}}(x_{i})|] \end{array} \right\} \end{aligned}$$
(3.7)

From the complement of GPFSS,

$$\begin{aligned} U(\Omega _{1}^{t}, \Omega _{2}^{t})= \frac{1}{5mn}\left\{ \begin{array}{ccl} \sum\limits_{i=1}^{m}\sum\limits_{j=1}^{n}[|u(\Omega _{2})_{ej}(x_{i})-u(\Omega _{1})_{ej}(x_{i})|+|v(\Omega _{2})_{ej}(x_{i})-v(\Omega _{1})_{ej}(x_{i})|\\ +|\pi (\Omega _{2})_{ej}(x_{i})-\pi (\Omega _{1})_{ej}(x_{i})|+|r(\Omega _{2})_{ej}(x_{i})-r(\Omega _{1})_{ej}(x_{i})|\\ +|d(\Omega _{2})_{ej}(x_{i})-d(\Omega _{1})_{ej}(x_{i})|+|u(g)_{e_{j}}(x_{i})-u(f)_{e_{j}}(x_{i})|\\ +|v(g)_{e_{j}}(x_{i})-v(f)_{e_{j}}(x_{i})|+|\pi (g)_{e_{j}}(x_{i})-\pi (f)_{e_{j}}(x_{i})|\\ +|r(g)_{e_{j}}(x_{i})-r(f)_{e_{j}}(x_{i})|+|d(g)_{e_{j}}(x_{i})-d(f)_{e_{j}}(x_{i})|] \end{array} \right\} \end{aligned}$$
(3.8)

the MS degrees of \(\Omega _{1}^{t}\) relate to the NMS degrees of \(\Omega _{2}^{t}\). Therefore, the distance between \(\Omega _{1}, \Omega _{2}\) is the same as the distance between \(\Omega _{1}^{t}, \Omega _{2}^{t}\). \(\square\)

Example 3.5

Let \(\Upsilon =\{x_{1}, x_{2}\}\), \(M=N=\{e_{1}, e_{2}, e_{3}\}\) be universe set and parameter sets, respectively. Choose two GPFFSs \(\Omega _{1}=(F, M, f)\), \(\Omega _{2}=(G, N, g)\) as follows:

$$\begin{aligned} F(e_{1})= & {} \{<x_{1}, 0.72, 0.55>,<x_{2}, 0.68, 0.61>\} \\ F(e_{2})= & {} \{<x_{1}, 0.84, 0.39>,<x_{2}, 0.75, 0.65>\} \\ F(e_{3})= & {} \{<x_{1}, 0.91, 0.27>,<x_{2}, 0.83, 0.45>\} \\ f= & {} \{<e_{1}, 0.48, 0.83>,<e_{2}, 0.91, 0.38>,\\< & {} e_{3}, 0.69, 0.54>\} \\ G(e_{1})= & {} \{<x_{1}, 0.66, 0.55>,<x_{2}, 0.76, 0.34>\} \\ G(e_{2})= & {} \{<x_{1}, 0.77, 0.52>,<x_{2}, 0.83, 0.25>\} \\ G(e_{3})= & {} \{<x_{1}, 0.80, 0.37>,<x_{2}, 0.74, 0.56>\} \\ g= & {} \{<e_{1}, 0.46, 0.63>,<e_{2}, 0.82, 0.56>,\\< & {} e_{3}, 0.67, 0.47>\} \end{aligned}$$

Hamming distance measure is calculated as follows:

$$\begin{aligned}&U_{H}(\Omega _{1}, \Omega _{2})\\&\quad = \frac{1}{5.3.2}\left[ \begin{array}{ccl} &{}&{}|0.72-.66|+|0.68-0.76|+|0.84-0.77|+|0.75-0.83|+|0.91-0.80|+|0.83-0.74|\\ &{}&{}+|0.4232-0.512|+|0.6772-0.56|+|0.38-0.37|+|0.1225-0.5|+|0.315-0.47|+|0.33-0.38|\\ &{}&{}+|0.906-0.86|+|0.736-0.56|+|0.93-0.93|+|0.992-0.87|+|0.95-0.88|+|0.994-0.93|\\ &{}&{}+|0.585-0.56|+|0.5331-0.73|+|0.7216-0.62||0.544-0.8|+|0.8-0.724|+|0.32-0.586|\\ &{}&{}+|0.48-0.46|+|0.91-0.82|+|0.69-0.67|+|0.83-0.63|+|0.38-0.56|+|0.54-0.47| \end{array} \right] . \end{aligned}$$

and \(U_{H}(\Omega _{1}, \Omega _{2})=0.11\). Similarly, Euclidean distance measure \(U_{E}(\Omega _{1}, \Omega _{2})=0.14\) and for \(\lambda =3\), generalized distance measure \(U_{G}(\Omega _{1}, \Omega _{2})=0.1646\).

4 The remoteness index-based VIKOR method for GPFSSs

The purpose of this section is to develop new remoteness index-based GPFSS-VIKOR methods that utilize the concepts of displaced and fixed remoteness indices based on new distance measures to address MCDA problems within the GPFSS environment. The definition of positive and negative ideal points for a solution is the basic idea of the VIKOR method. It focalizes on ranking and selecting from a finite set of possible alternatives in the presence of contradictory and non-measurable criteria. It interprets a multi-criteria ranking index based on the nearness to the ideal solution. By evaluating the alternatives according to each criterion, a compromise ranking can be obtained when the relative closeness measure is compared with the ideal alternative. Thus, the derived compromise solution is a feasible solution, which is the closest to the positive ideal solution and farthest from the negative ideal solution, and a compromise means an agreement established by mutual concessions made between the alternatives. The acquired CS can be admitted because it ensures a maximum group utility of the majority and a minimum of the individual regret of the opponent. The CSs can be the basis for negotiations, involving the decision maker’s preference on criteria weights. In the VIKOR method, the performance ratings of the alternatives according to a set of criteria are quantified as crisp values [11].

In this section, we define the displaced positive ideal PFV (DPI) and the displaced negative ideal PFV (DNI) which help to define the concept of displaced remoteness index. The fixed remoteness index is defined based on the fixed positive ideal PFS (FPI) and the fixed negative ideal PFV (FNI). The displaced and fixed group utility, individual regret and compromise indexes are the precise importance and Pythagorean fuzzy importance weights. The extra parameter parametric Pythagorean fuzzy set in GPFSS, which given by the head or director of decision-making committee, is used to define the precise importance and Pythagorean fuzzy importance weights. Each algorithm shows the complete procedure of the remoteness-based VIKOR method based on the displaced and fixed terminologies.

Let’s give the general framework of the MCDA problem about PF information:

Choose the set of candidate alternatives \(D=\{d_{1}, d_{2}, \cdots , d_{m}\}\), m \((m\ge 2)\) and the set of criteria \(K=\{k_{1}, k_{2}, \cdots , k_{n}\}\) n \((n\ge 2)\) as a discrete set of m and a finite set of n, respectively. Usually, the criteria set K is divided into two disjoint sets as \(K_{i}\) and \(K_{ii}\). \(K_{i}\) and \(K_{ii}\) refer to a collection of benefit criteria and cost criteria, respectively, which means a larger value better performance and a larger value worse performance, respectively. The assessment rating of alternative \(d_{i} \in D\) with respect to criterion \(k_{j} \in K\) is stated as a PF value \(p_{ij}=(u_{ij}, v_{ij})\). \(u_{ij}\) and \(v_{ij}\) indicate the degrees that alternative \(d_{i}\) satisfies and dissatisfies, respectively, criterion \(k_{j}\). The hesitation degree \(\pi _{ij}\) that corresponds to each \(p_{ij}\) is calculated as \(\pi _{ij}=\sqrt{1-(u_{ij})^{2}-(v_{ij})^{2}}\). Describing this situation, the PF matrix is as follows:

$$\begin{aligned} p_{ij}= \begin{bmatrix} (u_{11}, v_{11}) &{} (u_{12}, v_{12})&{} \cdots &{} (u_{1n}, v_{1n})\\ (u_{21}, v_{21}) &{} (u_{22}, v_{22})&{} \cdots &{} (u_{2n}, v_{2n})\\ \vdots &{} \vdots &{} \vdots &{} \vdots \\ (u_{m1}, v_{m1}) &{} (u_{m2}, v_{m2})&{} \cdots &{} (u_{mn}, v_{mn}) \end{bmatrix} \end{aligned}$$
(4.1)

In the PF environment, the importance weight of criterion \(k_{j} \in K\) is stated as a PF value \(w_{j}=({\overline{w}}_{j}, {\underline{w}}_{j})\), where \({\overline{w}}_{j}\) and \({\underline{w}}_{j}\) denote the degrees of importance and unimportance, respectively, of \(k_{j}\). Then the hesitation degree \(\tau _{j}=\sqrt{1-{\overline{w}}_{j}^{2}-{\underline{w}}_{j}^{2}}\).

4.1 The displaced and fixed ideal IF values

The appropriate reference points need to be designated for the PF decision process. These points will help define displaced and fixed remoteness indices. In other words, locating the PIPF and NIPF values from two different perspectives of displaced and fixed ideals is the main problem. Here, PIPF and NIPF values will be used as points of reference because these points can effectively form anchored judgments in subjective decision-making processes. To be more precise, for the anchored judgments, the displaced PIPF and NIPF values, which consist of all the best and worst criterion values, respectively, will be defined.

Firstly in light of anchored judgments with displaced ideals, the characterization of the PIPF and NIPF values with respect to each criterion can be identified according to all of the assessment ratings in the PF decision matrix p. This means that the PIPF and NIPF values can be often displaced because they are sensitive to changes in the available set of candidate alternatives.

Definition 4.1

Take a PF decision matrix \(p=[p_{ij}]_{(m\times n)}\). For the each criterion \(k_{j} \in K\) \((K=K_{i} \cup K_{ii})\), the displaced positive-ideal and negative-ideal PF values are defined as follows:

$$\begin{aligned} p(+)_{j}= & {} (u(+)_{j}, v(+)_{j})\nonumber \\= & {} \left\{ \begin{array}{ccl} \left( \max _{i=1}^{m}u_{ij}, \min _{i=1}^{m}v_{ij}\right), &{} if\; k_{j} \in K_{i}\\ \left( \min _{i=1}^{m}u_{ij}, \max _{i=1}^{m}v_{ij}\right), &{} if\; k_{j} \in K_{ii} \end{array}\right. \end{aligned}$$
(4.2)
$$\begin{aligned} p(-)_{j}= & {} (u(-)_{j}, v(-)_{j})\nonumber \\= & {} \left\{ \begin{array}{ccl} \left( \min _{i=1}^{m}u_{ij}, \max _{i=1}^{m}v_{ij}\right), &{} if\; k_{j} \in K_{i}\\ \left( \max _{i=1}^{m}u_{ij}, \min _{i=1}^{m}v_{ij}\right), &{} if\; k_{j} \in K_{ii} \end{array}\right. \end{aligned}$$
(4.3)

and also

$$\begin{aligned} \pi (+)_{j}= & {} \sqrt{1-(u(+)_{j})^{2}-(v(+)_{j})^{2}} \end{aligned}$$
(4.4)
$$\begin{aligned} \pi (-)_{j}= & {} \sqrt{1-(u(-)_{j})^{2}-(v(-)_{j})^{2}} \end{aligned}$$
(4.5)

Consider anchored judgments with fixed ideals. It is convenient to designate (1, 0) and (0, 1) as PIPF and NIPF values, respectively, with respect to benefit criteria, because (1, 0) and (0, 1) are the top and bottom elements, respectively, of the lattice \(({\mathscr {L}}, \le _{{\mathscr {L}}})\). Conversely, (0, 1) and (1, 0) can be considered the PIPF and NIPF values, respectively, with respect to cost criteria for simplicity. These ideal PF values are called fixed ideal PF values.

Definition 4.2

Take a PF decision matrix \(p=[p_{ij}]_{(m\times n)}\). For the each criterion \(k_{j} \in K\) \((K=K_{i} \cup K_{ii})\), the displaced PIPF and NIPF values are defined as follows:

$$\begin{aligned} p(\widehat{+})_{j}= & {} (u(\widehat{+})_{j}, v(\widehat{+})_{j})= \left\{ \begin{array}{ccl} (1,0) &, &{} if\; k_{j} \in K_{i}\\ (0,1) &, &{} if\; k_{j} \in K_{ii} \end{array}\right. \end{aligned}$$
(4.6)
$$\begin{aligned} p(\widehat{-})_{j}= & {} (u(\widehat{-})_{j}, v(\widehat{-})_{j})=\left\{ \begin{array}{ccl} (0,1) &{}, &\quad{} if\; k_{j} \in K_{i}\\ (1,0) &, &{} if\; k_{j} \in K_{ii} \end{array}\right. \end{aligned}$$
(4.7)

and also \(\pi (\widehat{+})_{j}=\pi (\widehat{-})_{j}=0\).

4.2 Remoteness index (RI)

In this subsection, it will be constructed the displaced and fixed RIs.

Usually, the larger and the smaller the assessment rating \(p_{ij}\) is, the greater the preference is for the criteria \(k_{j} \in K_{i}\) and \(k_{j} \in K_{ii}\), respectively. Therefore the displaced PIPF value \(p(+)_{j}\) and the fixed PIPF value \(p(\widehat{+})_{j}\) can be deemed the most favorable values for \(k_{j}\) on the basis of anchored judgments with displaced and fixed ideals, respectively. In this way, the characterization of the proposed displaced RI can be reasonably designated via the generalized distance measure \(U_{G}(p_{ij}, p(+)_{j})\). The smaller \(U_{G}(p_{ij}, p(+)_{j})\) value is, the closer to \(p(+)_{j}\) it is, and the greater the preference is for each assessment rating \(p_{ij}\). In a similar manner, the characterization of the proposed fixed RI can be defined using the generalized distance measure \(U_{G}(p_{ij}, p(\widehat{+})_{j})\). Moreover, the smaller \(U_{G}(p_{ij}, p(\widehat{+})_{j})\) value is, the closer \(p(\widehat{+})_{j}\) it is, and the greater the preference is for each assessment rating \(p_{ij}\).

However, it is worthwhile to mention that the \(U_{G}(p_{ij}, p(+)_{j})\) values for all \(k_{j} \in K\) do not have a consistent upper bound in most situations. The main reason is that the displaced PIPF values are most favorable PF values with all of the currently considered assessment rating with respect to each criterion; thus, they are usually different among criteria. As a result, the maximal possible value of the generalized PF distances, i.e., the \(U_{G}(p(-)_{j}, p(+)_{j})\) values, can be regarded as the upper bounds of the \(U_{G}(p_{ij}, p(+)_{j})\) values for all \(k_{j} \in K\).

Definition 4.3

Let \(p_{ij}, p(+)_{j}, p(-)_{j}\) denote an assessment rating, the displaced PIPF value, and the displaced NIPF value, respectively, in the PF decision matrix \(p=[p_{ij}]_{(m\times n)}\). Moreover, without loss of generality, assume that \(p(-)_{j}\ne p(+)_{j}\) for all \(k_{j} \in K\). The displaced RI \(RI(p_{ij})\) of \(p_{ij}\) based on the generalized distance measure \(U_{G}\) is defined as follows:

$$\begin{aligned} RI(p_{ij})=\frac{U_{G}(p_{ij}, p(+)_{j})}{U_{G}(p(-)_{j}, p(+)_{j})} \end{aligned}$$
(4.8)

Theorem 4.4

Choose the three assessment ratings \(p_{ij}, p_{i_{1}j}, p_{i_{2}j}\). The displaced RI satisfies the following properties:

  1. i.

    \(RI(p_{ij})=0\) \(\Leftrightarrow\) \(p_{ij}=p(+)_{j}\);

  2. ii.

    \(RI(p_{ij})=1\) \(\Leftrightarrow\) \(p_{ij}=p(-)_{j}\);

  3. iii.

    \(0 \le RI(p_{ij}) \le 1\);

  4. iv.

    For each \(k_{j} \in K_{i}\), \(RI(p_{i_{2}j})\le RI(p_{i_{1}j})\), if \(p_{i_{1}j} \le _{{\mathscr {L}}} p_{i_{2}j}\);

  5. v.

    For each \(k_{j} \in K_{ii}\), \(RI(p_{i_{1}j})\le RI(p_{i_{2}j})\), if \(p_{i_{1}j} \le _{{\mathscr {L}}} p_{i_{2}j}\).

Proof

  1. (i)

    \(\Rightarrow :\) If \(RI(p_{ij})=0\), then it implies that \(U_{G}(p_{ij}, p(+)_{j})=0\). According to the separability condition of Theorem 3.3, it is obtained that \(p_{ij}=p(+)_{j}\).

    \(\Leftarrow\) If \(p_{ij}=p(+)_{j}\), it is obvious that \(RI(p_{ij})=0\), because of the reflexivity condition \((U_{1})\) of Theorem 3.3.

  2. (ii)

    \(\Rightarrow :\) If \(RI(p_{ij})=1\), then it implies that \(U_{G}(p_{ij}, p(+)_{j})=U_{G}(p(-)_{j}, p(+)_{j})\). It follows that \(p_{ij}=p(-)_{j}\).

    \(\Leftarrow\) If \(p_{ij}=p(-)_{j}\), it is trivial to obtain that \(RI(p_{ij})=1\) by employing (4.8).

  3. (iii)

    Applying (4.2) and (4.3), it is obvious that \(u(-)_{j}\le u_{ij} \le u(+)_{j}\) and \(v(-)_{j}\ge v_{ij} \ge v(+)_{j}\) for each \(k_{j} \in K_{i}\). Moreover \(u(+)_{j}\le u_{ij} \le u(-)_{j}\) and \(v(+)_{j}\ge v_{ij} \ge v(-)_{j}\) for each \(k_{j} \in K_{ii}\). Thus, one can easily obtain \(p(-)_{j} \le _{{\mathscr {L}}} p_{ij} \le _{{\mathscr {L}}} p(+)_{j}\) and \(p(+)_{j} \le _{{\mathscr {L}}} p_{ij} \le _{{\mathscr {L}}} p(-)_{j}\) for \(k_{j} \in K_{i}\) and \(k_{j} \in K_{ii}\), respectively. According to \((U_{5})\) of Theorem 3.3, it is known that \(U_{G}(p_{ij}, p(+)_{j})\le (U_{G}(p(-)_{j}, p(+)_{j})\) and \(U_{G}(p(+)_{j}, p_{ij})\le (U_{G}(p(+)_{j}, p(-)_{j})\) for \(k_{i} \in K_{i}\) and \(k_{i} \in K_{ii}\), respectively. Next, combining the condition \(U_{G}(p_{ij}, p(+)_{j})\le (U_{G}(p(-)_{j}, p(+)_{j})\) and the boundedness property of Theorem 3.3, one can obtain \(0\le RI(p_{ij}) \le 1\).

  4. (iv)

    From the proof of (iii.), for each \(k_{i} \in K_{i}\), it is easily seen that \(p_{i_{1}j} \le _{{\mathscr {L}}} p_{i_{2}j}\le _{{\mathscr {L}}} p(+)_{j}\). According to the condition (v.) of Theorem 3.3, the condition \(p_{i_{1}j} \le _{{\mathscr {L}}} p_{i_{2}j}\le _{{\mathscr {L}}} p(+)_{j}\) implies \(U_{G}(p_{i_{2}j}, p(+)_{j})\le (U_{G}(p_{i_{1}j}, p(+)_{j})\). Note that in the denominator of (4.8), one has \(U_{G}(p(-)_{j}, p(+)_{j})>0\) because of the assumption of \(p(-)_{j} \ne p(+)_{j}\) for \(k_{i} \in K_{i}\). Thus, it is easily proven that \(RI(p_{i_{2}j})\le RI(p_{i_{1}j})\) for \(k_{i} \in K_{i}\).

The proof of condition (v) can be proven as similar to (iv). \(\square\)

Example 4.5

Take the set of candidate alternatives \(D=\{d_{1}, d_{2}, d_{3}\}\), and the criteria set \(K=\{k_{1}, k_{2}\}\), where \(K_{i}=\{k_{1}\}\) and \(K_{ii}=\{k_{2}\}\). The PF decision matrix can be given as follows:

$$\begin{aligned} p=[p_{ij}]_{3\times 2}= \begin{bmatrix} (0.54, 0.67) &{} (071, 0.43)\\ (0.88, 0.35) &{} (0.62, 0.49)\\ (0.50, 0.56) &{} (0.85, 0.36) \end{bmatrix} \end{aligned}$$

The displaced positive-ideal PF values are \(p(+)_{1}=(0.88, 0.35)\) and \(p(+)_{2}=(0.62, 0.49)\). Further, the displaced negative-ideal PF values are \(p(-)_{1}=(0.50, 0.67)\) and \(p(-)_{2}=(0.85, 0.36)\). Using the BPFSS and from 3.1, compute the Hamming distance (generalized distance measure for \(\lambda =1\)), which takes the form

$$\begin{aligned}&U_{H}(p_{ij}, p(+)_{j})\nonumber \\&\quad = \frac{1}{5}\left\{ \begin{array}{l} \bigg [|u_{(p_{ij})}-u_{(p(+)_{j})}|+|v_{(p_{ij})}-v_{(p(+)_{j})}|+|\pi _{(p_{ij})}-\pi _{(p(+)_{j})}|\\ +|r_{(p_{ij})}-r_{(p(+)_{j})}|+|d_{(p_{ij})}-d_{(p(+)_{j})}||\bigg ] \end{array} \right\} \end{aligned}$$
(4.9)

From (4.9), \(\Omega _{H}=(p(-)_{1}, p(+)_{1})=0.58\) and \(\Omega _{H}=(p(-)_{2}, p(+)_{2})=0.62\). \(RI(p_{11})=0.1205\), \(RI(p_{12})=0.176\), \(RI(p_{21})=0.101\), \(RI(p_{22})=0.13\), \(RI(p_{31})=0.155\), \(RI(p_{32})=0.127\).

Definition 4.6

Let \(p_{ij}, p(\widehat{+})_{j}, p(\widehat{-})_{j}\) denote an assessment rating, the fixed PIPF value, and the fixed NIPF value, respectively, in the PF decision matrix \(p=[p_{ij}]_{(m\times n)}\). The fixed RI \({\widehat{RI}}(p_{ij})\) of \(p_{ij}\) based on the generalized distance measure \(U_{G}\) is defined as follows:

$$\begin{aligned} {\widehat{RI}}(p_{ij})=\frac{U_{G}(p_{ij}, p(\widehat{+})_{j})}{U_{G}(p(\widehat{-})_{j}, p(\widehat{+})_{j})} \end{aligned}$$
(4.10)

Theorem 4.7

Choose the three assessment ratings \(p_{ij}, p_{i_{1}j}, p_{i_{2}j}\). The fixed remoteness index satisfies the following properties:

  1. i.

    \({\widehat{RI}}(p_{ij})=0\) \(\Leftrightarrow\) \(p_{ij}=p(\widehat{+})_{j}\);

  2. ii.

    \({\widehat{RI}}(p_{ij})=1\) \(\Leftrightarrow\) \(p_{ij}=p(\widehat{-})_{j}\);

  3. iii.

    \(0 \le RI({\widehat{d}})(p_{ij}) \le 1\);

  4. iv.

    For each \(k_{j} \in K_{i}\), \({\widehat{RI}}(p_{i_{2}j})\le {\widehat{RI}}(p_{i_{1}j})\), if \(p_{i_{1}j} \le _{{\mathscr {L}}} p_{i_{2}j}\);

  5. v.

    For each \(k_{j} \in K_{ii}\), \({\widehat{RI}}(p_{i_{1}j})\le {\widehat{RI}}(p_{i_{2}j})\), if \(p_{i_{1}j} \le _{{\mathscr {L}}} p_{i_{2}j}\).

This theorem can be proved similar to Theorem 4.4.

Example 4.8

Let DK and \(p=p_{ij}\) be as in Example 4.5. From Definition 4.2, the fixed positive-ideal PF values are \(p(\widehat{+})_{1}=(1,0)\) and \(p(\widehat{+})_{2}=(0,1)\). In the same way, the fixed negative-ideal PF values \(p(\widehat{-})_{1}=(0,1)\) and \(p(\widehat{-})_{2}=(1,0)\). Take \(\lambda =1\). Then, \(\Omega _{H}=(p(\widehat{-})_{1}, p(\widehat{+})_{1})=(p(\widehat{-})_{2}, p(\widehat{+})_{2})=0.6\). \({\widehat{RI}}(p_{11})=0.1033\), \({\widehat{RI}}(p_{12})=0.152\), \({\widehat{RI}}(p_{21})=0.098\), \({\widehat{RI}}(p_{22})=0.14\), \({\widehat{RI}}(p_{31})=0.1035\), \({\widehat{RI}}(p_{32})=0.15\).

As seen in this example, the oppositions between the remoteness indices become less discernible when anchoring subjective judgments with fixed ideal PF values. This situation is not so noticeable in Example 4.5.

4.3 Precise and PF importance weights

When decision-making in life, not all of the characteristics that affect these decisions are evaluated equally. Therefore, decision-making is determining alternatives according to the preferences and values of the decision-maker and choosing among them. Each situation has a weight for this choice. In DM operates, weights are assigned to the features. Two types of weights will be used in this study:

Example 4.9

  1. (i)

    The precise importance weights are obtained by expectation score function. The weight vector \(\omega\) calculated with \(\frac{\mathscr {ES}_{f(z)}}{\sum _{z}\mathscr {ES}_{f}(z)}\).

    Choose \(f=\{d_{1}, d_{2}, d_{3}\}=\{(k_{1}, 0.68, 0.33), (k_{2}, 0.75, 0.28), (k3, 0,84, 0.39)\}\). Using the (2.9), \(\mathscr {ES}_{f(d_{1})}=0.21435\), \(\mathscr {ES}_{f(d_{2})}=0.17955\), \(\mathscr {ES}_{f(d_{3})}=0.07115\) and \(\sum _{i=1}^{3}\mathscr {ES}_{f(d_{i})}=0.46505\). Then the weights are \(\omega _{1}=0.461\), \(\omega _{2}=0.386\) and \(\omega _{3}=0.153\) \((\sum _{i=1}^{3}\omega _{i}=1)\).

  2. (ii)

    The PF Importance Weights: The PF importance weights of criteria \(e_{j} \in \Sigma\) are displaced as the PFV \(\omega _{j}=({\overline{\omega }}_{j}, {\underline{\omega }}_{j})\), where \({\overline{\omega }}_{j}\) and \({\underline{\omega }}_{j}\) represent the importance and unimportance degrees of the criteria \(e_{j} \in \Sigma\), respectively. Throughout the DM process, the expert group examines the steps of the process and records their evaluations as a parametric Pythagorean fuzzy set (PPFS) at the end of the process. These PPFS built by experts will be used as PF importance weights in the DM process.

4.4 Remoteness index-based new VIKOR methods

New multiple criteria ranking indices for \(RI(p_{ij})\) and \({\widehat{RI}}(p_{ij})\) measures of remoteness to the displaced and fixed ideal PF values, respectively, will be offered. In this way, the remoteness index-based new VIKOR method will be obtained. These multiple criteria ranking indices are: remoteness-based group utility indices, individual regret indices, and compromise indices. They can provide a mechanism to trade off a maximum group utility for the majority and a minimum individual regret for the opponent and can be used as an aggregating function for a compromise ranking among alternatives.

The smaller the value obtained in the displaced distance index \(RI(p_{ij})\), the more priority the assessment rating of the \(p_{ij}\) will be. That is, the criterion-wise priority relationship among alternatives can be ordinarily identified according to ascending order of each \(RI(p_{ij})\) value. This is because displaced ideal PF values are to facilitate anchored judgments. In this study, a displaced remoteness-based compromise index by combining the suggested new definitions of a group utility index and an individual regret index is offered. The reason for this is to synthesize criteria-based priority relationships to obtain a compromise ranking on all criteria.

Definition 4.10

The displaced remoteness-based group utility index, the displaced remoteness-based individual regret index and the displaced remoteness-based compromise index are defined as follows:

$$\begin{aligned} \mathbb {\widetilde{S}}(d_{i})= & {} \sum _{j=1}^{n}\left( RI(p_{ij}).\omega _{j}\right) , (k_{j} \in K), \end{aligned}$$
(4.11)
$$\begin{aligned} \mathbb {\widetilde{R}}(d_{i})= & {} \max _{j=1}^{n}\left( RI(p_{ij}).\omega _{j}\right) , (k_{j} \in K), \end{aligned}$$
(4.12)
$$\begin{aligned} \mathbb {\widetilde{Q}}(d_{i})= & {} \delta .\frac{\mathbb {S}(d_{i})-\min _{i^{'}=1}^{m}\mathbb {S}(d_{i^{'}})}{\max _{i^{'}=1}^{m}\mathbb {S}(d_{i^{'}})-\min _{i^{'}=1}^{m}\mathbb {S}(d_{i^{'}})}\nonumber \\&+(1-\delta ).\frac{\mathbb {R}(d_{i})-\min _{i^{'}=1}^{m}\mathbb {R}(d_{i^{'}})}{\max _{i^{'}=1}^{m}\mathbb {R}(d_{i^{'}})-\min _{i^{'}=1}^{m}\mathbb {R}(d_{i^{'}})} \end{aligned}$$
(4.13)

where \(p=[p_{ij}]_{m\times n}\) and \(\omega _{j}\) are PF decision matrix and a set of precise importance weights, respectively, and \(\delta \in [0,1]\). Further, the conditions \(0\le \omega _{j} \le 1\) and \(\sum _{j=1}^{n}\omega _{j}=1\) are the normalization conditions of precise weights.

Theorem 4.11

The ranking indices \(\mathbb {\widetilde{S}}(d_{i}), \mathbb {\widetilde{R}}(d_{i}), \mathbb {\widetilde{Q}}(d_{i})\) satisfy the following properties:

  1. i.

    \(0\le \mathbb {\widetilde{S}}(d_{i}) \le 1\), \(0\le \mathbb {\widetilde{R}}(d_{i}) \le 1\) and \(0\le \mathbb {\widetilde{Q}}(d_{i}) \le 1\);

  2. ii.

    If \(p_{i_{2}j} \le _{{\mathscr {L}}} p_{i_{1}j}\) and \(p_{i_{1}j} \le _{{\mathscr {L}}} p_{i_{2}j}\) for all \(k_{i} \in K_{i}\) and \(k_{i} \in K_{ii}\), then \(\mathbb {\widetilde{S}}(d_{i_{1}}) \le \mathbb {\widetilde{S}}(d_{i_{2}})\), \(\mathbb {\widetilde{R}}(d_{i_{1}}) \le \mathbb {\widetilde{R}}(d_{i_{2}})\) and \(\mathbb {\widetilde{Q}}(d_{i_{1}}) \le \mathbb {\widetilde{Q}}(d_{i_{2}})\), respectively.

Proof

  1. (i)

    Since \(0 \le RI(p_{ij}) \le 1\)(the condition of (iii) of Theorem 4.4) and normalization conditions of precise weights, \(0\le \mathbb {\widetilde{S}}(d_{i}) \le 1\), \(0\le \mathbb {\widetilde{R}}(d_{i}) \le 1\) hold. Additionally, since \(\mathbb {\widetilde{Q}}(d_{i})\) is merely a linear combination of the normalized values of \(\mathbb {\widetilde{S}}(d_{i}), \mathbb {\widetilde{R}}(d_{i})\), one can easily infer that \(0\le \mathbb {\widetilde{Q}}(d_{i}) \le 1\).

  2. (ii)

    From the condition of (iv) of Theorem 4.4, since \(p_{i_{2}j} \le _{{\mathscr {L}}} p_{i_{1}j}\) for each \(k_{i} \in K_{i}\), then \(\mathbb {R}(p_{i_{1}j}) \le \mathbb {R}(p_{i_{2}j})\). Since \(p_{i_{1}j} \le _{{\mathscr {L}}} p_{i_{2}j}\) for each \(k_{i} \in K_{ii}\), \(\mathbb {R}(p_{i_{1}j}) \le \mathbb {R}(p_{i_{2}j})\) holds from the condition (v) of Theorem 4.4. That is, \(\mathbb {R}(p_{i_{1}j}) \le \mathbb {R}(p_{i_{2}j})\) is fulfilled for \(k_{i} \in K\). So, the inequalities \(\mathbb {\widetilde{S}}(d_{i_{1}}) \le \mathbb {\widetilde{S}}(d_{i_{2}})\) and \(\mathbb {\widetilde{R}}(d_{i_{1}}) \le \mathbb {\widetilde{R}}(d_{i_{2}})\) and \(\mathbb {\widetilde{Q}}(d_{i_{1}}) \le \mathbb {\widetilde{Q}}(d_{i_{2}})\) are valid.

\(\square\)

Next, consider an MCDA problem involving the PF decision matrix and PF importance weights. In a similar manner, this paper defines the multiple criteria ranking indices \(\mathbb {S}(d_{i}), \mathbb {R}(d_{i}), \mathbb {Q}(d_{i})\) under the condition of PF preference information. Analogous to Definition 4.10 , the product of \(RI(p_{ij})\) multiplied by \(\omega _{j}\) must be calculated before determining these multiple criteria ranking indices. However, a difficulty in comparing the multiple criteria ranking indices among alternatives is encountered because \(RI(p_{ij}).\omega _{j}\) s a PF value when the PF importance weight \(\omega _{j}\) is taken into account. Here, a score function approach to acquire a comparable value of \(RI(p_{ij}).\omega _{j}\) and then to identify the multiple criteria ranking indices is offered. It is worth noting that the score function approach leads to different ranges among \(\mathbb {S}(d_{i}), \mathbb {R}(d_{i}), \mathbb {Q}(d_{i})\) because the score function defined as in (2.4) is between \(-1\) and 1.

Definition 4.12

The displaced remoteness-based group utility index, the displaced remoteness-based individual regret index and the displaced remoteness-based compromise index are defined as follows:

$$\begin{aligned}&\mathbb {S}(d_{i})=\sum _{j=1}^{n}\mathscr {ES}\left( RI(p_{ij}).\omega _{j}\right) \nonumber \\&\quad =\sum _{j=1}^{n}\left\{ \frac{1}{2}\left[ 2-\left( 1-({\underline{\omega }}_{j})^{2}\right) ^{RI(p_{ij})}-({\overline{\omega }}_{j})^{2.RI(p_{ij})}\right] \right\} , (k_{j} \in K), \end{aligned}$$
(4.14)
$$\begin{aligned}&\mathbb {R}(d_{i})=\max _{j=1}^{n}\mathscr {ES}\left( RI(p_{ij}).\omega _{j}\right) \nonumber \\&\quad =\max _{j=1}^{n}\left\{ \frac{1}{2}\left[ 2-\left( 1-({\underline{\omega }}_{j})^{2}\right) ^{RI(p_{ij})}-({\overline{\omega }}_{j})^{2.RI(p_{ij})}\right] \right\} , (k_{j} \in K), \end{aligned}$$
(4.15)
$$\begin{aligned}&\mathbb {Q}(d_{i})=\delta .\frac{\mathbb {S}(d_{i})-\min _{i^{'}=1}^{m}\mathbb {S}(d_{i^{'}})}{\max _{i^{'}=1}^{m}\mathbb {S}(d_{i^{'}})-\min _{i^{'}=1}^{m}\mathbb {S}(d_{i^{'}})}\nonumber \\&\qquad\qquad +(1-\delta ).\frac{\mathbb {R}(d_{i})-\min _{i^{'}=1}^{m}\mathbb {R}(d_{i^{'}})}{\max _{i^{'}=1}^{m}\mathbb {R}(d_{i^{'}})-\min _{i^{'}=1}^{m}\mathbb {R}(d_{i^{'}})} \end{aligned}$$
(4.16)

where \(p=[p_{ij}]_{m\times n}\) and \(\omega _{j}\) are PF decision matrix and a set of PF importance weights, respectively, and \(\delta \in [0,1]\).

Theorem 4.13

The ranking indices \(\mathbb {S}(d_{i}), \mathbb {R}(d_{i}), \mathbb {Q}(d_{i})\) satisfy the following properties:

  1. i.

    \(-n\le \mathbb {S}(d_{i}) \le n\), \(-1\le \mathbb {R}(d_{i}) \le 1\) and \(0\le \mathbb {Q}(d_{i}) \le 1\);

  2. ii.

    If \(p_{i_{2}j} \le _{{\mathscr {L}}} p_{i_{1}j}\) and \(p_{i_{1}j} \le _{{\mathscr {L}}} p_{i_{2}j}\) for all \(k_{i} \in K_{i}\) and \(k_{i} \in K_{ii}\), then \(\mathbb {S}(d_{i_{1}}) \le \mathbb {S}(d_{i_{2}})\), \(\mathbb {R}(d_{i_{1}}) \le \mathbb {R}(d_{i_{2}})\) and \(\mathbb {Q}(d_{i_{1}}) \le \mathbb {Q}(d_{i_{2}})\), respectively.

Proof

According to the (2.9), the inequality \(-1\le \mathscr {ES}(RI(p_{ij}).\omega _{j}) \le 1\) \((\forall k_{j}\in K)\) holds. Then, \(-1\le \max _{j=1}^{n}\mathscr {ES}(RI(p_{ij}).\omega _{j}) \le 1\) is valid. From here, \(-1\le \mathbb {R}(d_{i}) \le 1\) is obtained. At the same time, \(-n\le \sum _{j=1}^{n}\mathscr {ES}(RI(p_{ij}).\omega _{j}) \le n \Rightarrow -n\le \mathbb {R}(d_{i}) \le n\). Consider (4.16). Since \(\mathbb {Q}(d_{i})\) is merely a linear combination of the normalized values \(\mathbb {S}(d_{i})\), and \(\mathbb {R}(d_{i})\), then, it is obvious that it will be \(0\le \mathbb {Q}(d_{i}) \le 1\).

The proof of the condition (ii.) can be easily get from part (ii.) of Theorem 4.11\(\square\)

For fixed ideal PF values, new definitions are given as follows:

Definition 4.14

The fixed remoteness-based group utility index, the fixed remoteness-based individual regret index and the fixed remoteness-based compromise index are defined as follows:

$$\begin{aligned} \mathbb {\widetilde{SF}}(d_{i})= & {} \sum _{j=1}^{n}\left( {\widehat{RI}}(p_{ij}).\omega _{j}\right) , (k_{j} \in K), \end{aligned}$$
(4.17)
$$\begin{aligned} \mathbb {\widetilde{RF}}(d_{i})= & {} \max _{j=1}^{n}\left( {\widehat{RI}}(p_{ij}).\omega _{j}\right) , (k_{j} \in K), \end{aligned}$$
(4.18)
$$\begin{aligned} \mathbb {\widetilde{QF}}(d_{i})= & {} \delta .\frac{\mathbb {\widetilde{SF}}(d_{i})-\min _{i^{'}=1}^{m}\mathbb {\widetilde{SF}}(d_{i^{'}})}{\max _{i^{'}=1}^{m}\mathbb {\widetilde{SF}}(d_{i^{'}})-\min _{i^{'}=1}^{m}\mathbb {\widetilde{SF}}(d_{i^{'}})}\nonumber \\&+(1-\delta ).\frac{\mathbb {\widetilde{RF}}(d_{i})-\min _{i^{'}=1}^{m}\mathbb {\widetilde{RF}}(d_{i^{'}})}{\max _{i^{'}=1}^{m}\mathbb {\widetilde{RF}}(d_{i^{'}})-\min _{i^{'}=1}^{m}\mathbb {\widetilde{RF}}(d_{i^{'}})} \end{aligned}$$
(4.19)

where \(p=[p_{ij}]_{m\times n}\) and \(\omega _{j}\) are PF decision matrix and a set of precise importance weights, respectively, and \(\delta \in [0,1]\). Further, the conditions \(0\le \omega _{j} \le 1\) and \(\sum _{j=1}^{n}\omega _{j}=1\) are the normalization conditions of precise weights.

Definition 4.15

The fixed remoteness-based group utility index, the fixed remoteness-based individual regret index and the fixed remoteness-based compromise index are defined as follows:

$$\begin{aligned}&\mathbb {SF}(d_{i})=\sum _{j=1}^{n}\mathscr {ES}\left( {\widehat{RI}}(p_{ij}).\omega _{j}\right) \nonumber \\&\quad =\sum _{j=1}^{n}\left\{ \frac{1}{2}\left[ 2-\left( 1-({\underline{\omega }}_{j})^{2}\right) ^{{\widehat{RI}}(p_{ij})}-({\overline{\omega }}_{j})^{2.{\widehat{RI}}(p_{ij})}\right] \right\} , (k_{j} \in K), \end{aligned}$$
(4.20)
$$\begin{aligned}&\mathbb {RF}(d_{i})=\max _{j=1}^{n}\mathscr {ES}\left( {\widehat{RI}}(p_{ij}).\omega _{j}\right) \nonumber \\&\quad =\max _{j=1}^{n}\left\{ \frac{1}{2}\left[ 2-\left( 1-({\underline{\omega }}_{j})^{2}\right) ^{{\widehat{RI}}(p_{ij})}-({\overline{\omega }}_{j})^{2.{\widehat{RI}}(p_{ij})}\right] \right\} , (k_{j} \in K), \end{aligned}$$
(4.21)
$$\begin{aligned}&\mathbb {QF}(d_{i})=\delta .\frac{\mathbb {SF}(d_{i})-\min _{i^{'}=1}^{m}\mathbb {SF}(d_{i^{'}})}{\max _{i^{'}=1}^{m}\mathbb {SF}(d_{i^{'}})-\min _{i^{'}=1}^{m}\mathbb {SF}(d_{i^{'}})}\nonumber \\&\qquad +(1-\delta ).\frac{\mathbb {RF}(d_{i})-\min _{i^{'}=1}^{m}\mathbb {RF}(d_{i^{'}})}{\max _{i^{'}=1}^{m}\mathbb {RF}(d_{i^{'}})-\min _{i^{'}=1}^{m}\mathbb {RF}(d_{i^{'}})} \end{aligned}$$
(4.22)

where \(p=[p_{ij}]_{m\times n}\) and \(\omega _{j}\) are PF decision matrix and a set of PF importance weights, respectively, and \(\delta \in [0,1]\).

Theorem 4.16

The ranking indices \(\mathbb {\widetilde{SF}}(d_{i}), \mathbb {\widetilde{RF}}(d_{i}), \mathbb {\widetilde{QF}}(d_{i})\) satisfy the following properties:

  1. i.

    \(0\le \mathbb {\widetilde{SF}}(d_{i}) \le 1\), \(0\le \mathbb {\widetilde{RF}}(d_{i}) \le 1\) and \(0\le \mathbb {\widetilde{QF}}(d_{i}) \le 1\);

  2. ii.

    If \(p_{i_{2}j} \le _{{\mathscr {L}}} p_{i_{1}j}\) and \(p_{i_{1}j} \le _{{\mathscr {L}}} p_{i_{2}j}\) for all \(k_{i} \in K_{i}\) and \(k_{i} \in K_{ii}\), then \(\mathbb {\widetilde{SF}}(d_{i_{1}}) \le \mathbb {\widetilde{SF}}(d_{i_{2}})\), \(\mathbb {\widetilde{RF}}(d_{i_{1}}) \le \mathbb {\widetilde{RF}}(d_{i_{2}})\) and \(\mathbb {\widetilde{QF}}(d_{i_{1}}) \le \mathbb {\widetilde{QF}}(d_{i_{2}})\), respectively.

Theorem 4.17

The ranking indices \(\mathbb {SF}(d_{i}), \mathbb {RF}(d_{i}), \mathbb {QF}(d_{i})\) satisfy the following properties:

  1. i.

    \(-n\le \mathbb {SF}(d_{i}) \le n\), \(-1\le \mathbb {RF}(d_{i}) \le 1\) and \(0\le \mathbb {QF}(d_{i}) \le 1\);

  2. ii.

    If \(p_{i_{2}j} \le _{{\mathscr {L}}} p_{i_{1}j}\) and \(p_{i_{1}j} \le _{{\mathscr {L}}} p_{i_{2}j}\) for all \(k_{i} \in K_{i}\) and \(k_{i} \in K_{ii}\), then \(\mathbb {SF}(d_{i_{1}}) \le \mathbb {SF}(d_{i_{2}})\), \(\mathbb {RF}(d_{i_{1}}) \le \mathbb {RF}(d_{i_{2}})\) and \(\mathbb {QF}(d_{i_{1}}) \le \mathbb {QF}(d_{i_{2}})\), respectively.

These theorems can be proved as Theorems 4.11 and 4.13.

The decision process coefficient is shown by \(\delta\). One can modify decision making strategy by changing the value of \(\delta\). The value of parameter \(\delta\) indicates the importance of maximum group utility, while \(1-\delta\) indicates the importance of individual regrets. In the classical VIKOR method, the higher the value of \(\delta\) (when \(\delta >0.5\)), the compromise ranking procedure is categorized as the procedure with voting by majority. The compromise ranking procedure is categorized as the procedure with veto when \(\delta <0.5\). The consensus is achieved in the compromise ranking procedure at \(\delta =0.5\).

Example 4.18

Consider the (FMf) of Example 3.5. Here, Let \(M_{1}=\{e_{1}, e_{2}\}\) and \(M_{2}=\{e_{3}\}\) (\(M=M_{1} \cup M_{2}\)). For \(\delta =0.5\) and \(\lambda =1\),

(i).:

Using the expectation score function, calculate the precise importance weights: \(\omega _{1}=0.158811626\), \(\omega _{2}=0.493797108\), \(\omega _{3}=0.347391266\). Choose \(\lambda =1\) and \(\delta =0.5\). From Definition 4.10, results are obtained: \(RI(p_{11})=0\), \(RI(p_{12})=0\), \(RI(p_{13})=0.911\), \(RI(p_{21})=0.0406\), \(RI(p_{22})=0.98\), \(RI(p_{23})=0\) and

$$\begin{aligned} \mathbb {\widetilde{S}}(x_{1})= & {} RI(p_{11}).\omega _{1}+RI(p_{12}).\omega _{2}+RI(p_{13}).\omega _{3}=0.32,\\ \mathbb {\widetilde{S}}(x_{2})= & {} RI(p_{21}).\omega _{1}+RI(p_{22}).\omega _{2}+RI(p_{23}).\omega _{3}=0.659 \end{aligned}$$

and \(\mathbb {\widetilde{R}}(x_{1})=0.32\), \(\mathbb {\widetilde{R}}(x_{2})=0.484\). Then, \(\mathbb {\widetilde{Q}}(x_{1})=0\) and \(\mathbb {\widetilde{Q}}(x_{2})=1\). Therefore, the rankings are obtained: \(x_{1}<x_{2}\) for \(\mathbb {\widetilde{S}}(x_{i})\), \(x_{1}<x_{2}\) for \(\mathbb {\widetilde{R}}(x_{i})\),and \(x_{1}<x_{2}\) for \(\mathbb {\widetilde{Q}}(x_{i})\). Similar calculations can be made for Definition 4.14.

(ii).:

The PF importance weights are given as \(\omega _{1}=(\underline{\omega _{1}}, \overline{\omega _{1}})=(0.65, 0.35)\) and \(\omega _{2}=(0.55, 0.50)\), \(\omega _{3}=(0.60, 0.60)\). From Definition 4.12, results are obtained: \(\widetilde{S}(x_{1})=0.467\), \(\widetilde{S}(x_{2})=1.195\), \(\widetilde{R}(x_{1})=0.467\), \(\widetilde{R}(x_{2})=0.675\) and \(\widetilde{Q}(x_{1})=0\), \(\widetilde{Q}(x_{2})=1\). Therefore, the rankings are obtained: \(x_{1}<x_{2}\) for \(\widetilde{S}(x_{i})\), \(x_{1}<x_{2}\) for \(\widetilde{R}(x_{i})\),and \(x_{1}<x_{2}\) for \(\widetilde{Q}(x_{i})\). Similar calculations can be made for Definition 4.15.

5 Quarantine situations application

In this section, firstly, using newly defined distance measurements, weights, displaced and fixed ideals, and displaced and fixed remoteness indices novel algorithms will be given. Later, these algorithms will be applied to the cognitive and behavioral development of early childhood children staying at home in the COVID-19 quarantine.

5.1 Algorithms

For the four different scenarios, algorithms are given as follows:

figure a
figure b
figure c
figure d

5.2 Quarantine application for early childhood children

In natural disasters and especially in epidemic diseases, some measures are taken to protect people from the negative effects of the situation. One of the measures that can be taken is quarantine. The quarantine is used to indicate restrictions on the activities of people or animals exposed to infectious diseases during the infectious period. Children, who are members of society and cannot be isolated from society, should be informed correctly and sufficiently to prevent them from being affected by both the biological effect and the psychological effect of the epidemic. It is extremely important to inform and raise awareness of children beforehand in order to prevent a pandemic because children may face troubles due to the long duration of natural disasters and measures such as quarantine restricting people. Children may face personal losses, collective deaths, and discomfort caused by the diseases caught in natural disasters and outbreaks. These situations can cause adversities such as stress, anxiety, depression, and behavioral disorders in children.

Children learn a lot of the information they learn through environmental stimuli. The interaction of the child with his/her environment, social relationships, other people, especially adults, plays a very important role in cognitive development. The stimuli that it is exposed to in the pandemic process direct the perception of children to the pandemic. In this case, it is clear that children will pay more attention to the pandemic, quarantine, and related stimuli. In the process, the vast majority of stimuli around children, including parents and digital media, lead their perception of COVID-19. If this perception cannot be controlled properly, a false cognition and belief in children will be inevitable.

In this study, the survey model was used. The survey was prepared to be answered on the internet. Survey questions were asked to children aged 5-6 and their families. The linguistic terms and their numeric labels are: For Questions to be asked to the child: Yes(1), maybe/some (2), no(3). For Questions to be asked to parents: too much(1), much (2), some (3), too little(4), none (5).

The survey included the following questions:

Questions to be asked to the child:

  • Do you know Corona-virus?

  • Are you afraid of Corona-virus?

  • Does Corona-virus harm people?

  • Does Corona-virus harm animals?

  • Can Corona-virus be prevented?

  • Do you think it’s nice not to go to school?

  • Are you upset that you can’t go to school?

  • Is the obligation to stay home boring?

  • Can we be protected from Corona-virus by staying at home?

  • Do you think you can go to school from now on?

Questions to be asked to parents:

  • Does your child behave anxiously after Corona-virus?

  • Is your child afraid when a conversation about Corona-virus has passed?

  • Does your child ask about Corona-virus?

  • Does your child pay more attention to cleaning after Corona-virus?

  • Has your child’s sleep pattern been impaired after Corona-virus?

  • Have there been changes in your child’s nutritional habits after Corona-virus?

  • Did your child develop undesirable behaviour after Corona-virus?

  • Is your child happy because she/he can’t go to school?

  • Has the time your child spent on the Internet after Corona-virus increased?

  • Has the time your child spent in front of the TV increased after Corona-virus?

The cognitive and behavioural distributions of questions are as follows:

For children’s cognition;

\(c_{1}\)::

Do children know about the current situation? (4 questions)

\(c_{2}\)::

Does the current situation affect children’s emotions? (4 questions)

\(c_{3}\)::

Does the current situation affect children’s thoughts? (2 questions)

For children’s behavioural;

\(d_{1}\)::

Has Corona-virus changed the basic habits of children? (3 questions)

\(d_{2}\)::

Did behaviour change occur in children after quarantine? (5 questions)

\(d_{3}\)::

Did children’s behaviour regarding information technologies increase after quarantine? (2 questions)

The target audience is children aged 5-6 in early childhood. Children of this age group are in the process of gaining skills in expressing their feelings during this period. In addition, the emotional responses of these children can be noticed by a careful observer or even an expert. This section is to investigate the effects of quarantine status due to the COVID-19 pandemic on the cognitions and behaviors of children who stay at home. For this, questions were asked to the children and their parents. Opinions of each child and each parent about the questions asked were got. The effect of quarantine on their own cognition in line with the answers given by the children and the effect of the behavior of their children in line with the observations of the parents have been revealed.

Let \(C=\{c_{1}, c_{2}, c_{3}\}\) and \(D=\{d_{1}, d_{2}, d_{3}\}\) represent the sets of the questions related to children’s cognition and the questions related to children’s behavioral, respectively. Let the sets \(M=\{m_{1}, m_{2}, m_{3}\}=\{yes, maybe, no\}\) and \(N=\{n_{1}, n_{2}, n_{3}\}=\{much, some, none\}\) represent the answer given by the children to the questions asked and answer given by the parents to the questions asked, respectively.

An expert group was formed to ask questions to children and evaluate their answers. The experts assessed the answers according to criteria. The experts gave their judgments based on their knowledge and previous statistical measures. The committee preferences are stored in the form of GPFSS. These results are shown in Tables 2 and 3. f and g in these tables represent expert opinions. Algorithm-based solutions will be made for both cognitive and behavioral developments.

Table 2 Experts’ opinion related to the cognitive development
Table 3 Experts’ opinion related to the behavioral development

As demonstrated in Theorem 3.3 and Note 1, the Hamming distance \(\Omega _{H}\) is a metric because of its reflexivity, symmetry, separability and triangle inequality. Thus, when computing the generalized distance measure \(\Omega _{G}\), this paper sets the distance parameter \(\lambda =1\) for the application. Moreover, assume that the decision-maker would like to achieve the goal of “consensus,” and thus, set \(\delta =0,5\).


Algorithm 1 for (FMf):

A simple MCDA problem concerning three alternatives and three benefit criteria for evaluating the alternatives, where \(K=K_{i}=\{m_{1}, m_{2}, m_{3}\}\) and \(K_{ii}=\emptyset\).

The values \(p(+)_{j}\) and \(p(-)_{j}\) for (FMf) are calculated by using (4.2) and (4.3):

$$\begin{aligned} p(+)_{j}= & {} \left\{ (0.84, 0.34), (0.89, 0.19), (0.87, 0.27)\right\} \\ p(-)_{j}= & {} \left\{ (0.69, 0.54), (0.72, 0.41), (0.56, 0.59)\right\} \end{aligned}$$

Using (4.9), the distance measures are computed for \(\Omega _{H}(p(-)_{j}, p(+)_{j})\) and \(\Omega _{H}(p_{ij}, p(+)_{j})\): \(\Omega _{H}(p(-)_{1}, p(+)_{1})=0.124\), \(\Omega _{H}(p(-)_{2}, p(+)_{2})=0.164\), \(\Omega _{H}(p(-)_{3}, p(+)_{3})=0.18\).

$$\begin{aligned} \Omega _{H}(p_{ij}, p(+)_{j})= \begin{bmatrix} 0.106 &{} 0.136 &{} 0\\ 0.096 &{} 0 &{} 0.19\\ 0 &{} 0.162 &{} 0.222 \end{bmatrix} \end{aligned}$$

The displaced RIs are calculated as:

$$\begin{aligned} RI(p_{ij})= \begin{bmatrix} 0.855 &{} 0.83 &{} 0\\ 0.774 &{} 0 &{} 1.056\\ 0 &{} 0.99 &{} 1.23 \end{bmatrix} \end{aligned}$$
(5.13)

The PPFS values of the PF Importance weights are established by the expert group by evaluating the process: For \(\omega _{j}=({\underline{\omega }}_{j}, {\overline{\omega }}_{j})\), \(\omega _{1}=(0.8256, 0.1935), \omega _{2}=(0.7756,0.2559), \omega _{3}=(0.7226,0.2234)\).

For compromise ranking of the alternatives, the values \(\widetilde{{\mathbb {S}}}(c_{i})\), \(\widetilde{{\mathbb {R}}}(c_{i})\) and \(\widetilde{{\mathbb {Q}}}(c_{i})\) are as in Table 4:

Table 4 Ranking with respect to the cognitive development for Algorithm 1

Now, \(\widetilde{{\mathbb {Q}}}^(c_{2})-\widetilde{{\mathbb {Q}}}^(c_{1})=0.49-0.222=0.268< \frac{1}{3-1}=0.5\). The first condition in S 8 of Algorithm 1 is not satisfied. Then, the ultimate compromise solution is proposed. Therefore the order is \(\{c_{1}, c_{2}\}<c_{3}\).


Algorithm 1 for (FNg):

Consider \(K=K_{i}=\{n_{1}, n_{2}, n_{3}\}\) and \(K_{ii}=\emptyset\). The values \(p(+)_{j}\) and \(p(-)_{j}\):

$$\begin{aligned} p(+)_{j}= & {} \left\{ (0.73, 0.42), (0.95, 0.09), (0.77, 0.37)\right\} \\ p(-)_{j}= & {} \left\{ (0.61, 0.64), (0.62, 0.51), (0.51, 0.54)\right\} \end{aligned}$$

The distance measures are computed as \(\Omega _{H}(p(-)_{1}, p(+)_{1})=0.126\), \(\Omega _{H}(p(-)_{2}, p(+)_{2})=0.3162\), \(\Omega _{H}(p(-)_{3}, p(+)_{3})=0.183\).

$$\begin{aligned} \Omega _{H}(p_{ij}, p(+)_{j})= \begin{bmatrix} 0.1276 &{} 0.322 &{} 0\\ 0.0364 &{} 0 &{} 0.06\\ 0 &{} 0.3832 &{} 0.16 \end{bmatrix} \end{aligned}$$

The displaced RIs are calculated as:

$$\begin{aligned} RI(p_{ij})= \begin{bmatrix} 1.013 &{} 1.018 &{} 0\\ 0.289 &{} 0 &{} 0.33\\ 0 &{} 1.212 &{} 0.874 \end{bmatrix} \end{aligned}$$
(5.14)

The PPFS values of the PF Importance weights: \(\omega _{1}=(0.7139, 0.2821), \omega _{2}=(0.6983,0.3467), \omega _{3}=(0.6762,0.3345)\).

For compromise ranking of the alternatives, the values \(\widetilde{{\mathbb {S}}}(c_{i})\), \(\widetilde{{\mathbb {R}}}(c_{i})\) and \(\widetilde{{\mathbb {Q}}}(c_{i})\) are as in Table 5:

Table 5 Ranking with respect to the behavioral development for Algorithm 1

\(\widetilde{{\mathbb {Q}}}^(d_{1})-\widetilde{{\mathbb {Q}}}^(d_{2})=0.313-0.000=0.313< \frac{1}{3-1}=0.5\). The first condition in S 8 of Algorithm 1 is not satisfied. Then, the ultimate compromise solution is proposed. Therefore the order is \(\{d_{1}, d_{2}\}<d_{3}\).


Algorithm 2 for (FMf):

The values \(\mathbb {SF}(c_{i})\), \(\mathbb {RF}(c_{i})\) and \(\mathbb {QF}(c_{i})\) as in Table 6:

Table 6 Ranking with respect to the cognitive development for Algorithm 2

\(\mathbb {QF}(c_{2})-\mathbb {QF}(c_{1})=0.58 - 0 = 0.58 > \frac{1}{3-1}=0.5\). Then, both of the conditions in S 8 of Algorithm 2 are satisfied. Thus, the order is \(c_{1}<c_{2}<c_{3}\).


Algorithm 2 for (FNg) :

The values \(\mathbb {SF}(d_{i})\), \(\mathbb {RF}(d_{i})\) and \(\mathbb {QF}(d_{i})\) as in Table 7:

Table 7 Ranking with respect to the behavioral development for Algorithm 2

\(\mathbb {QF}(d_{2})-\mathbb {QF}(d_{1})=0.87 - 0 = 0.87 > \frac{1}{3-1}=0.5\). Then, both of the conditions in S 8 of Algorithm 2 are satisfied. Thus, the order is \(d_{2}<d_{1}<d_{3}\).


Algorithm 3 for (FMf):

The values \(p(+)_{j}\) and \(p(-)_{j}\) for (FMf) as follows:

$$\begin{aligned} p(+)_{j}= & {} \left\{ (0.84, 0.34), (0.89, 0.19), (0.87, 0.27)\right\} \\ p(-)_{j}= & {} \left\{ (0.69, 0.54), (0.72, 0.41), (0.56, 0.59)\right\} . \end{aligned}$$

The distance measures are computed for \(\Omega _{H}(p(-)_{j}, p(+)_{j})\) and \(\Omega _{H}(p_{ij}, p(+)_{j})\): \(\Omega _{H}(p(-)_{1}, p(+)_{1})=0.124\), \(\Omega _{H}(p(-)_{2}, p(+)_{2})=0.164\), \(\Omega _{H}(p(-)_{3}, p(+)_{3})=0.18\) and \(\Omega _{H}(p_{11}, p(+)_{1})=0.106\), \(\Omega _{H}(p_{12}, p(+)_{2})=0.136\), \(\Omega _{H}(p_{13}, p(+)_{3})=0\), \(\Omega _{H}(p_{21}, p(+)_{1})=0.096\), \(\Omega _{H}(p_{22}, p(+)_{2})=0\), \(\Omega _{H}(p_{23}, p(+)_{3})=0.19\), \(\Omega _{H}(p_{31}, p(+)_{1})=0\), \(\Omega _{H}(p_{32}, p(+)_{2})=0.162\), \(\Omega _{H}(p_{33}, p(+)_{3})=0.222\). Further, consider the RI values in (5.13).

Using the expectation score function, calculate the weights: \(\mathscr {ES}_{1}=0.6024\), \(\mathscr {ES}_{2}=0.5828\), \(\mathscr {ES}_{3}=0.7898\) and \(\omega _{1}=0.305012658\), \(\omega _{2}=0.295088608\), \(\omega _{3}=0.399898734\) (\(\sum _{j=1}^{3}\omega _{j}=1\)). The RIs values are multiplied by precise weights \(\omega _{j}\):

$$\begin{aligned} RI(p_{ij}).\omega _{j}= \begin{bmatrix} 0.261 &{} 0.245 &{} 0\\ 0.236 &{} 0 &{} 0.422\\ 0 &{} 0.292 &{} 0.492 \end{bmatrix} \end{aligned}$$
Table 8 Ranking with respect to the cognitive development for Algorithm 3

\(\widetilde{{\mathbb {Q}}}(d_{2})-\widetilde{{\mathbb {Q}}}(d_{1})=0.6234 - 0 = 0.6234 > \frac{1}{3-1}=0.5\). Then, both of the conditions in S 8 of Algorithm 3 are satisfied. Thus, the order is \(c_{1}<c_{2}<c_{3}\).


Algorithm 3 for (FNg):

$$\begin{aligned} p(+)_{j}= & {} \left\{ (0.73, 0.42), (0.95, 0.09), (0.77, 0.37)\right\} \\ p(-)_{j}= & {} \left\{ (0.61, 0.64), (0.62, 0.51), (0.51, 0.54)\right\} \end{aligned}$$

The distance measures are computed as \(\Omega _{H}(p(-)_{1}, p(+)_{1})=0.126\), \(\Omega _{H}(p(-)_{2}, p(+)_{2})=0.3162\), \(\Omega _{H}(p(-)_{3}, p(+)_{3})=0.183\) and \(\Omega _{H}(p_{11}, p(+)_{1})=0.1276\), \(\Omega _{H}(p_{12}, p(+)_{2})=0.322\), \(\Omega _{H}(p_{13}, p(+)_{3})=0\), \(\Omega _{H}(p_{21}, p(+)_{1})=0.0364\), \(\Omega _{H}(p_{22}, p(+)_{2})=0\), \(\Omega _{H}(p_{23}, p(+)_{3})=0.06\), \(\Omega _{H}(p_{31}, p(+)_{1})=0\), \(\Omega _{H}(p_{32}, p(+)_{2})=0.3832\), \(\Omega _{H}(p_{33}, p(+)_{3})=0.16\). Further, consider the RI values in (5.14).

Using the expectation score function, calculate the weights: \(\mathscr {ES}_{1}=0.5812\), \(\mathscr {ES}_{2}=0.8186\), \(\mathscr {ES}_{3}=0.692\) and \(\omega _{1}=0.27784683\), \(\omega _{2}=0.391337604\), \(\omega _{3}=0.330815566\) (\(\sum _{j=1}^{3}\omega _{j}=1\)). The RIs values are multiplying by precise weights \(\omega _{j}\):

$$\begin{aligned} RI(p_{ij}).\omega _{j}= \begin{bmatrix} 0.2815 &{} 0.4 &{} 0\\ 0.803 &{} 0 &{} 0.11\\ 0 &{} 0.474 &{} 0.29 \end{bmatrix} \end{aligned}$$
Table 9 Ranking with respect to the behavioral development for Algorithm 3

\(\widetilde{{\mathbb {Q}}}(d_{2})-\widetilde{{\mathbb {Q}}}(d_{1})=0.28 - 0 = 0.28 < \frac{1}{3-1}=0.5\). Then, both of the conditions in S 8 of Algorithm 3 are not satisfied. Thus, the order is \(\{d_{1}, d_{2}\}<d_{3}\).


Algorithm 4 for (FMf):

The values \(\widetilde{\mathbb {SF}}(c_{i})\), \(\widetilde{\mathbb {RF}}(c_{i})\) and \(\widetilde{\mathbb {QF}}(c_{i})\) are as in Table 10:

Table 10 Ranking with respect to the behavioral development for Algorithm 4

\(\widetilde{\mathbb {QF}}(c_{2})-\widetilde{\mathbb {QF}}(c_{1})=1 - 0.17 = 0.83 > \frac{1}{3-1}=0.5\). Then, both of the conditions in S 8 of Algorithm 4 are satisfied. Thus, the order is \(c_{1}<c_{3}<c_{2}\).


Algorithm 4 for (FNg) :

The values \(\widetilde{\mathbb {SF}}(d_{i})\), \(\widetilde{\mathbb {RF}}(d_{i})\) and \(\widetilde{\mathbb {QF}}(d_{i})\) as in Table 11:

Table 11 Ranking with respect to the behavioral development for Algorithm 4

\(\widetilde{\mathbb {QF}}(d_{2})-\widetilde{\mathbb {QF}}(d_{1})=0.242 - 0 = 0.242 < \frac{1}{3-1}=0.5\). Then, both of the conditions in S 8 of Algorithm 4 are not satisfied. Thus, the order is \(\{d_{1}, d_{2}\}<d_{3}\).

Table 12 Results of Solutions for \(\lambda =1\) and \(\delta =0.5\)

It is obvious that the remoteness index-based GPFSS-VIKOR algorithms are simple and easy to perform. Besides facilitating operation and computational efficiency, it has the ability to generate admissibly and reasonable outcomes for the final preference rankings. Moreover, the algorithmic process corresponding to the suggested methods can ensure flexibility to adapt to the particularities in real-world decision cases, such as different preference information of criterion significance and distinct points of reference from displaced and fixed ideal views. As a whole, the practicability and applicability of the developed methods in the GPFSS have been proved via real-world practice to the cognitive and behavioral development of early childhood children. Furthermore, as opposed to the available VIKOR approaches the suggested methods can sufficiently describe more uncertainty and ambiguity inherent in the subjective DM process by employing the GPFSS information.

5.3 Stability analysis

In the application for algorithms, \(\lambda = 1\) and \(\delta =0.5\) were chosen. Solution results for \(\lambda = 1\) and \(\delta =0.5\) are shown in Table 12. When the value of parameter \(\delta\) decreases from 1 to 0, the individual regret gets importance. The ranking of alternatives varying when the values of \(\delta\) change, but the alternative y4 remains the best option from alternatives. Also, different values of the parameter p are used for Algorithms 1-4. Algorithms 1 and 3 used the displaced terminologies and concepts, while Algorithms 2 and 4 are based on the fixed terminologies and concepts. Algorithms 1 and 3 are more sensitive as compared to Algorithms 2 and 4, i.e., the ranking of the alternatives changes rapidly by changing the values of the parameter \(\delta\). For more details, we refer to Tables 13 and 14.

Table 13 Stability analysis for the GPFSS (FMf)
Table 14 Stability analysis for the GPFSS (FNg)

5.4 Comparison with some methods

If there are conflicting criteria among the alternatives, the sorting and selection procedures for these criteria can be carried out. These are the main features of the VIKOR method. In this method, majority and opponent describe maximum group utility and minimum individual regret, respectively. Providing these situations with the VIKOR approach helps decision-makers to obtain a compromise solution.

While the VIKOR method expresses the aggregation operator representing the closest to the ideal, TOPSIS determines a solution with the shortest distance to the ideal solution and the furthest distance to the negative ideal solution. Both methods consider that there exists a performance matrix obtained by the evaluation of all the alternatives in terms of each criterion. Normalization is used to eliminate the units of criterion values. The VIKOR method was created with the idea of nearness to the ideal and emerged from an aggregation operator using linear normalization. TOPSIS method uses vector normalization as normalization. This normalization does not take into account the relative importance of the distances from these points while using two reference points. The main difference appears in the aggregation approaches. An aggregating function is formulated, and it is used as a ranking index. In addition to ranking, the VIKOR method proposes a compromise solution with an advantage rate.

Ranking by PROMETHEE [9, 10], with a linear preference function, gives the same outcomes as ranking by VIKOR, with measures symbolizing group utility. The different weights are less important in the PROMETHEE method than in the VIKOR method. Rankings obtained using TOPSIS and PROMETHEE usual methods have a strong relationship between themselves. Although the PROMETHEE method does not use normalization, it gives quite different results depending on the parameter set used. This indicates that method selection and configuration are important for the DM problem.

ELECTRE is an outranking method that is simple but provides partial ranking [8, 15]. The ELECTRE method is the concept of elimination from a normalized result through a function that is further reduced to concordance and discordance to rank. The ELECTRE method has a basic understanding of the same data as a previously normalized VIKOR with rules. The largest value is the best value, or the smallest value is the best value. The end result can be searched by mathematical deduction of concordance and discordance resulting in a ranking of alternatives that can be decided.

VIKOR technique provides an effective DM method in domains in which the selection of the best alternative is a highly complex environment. Therefore, the VIKOR method may give better results than the DM methods discussed here.

6 Discussion and conclusion

The concept of GPFFS has defined by Kirisci [27]. This new definition adds the generalization parameter to the pool of PFNs and extends PFSS to GPFSS. The concept of GPFSS, which emerged from PFSS, functions to evaluate the information obtained with a parameter that reflects the views of an expert. It is important that the uncertainty associated with an observer be verified by an expert in the choice of degree of membership, and this is the main motivation for developing the GPFSS concept. The Pythagorean fuzzy soft sets are generalized by introducing the generalization parameter given by an expert to validate the original Pythagorean fuzzy values. The role of the expert is to refine PFNs with his domain-specific knowledge. It may be noted that information of any sort often gets misinterpreted during its presentation. This usually happens as the presenter has no domain-specific knowledge or lacks the standard terminologies, and this is sought to be corrected by the proposed GPFSS. An expert can be viewed as a domain professional having years of experience. The expert’s opinion indicates the credibility of the evaluation of the alternative by the experts. This opinion coming from the expert constitutes the generalization parameter. The severity of symptoms as reported by the patients in the linguistic form as membership grades can be used in the PFSS.

In this study, new Euclidean, Hamming, and generalized distance measures were given for GPFSSs according to the new PFNs \(P=(r_{P},d_{P})\), and their basic properties were examined. Pythagorean membership grade was defined by Yager [45] as a new pair of PFNs \((r_{P},d_{P})\). Here, \(r_{P}\) is stated as the strength of commitment and \(d_{P}\) as the direction of commitment. \(r_{p}\) and \(d_{p}\) are characterized by \(u_{p}=r_{p}\cos (\theta _{p})\), \(v_{p}=r_{p}\sin (\theta _{p})\). These values are related to MS and NMS grades. The new distance measures have been used to build useful concepts of remoteness indices, consisting of displaced positive- and negative-ideal remoteness indices as well as fixed positive- and negative-ideal remoteness indices. After distance measurements, new remoteness index-based GPFSS-VIKOR methods have been developed to address MCDA problems in the GPFSS environment, using the concepts of displaced and fixed remoteness indices based on new distance measures. These VIKOR methods are a generalization of the standard VIKOR. This generalization is a new compromising decision-making method based on the proposed generalized distance measure and the displaced and fixed remoteness index. In the DM process, appropriate reference points need to be determined to facilitate the definition of displaced and fixed remoteness indices because these points can effectively form anchored judgments in subjective decision-making processes. Finding these points is the first step for remoteness index-based GPFSS-VIKOR. Then, the displaced and fixed Remoteness Indexes were built. The displaced RI and fixed RI were defined, proved the properties, and given the relevant examples. It has been seen from the examples that the contrasts between the remoteness indices become less noticeable when anchoring subjective judgments with fixed ideal PF values. However, for the displaced ideal PF values, this situation is not very clear.

In decision-making processes, the criteria that determine the decisions are not of equal importance. Therefore, the criteria are weighted when making a decision. These weights are used to determine alternatives and make choices according to the values of subjective evaluations. Two types of weights were used in this study. The first of these weights is calculated with the expectation score function. The second is the weights obtained as a result of recording the evaluations of the experts as PFSS.

In this study, a multi-criteria ranking was created with remoteness indices based on displaced and fixed ideal PF values. Thus, remoteness-based VIKOR methods were built. The smaller the displaced remoteness index is, the more prior the evaluative rating is. In the new VIKOR methods, the priority relationship between the alternatives is evaluated on the basis of criteria. Thus, this priority relationship is defined in ascending order according to each displaced remoteness index because the placed ideal PF values have properties that facilitate anchored judgments. Since a compromise ranking is expected on all criteria, new definitions of a group utility index and an individual regret index are given to synthesize criterion-based priority relationships. As it can be understood from these explanations, in this study, the basic structure of VIKOR has been expanded and a new compromising decision-making method has been established with the proposed distance measures, the displaced remoteness index, and the fixed remoteness index. The method proposed by the study can be easily used to solve MCDA problems where data describing the performance of alternatives are characterized by PF values.

In this study, a new way is proposed to evaluate an aggregated result for each alternative and to obtain a compromise solution. The remoteness indexes obtained based on the new distance measures have the function of formulating the remoteness based on multiple criteria ranking. The computational processes of the four algorithms presented here suggest various preference structures and are capable of fully exploiting the complex information contained in the PF decision matrix. Considering the VIKOR methods in the PF environment, the remoteness-based VIKOR methods proposed in this study are extremely flexible since they detect the location of the displaced and fixed reference points. These new VIKOR methods can be easily applied to real-life problems and solve many challenging decision-making problems. As seen in the real-life example given here, new remoteness-based VIKOR methods can be easily used in very difficult processes such as monitoring and evaluating cognitive and behavioral development. The target audience of the sample discussed is children aged 5-6 years. The effect of this age group’s mandatory stay at home in COVID-19 quarantines on their development has been examined. The factors affecting the developmental processes of children aged 5-6 years were evaluated with the newly given remoteness-based VIKOR methods, and the factors affecting development were listed.