Introduction

With the rapid development of emerging technologies and continuous upgrade of municipal and rural tourism consumption and services, sharing accommodation has gradually gone into people’s sight. Figure 1 indicates the high growth of sharing accommodation orders in some Chinese cities. Sharing accommodation has characteristics of diversified supply subjects and service contents and socialized user experience, which can meet diverse accommodation need. At the same time, sharing accommodation can reduce the information asymmetry and transaction risk between the landlords and the tenants. It offers a better service experience to landlords [1]. Not only that the economic development of sharing accommodation plays a positive role in expanding employment and entrepreneurship channels, increasing the income of urban and rural residents, and helping rural revitalization [2]. In 2017, Airbnb created over 20,000 employment opportunities in Beijing and Shanghai. According to a survey of Airbnb’s domestic landlords, landlords can significantly improve their daily living standards by house sharing income. More than that, Airbnb participates in the innovative tourism targeted poverty alleviation project, explores new models for rural housing poverty relief, and helps rural poor families get rid of poverty. Sharing accommodation has made great contributions to social development.

Fig. 1
figure 1

The growth of sharing accommodation orders in some Chinese cities

The report on development of China’s sharing accommodation shows that the market scale of sharing accommodation industry will maintain a growth rate of 50% in the next 3 years, as well as the number of online housing resource will get an increase [1]. In such a trend, it is more difficult to obtain the required housing resource information by users themselves. This situation is bad for the development of sharing accommodation. Therefore, it is worth thinking about how to response to this trend and provide accurate accommodation information for users.

Recommender systems are used as decision support systems to overcome information overload in online environments [3]. Collaborative filtering (CF) [4] recommender system is one of the most popular types of recommender systems [5]. It is mainly used as information retrieval system based on past users’ activities and ratings on the items. Content-based recommendation is the continuation and development of CF recommendation. It does not rely on users’ historical activities and ratings on the items, but calculates the similarity degree between users based on product content information [6]. Besides, according to online review information, scholars have put forward some recommendation systems based on users’ historical information [7, 8]. By extracting the online review information, the text information needs to be transformed into computable data, and it will inevitably lead to the loss of information. Therefore, it is necessary to propose a new recommendation method to solve the existed problems.

Without considering the content characteristics of users and products, the network-based recommendation algorithm directly regards them as abstract nodes. Owing to the relatively clear data structure, bipartite networks have attracted widespread attention. They have been successively applied to cluster analysis [9,10,11], pattern recognition [12,13,14], recommender systems [15,16,17,18], and other fields. Zhou et al. presented a bipartite network algorithm based on resource-allocation dynamics [19]. The similarity degree among homogeneous nodes is utilized to implement personal recommendation for users. Subsequently, Liu et al. analyzed the similarities and differences between weighted and unweighted bipartite graph networks by combining collaborative filtering with network reasoning [20]. It is concluded that recommendation accuracy of the weighted bipartite graph network is higher than that of the unweighted one under the same amount of calculations. The attendant problem is how to determine the weights of adjacent edge for improving algorithm accuracy. Wang et al. introduced a half cumulative distribution method to normalize the weights of adjacent edges, and mapped the original scores to the probability of users’ preference [21]. Pan et al. defined the weights of edges based on the distribution of degree between user and object [22]. Then, they applied the weights of edges in the quality diffusion to improve the accuracy of recommendation algorithm. Li et al. used monotonic saturation function to calculate the weights of edges [23]. An approach to calculate the weights which is based on balance factor was proposed by Song et al. [24]. In the existed methods, the connection degree and the distribution of two kinds of nodes are two commonly used tools to calculate the weights of adjacent edges, which are based on users’ ratings or historical usage. As the amount of information increases, users are accustomed to utilize key words or attributes to filter information and services. Besides, when we determine the recommended content based on the users’ historical data, it will not only focus on the direct evaluation on object, but also analyze the root cause of users’ choices, because users pay attention to the characteristics of the content, which is one of the research contents in this paper.

Users often passively accept recommendations rather than making their own decisions; that is to say, the information service platform recommends information or services to users depending on their historical data. In some cases, decision-makers are not always able to obtain “certain information” when evaluating alternatives because of the diversification of users’ demands. How to depict the evaluation information more reasonably is another urgent problem. To describe the uncertainty, Zadeh proposed the concept of fuzzy sets [25, 26]. With the continuous process of human practice and increasingly complicated decision-making environment, some extended forms of fuzzy sets have been proposed to model uncertainty, such as the interval-valued fuzzy sets [27], the intuitionistic fuzzy sets [28,29,30], the type-2 fuzzy sets [31, 32], the fuzzy multi-sets [33], the Pythagorean fuzzy sets [34, 35], the hesitant fuzzy sets (HFSs) [36, 37], the dual hesitant fuzzy soft sets [38], the probabilistic hesitant fuzzy sets [39], and so on. Since the hesitant fuzzy set was proposed, it has been widely applied in many fields, such as decision-making [40,41,42,43], emergency management [44, 45], online reviews [46], supply chain management [47], and so on. When defining the membership of an element, the difficulty of establishing the membership degree is not because we have a margin of error [as in intuitionistic fuzzy sets (IFSs)], but because we have a set of possible values [36]. For example, a decision organization that contains ten experts is authorized to assess the degree that an alternative should satisfy an attribute. Some of them think that 0.9 is most appropriate and the rest insist that it is 0.5, and these two different groups of experts cannot persuade each other. It is noted that the HFE {0.5, 0.9} can describe the above situation more objectively than the fuzzy number 0.5 (or 0.9) or the interval-value fuzzy set [0.5, 0.9]. Therefore, HFSs can express the decision-makers’ judgments more informative and reliable in this case. Probabilistic hesitant fuzzy set (PHFS), as an extension of HFS, not only contains different membership degrees, but also gives the occurrence probabilities associated with membership degrees, which reserves more original information than HFSs. Therefore, it can fully describe the uncertainty of decision-making information, which has attracted more and more scholars’ attention. Zhou and Xu defined the basic operation of PHFSs [48]. Furthermore, the operation rules and aggregation operators of PHFSs are defined by Zhang et al. [49]. Subsequently, it was applied to various fields such as decision-making [50, 51], supply chain management [52], project investment [53], and so on. Due to the uncertainty of users’ recommendation content, PHFS which can be thought as an effective tool for modeling the uncertainty is utilized to present the evaluation information in this paper. Considering those advantages of PHFS, we introduce it into the bipartite network projection, and establish a novel bipartite graph network which integrates users’ decision information and experts’ evaluation information simultaneously.

