arxiv: v1 [cs.gt] 15 Oct 2017

Transcription

1 Stable Matchings in Metric Spaces: Modeling Real-World Preferences using Proximity Hossein Karkeh Abadi Stanford University Balaji Prabhakar Stanford University arxiv: v1 [cs.gt] 15 Oct 2017 Suppose each of n men and n women is located at a point in a metric space. A woman ranks the men in order of their distance to her from closest to farthest, breaking ties at random. The men rank the women similarly. An interesting problem is to use these ranking lists and find a stable matching in the sense of Gale and Shapley. This problem formulation naturally models preferences in several real world applications; for example, dating sites, room renting/letting, ride hailing and labor markets. Two key questions that arise in this setting are: (a) When is the stable matching unique without resorting to tie breaks? (b) If X is the distance between a randomly chosen stable pair, what is the distribution of X and what is E(X)? These questions address conditions under which it is possible to find a unique (stable) partner, and the quality of the stable matching in terms of the rank or the proximity of the partner. We study dating sites and ride hailing as prototypical examples of stable matchings in discrete and continuous metric spaces, respectively. In the dating site model, each man/woman is assigned to a point on the k-dimensional hypercube based on their answers to a set of k questions with binary answers (e.g., like/dislike). We consider two different metrics on the hypercube: Hamming and Weighted Hamming (in which the answers to some questions carry more weight). Under both metrics, there are exponentially many stable matchings when k = log n. There is a unique stable matching, with high probability, under the Hamming distance when k = Ω(n 6 ), and under the Weighted Hamming distance when k > (2 + ϵ) log n for some ϵ > 0. Furthermore, under the Weighted Hamming distance, we show that log(x )/log(n) 1, as n, when k > (1 + ϵ) log n for some ϵ > 0. In the ride hailing model, passengers and cabs are modeled as points on the line and matched based on Euclidean distance (a proxy for pickup time). Assuming the locations of the passengers and cabs are independent Poisson processes of different intensities, we derive bounds on the distribution of X in terms of busy periods at a last-come-first-served preemptive-resume (LCFS-PR) queue. We also get bounds on E(X ) using combinatorial arguments.

2 1 INTRODUCTION The stable marriage problem was first introduced by Gale and Shapley [1962] as a way of modeling the college admissions process, in which students are matched with colleges, and the process of courtship leading to marriage, in which women and men are matched. They introduced two key properties of matchings: stability and optimality. These properties are quite well-known and we will recall them formally later; for now, we proceed informally. Stability captures the requirement that a matching should not pair a man M and a woman W with partners whom they both prefer less than each other. Should this happen, M and W are both incentivized to break up with their assigned partners and match with each other. Gale and Shapley [1962] show that there is always at least one stable matching and present the deferred-acceptance algorithm for finding it. Optimality refers to the quality of a matching in terms of the rank of men in their partners preference lists, and vice versa. The stable marriage problem has also been studied in several other real world settings. One famous example is the National Resident Matching Program (NRMP) [Roth, 1984, 1996] where medical school students are matched to residency programs through a centralized stable matching mechanism. Other examples include online dating [Hitsch et al., 2010], sorority rush [Mongell and Roth, 1991], and school choice [Abdulkadiroglu and Sönmez, 2003]. Knuth [1976] initiated the theoretical analysis of large-scale instances of the stable marriage problem under the random preference list assumption, where the preference lists of each man and woman are drawn independently and uniformly from the set of all permutations. Knuth poses the question of estimating the average number of stable matchings when n, the number of men (equal to the number of women) grows large, and provides an integral formula for the probability that a given matching is stable. Pittel [1989, 1992] evaluated this integral and showed that the average number of stable matchings is asymptotic to e 1 n ln n as n, and that any given woman (or man) has Θ(log n) stable partners, on average. We mention a few other results under the random preference list assumption relevant to our work: Immorlica and Mahdian [2005] proved that if the preference list of each woman has only a constant number of entries, then the number of people with multiple stable partners is vanishingly small. 1 Ashlagi et al. [2015] studied the unbalanced case when there are n men and n +k women, for k 1. They show that, with high probability, 2 the fraction of men and women with multiple stable partners tends to zero as n. This line of work is theoretically very interesting, but preference lists in the real world are rarely drawn at random there can exist a significant correlation in the choices people and organizations make. For example, Roth and Peranson [1999] empirically observed correlations in the NRMP preference lists; the applicants largely prefer the same programs and the programs tend to rank the applicants similarly (i.e., a top applicant in one program was also top-ranked in other programs). They note that these correlations can result in a small set of stable matchings. Holzman and Samet [2014] make the previous observation mathematically precise by assuming each participant picks their preference list from a small set of permutations. While the above correlations capture a sameness in the preferences of people and organizations, in this paper we consider correlations due to proximity. Proximity can arise from a coincidence of likes and interests between members of the two sides of a matching market. For example, each member of a matching market answers a questionnaire describing their likes, dislikes, interests or requirements. The questionnaire can either be the same for both sides of the matching market (e.g., online dating) or different (e.g., renters answer questions describing their preferred properties while lessors describe attributes of their preferred renters). The vector of answers can be viewed as points in a metric space and proximity is equated with distance in the metric space. Each participant in the market ranks members of the other side based on their proximity to the participant, from closest to farthest. Distance also arises naturally in the case of ride 2 1 Roth and Peranson [1999] also empirically observe this phenomenon in the context of candidates interviewing for jobs. 2 We say a sequence of events E n occur with high probability if lim n P(E n ) = 1.