In addition, the satisfaction degree of PHFE is proposed, and the weights of adjacent edges are determined in this paper. We combine bipartite network with PHFS to address complex decision-making problems in recommendation system. The main contributions and innovation of this paper are mainly reflected in the following aspects:

  1. (1)

    Recommendation system is one of the most effective ways to overcome information overload. A novel recommender system is proposed to aid decision-making in this paper. It could assist us in making decisions with huge information.

  2. (2)

    Based on bipartite network, a recommendation decision-making method is proposed. To measure the degree of diversification for users’ demands, probabilistic hesitant fuzzy set is used to model the uncertainty, and thus, a probabilistic hesitant fuzzy bipartite network is established in this paper to aid decision-making.

  3. (3)

    The weights of adjacent edges play a vital role in bipartite network projection. Then, a probabilistic hesitant fuzzy satisfaction degree is presented to calculate the weight of adjacent edges. Our method expresses the uncertainty flexibly, and, at the same time, greatly improves the accuracy of recommendation system.

  4. (4)

    The proposed method can address a series of problems such as information overloading, ambiguity of user’s requirement, and so on, and, thus, provide effective measures for the development of sharing economy.

The remaining sections of this paper are organized as follows. Some basic concepts related to PHFS and bipartite network projection are introduced in Sect. 1. In Sect. 2, we establish a bipartite graph connecting users and alternatives. Then, the satisfaction degree of probabilistic hesitant fuzzy element is defined. In Sect. 3, an illustrative example of recommending a hotel for Airbnb user is used to explain the feasibility and effectiveness of the proposed method. And conclusions are drawn in Sect. 4.

Preliminaries

Some basic concepts related to PHFS and the preliminaries used throughout this paper are introduced in this section.

PHFSs

To model preferences of decision-makers, Zhu defined PHFSs as follows [39].

Definition 1

Let \( X = \left\{ {x_{1} ,x_{2} , \ldots ,x_{n} } \right\} \) be a fixed set, a PHFS on \( X \) is defined in terms of a function \( h_{x} \left( {p_{x} } \right) \) that when applied to \( X \) returns a subset of [0, 1], which is expressed as:

$$ H = \left\{ {\left\langle {x,h_{x} \left( {p_{x} } \right)} \right\rangle \left| {x \in X} \right.} \right\}, $$
(1)

where the function \( h_{x} \left( {p_{x} } \right) \) is a set of some values in [0, 1]. \( h_{x} \) denotes the possible membership degrees of the element \( x \) in \( X \) to the set \( H \), and \( p_{x} \) is a set of probabilities associated with \( h_{x} \). \( h_{x} \left( {p_{x} } \right) \) is called probabilistic hesitant fuzzy element (PHFE). For convenience, we express \( h_{x} \left( {p_{x} } \right) \) as \( h\left( p \right) \):

$$ h\left( p \right) = h_{x} \left( {p_{x} } \right) = \left\{ {\gamma_{l} \left( {p_{l} } \right),l = 1,2, \ldots ,\left| {h\left( p \right)} \right|} \right\}, $$
(2)

where \( p_{l} \) is the probability of the membership degree of \( \gamma_{l} \), \( \gamma_{l} \left( {p_{l} } \right) \) is called PHFE, \( \left| {h\left( p \right)} \right| \) is the number of all different membership degrees, and \( \sum\nolimits_{l = 1}^{{\left| {h\left( p \right)} \right|}} {p_{l} } = 1 \).

To compare PHFEs, the score function and the deviation function of PHFEs are defined as follows.

Definition 2

Let \( h\left( p \right) \) be a PHFE. Then:

$$ E\left( {h\left( p \right)} \right) = \sum\limits_{l = 1}^{{\left| {h\left( p \right)} \right|}} {\gamma_{l} \times p_{l} } , $$
(3)

is called the score function of the \( h\left( p \right) \).

Definition 3

Let \( h\left( p \right) \) be a PHFE. Then,

$$ D\left( {h\left( p \right)} \right) = \sum\limits_{l = 1}^{{\left| {h\left( p \right)} \right|}} {\left( {p_{l} \gamma_{l} - E\left( {h\left( p \right)} \right)} \right)^{2} } , $$
(4)

is called the deviation function of \( h\left( p \right) \).

Bipartite network projection and personality recommendation

In the bipartite graph network, there are two types of node sets: \( X \) nodes and \( Y \) nodes. Arbitrary two nodes in the same set are not directly connected. If there is a connection between two nodes in different sets, there must be an association or relationship between the two nodes. For different models, the nature, the content, and the sense of relationships between nodes are different. The bipartite graph network contains a wealth of information, which has great significance both in theoretical research and practical applications. It has been used to investigate the relationship between nodes in the same set by projection compression.

Bipartite network projection

Based on resource allocation, Zhou et al. [19] proposed a bipartite projection network, and analyzed the inherent relationship between nodes in the same set by correlating the nodes in two different sets.

For a given bipartite graph network, as shown in Fig. 2a, the internal relationship between nodes can be studied by two-step projection compression: first, the resources flow from \( X \) nodes to \( Y \) nodes. The resources obtained by \( Y \) nodes are shown in Fig. 2b. Then, resources flow back to \( X \) nodes. The resources on \( X \) nodes are shown in Fig. 2c. In this case, the resources of nodes are equally allocated to their adjacent edges.

Fig. 2
figure 2

Bipartite projection network

The final resources located in those three nodes are denoted by \( x^{\prime},y^{\prime},\,{\text{and}}\,z^{\prime} \), which can be expressed by:

$$ \left( \begin{aligned} {x^{\prime}} \hfill \\ {y^{\prime}} \hfill \\ {z^{\prime}} \hfill \\ \end{aligned} \right) = \left( {\begin{array}{*{20}c} {11/18} & {1/6} & {5/18} \\ {1/9} & {5/12} & {5/18} \\ {5/18} & {5/12} & {4/9} \\ \end{array} } \right)\left( {\begin{array}{*{20}c} x \\ y \\ z \\ \end{array} } \right). $$
(5)

Here, the element in the \( i{\text{th}} \) row and \( j{\text{th}} \) column represents the fraction of resource that the \( j{\text{th}} \) node of X transfers to the \( i{\text{th}} \) node. The larger the value is, the closer the relationship between the two nodes is.

However, in real life, the attribute weights are not always the same. The definition of bipartite projection with general weight is given as below [19].

Suppose \( G\left( {X,Y,\varGamma } \right) \) is a general bipartite network, where \( \varGamma \) is the set of edges. \( X \) is the set of upper nodes, \( Y \) is the set of lower nodes, and the nodes in \( X \) and \( Y \) are denoted by \( x_{1} ,x_{2} , \ldots ,x_{n} \) and \( y_{1} ,y_{2} , \ldots ,y_{m} \), respectively. We first assume that the connection of nodes only exists between the different sets, and it does not exist in the same set. The initial resource located on \( x_{i} \) node meets \( f\left( {x_{i} } \right) \ge 0 \).

All the resources in \( X \) flow to \( Y \), and the resource of the \( y_{l} \) node is denoted as:

$$ f\left( {y_{l} } \right) = \sum\limits_{i = 1}^{n} {\frac{{a_{il} f\left( {x_{i} } \right)}}{{k\left( {x_{i} } \right)}}} , $$
(6)

where \( k\left( {x_{i} } \right) \) denotes the degree of \( x_{i} \) and \( a_{il} \) is a characteristic function of set \( \varGamma = \left\{ {o\left( {i,l} \right)\left| {x_{i} \in X,\quad y_{l} \in Y} \right.} \right\} \):