3 hailing, where it is desirable to match a hailer with the closest available car. Thus, a wide variety of real world applications can be modeled in this framework; for example, dating sites, 3 renting/letting, 4 labor markets, 5 and ride hailing. 6 Our results. We analyze stable matchings in discrete and continuous metric spaces as the number of participants grows large. We make distributional assumptions on the distances between the participants (hence on the preference lists) and analyze the number and quality of stable matchings. The quality of a stable matching is captured by how small the distances are between stable partners in the matching. When the metric space is continuous, the stable matching is almost surely unique under very mild and natural distributional assumptions. However, this is not necessarily true in discrete metric spaces. An interesting finding of our work is that a participant (on either side of the market) is at the same distance from their partner in all stable matchings. Thus, it makes sense to consider X, the distance between a randomly chosen stable pair (regardless of which stable matching they re picked from, should there be more than one stable matching). We are interested in the distribution of X and E(X ) as the number of participants grows large. We explore these quantities in the dating sites and ride hailing settings. Dating sites. Suppose the men and women of a community are seeking to get matched to a partner in a dating site. At the time of signing up, participants are usually asked to answer a fixed set of k yes/no questions about their preferences, (e.g., Do you like pets?, Are you a morning person? ). We call the k-bit vector representing a participant s answers to these questions the participant s profile. Each profile can be modeled as a point on the k-dimensional hypercube, Q k. The aim is to match a woman to a man whose profile is closest or most similar to hers. We consider two different metrics on Q k for measuring this similarity: the Hamming distance and the Weighted Hamming distance. The Hamming distance between two profiles is equal to the number of entries at which they disagree. The Weighted Hamming distance weighs some disagreements more; the details are in Section 3. Since the distances are not necessarily distinct, we also assume that each person has a tie-breaking preference list for ranking members of the other side and uses this to break ties. One way to think of the actual preference list of a woman is that it ranks the men by distance, closest first. Men at the same distance are ranked according to her tie-breaking preference list. The men form their preference lists similarly. 7 We consider the setting in which profiles are picked independently and uniformly at random from Q k, and the tie-breaking preference lists are chosen independently and uniformly from the set of all permutations. Let ϵ > 0 be an arbitrary positive number. We shall prove that under both the Hamming and the Weighted Hamming distances, for k < (1 ϵ) log n, the fraction of people with multiple stable partners tends to zero, with high probability, as n. However, if k = log n, there are exponentially many stable matchings. We show that, with high probability, the stable matching is unique under the Hamming distance for k = Ω(n 6 ), and it is unique under the Weighted Hamming distance for k > (2+ϵ) log n, without resorting to tie breaks. 8 We derive a lower bound on X under the Hamming distance. Under the Weighted Hamming distance, we prove that if k > (1 + ϵ) log n, then log X /log n 1 in probability. 3 Tinder ( Zoosk ( 4 Airbnb ( Zillow ( 5 LinkedIn ( 6 Uber ( Lyft ( 7 One way to generalize this model to matching markets with two different questionnaires (one for each side of the market) is to ask each participant to answer their questionnaire and also to indicate their best answers from participants on the other side of the market (e.g., renters and lessors answer their questions and that of an ideal response from the other side). The overall profile is then formed by concatenating the answers to both questionnaires. 8 Tie-breaking represents chance, which, in the context of dating, could reasonably be thought of as being less preferable to choice. In other words, a participant would prefer to find his/her partner from their profile rather than through a process involving a coin flip. 3

4 Ride hailing. Consider the problem of matching passengers and cabs on a street. Let blue and red points on the real line represent the location of passengers and cabs, respectively. Suppose the blue and red points occur according to two independent Poisson processes with respective intensities λ and µ. Each point forms its preference list by ranking points of the other color in an increasing order of their Euclidean distance to it. Holroyd et al. [2009] studied translation-invariant matchings between the points of two d-dimensional Poisson processes with the same intensities (λ = µ). They show the natural algorithm of matching mutually closest pairs of points iteratively yields an almost surely unique stable matching. They analyze the tail behavior of X, the distance between a typical pair of stable partners. In the 1-dimensional case, they derive power law upper and lower bounds for the tail distribution of X. In this paper, we study the stable matching problem between two Poisson processes on the real line in the unbalanced case where λ < µ. We derive bounds on the distribution of X in terms of the busy cycles of a last-come-first-served preemptive-resume (LCFS-PR) queue. 9 Using combinatorial arguments, we prove that E(X) ( 1 + ln µ+λ ) µ λ /(µ λ). The rest of the paper is organized as follows. In Section 2 we define the stable matching problem, introduce relevant notation, and state some known results. In Section 3 we describe the stable matching problem on hypercubes and present our results in this model. In Section 4 we analyze the stable matching problem on the real line. Section 5 concludes the paper. 2 BACKGROUND AND PREVIOUS WORK A community of n men and n women is represented by sets M and W, respectively. Suppose each person x in the community has a strict preference list, x, which ranks members of the opposite gender. Thus, y 1 x y 2 means x prefers y 1 to y 2. A matching µ is a mapping from M W to itself, such that for each man m, µ(m) W {m}, for each woman w, µ(w) M {w}, and for any m,w M W, µ(m) = w implies µ(w) = m. A man or woman x is unmatched under µ if µ(x) = x. A pair (m,w) M W is called a blocking pair for a matching µ if w m µ(m) and m w µ(w). A matching is called stable if it does not have any blocking pairs. If a man m and a woman w are matched to each other in a stable matching, we say w and m are a stable partner of each other. The problem of stable matching was first introduced by Gale and Shapley [1962]. They proved that there always exists a stable matching, which can be found using an iterative algorithm called the deferredacceptance algorithm. This algorithm proceeds in a series of proposals and tentative approvals until there is a one-to-one matching between the men and women. When the women propose, they each end up with the best stable partner they can have in any stable matching. This matching, often called woman-optimal, also pairs each man with his lowest-ranked stable partner. The man-optimal stable matching, which results when the men do the proposing, may be distinct from the woman-optimal stable matching; thus, there may be many stable matchings. Under the random preference list assumption, Pittel [1989, 1992] proved that the average number of stable matchings is asymptotic to e 1 n ln n as n, and each person has Θ(log n) stable partners, on average. The stable marriage problem can be extended to the unbalanced case where the number of men and women is not equal. It is clear that for any stable matching in the unbalanced case, there are some people who remain unmatched. This may also happen in the balanced case if the preference lists of some men or women are not complete. We state the following theorems for ready reference. Theorem 2.1 (Rural Hospital). [McVitie and Wilson, 1970, Roth, 1986] The set of men and women who are not matched is the same for all stable matchings. Theorem 2.2. [Immorlica and Mahdian, 2005] Consider the stable marriage problem with n men and n women. Suppose the preference lists of the women are drawn independently and uniformly at random from the 9 Such a queue is also called a stack [Kelly and Yudovina, 2014]. 4

5 5 set of all orderings of men. For a fixed k 1, let the preference lists of the men be drawn independently and uniformly at random from the set of all ordered lists of any k women. (The k women on two different men s preference lists may be different.) In this setting, the expected number of women who have multiple stable partners is o(n). Theorem 2.3. [Ashlagi et al., 2015] Consider a stable marriage problem with n men and n + k women, for arbitrary k = k(n) 1. Suppose the preference lists of women are drawn independently and uniformly at random from the set of all orderings of men, and the preference lists of men are drawn independently and uniformly at random from the set of all orderings of women. The fractions of men and women who have multiple stable partners tends to zero, with high probability, as n. The independence of the randomly drawn preference lists is the key assumption in the analysis of both Theorem 2.2 and Theorem 2.3. Under this assumption, Theorem 2.2 shows that if the preference lists of one side of the market is limited to a fixed k 1 entries, the fraction of men and women with multiple stable partners is vanishingly small. Theorem 2.3 proves the same result for unbalanced markets where there is a size k 1 discrepancy between the number of men and women. In the following section, we derive similar results for the matching markets with correlated preference lists where each person reveals k 1 bits of information about their preference by answering k yes/no questions. 3 STABLE MATCHING ON HYPERCUBES Consider a dating site with n men and n women, represented by sets M and W. Let S = M W and let k be a positive integer. For each x S, let the k-bit vector representing their profile be denoted by a k (x) = ( a 1 (x),..., a k (x) ) {0, 1} k, where a i (x) = 0 if x s answer to the i th question is no, and a i (x) = 1 otherwise. Thus, each profile is a point on the k-dimensional hypercube, Q k = {0, 1} k. For simplicity, we shall suppress the subscript k from a k whenever it can be inferred. In this setting, participants prefer to be matched to someone with a similar profile. Similarity is measured using two metrics on Q k : The Hamming distance and the Weighted Hamming distance. The Hamming distance d h (a, a ) between a and a equals d h (a, a ) k 1(a i a i ). i=1 The Hamming distance assumes that all questions have the same weight. However, some questions may have higher importance than others. For example, Are you allergic to cats? will likely outweigh Do you like caramel?. The Weighted Hamming distance, d w (a, a ) k 2 i 1(a i a i ), addresses this by assigning different weights to different questions. i=1 Remark. Our results for the Weighted Hamming distance (Theorem 3.4) can be extended to any exponentially decaying weights. Remark. When making statements which apply to both metrics we shall use the notation d(.,.). We shall use d(x,y) to denote the distance between the profiles of participants x and y. The preference list of x is arranged according to distance, as follows: for x,y,y S, y x y d(x,y) d(x,y ).

6 6 Since distances are not necessarily distinct, a tie-breaking rule is needed to strictly order preference lists. As mentioned in the Introduction, participant x uses their tie-breaking list, T x, to break ties. Thus, each woman w, ranks men in increasing order of their distance to her and arranges men at the same distance according to their order in her tie-breaking list, T w. 10 For any x and y in S, T x is not necessarily equal to T y. Let the final strict preference list of user x be denoted by P x. We shall use x to indicate ordering in this list. We are now ready to state The Profile Matching Problem (PMP). Given n men and n women and their strict preference lists, the profile matching problem seeks to find a stable matching between the men and the women. A priori, it seems there may be many stable matchings and multiple stable partners for some women and men. However, we shall see in Lemma 3.1 that the multiple stable matchings, should they exist, are all essentially equal in quality. Suppose µ is a stable matching for the PMP. Let d µ (x) be the distance between x and µ(x), d µ (x) d ( x, µ(x) ). Lemma 3.1. Let µ 1 and µ 2 be two stable matchings for the Profile Matching Problem. Then d µ1 (x) = d µ2 (x) for every x S. Proof. See appendix A.1. According to Lemma 3.1, d µ (x) does not depend on µ. Hence, we shall simply denote d µ (x) by d(x) and call d(x) the matching distance of x. Let the random variable X denote the matching distance of a randomly chosen participant x. We analyze X in the following section. 3.1 The Random Profile Matching Problem (RPMP) We now analyze the PMP under certain distributional assumptions of preference lists and profiles when the number of participants grows large. Our main goals are to understand the following questions: How many questions are needed to find a unique partner for each participant without resorting to tie-breaking? What is the matching distance, X? These questions will be answered under the Hamming and the Weighted Hamming metrics. Probabilistic assumptions. We assume each participant answers each of the k questions equally likely with a yes or a no. Further, the answers to all questions by all the participants are independent. Geometrically, this assumption places the k-bit profile vector of each participant (or, equivalently, the participant) at one of the 2 k vertices of Q k, independently and uniformly at random. The preference lists are then generated based on the distances induced by the above placement and the tie-breaking lists T x, x S. We assume each T x is generated independently and uniformly at random from the set of all orderings of men (or women, depending on x). The RPMP-k. Given n men and n women, each of whose preference lists are generated according to the above probabilistic assumptions, the RPMP-k aims to find a stable matching between the men and the women. Remark. Note that RPMP-0 is equivalent to the standard stable matching problem with randomly generated preference lists. 10 If w breaks ties at random, then T w is a random ordering of all the men.

7 3.2 Our results In this section we present our main results for the RPMP-k. Theorem 3.2 considers the case where k log n and Theorem 3.3 and Theorem 3.4 study larger values of k. Due to page limitation we moved all the proofs to appendix. Theorem 3.2. Consider the RPMP-k for k 1. Fix ϵ > 0. Under any metric on Q k, the following statements hold with high probability: (i) if k < (1 ϵ) log n, the fraction of users with multiple stable partners tends to zero as n, so long as tie-breaking is used; and (ii) if k = log n, there are O(n) users with multiple stable partners and there are exponentially many stable matchings. Theorem 3.3. Under the Hamming distance, with high probability, we have the following: (i) if k(n) = Ω(n 3 ), the fraction of users with multiple stable partners tends to zero as n ; (ii) if k(n) = Ω(n 6 ), the stable matching is unique without resorting to tie-breaking; and (iii) for any β > 1, ( P X < k 2 ) βk log n n 1 β. Theorem 3.4. Fix ϵ > 0. Under the Weighted Hamming distance, with high probability, we have the following: (i) if k(n) > (1 + ϵ) log n, the fraction of users with multiple stable partners tends to zero as n. Moreover, log X log n p 1, where p represents convergence in probability; and (ii) if k(n) > (2 + ϵ) log n, the stable matching is unique, without resorting to tie-breaking. According to Theorem 3.2, in large instances of the RPMP-k, if users answer even one question (k = 1), the preference lists become skewed so that, with high probability, any given participant has a unique stable partner. This contrasts starkly with the case k = 0, where Pittel [1992] showed that each participant has, on average, Θ(log n) stable partners. In Theorem 3.3 and 3.4 we distinguish the statements the fraction of participants with a unique stable partner goes to 1 with high probability from the statement there is a unique stable matching, since the former does not imply the latter. Moreover, our method of proving the latter consists of proving the following two steps: (i) if the distances of each man from a given woman are distinct, then she will have a unique stable partner (see Lemma A.4); and (ii) if this holds for all the women (or all the men), then the stable matching is unique. From a market design perspective the uniqueness of the stable matching is important to achieve a shape prediction of the market. Theorem 3.3 shows that under the Hamming distance, if k(n) = Ω(n 6 ), with high probability, there exists a unique stable matching without resorting to tie-breaking. However, asking that many questions from users is not feasible. On the bright side, Theorem 3.4 shows that if the answers to questions carry different weights, we can achieve a unique stable matching with k(n) = O(log n) questions. These theorems also study the matching distance, X. It will be clear from the proof of Theorem 3.2 that X = 0, with high probability, when k < (1 ϵ) log n. Theorem 3.3 establishes an upper bound on the matching distance X. Theorem 3.4 covers the case of the Weighted Hamming metric. Remark. All above theorems can be extended to unbalanced markets with n men and n + r women. 7

8 4 STABLE MATCHING ON THE LINE Consider the problem of matching passengers and cabs on a street. Suppose the passengers and cabs are represented as blue and red points, respectively, on R. Let B and R denote the set of blue and red points, respectively. Let S = B R. A matching between B and R is a mapping M from S to S { }, such that for every red point r, M(r) B { }, for every blue point b, M(b) R { }, and for every b, r S, M(r) = b implies M(b) = r. A point x S is unmatched if M(x) =. The preference list of each point is based on its Euclidean distance to the points with a different color, closest first. A matching M is stable if there is no pair (b, r) B R such that b M(r) and r b < min ( r M(r), b M(b) ). For any matching M and any point x S, let I M (x) R denote the open interval which has x and M(x) at its end-points, and let d M (x) represent the length of I M (x), i.e., d M (x) = x M(x). We call I M (x) the matching segment of x, and d M (x) the matching distance of x in M. With the above definitions, suppose that points in B and R occur according to independent Poisson processes with rates µ and λ, respectively, where λ µ. We call the matching problem defined above as the Poisson Matching problem and denote it by PM(λ, µ). As mentioned in the Introduction, Holroyd et al. [2009] studied translation-invariant matchings between two d-dimensional Poisson processes with the same intensities; in particular, they studied stable matchings. They showed that the following algorithm finds a unique stable matching: Each blue point simultaneously emits two rays, one in each direction, such that at any time t, each ray is at distance t from its emitter. Once a ray hits an unmatched red point r, the emitter b will be matched to r, and both points leave the system. Denote the unique stable matching by M s and let x B be an arbitrary blue point. Define the random variable X to be x s matching distance in M s, i.e., X d Ms (x). Holroyd et al. [2009] proved that if µ = λ, E(X 1/2 ) =. Theorem 4.1. [Holroyd et al., 2009] Let B and R be independent 1-dimensional Poisson processes of intensity 1, and let X represent the matching distance of an arbitrary point in the stable matching between B and R. We have, for some constant C (0, ). E(X 1/2 ) = and P(X > r) Cr 1/2 r > 0, In this section we analyze the 1-dimensional PM(λ, µ) problem for λ < µ. This models the situation in which there are fewer passengers than cabs and sheds light on the time it would take for a passenger to be picked up by the nearest cab that is assigned to pick up the passenger. 11 Thus, we shall be interested in the distribution (Theorem 4.5) and the expected value (Theorem 4.7) of X. However, in order to get at these quantities, we need to introduce various ideas such as the relationship among PM(λ, µ), last-come-firstserved preemptive-resume (LCFS-PR) queue, and nested matchings. We believe these ideas are interesting in their own right. 4.1 Queue Matching Red partners in a stable matching may be either to the left or to the right of the corresponding blue points. However, in queue matchings they are either only on the left or only on the right. Consider PM(λ, µ) with the constraint that each blue point can only be matched to red points that are on its right. In the passenger-cab scenario, this constraint can be the result of having a one-way street or a road divider, where each cab can only pick up passengers on its left. In order to find the stable matching, all the blue points 11 Note the nearest cab may not be able to pick up a passenger since it may be assigned to pick up another passenger who is nearer to the cab than the first passenger. Hence, stable matchings are quite natural in this setting. 8

9 9 Fig. 1. Different matchings between two point processes. Top: forward queue matching M + s, Middle: backward queue matching M s, Bottom: stable matching M s. simultaneously emit a ray to their right at time 0. Once a ray hits an unmatched red point r, the emitter b will be matched to r. It is clear that this algorithm is equivalent to running an LCFS-PR queue where the time of job arrivals and departures in this queue are represented as blue points and red points, respectively. The arrival rate is λ and the service rate is µ (the service times are i.i.d. exponentials of rate µ). We call the resulting stable matching, M + s, the forward queue matching, corresponding to running the queue forward in time. Similarly, we can define a backward queue matching, M q, where each blue point is matched to a red point on its left, and can be found by running the LCFS-PR queue backward in time. Figure 1 shows M + s, M s, and M s for an instance of the problem. The following are well-known facts about LCFS-PR queues with rate λ Poisson arrivals and rate µ > λ i.i.d. exponential service times which are independent of the arrival process. Since λ < µ, the queue is stable and each blue point in B almost surely has a partner in R. Let x B be an arbitrary blue point and let X + be x s matching distance in M + s. It is clear that X + has the same distribution as the busy cycle in the corresponding LCFS-PR queue, where the busy cycle is the duration of time from the arrival of a job at an empty queue to the time the job leaves the queue. It is known [Gross and Harris, 1998] that the probability density function of the busy cycle is given by f λ, µ (t) = 1 t ρ e(λ+µ)t I 1 (2t λµ), t > 0, where ρ = λ µ, and I 1 is the modified Bessel function of the first kind. Let B(λ, µ) represent this distribution. The average busy cycle duration is E(X + ) = 1/(µ λ). In the following section we introduce a class of matchings which includes both stable and queue matchings. 4.2 Nested Matching For any interval I R, represent its closure by Ī. A matching M is said to be nested if for any x,y S, x I M (y) implies M(x) I M (y). Therefore, in any nested matching if I M (x) I M (y), then one of the matching segments is nested inside the other one. Remark. Since the matching segment of an unmatched point x is (x, ), there is no matching segment of a matched point in a nested matching which contains an unmatched point. From the discussion in the previous section it is easy to see that any queue matching is nested. The following lemma proves that the stable matching M s is also nested. Lemma 4.2. The stable matching M s is nested. Proof. See appendix A.2.

10 Let A be the set of all nested matchings between points in B and R. We say a red point r is a potential match for a blue point b, if there exists a nested matching in which b is matched to r. For any blue point b B define P(b) to be the set of all potential matches of b, P(b) = {r R : M A s.t. M(b) = r }. The following lemma shows that the set of potential matches of any two blue points are either disjoint or the same. Lemma 4.3. For any b 1,b 2 B, P(b 1 ) P(b 2 ) implies P(b 1 ) = P(b 2 ). 10 Proof. See appendix A.2. Now define the relation on points in B as follow: b 1 b 2 P(b 1 ) P(b 2 ). According to Lemma 4.3, for any b 1,b 2 B, if P(b 1 ) and P(b 2 ) are not disjoint, then they are the same. Therefore, is an equivalence relation on B. For any blue point b B, define [b] to be b s equivalence class in B, i.e., [b] {b B : P(b ) = P(b)}. Let N (b) represent the size of b s equivalence class, i.e., N (b) = [b]. Also define N + (b) = {b [b],b b} and N (b) = {b [b],b b}. It is clear that N (b) = N + (b) + N (b) 1. In the following lemma we prove some facts about the structure of the equivalence classes. Lemma 4.4. Suppose λ < µ and let B and R represent the set of blue and red points in a Poisson matching problem PM(λ, µ), respectively. For any given blue point b B we have (i) there exist r 1, r 2 P(b) such that r 1 > b and r 2 < b, almost surely; (ii) there exists exactly one potential red point r P(b) between every two consecutive blue points in [b]; (iii) P(b) = N (b) + 1 on { P(b) < }; and (iv) N + (b) and N (b) are independent geometric random variables with parameter 1 λ/µ. Proof. See appendix A.2. Let b B be an arbitrary blue point. Since λ < µ, then almost surely P(b) <. From Lemma 4.4, we can conclude that blue and red points in [b] P(b) form a finite sequence {w i }, for 2N (b)+1 i 2N + (b) 1, where b 0 = b, [b] = {w i : i is even}, and P(b) = {w i : i is odd}. In other words, this sequence starts with a potential red point, alternates between points in [b] and P(b), and ends with another potential red point. We call the sequence {w i } b s potential wave and denote it by W(b). Figure 2 shows potential waves of an instance of PM(λ, µ). 12 A key observation here is that in any nested matching, any blue point in [b] should be matched to a red point in P(b). Therefore, a nested matching first partitions S into potential waves and then matches points within each wave, separately. In the following section we present our results on the analysis of the matching distance X in the stable matching M s. 12 This instance is the same as the instance in Figure 1. As we can see in all the matchings shown in Figure 1, each point is matched to a point within the same potential wave shown in Figure 2.

11 11 Fig. 2. Different potential waves of a sample problem. 4.3 Matching distance, X The following theorem, proves bounds on the distribution of X, in terms of busy cycles. Theorem 4.5. Consider an instance of a Poisson matching problem PM(λ, µ), where λ < µ. Then we have { 2N 1 } 2N 1 min{b 1, B 1 } < X < max B i, where B i and B i are i.i.d. random variables with distribution B(λ, µ), and N and N are i.i.d. geometric random variables with parameter 1 λ/µ. i=1 i=1 B i Proof. See appendix A.2. Using Theorem 4.5 we can find the following upper bound for the expected matching distance E(X ). Corollary 4.6. Proof. See appendix A.2 E(X ) ( 1 + µ + λ ) µ λ 1 µ λ. Remark. Note that the results of Theorem 4.5 and Corollary 4.6 also hold if X is the matching distance in any nested matching. In the next theorem we improve the upper bound given in Corollary 4.6 for the expected matching distance E(X). The proof of this theorem is extensive and requires some detailed combinatorial arguments. For more details see appendix A.3. Theorem 4.7. For the stable matching M s, we have ( E(X) Proof. See appendix A ln ( µ + λ µ λ ) ) 1 µ λ. In order to evaluate the goodness of the bound in Theorem 4.7, note that for large values of µ (µ λ), with a high probability, each blue point will be matched to the closest red point to it. Therefore, as µ/λ, X converges to an exponential distribution with rate 2µ (minimum of two i.i.d. exponentials with rate µ) and E(X ) 1 2µ. However, from Theorem 4.7, in the limit as µ/λ, E(X) is upper bounded by 1 µ.

12 5 CONCLUSION This paper introduced a model for studying matching markets in which preference lists are drawn according to distances in appropriate metric spaces, either between the profiles of participants or between the participants themselves. The model naturally captures several aspects of real world matching markets. Various results regarding the uniqueness and quality of stable matchings were obtained. Specifically, for matchings on the hypercube under the Hamming and Weighted Hamming distances, lower and upper bounds were obtained on the dimension of the hypercube (equal to the number of questions a participant in a dating site needs to answer) so as to obtain unique stable partners or stable matchings. Furthermore, bounds on the distribution and the average value of the matching distance of a typical participant (a measure of the quality of the stable matching) were obtained for stable matchings on the hypercube and on the real line. We view this work as a first step in studying matching markets in the metric space setting. Several obvious next steps suggest themselves, notably studying the problem under dynamic inputs; i.e., as participants arrive and depart. A APPENDIX A.1 Proofs omitted from section 3 Proof of Lemma 3.1: Assume, by contradiction, that there exist stable matchings µ 1 and µ 2 so that for some z S, r 1 = d µ1 (z) > d µ2 (z) = r 2. Let r = (r 1 + r 2 )/2. For each x S, let P x = {y P x : d(x,y) r } be the preference list P x truncated to contain only those participants who are at a distance less or equal to r from x. Let the ordering in the truncated list P x be denoted by x. Call the PMP restricted to the truncated preference lists as the truncated matching problem. In the truncated matching problem, each person prefers to remain unmatched than to match with a person at a distance greater than r from them. Let µ be a stable matching for the PMP which has stable partners with a matching distance greater than r. Construct the partial matching µ from µ by removing all pairs with a distance greater than r. We show that µ is a stable matching for the truncated matching problem. Suppose m and w are not matched to each other in µ. If d(m,w) > r, then clearly (m,w) cannot form a blocking pair for the truncated matching problem. Suppose d(m,w) r. Since m and w are not matched to each other in µ, they cannot be matched to each other in µ. Moreover, since µ is stable, either µ(m) m w, or µ(w) w m. Without loss of generality, assume µ(m) m w. Therefore, d ( m, µ(m) ) d(m,w) r µ(m) m w. This implies m is also matched to µ(m) in µ ; i.e., µ (m) = µ(m). Therefore, µ (m) m w and (m,w) cannot be a blocking pair for µ. This proves that µ is a stable matching for the truncated matching problem. Now define S r (µ) to be the set of all users who are matched to someone at a distance greater than r in µ, S r (µ) = { x S : d µ (x) > r }. It is clear that S r (µ) is the same set of users who are not matched in µ. By the Rural Hospital Theorem, the set of unmatched men and women in the truncated matching problem is the same in all stable matchings. This implies S r (µ) does not depend on µ. This contradicts our initial assumption, since z S r (µ 2 ) but z S r (µ 1 ), proving the lemma. 12 Consider the RPMP-k. Let M and W represent the set of n men and n women, respectively. For any profile a = (a 1,..., a k ) Q k, let M a and W a be the sets of all men and women whose profiles equal a, respectively. Define S a = M a W a.

13 Lemma A.1. Fix a Q k and without loss of generality assume M a W a. We claim that in every stable matching, each man in M a will be matched to a woman in W a. Proof. Suppose to the contrary that there is a stable matching µ and an m M a such that µ(m) W a. Since M a W a, there should also exist a woman w W a such that µ(w) M a. However, since d(m,w) = 0, w m µ(m) and m w µ(w). Therefore, (m,w) forms a blocking pair for µ, which is a contradiction. Thus, for any a Q k, every stable matching should first try to match men in M a with women in W a according to their tie-breaking preference lists. Any one unmatched woman in W a will be matched to someone at a further distance. Define U a to be the set of all users with profile a, which are matched to someone with a profile different from a. Note that according to the Rural Hospital Theorem, U a is the same for all stable matchings and U a = Ma W a. Define Sa = M a W a. The following lemma shows that if k is constant, then for any a Q k, U a = O( n). Lemma A.2. For any arbitrary profile a Q k, as n, M a W a 2p(1 p)n d Z 1 and S a 2np 2p(1 p)n d Z 2, where p = 2 k, Z 1 and Z 2 are independent standard normal N(0, 1) random variables, and d represents convergence in distribution. Proof. First note that M a and W a are i.i.d. with a Binomial(n,p) distribution. From Central Limit Theorem (CLT) we have that as n, M a np np(1 p) d N 1, W a np np(1 p) d N 2, where N 1 and N 2 are two independent random variables with a standard normal distribution, i.e., N 1, N 2 N(0, 1). Therefore, as n, M a W a 2p(1 p)n d N 1 N 2 2, S a 2np 2p(1 p)n d N 1 + N 2 2. Define Z 1 = (N 1 N 2 )/ 2 and Z 2 = (N 1 + N 2 )/ 2. It is clear that Z 1, Z 2 N(0, 1). Moreover, since N 1 and N 2 are independent, Z 1 and Z 2 are also independent. This completes the proof. Since Lemma A.1 requires each stable matching µ to first match men and women in S a using their tie-breaking preference lists, the O( n) discrepancy between the number of men and women in S a makes this sub-problem significantly unbalanced. Using the approach of Ashlagi et al. [2015], we prove some useful bounds on the number of stable partners in unbalanced matching problems which is true for every n. Lemma A.3. Let r 1 and consider an unbalanced two-sided matching problem with n men and n + r women represented by M and W, respectively. Suppose the men s preference lists are generated independently and uniformly at random from the set of all orderings of women in W. Similarly, suppose the women s preference lists are generated independently and uniformly at random from the set of all orderings of men in M. For any given x M W, let N (x) represent the number of x s stable partners. We have that P ( N (x) > 1 ) 1 r + 1 and E ( N (x) ) r. 13

14 Proof. Let µ M represent the men-optimal stable matching found by running the men-proposing deferred acceptance algorithm, and let U be the set of all women who are not matched in µ M. According to the Rural Hospital Theorem, the set of women who are unmatched is the same as U for all stable matchings. Let w W\U be an arbitrary woman. In order to find all the stable partners of w, we employ the same algorithm that is used in McVitie and Wilson [1970], Immorlica and Mahdian [2005], and Ashlagi et al. [2015]. It has been proved by Immorlica and Mahdian [2005] that the following algorithm outputs all the stable partners of w. Algorithm I (1) Run the men-proposing algorithm to find the men-optimal stable matching µ M. If w is unmatched in µ M, output. Initialze µ = µ M. (2) Set m = µ(w) and output m as one of the stable partners of w. Then have w reject m and remove the pair (m, w) from µ. Set u = m. (u represents the current unmatched man.) (3) Let w be the next woman in u s preference list whom he has not proposed to yet. If w is unmatched in µ M, terminate the algorithm. (4) (a) If w has already received a proposal from someone better than, she simply rejects u and the algorithm continues to step 3. (b) If not, w accepts u s proposal. If w = w, the algorithm continues to step 2. Otherwise, set u = µ(m ) and the algorithm continues to step 3. In order to analyze algorithm I, we use the principle of deferred decision which assumes that the random preference lists are not known in advance and rather unfold step by step in the algorithms when a proposal/rejection happens. Let t i be the time of the i th visit of the algorithm at step 3, and define u i and w i to be the unmatched man and the next woman who u i wants to propose to at time t i. Also define X i to be the set of all women who u i has not proposed to yet at time t i. Since we are using the principle of deferred decision, at any time t i, rankings of women in X i are not yet unfolded in u i s preference list. Therefore, at any time t i, every woman in X i has the same chance of 1/ X i to receive the next proposal from u i. Define the events E i = { w i {w} U }. Since the algorithm has not been terminated by time t i, U X i. Therefore, given E i, the probability that u i proposes to w is at most 1/(r + 1), and the probability that the algorithm terminates is at least r/(r + 1), i.e., and P(w i = w E i) 1 r + 1, P(The algorithm terminates at t i E i ) r r + 1. As the algorithm progresses, woman w finds a new stable partner only if she receives a proposal from an unmatched man at step 3 of the algorithm. Let V i be the total number of proposals received by woman w from time t 1 to time t i. If E i does not occur then V i+1 = V i, and if E i occurs then V i+1 = V i + 1 with a probability of at most 1/(r + 1) and the algorithm terminates with a probability at least r/(r + 1). Therefore, if V represents the total number of proposals received by w after time t 1, V is stochastically dominated by a geometric random variable with rate p = r/(r + 1). Thus, P(V > 0) 1 r + 1, 14

15 and E(V ) 1 p p Since N(w) 1 + V, the proof is complete for any w W. It remains to prove the inequalities for x M. Fix x M. Note that the two events {N(x) > 1} and {N(µ M (x)) > 1} are equivalent. Therefore, Since µ M (x) W, P ( N(x) > 1 ) = P ( N(µ M (x)) > 1 ) 1 r + 1. Moreover, since m M N(m) = w W\U N(w) (both are equal to the total number of stable partner pairs), from symmetry we have, = 1 r. 15 E (N (x)) = 1 n E (N(m)) = 1 n m M E (N(w)) r. w W\U We now prove Theorem 3.2 by using Lemmas A.2 and A.3. Proof of Theorem 3.2: Part (i). We prove this part of the theorem only for constant profile size k 1. The proof for arbitrary profile size k (1 ϵ) log n is similar. Fix n and consider an instance of the random profile matching problem with n men and n women represented by M (n) and W (n), respectively. Let x M (n) W (n) be an arbitrary user and let a Q k represent his/her profile. Define M a (n) S a (n) as before. Also define Z (n) 1 = ( S a (n) 2np)/ 2np(1 p) and Z (n) 2 = ( W a (n) M a (n) According to Lemma A.2, as n, Z (n) 1 Z 1 and Z (n) 2 are independent. For any ϵ > 0 define subsets A ϵ, B ϵ R as follows, d, W a (n), and )/ 2np(1 p). d Z 2, where Z 1, Z 2 N(0, 1), and Z 1 and Z 2 A ϵ = [, 1 ϵ ] [1 ϵ, + ], B ϵ = [0, ϵ] [ 1 ϵ, + ]. Let δ > 0 be an arbitrary positive number. Choose ϵ > 0 small enough to have, P (Z 1 A ϵ ) δ 8 and P ( Z 2 B ϵ ) δ (n) 8. Since Z 1 and Z (n) 2 converge in distribution to Z 1 and Z 2, respectively, there exists a large number N 1 such that for any n > N 1, P(Z1 A ϵ ) P(Z (n) 1 A ϵ ) δ < 8, P(Z2 B ϵ ) P(Z (n) 2 B ϵ ) δ < 8. Therefore, for any n > N 1 we have, ( ) P Z (n) 1 A ϵ or Z (n) 2 B ϵ P(Z (n) 1 A ϵ ) + P(Z (n) 2 B ϵ ) P(Z 1 A ϵ ) + P(Z (n) 1 A ϵ ) P(Z 1 A ϵ ) + P(Z 2 B ϵ ) + P(Z (n) 2 B ϵ ) P(Z 2 B ϵ ) δ 2. Therefore, with a probability of at least 1 δ/2, the following event occurs: { E = S (n) a [2np C 1 n, 2np + C1 n] and W (n) a M a (n) [C2 n,c3 n] },

16 16 where C 1 = 2p(1 p) ϵ, C 2 = ϵ p(1 p), and C 3 = p(1 p) ϵ. Without loss of generality assume M a (n) W a (n) and define r = W a (n) M a (n). The problem of matching men in M a (n) and women in W a (n) according to preference lists given by P is an unbalanced matching problem with a discrepancy equal to r between the number of men and the number of women. Let U a (n) represent the set of unmatched women in the men-optimal stable matching for this unbalanced matching problem. Therefore, if N (n) (x) represents the number of x s stable partners, we have, P ( N (n) (x) > 1 ) P(E o ) + P ( N (n) (x) > 1 ) E = P(E o ) + P ( x U (n) ) ( E + P N (n) (x) > 1 E, x U (n)) P(E o ) + r S a (n) + 1 r + 1 δ 2 + C 2 n 1 +, 2np C 1 n C 1 n + 1 a where in the first inequality we used the results of the Lemma A.3, and in the last inequality we used the bounds on r and S a (n) given by the event E. Pick N 2 large enough to have, C 2 n 1 + < δ/2, n N 2. 2np C 1 n C 1 n + 1 Define N = max{n 1, N 2 }. Therefore, for any n > N, P ( N (n) (x) > 1 ) δ. Since δ > 0 is arbitrary, P ( N (n) (x) = 1 ) 1 as n. This implies that with high probability the fraction of users with multiple stable partners tends to zero as n. Part (ii). In order to prove the second part of the theorem, note that since M a and W a are Binomial random variables with parameters n and p = 2 k = 1/n, according to the well-known Poisson limit theorem, both converge to the Poisson(1) distribution as n goes to infinity. Therefore, in the limit, with a positive probability of c = e 2 /4 there are exactly two men and two women whose profiles are equal to a. On the other hand, it is easy to see that in a random stable matching problem with two men and two women, the probability of having exactly two stable matchings is equal to 1/8. Therefore, for any given profile a Q k, with a positive probability of β = c/8 > 0, there are exactly two men and two women with profile a who have multiple stable partners. This proves that the expected number of users with multiple stable partners is O(n). Moreover, since the number of such profiles is O(n), in expectation there are exponentially many stable matchings. The following lemma shows that if the preference list of a user is uniquely identified by profile distances and no further tie-breaking is required, then he/she has a unique stable partner. Lemma A.4. In a profile matching problem, if the distances of a given user x from all the members of the opposite sex are distinct, then x has a unique stable partner. Proof. By contradiction, suppose x has two different stable partners y 1 and y 2. According to Lemma 3.1, y 1 and y 2 should be at the same distance from x. But, this contradicts with the assumption that x has different distances from y 1 and y 2. Therefore, x has a unique stable partner. In order to apply Lemma A.4, k should be large enough to have a unique stable matching without resorting to tie-breaks. Now we prove Theorems 3.3 and 3.4. a

17 Proof of Theorem 3.3: Part (i). Let x be an arbitrary user and without loss of generality, assume x W. Suppose x has multiple stable partners and let y and y be two different stable partners of x. Since x has multiple sable partners, y also has another stable partner x (different from x). According to Lemma 3.1, d h (x,y ) = d h (x,y) = d h (x,y). For any z M define the following event E z = { x W\{x}, z M\{z};d h (x, z ) = d h (x, z) = d h (x, z) }. Using the union bound we have, ( ) ( ) P(E z ) = P x W\{x};d h (x, z) = d h (x, z) P z M\{z};d h (x, z) = d h (x, z ) ((n 1) P ( d h (x,m) = d h (x, z) )) ( (n 1) P ( d h (w, z) = d h (x, z) )) n 2 P ( d h (x,m) = d h (x, z) ) 2 where m and w are a man and a woman who are chosen randomly from M and W, respectively. Note that in the last inequality we used the existing symmetry in the problem. Since d h (x,m) has a binomial distribution (as a function of the random variable m), the maximum value of P ( d h (x,m) = d h (x, z) ) is at d h (x, z) = k/2. Using the Sterling approximation we have: Therefore, P ( d h (x,m) = k ) k! = 2 k 2! k 2!2 k 2πk( k e )k πk( k 2e )k 2 k = 2 π 1 k. P(x has multiple stable partners) P ( z M E z ) n 3 2 πk. Since k = Ω(n 3 ), the right hand side of the above inequality tends to zero as n goes to infinity. This implies that with high probability the fraction of users with multiple stable partners tends to zero as n goes to infinity. Note that using the union bound, we can conclude that if k = Ω(n 4 ), with high probability, there exists a unique stable matching. Part (ii). Fix a woman x W. For any y M, define the event E y = { y M\{y};d h (x,y) = d h (x,y )}. Also define A x to represent the event that the distances of x from all men in M are distinct. Similar to part (i) we have P(A c x ) P ( ) y M E y n 2 2 kπ. According to Lemma A.4, if the event A x happens for every x W, there is a unique stable matching without resorting to tie-breaking. Therefore, P(There are multiple stable matchings) P ( x W A c ) x n 3 2 kπ. Since k = Ω(n 6 ), the probability that there are multiple stable matchings goes to zero as n goes to infinity. Part (iii). Without loss of generality assume x W and let X i represent the distance of x from man m i, i.e., d i = d h (x,m i ). Clearly, X i s are i.i.d. with Binomial distribution with parameters k and 1/2. Define Z = min i X i. Clearly d h (x) Z. Therefore, for any positive number r > 0, P(d h (x) r) P(Z r) = P(min i X i r) = n P(X i r) = P(X 1 r) n. i=1 17