$$ a_{il} = \left\{ {\begin{array}{*{20}l} {1,} \hfill & {o\left( {i,l} \right) \in \varGamma } \hfill \\ {0,} \hfill & {\text{otherwise}} \hfill \\ \end{array} } \right.. $$
(7)

Then, all the resources flow back to \( X \), and the final resource of node \( x_{i} \) is:

$$ f^{\prime}\left( {x_{i} } \right) = \sum\limits_{l = 1}^{m} {a_{il} } f\left( {y_{l} } \right)/k\left( {y_{l} } \right) = \sum\limits_{l = 1}^{m} {\frac{{a_{il} }}{{k\left( {y_{l} } \right)}}} \sum\limits_{j = 1}^{n} {\frac{{a_{jl} f\left( {x_{j} } \right)}}{{k\left( {x_{j} } \right)}}} . $$
(8)

This can be expressed as:

$$ f^{\prime}\left( {x_{i} } \right) = \sum\limits_{j = 1}^{n} {w_{ij} f\left( {x_{j} } \right)} , $$
(9)

where

$$ w_{ij} = \frac{1}{{k\left( {x_{j} } \right)}}\sum\limits_{l = 1}^{m} {\frac{{a_{il} a_{jl} }}{{k\left( {y_{l} } \right)}}} . $$
(10)

Here, \( k\left( {x_{i} } \right) \) denotes the degree of \( x_{i} \) and \( k\left( {y_{l} } \right) \) denotes the degree of \( y_{l} \).

Personal recommendation based on weighted bipartite network projection

Wang et al. [21] proposed a weighted bipartite network projection for personalized recommendations.

In a bipartite network \( G\left( {U,A,\varGamma } \right) \), \( \varGamma \) is the set of edges, \( U \) is the upper nodes set for users, and \( A \) is the lower nodes set for alternatives. The nodes in \( U \) and \( A \) are denoted by \( U_{1} ,U_{2} , \ldots ,U_{n} \) and \( A_{1} ,A_{2} , \ldots ,A_{m} \), respectively. The connection exists only between \( U \) and \( A \), and it is not directly connected in the same set. For user \( U_{i} \), the weighted network-based method starts by assigning the initial resource for alternatives. If the alternative connects with the user, it would be assigned a unit resource as its initial resource, otherwise zero. The resource-allocation process includes two stages:

  1. (1)

    The resource flows from \( A \) to \( U \).

    The resource of alternative \( A_{i} \) is assigned to its neighbor users according to the ratio of the edge weights. The total resources of user \( U_{k} \) are as follows:

    $$ g\left( {U_{k} } \right) = w_{ki} f\left( {A_{i} } \right). $$
    (11)

    Here, \( g\left( {U_{k} } \right) \) is the resource that user \( U_{k} \) would obtain from its neighbor alternatives; \( f\left( {A_{i} } \right) \) is the initial resource for alternative \( A_{i} \).

  2. (2)

    The resource flows from \( {\text{U}} \) to \( {\rm A} \).

    Similarly, the finally resources that alternative \( A_{i} \) obtained from its connected users across the process of allocation are as follows:

    $$ f^{\prime}\left( {A_{i} } \right) = {{w_{ki} g\left( {U_{k} } \right)} \mathord{\left/ {\vphantom {{w_{ki} g\left( {U_{k} } \right)} {\sum\limits_{k = 1}^{s} {w_{ki} } }}} \right. \kern-0pt} {\sum\limits_{k = 1}^{s} {w_{ki} } }}. $$
    (12)

    By plugging Eqs. (11) into (12), \( f^{\prime}\left( {A_{i} } \right) \) can be expressed by \( f^{\prime}\left( {A_{i} } \right) = \rho_{ki} f\left( {A_{i} } \right) \), where \( \rho_{ki} = w_{ki}^{2} /\sum\nolimits_{k = 1}^{s} {w_{ki} } \), and \( \rho_{ki} \) denotes the similarity coefficient of user \( U_{k} \) and user \( U_{i} \). The greater the similarity coefficient between user \( U_{k} \) and user \( U_{i} \) is, the more possibility they choose the same alternatives. Then, the alternative can be recommended for the user \( U_{k} \) based on the choice of user \( U_{i} \).

Recommendation decision-making algorithm based on PHFSs and bipartite network projection

The probabilistic hesitant fuzzy satisfaction degree

For a probabilistic hesitant fuzzy multi-attribute decision-making problem, there is a set of \( n \) alternatives \( A = \left\{ {A_{1} ,A_{2} , \ldots ,A_{n} } \right\} \), and a set of \( m \) attributes \( C = \left\{ {C_{1} ,C_{2} , \ldots ,C_{m} } \right\} \). Experts evaluate the alternatives for each attribute, and the set of all possible evaluations for an alternative \( A_{i} \left( {i = 1,2, \ldots ,n} \right) \) under the attribute \( C_{j} \left( {j = 1,2, \ldots ,m} \right) \) can be considered as a PHFE \( h_{ij} \left( {p_{ij} } \right) \). Thus, based on the above description, we can construct a probabilistic hesitant fuzzy decision matrix \( \varTheta \) as follows:

$$ \varTheta = \left[ {\begin{array}{*{20}c} {h_{11} \left( {p_{11} } \right)} & {h_{12} \left( {p_{12} } \right)} & \cdots & {h_{1m} \left( {p_{1m} } \right)} \\ {h_{21} \left( {p_{21} } \right)} & {h_{22} \left( {p_{22} } \right)} & \cdots & {h_{2m} \left( {p_{2m} } \right)} \\ \vdots & \vdots & \ddots & \vdots \\ {h_{n1} \left( {p_{n1} } \right)} & {h_{n2} \left( {p_{n2} } \right)} & \cdots & {h_{nm} \left( {p_{nm} } \right)} \\ \end{array} } \right]. $$
(13)

To measure the satisfaction degree of alternatives better and understand user’s requirement, we define the satisfaction degree of PHFEs which is inspired by the idea of Liu et al. [54].

Definition 4

Let \( h\left( p \right) = \{ \gamma_{l} \left( {p_{l} } \right),l = 1,2, \ldots ,\left| {h\left( p \right)} \right|\,{\text{and}}\,\sum\nolimits_{l = 1}^{{\left| {h\left( p \right)} \right|}} {p_{l} = 1} \} \) be a PHFE. We define the satisfaction degree of \( h\left( p \right) \) as:

$$ \varphi \left( {h\left( p \right)} \right) = \frac{{E\left( {h\left( p \right)} \right)}}{{1 + D\left( {h\left( p \right)} \right)}} = \frac{{E\left( {h\left( p \right)} \right)}}{{1 + \sum\limits_{l = 1}^{{\left| {h\left( p \right)} \right|}} {\left( {p_{l} \gamma_{l} - E\left( {h\left( p \right)} \right)} \right)^{2} } }}. $$
(14)

Here, \( E\left( {h\left( p \right)} \right) \) is the score function and \( D\left( {h\left( p \right)} \right) \) is the deviation function of \( h\left( p \right) \), which can be calculated by Eqs. (2) and (3), respectively.

The decision-making information can be utilized to measure the satisfaction of alternatives more completely by considering the score function and deviation function simultaneously. The larger the \( E\left( {h\left( p \right)} \right) \) is, the higher the satisfaction is. While the deviation function can reflect the level of disagreement among the decision-makers. Intuitively, the smaller the deviation \( D\left( {h\left( p \right)} \right) \) is, the higher the satisfaction degree is.

The recommendation decision-making method

In the traditional bipartite graph recommendation system, the evaluation value of user is represented as crisp number. However, in the actual decision-making, users are often hesitant to express their preferences for alternatives. To make the recommendation content closer to the users’ preferences, a bipartite network which considers the user’s demand for product attributes and the satisfaction degree of the recommendation content is established. The satisfaction degree is used to determine the adjacent edge weight of the bipartite network. Combining probabilistic hesitant fuzzy information with the bipartite graph network, a recommendation decision-making algorithm is proposed.

Assume that \( A = \left\{ {A_{1} ,A_{2} , \ldots ,A_{n} } \right\} \) is the set of alternative, \( C = \left\{ {C_{1} ,C_{2} , \ldots ,C_{m} } \right\} \) is the set of attribute, \( U = \left\{ {U_{1} ,U_{2} , \ldots ,U_{\text{s}} } \right\} \) is the set of user, and \( \varGamma \) is the set of adjacent edge. According to users’ requirement, the evaluation values for alternatives take the form of PHFEs. Then, the probabilistic hesitant fuzzy decision-making matrix is expressed as follows:

$$ \varTheta = \left[ {\begin{array}{*{20}c} {h_{11} \left( {p_{11} } \right)} & {h_{12} \left( {p_{12} } \right)} & \cdots & {h_{1m} \left( {p_{1m} } \right)} \\ {h_{21} \left( {p_{21} } \right)} & {h_{22} \left( {p_{22} } \right)} & \cdots & {h_{2m} \left( {p_{2m} } \right)} \\ \vdots & \vdots & \ddots & \vdots \\ {h_{n1} \left( {p_{n1} } \right)} & {h_{n2} \left( {p_{n2} } \right)} & \cdots & {h_{nm} \left( {p_{nm} } \right)} \\ \end{array} } \right]. $$
(15)

The satisfaction degree of PHFE can be calculated by Eq. (14), and then, the probabilistic hesitant fuzzy satisfaction degree matrix \( \varPsi \) can be obtained:

$$ \varPsi = \left[ {\begin{array}{*{20}c} {\varphi \left( {h_{11} \left( {p_{11} } \right)} \right)} & {\varphi \left( {h_{12} \left( {p_{12} } \right)} \right)} & \cdots & {\varphi \left( {h_{1m} \left( {p_{1m} } \right)} \right)} \\ {\varphi \left( {h_{21} \left( {p_{21} } \right)} \right)} & {\varphi \left( {h_{22} \left( {p_{22} } \right)} \right)} & \cdots & {\varphi \left( {h_{2m} \left( {p_{2m} } \right)} \right)} \\ \vdots & \vdots & \ddots & \vdots \\ {\varphi \left( {h_{n1} \left( {p_{n1} } \right)} \right)} & {\varphi \left( {h_{n2} \left( {p_{n2} } \right)} \right)} & \cdots & {\varphi \left( {h_{nm} \left( {p_{nm} } \right)} \right)} \\ \end{array} } \right]. $$
(16)

The standard frequency matrix \( \varOmega \) of the user’s requirement for attributes can be calculated as below:

$$ \varOmega = \left[ {\begin{array}{*{20}c} {t_{11} } & {t_{12} } & \cdots & {t_{1m} } \\ {t_{21} } & {t_{22} } & \cdots & {t_{2m} } \\ \vdots & \vdots & \ddots & \vdots \\ {t_{s1} } & {t_{s2} } & \cdots & {t_{sm} } \\ \end{array} } \right], $$
(17)

where \( \sum\nolimits_{j = 1}^{m} {t_{kj} = 1,(k = 1,2, \ldots ,s)} \).

Thus, combining matrix \( \varPsi \) and matrix \( \varOmega \), we could acquire the set of adjacent edge \( o\left( {i,k} \right) \) between user \( U_{k} \left( {k = 1,2 \ldots ,s} \right) \) and alternative \( A_{i} \left( {i = 1,2, \ldots ,n} \right) \):

$$ \varGamma = \left\{ {o\left( {i,k} \right)\left| {{\text{for}}\,{\text{every }}t_{jk} \ne 0,\,{\text{we}}\,{\text{have}}\,\varphi_{ij} \ne } \right.0,j = 1,2, \ldots m} \right\}. $$
(18)

Then, the weight of adjacent edge \( \phi_{ki}^{{}} \) between alternative \( A_{i} \left( {i = 1,2, \ldots ,n} \right) \) and user \( U_{k} \left( {k = 1,2 \ldots ,s} \right) \) can be calculated by the following equation:

$$ \phi_{ki} = \sum\limits_{j = 1}^{m} {t_{kj} \varphi_{ij} } . $$
(19)

In the recommendation decision-making process, the recommended user \( U_{\text{r}} \) is set to gray, while the initial resource of the alternative \( A_{\text{c}} \) connected to the recommended user is set to 1, otherwise 0. The projection steps are shown as follows:

(1) The resource flows from \( A \) to \( U \), and the resource of user \( U_{k} \) is:

$$ g\left( {U_{k} } \right) = w_{ki} f\left( {A_{i} } \right). $$
(20)

Let

$$ g\left( {U_{h} } \right) = { \hbox{max} }g\left( {U_{k} } \right),(k = 1,2, \ldots ,s{\text{ and }}k \ne r), $$
(21)

where \( w_{ki} = {{\phi_{ki}^{{}} } \mathord{\left/ {\vphantom {{\phi_{ki}^{{}} } {\sum\limits_{i = 1}^{n} {\phi_{ki}^{{}} } }}} \right. \kern-0pt} {\sum\nolimits_{i = 1}^{n} {\phi_{ki}^{{}} } }} \), and \( f\left( {A_{i} } \right) \) is the initial resource of alternative \( A_{i} \). The resource of user \( U_{k} \) is obtained from its neighbor edge.

(2) The resource flows from \( U \) to \( A \), and the resources owned by alternative \( A_{i} \) are \( f^{\prime}\left( {A_{i} } \right) \).

$$ f^{\prime}\left( {A_{i} } \right) = \rho_{ki} f\left( {A_{i} } \right). $$
(22)

Let

$$ F\left( {A_{o} } \right) = \hbox{max} \, f^{\prime}\left( {A_{i} } \right),\left( {i = 1,2, \ldots ,n\,{\text{and}}\,i \ne c} \right), $$
(23)

where \( \rho_{ki} = w_{ki}^{2} /\sum\nolimits_{k = 1}^{s} {w_{ki} } \), and \( \rho_{ki} \) denotes the similarity degree between user \( U_{k} \) and user \( U_{i} \). The higher the similarity degree is, the more similar the preferences for alternatives are. Thus, we could recommend alternative \( A_{o} \) for user \( U_{r} \) by analyzing the preference of user \( U_{h} \).

The steps of the proposed recommendation decision-making algorithm

Based on the above analysis, a recommendation decision-making algorithm based on probabilistic hesitant fuzzy bipartite graph network is summarized as below:

Step 1. According to users’ requirements for each attribute \( C_{j} (i = 1, \ldots ,m) \), we can obtain the corresponding frequency matrix \( \varOmega \). Experts evaluate alternatives \( A_{i} \left( {i = 1,2, \ldots ,n} \right) \), and the corresponding probability hesitant fuzzy evaluation matrix \( \varTheta \) can be obtained.

Step 2. Calculate the satisfaction degree of PHFE in \( \varTheta \) by Eq. (14) and the corresponding probabilistic hesitant fuzzy satisfaction degree matrix \( \varPsi \) can be derived.

Step 3. Combining matrix \( \varPsi \) and matrix \( \varOmega \), the set \( \varGamma \) of adjacent edge \( o\left( {i,k} \right) \) between user \( U_{k} \left( {k = 1,2 \ldots ,s} \right) \) and alternative \( A_{i} \left( {i = 1,2, \ldots ,n} \right) \) can be obtained by Eq. (18). The adjacent edge weights \( \phi_{ki} (k = 1, \ldots ,s,i = i, \ldots ,n) \) of bipartite networks are calculated by Eq. (19).

Step 4. The recommended user \( U_{r} \) is set to gray, while the initial resource of the alternative \( A_{\text{c}} \) connected to the recommended user is set to 1, otherwise 0. Calculate \( g\left( {U_{h} } \right) \) and \( F\left( {A_{o} } \right) \) by Eqs. (21) and (23), respectively.

Step 5. Output the recommended alternative \( A_{o} \). Therefore, it is prior to recommend alternative \( A_{o} \) for user \( U_{\text{r}} \) by analyzing the preference of user \( U_{h} \).

Step 6. End.

To facilitate the reader’s understanding, the flowchart of the algorithm is shown in Fig. 3.

Fig. 3
figure 3

The flowchart of the algorithm

Illustrative example

Airbnb is one of the leading rental communities, where users can post or search vacation rental information and complete online booking programs via its web or mobile apps. It is convenient to obtain house information on Airbnb mobile application. Then, an eligible list will be shown for users on the website after inputting destination, dates, guests’ number, etc. It can be further set property prices, community facilities, and other search results for secondary screening. By clicking on a link in the message, user could obtain more detailed information, including description, photos, and reviews of previous guests. The search page of Airbnb mobile software is shown in Fig. 4.

Fig. 4
figure 4

The search page of Airbnb mobile application

In 2015, Airbnb settled in China. Due to the fierce competition in the Chinese market, it did not get the expected growth in time. To get a breakthrough in the Chinese market, Airbnb needs to carry out adequate public relations activities and adopts precise market positioning and marketing strategies [55]. According to the report released by Airbnb in 2019 [56], the disposable income of residents has been continuously increasing in recent years and the consumption structure of residents has been gradually optimized. Moreover, the consumption of education and entertainment continues to grow. In recent years, tourism market presents sound and rapid development momentum, with 4.0 domestic per capita trips and 10.7 outbound trips per 100 people. With comprehensive penetration of the Internet, people’s lifestyle and consumption habits are constantly changing. At the same time, the tourism market is increasingly delivered over the Internet.

The millennials (who were born in 1982–2000) constitute nearly 80% of online travelers, and become the main force of tourism market. They pursue a high-quality and fashion lifestyle, but also explore the unknowns without neglecting to interact with family members. Nearly half of millennials have a higher expense on accommodation than before. In addition, the pursuit of comfortable sleep and the general relaxation environment like home and a colorful accommodation experience will greatly enhance travel happiness. With the development of the sharing economy and the popularity of sharing accommodation, the specialty home stays and the short-term rentals are ever more popular among millennial travelers. However, it is hard to obtain an accurate house information in travel. Hence, we propose a method to recommend accommodation for Airbnb users in this paper.

A handful of studies on tourists’ motivations for choosing Airbnb have been investigated. Quinby and Gasdia found that home amenities and space were the two main motivations for users to choose sharing accommodation [57]. Several studies show that price (or economic benefits) has been recognized sometimes as the most important factor [57,58,59]. Nowak et al. considered location to be the second most important factor [59]. Besides, some scholars analyzed the user’s intention of using Airbnb [60,61,62]. According to the above literatures and attributes collected from Airbnb community (https://www.airbnb.cn), factors influencing user’s decision are shown in Table 1.

Table 1 Influence factor

Five different houses near Nanjing Confucius Temple are to be evaluated (as shown in Table 2). Four users who have stayed at the five houses are randomly selected. According to their habits of booking rooms at Airbnb and taking those houses as an example, we illustrate how to select a recommendation alternative using the proposed method in this paper.

Table 2 House information

Step 1. According to requirements of users, the absolute frequency matrix is shown in Table 3.

Table 3 The absolute frequency matrix

The corresponding requirement relative frequency matrix is shown in Table 4.

Table 4 The relative frequency matrix

The alternative is evaluated by experts based on the required attributes of users (https://share.weiyun.com/5J3KcW8), and the corresponding probabilistic hesitant fuzzy evaluation matrix is obtained which is shown in Table 5.

Table 5 Probabilistic hesitant fuzzy evaluation matrix

Step 2. Calculate the satisfaction degree of PHFE according to Eq. (14), and they are shown in Table 6.

Table 6 The satisfaction degree of PHFE

Step 3. We establish bipartite network between users and alternatives based on PHFEs. Users’ requirements for attributes and the satisfaction degree of PHFEs are considered simultaneously. Specific steps of recommendation method are as follows:

First, determine the set of adjacent edges according to Eq. (18):

$$ \begin{aligned}\varGamma &= \left\{ o\left( {1,1} \right),o\left( {1,4} \right),o\left( {1,5} \right),o\left( {2,1} \right),o\left( {2,2} \right),o\left( {2,3} \right),\right.\\ &\quad \left. o\left( {2,4} \right),o\left( {3,3} \right),o\left( {4,1} \right),o\left( {4,2} \right),o\left( {4,3} \right),o\left( {4,4} \right) \right\}. \end{aligned}$$

Then, calculate the weight of adjacent edge according to Eq. (19) as shown in Table 7.

Table 7 The weight matrix of adjacent edges

Step 4. In the established bipartite graph network, the recommended user \( U_{ 3} \) was filled with gray. If the alternatives have been connected with user \( U_{ 3} \), its initial resource would be assigned a unit resource, otherwise zero. As shown in Fig. 5a, initial resource of alternative \( A_{3} \) is 1, and initial resource of the other alternatives is 0. First, the resource flows from alternative to user, and the resource allocation is shown in Fig. 5b. Then, resource flows back to alternative, and the resource allocation is shown in Fig. 5c.

Fig. 5
figure 5

Bipartite network projection based on satisfaction degree of PHFEs

It is easy to see from Fig. 5 that \( g\left( {U_{h} } \right) = g\left( {U_{4} } \right) \) and \( F\left( {A_{o} } \right) = f^{\prime}(A_{2} ) \).

Step 5. Finally, it is prior to recommend alternative \( A_{2} \) for user \( U_{ 3} \). Because user \( U_{ 3} \) is similar to user \( U_{ 4} \) and user \( U_{ 4} \) gives a high evaluation for alternative \( A_{2} \).

Step 6. End.

Comparative analysis

Comparison analysis with the bipartite network projection recommendation method

Personalized recommendation method based on resource allocation was first introduced by Zhou et al. [19]. They proposed recommending objects to users in accordance with resource allocation. In what follows, we use bipartite projection network personal recommendation method to recommend a suitable accommodation for user \( U_{ 3} \). The bipartite network between users and alternatives is established and shown in Fig. 6.

Fig. 6
figure 6

Bipartite network projection personal recommendation

As shown in Fig. 6, the final resources obtained by each program are:

$$ A_{1} = 4/27,A_{2} = 2/9,A_{3} = 13/27,A_{4} = 4/27,A_{5} = 0. $$

Obviously, alternative \( A_{2} \) obtained maximum resources, so alternative \( A_{2} \) is the prior alternative which is recommended to users \( U_{ 3} \). However, it is worth noting that if the recommended content is determined according to the final resource allocation, and the resources obtained by the alternative \( A_{1} \) and alternative \( A_{4} \) are equal. In fact, from Table 3, it can be seen that the user \( U_{ 3} \) have a requirement of attributes \( C_{1} \), \( C_{3} \), \( C_{5} \), and \( C_{6} \). Obviously, alternative \( A_{1} \) and alternative \( A_{4} \) have different degree of satisfying those attributes, which implies that the resources obtained by the two alternatives should be different. The main reason is that the calculating method for the edge weight is not reasonable, and thus, it cannot be used to determine the edge weight.

Comparison analysis with the method based on score values of PHFEs

Pan et al. [22] defined the weight of edges based on the distribution of degree between user and object. In most cases, users do not directly express their preferences for objects, but give their demand. Our method can deal with this situation effectively. To satisfy user’s requirement, we use the evaluation information of experts to determine the possible election of users. Not only that a novel edge weight calculation method based on satisfaction degree of alternatives is proposed in this paper. This method takes uncertainty of recommendation decision information into consideration. To prove the validity of the proposed method, a comparison with the edge weight calculating method based on score function of alternatives is conducted in the following.

Step 1. According to requirements of users, the absolute frequency matrix is obtained and shown in Table 3.

Step 2. We calculate the score values of PHFEs in Table 5, and the results are shown in Table 8.

Table 8 The score value matrix of alternatives

The value in \( i{\text{th}} \) row and \( j{\text{th}} \) column is the score of PHFE, where the score can be calculated by Eq. (3).

Step 3. First, determine the set of adjacent edges according to Eq. (18):

$$ \begin{aligned}\varGamma &= \left\{ o\left( {1,1} \right),o\left( {1,4} \right),o\left( {1,5} \right),o\left( {2,1} \right),o\left( {2,2} \right),o\left( {2,3} \right),\right.\\ &\quad \left.o\left( {2,4} \right),o\left( {3,3} \right),o\left( {4,1} \right),o\left( {4,2} \right),o\left( {4,3} \right),o\left( {4,4} \right) \right\}.\end{aligned} $$

Then, the weight of adjacent edge \( \phi_{ki}^{{}} \) between alternative \( A_{i} \left( {i = 1,2, \ldots ,n} \right) \) and user \( U_{k} \left( {k = 1,2 \ldots ,s} \right) \) can be calculated by \( \phi_{ki} = \sum\nolimits_{j = 1}^{m} {E_{kj} \varphi_{ij} } \), where \( E_{kj} \) is the score value in Table 8. The weights of adjacent edge are shown in Table 9.

Table 9 The weight matrix of adjacent edges

Step 4. We establish bipartite network between users and alternatives based on PHFEs. According to Eqs. (20) and (21), the final projection result is shown in Fig. 7.

Fig. 7
figure 7

Bipartite network projection based on score values of PHFEs

It is easy to see from Fig. 7 that \( g\left( {U_{h} } \right) = g\left( {U_{4} } \right) \) and \( F\left( {A_{o} } \right) = f^{\prime}(A_{2} ) \).

Step 5. It is prior to recommend alternative \( A_{2} \) for user \( U_{ 3} \). Because user \( U_{ 3} \) is similar to user \( U_{ 4} \) and user \( U_{ 4} \) gives a high evaluation for alternative \( A_{2} \).

Step 6. End.

We can find that the result is the same as that by our proposed method, which implies that the proposed method is reasonable. Different from the edge weight determination method based on score values of PHFE, the proposed method takes the deviation degree of decision-makers into account, which makes the results more comprehensive and reliable.

Comparison analysis with recommendation decision-making algorithm based on HFSs and bipartite network projection

To further explicate the superiority of using PHFSs, we make a comparison with the algorithm based on HFSs and bipartite network projection. Correspondingly, the hesitant fuzzy satisfaction degree in reference [54] is used to calculate the adjacent weights.

Step 1. According to requirements of users, the absolute frequency matrix is shown in Table 3. The corresponding requirement relative frequency matrix is shown in Table 4. The alternative is evaluated by experts based on the required attributes of users, and the corresponding hesitant fuzzy evaluation matrix is obtained which is shown in Table 10.

Table 10 Hesitant fuzzy evaluation matrix

Step 2. Calculate the satisfaction degree of HFE, and the satisfaction degree matrix can be derived, which is shown in Table 11.

Table 11 The satisfaction degree matrix of evaluation information

Step 3. We establish bipartite network between users and alternatives based on HFEs. Users’ requirements for attributes and the satisfaction degree of HFEs are considered simultaneously. Specific steps of recommendation method are as follows:

First, determine the set of adjacent edges according to Eq. (18):

$$\begin{aligned} \varGamma &= \left\{ o\left( {1,1} \right),o\left( {1,4} \right),o\left( {1,5} \right),o\left( {2,1} \right),o\left( {2,2} \right),o\left( {2,3} \right),\right.\\ &\quad\left. o\left( {2,4} \right),o\left( {3,3} \right),o\left( {4,1} \right),o\left( {4,2} \right),o\left( {4,3} \right),o\left( {4,4} \right) \right\}.\end{aligned} $$

Then, calculate the weight of adjacent edge according to Eq. (19) as shown in Table 12.

Table 12 The weight matrix of adjacent edges

Step 4. In the established bipartite graph network, the recommended user \( U_{ 3} \) was filled with gray. If the alternatives have been connected with user \( U_{ 3} \), its initial resource would be assigned a unit resource, otherwise zero. As shown in Fig. 8a, initial resource of alternative \( A_{3} \) is 1, and initial resource of the other alternatives is 0. First, the resource flows from alternative to user, and the resource allocation is shown in Fig. 8b. Then, resource flows back to alternative, and the resource allocation is shown in Fig. 8c.

Fig. 8
figure 8

Bipartite network projection based on satisfaction degree of HFEs

It is easy to see from Fig. 8 that \( g\left( {U_{h} } \right) = g\left( {U_{4} } \right) \) and \( F\left( {A_{o} } \right) = f^{\prime}(A_{2} ) \).

Step 5. It is prior to recommend alternative \( A_{2} \) for user \( U_{ 3} \). Because user \( U_{ 3} \) is similar to user \( U_{ 4} \) and user \( U_{ 4} \) gives a high evaluation for alternative \( A_{2} \).

Step 6. End.

We can find that the result is the same as that by the proposed method, which implies that the proposed method is effective and reasonable. However, as an extension of HFSs, PHFSs contains not only the different membership degrees but also the occurrence probabilities associated with the membership degrees, which reserves more original information than HFSs. Therefore, our method can describe uncertainty of decision-making information better.

Conclusion

Based on users’ decision-making information and experts’ evaluation information, a bipartite network between users and alternatives is established in this paper. And a probabilistic hesitant fuzzy recommendation decision-making method based on bipartite network projection is proposed. The proposed method is utilized to recommend accommodation for Airbnb users in this paper. It is worth mentioning that the satisfaction degree of PHFE is applied to calculate the edge weight of bipartite network. And the priority recommendation of alternative is determined by the resource allocation. We consider the uncertainties of recommendation decision-making information and catch the root cause of users’ choice. Thus, a more reasonable recommendation result could be obtained by our method.

In future research, we will continue to focus on bipartite network projection, and expand the application fields of the proposed method, such as online book recommendation, gourmet recommendation, and so on. In addition, we will extend the proposed method to accommodate probabilistic dual hesitant fuzzy or probabilistic linguistic environments, and apply the proposed method to pattern recognition and online reviews.