This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

11institutetext: UNSW Sydney, Australia
22institutetext: Data 61, CSIRO, Australia 33institutetext: TU Berlin, Germany

From Matching with Diversity Constraints to Matching with Regional Quotas

Haris Aziz 1122    Serge Gaspers Zhaohong Sun 11221122    Toby Walsh 112233
Abstract

In the past few years, several new matching models have been proposed and studied that take into account complex distributional constraints. Relevant lines of work include (1) school choice with diversity constraints where students have (possibly overlapping) types and (2) hospital-doctor matching where various regional quotas are imposed. In this paper, we present a polynomial-time reduction to transform an instance of (1) to an instance of (2) and we show how the feasibility and stability of corresponding matchings are preserved under the reduction. Our reduction provides a formal connection between two important strands of work on matching with distributional constraints. We then apply the reduction in two ways. Firstly, we show that it is NP-complete to check whether a feasible and stable outcome for (1) exists. Due to our reduction, these NP-completeness results carry over to setting (2). In view of this, we help unify some of the results that have been presented in the literature. Secondly, if we have positive results for (2), then we have corresponding results for (1). One key conclusion of our results is that further developments on axiomatic and algorithmic aspects of hospital-doctor matching with regional quotas will result in corresponding results for school choice with diversity constraints.

1 Introduction

Real-life matching markets are often associated with various distributional constraints. In view of these constraints, there is a growing literature on matching markets that models and deals with such constraints. There are at least two distinct research directions in this growing literature.

The first one is school choice with diversity constraints, studied intensely in the controlled school choice problem, in which students have types such as race, gender, or socio-economic status. Each school is endowed with a lower and an upper quota for each distinct type. Such type-specific quotas are taken into account while determining the outcome. For example, a school may impose a target lower quota for accepting students from some disadvantaged group. The seeds for considering models where students may be of different types were already sown in the seminal paper on school choice (Abdulkadiroğlu and Sönmez, 2003). Since the publication of the paper, there have been significant developments on work concerning fairness requirement and algorithm design (Abdulkadiroğlu, 2005; Kojima, 2012; Hafalir et al., 2013; Ehlers et al., 2014). The most general model in the line of work is school choice with overlapping types where students can belong to multiple types (Kurata et al., 2017). To overcome the non-existence of feasible and stable outcomes, type-specific quotas of schools may be viewed as soft requirements (Ehlers et al., 2014; Kurata et al., 2017).

Another research direction arises in the context of hospitals and doctors matching with a restriction on the number of doctors that are allowed to be matched to certain subsets of hospitals. This form of distributional constraints can be modeled as hospital-doctor matching with regional quotas, in which doctors are matched to hospitals, hospitals are associated with regions, and both hospitals and regions are subject to quotas. For example, an upper quota may be imposed on urban regions with several hospitals to ensure that enough doctors are hired in rural regions (Kamada and Kojima, 2015, 2017a). Individual minimum quotas are studied in school admissions, motivated by the fact that each school may require a minimum number of students to operate (Biró et al., 2010; Fragiadakis et al., 2016). Under general regional quotas, the set of stable outcomes may be empty (Kamada and Kojima, 2017a, b) and it is NP-complete to check whether there exists a feasible outcome (Goto et al., 2014, 2016). Due to these negative results, most work concentrate on special cases with restrictions on the structure of regions for which they proposed algorithms (Kamada and Kojima, 2012, 2015; Goto et al., 2014, 2015, 2016).

Although both lines of work have progressed in the past few years, their development has been generally distinct from each other. Since each of the lines of work stems from different real-life requirements, there has not been much work on identifying formal connections between different new models. In particular, several influential papers mention that one setting is different from the other. For example, Hafalir et al. (2013) note in their seminal paper on school choice with diversity constraints that

“Kamada and Kojima (2011) study the Japanese Residency Matching Program, where there are quotas (regional caps) on the number of residents that each region can admit. […] Although the idea of their paper is similar to ours, the setups are completely different (for instance, there are no doctor types in their model) as are the suggested solutions.”

And Goto et al. (2016) remark in their paper on matching with regional minimum and maximum quotas that

“However, models and theoretical properties in a controlled school choice program setting are quite different from the setting used in our paper.”

The remarks were made because the two models address different concerns and intuitively appear different as well. In this paper, however, we demonstrate that although the two setups discussed above seem different, there are strong mathematical connections between them. Identifying formal connections between matching models have several advantages (1) they help unify the literature, and (2) they provide an efficient route to translate results from one model to another. In fact, one of the major success stories of matching markets has been the identification of general structure over the preferences of hospitals that guarantees the existence of stable matchings (Hatfield and Milgrom, 2005; Hatfield and Kojima, 2008).

\Longstack[c]School Choice with
Diversity Constraints
\Longstack[c]Hospital-doctor Matching
with Regional Quotas
ReductionNP-completenessAlgorithm
Figure 1: Implications of the reduction from school choice with diversity constraints to hospital-doctor matching with regional quotas

Contributions

In this paper, we present a polynomial-time reduction to transform an instance of (1) school choice with diversity constraints into an instance of (2) hospital-doctor matching with regional quotas. We show how the feasibility and stability of corresponding matchings are preserved under the reduction. Then we apply the reduction in two ways as described in Figure 1. First, we study the complexity issues on computing a feasible and stable outcome. We prove that it is NP-complete to check the existence of feasible and stable outcomes for (1). Our reduction implies that these complexity results hold for (2) as well. In view of this, we help unify some of the results that have been presented in the literature. Second, if we have positive results, such as polynomial-time algorithms that guarantee the existence of some weakly stable outcomes for the model with regional quotas, then we have corresponding results for school choice with diversity constraints. One key conclusion of our results is that further developments on axiomatic and algorithmic aspects of hospital-doctor matching with regional quotas will result in corresponding results in school choice with diversity constraints. In addition, we consider how to convert regional minimum quotas into regional maximum quotas and we show the difference between regional minimum quotas and regional maximum quotas.

2 Model

School choice

An instance ISI^{S} of the basic school choice problem consists of a tuple (S,C,qC,𝒳,S,C)(S,C,q_{C},\mathcal{X},\succ_{S},\succ_{C}).

There is a set of students S={s1,s2,,sn}S=\{s_{1},s_{2},...,s_{n}\} and a set of schools C={c1,c2,,cm}C=\{c_{1},c_{2},...,c_{m}\}. Each school cCc\in C has a capacity qcq_{c} and let qC=(qc)cCq_{C}=(q_{c})_{c\in C} be a capacity vector consisting of all schools’ capacities.

Each contract x=(s,c)x=(s,c) is a student-school pair indicating that student ss is matched with school cc. Let 𝒳S×C\mathcal{X}\subseteq S\times C denote the set of available contracts. For any X𝒳X\subseteq\mathcal{X}, denote Xs=X_{s}= {(s,c)X|cC}\{(s,c)\in X|c\in C\} as the set of contracts involving student ss and Xc={(s,c)X|sS}X_{c}=\{(s,c)\in X|s\in S\} as the set of contracts involving school cc in XX.

Each student ss has a strict preference ordering s\succ_{s} over 𝒳s{(s,)}\mathcal{X}_{s}\cup\{(s,\emptyset)\} where (s,)(s,\emptyset) denotes the option of being unmatched for student ss. A contract (s,c)(s,c) is acceptable to student ss if (s,c)(s,c) s\succ_{s} (s,)(s,\emptyset) holds. The preference profile of all students is denoted by S={s1,,sn}\succ_{S}=\{\succ_{s_{1}},...,\succ_{s_{n}}\}. Each school cc has a strict priority ordering c\succ_{c} over 𝒳c{(,c)}\mathcal{X}_{c}\cup\{(\emptyset,c)\}, where (,c)(\emptyset,c) represents the option of leaving a seat vacant for school cc. A contract (s,c)(s,c) is acceptable to school cc if (s,c)c(,c)(s,c)\succ_{c}(\emptyset,c) holds. Let C={c1,,cm}\succ_{C}=\{\succ_{c_{1}},...,\succ_{c_{m}}\} denote the priority profile of all schools. Given any two preference (or priority) orderings p\succ_{p} and q\succ_{q}, we say preference (or priority) ordering p\succ_{p} is consistent with preference (or priority) ordering q\succ_{q} if for any two contracts x,yx,y, when xpyx\succ_{p}y holds, it implies xqyx\succ_{q}y.

An outcome (or a matching) XX is a subset of 𝒳\mathcal{X}. Denote Sc(X)={sS|(s,c)X}S_{c}(X)=\{s\in S|(s,c)\in X\} as the set of students matched to school cc and Cs(X)={cC|(s,c)X}C_{s}(X)=\{c\in C|(s,c)\in X\} as the set of schools matched to student ss in the outcome XX. An outcome XX is feasible for ISI^{S} if i) for each student s,|Xs|1s,|X_{s}|\leq 1, ii) for each school c,|Xc|qcc,|X_{c}|\leq q_{c}. A feasible outcome XX is individually rational if each contract (s,c)X(s,c)\in X is acceptable to both student ss and school cc.

A mechanism ϕ\phi is a function that takes an instance as input and returns a matching as an outcome. A mechanism ϕ\phi is strategy-proof for students if any student sSs\in S cannot be admitted to a better school by misreporting his preference.

School choice with diversity constraints

An instance ID\mathit{I^{D}} of school choice with diversity constraints is an extension of school choice problem, denoted by a tuple (S(S,CC,qCq_{C},TT,τ\tau,η¯\overline{\eta},η¯\underline{\eta}, 𝒳\mathcal{X},S\succ_{S},C)\succ_{C}).

Let T={t1,t2,,tk}T=\{t_{1},t_{2},...,t_{k}\} be the type space with |T||S||T|\leq|S|. A type vector τs=(τst)tT\tau_{s}=(\tau_{s}^{t})_{t\in T} of student ss consists of 11’s and 0’s such that τst=1\tau_{s}^{t}=1 if student ss belongs to type tt and τst=0\tau_{s}^{t}=0 otherwise. Let τ\tau be the type matrix of all students’ type vectors. Let η¯c=(η¯ct)tT\overline{\eta}_{c}=(\overline{\eta}_{c}^{t})_{t\in T} be a vector of school cc’s type-specific maximum quotas where η¯ct\overline{\eta}_{c}^{t} is school cc’s maximum quota for type tt. Similarly, η¯c=(η¯ct)tT\underline{\eta}_{c}=(\underline{\eta}_{c}^{t})_{t\in T} is a vector of school cc’s type-specific minimum quotas. Let η¯\overline{\eta} and η¯\underline{\eta} be two matrices consisting of all schools’ type-specific maximum vectors and minimum vectors respectively.

For any two vectors consisting of non-negative integers ω=(ω1,,ωk)\omega=(\omega_{1},...,\omega_{k}) and ω=(ω1,,ωk),\omega^{\prime}=(\omega_{1}^{\prime},...,\omega_{k}^{\prime}), we compare them in the following way: i) ωω\omega\leq\omega^{\prime} if for each i[1,k],ωiωii\in[1,k],\omega_{i}\leq\omega_{i}^{\prime} and ii) ω<ω\omega<\omega^{\prime} if for each i[1,k],ωi<ωii\in[1,k],\omega_{i}<\omega_{i}^{\prime}. An outcome X𝒳X\subseteq\mathcal{X} is feasible for IDI^{D} with diversity constraints if it is feasible for ISI^{S} and it respects diversity constraints, i.e., for each school cCc\in C, we have ηc¯sSc(X)τsηc¯\underline{\eta_{c}}\leq\sum_{s\in S_{c}(X)}\tau_{s}\leq\overline{\eta_{c}}.

The following definition is a natural extension of classical stability concept by Roth (1985) to the setting of diversity constraints.

Definition 1 (Stability)

Given a feasible outcome XX for instance IDI^{D} with diversity constraints, a student sSs\in S and a school cCc\in C with (s,c)X(s,c)\notin X will form a blocking pair if (s,c)sXs(s,c)\succ_{s}X_{s} and there exists a set of students SSc(X)S^{\prime}\subseteq S_{c}(X) such that i) for each student sSs^{\prime}\in S^{\prime}, we have (s,c)(s,c) c\succ_{c} (s,c)(s^{\prime},c) and ii) the new outcome X{(s,c)}(sS(s,c){Xs})X\cup\{(s,c)\}\setminus(\bigcup_{s^{\prime}\in S^{\prime}}(s^{\prime},c)\cup\{X_{s}\}) is feasible for instance IDI^{D}. A feasible outcome XX is stable if it is individually rational and it admits no blocking pair.

Definition 1 states that given a feasible outcome XX, student ss and school cc will form a blocking pair if student ss prefers school cc to his current assignment XsX_{s} (which could be empty if Xs=(s,)X_{s}=(s,\emptyset)), and there exists a set of students SS^{\prime} (which could be empty) that are matched to school cc in the outcome XX such that i) each student sSs^{\prime}\in S^{\prime} has lower priority than student ss and ii) school cc could admit student ss by (possibly) removing the set of students SS^{\prime}.

We can decompose the blocking pair in Definition 1 into two cases: when the set of students SS^{\prime} is non-empty, we say student ss has justified envy towards students SS^{\prime} and a feasible outcome XX is fair if XX admits no justified envy. If the set of students SS^{\prime} is empty, then we say the outcome XX is wasteful. Alternatively, a feasible outcome XX is stable if it is individually rational, fair and non-wasteful.

Hospital-doctor matching

An instance IH\mathit{I^{H}} of hospital-doctor matching is isomorphic to an instance IS\mathit{I^{S}} of basic school choice problem, denoted by a tuple (D,H,qH,𝒴,D,H)(D,H,q_{H},\mathcal{Y},\succ_{D},\succ_{H}).

Let D={d1,,dn}D=\{d_{1},...,d_{n}\} denote a set of doctors and H={h1,,hm}H=\{h_{1},...,h_{m}\} denote a set of hospitals. Let qH=(qh)hHq_{H}=(q_{h})_{h\in H} be the vector consisting of each hospital’s capacity where qhq_{h} is the capacity of hospital hh.

A contract y=(d,h)y=(d,h) is a doctor-hospital pair such that dd is matched with hh. Let 𝒴D×H\mathcal{Y}\subseteq D\times H denote a finite set of available contracts. For any Y𝒴,Y\subseteq\mathcal{Y}, let Yd={(d,h)Y|hH}Y_{d}=\{(d,h)\in Y|h\in H\} be the set of contracts involving doctor dd and Yh={(d,h)Y|dD}Y_{h}=\{(d,h)\in Y|d\in D\} be the set of contracts involving hospital hh.

Each doctor dd has a strict preference ordering d\succ_{d} over 𝒴d{(d,)}\mathcal{Y}_{d}\cup\{(d,\emptyset)\} and a contract (d,h)(d,h) is acceptable to doctor dd if (d,h)d(d,)(d,h)\succ_{d}(d,\emptyset). Each hospital hh has a strict priority ordering over 𝒴h{(,h)}\mathcal{Y}_{h}\cup\{(\emptyset,h)\} and a contract (d,h)(d,h) is acceptable to hospital hh if (d,h)h(,h)(d,h)\succ_{h}(\emptyset,h). Let D\succ_{D} and H\succ_{H} denote the preference and priority profiles of DD and HH respectively.

An outcome (or a matching) YY is a subset of 𝒴\mathcal{Y}. Denote Dh(Y)={dD|(d,h)Y}D_{h}(Y)=\{d\in D|(d,h)\in Y\} as the set of doctors matched to hospital hh and Hd(Y)={hH|(d,h)Y}H_{d}(Y)=\{h\in H|(d,h)\in Y\} as the set of hospitals matched to doctor dd in the outcome YY. An outcome Y𝒴Y\subseteq\mathcal{Y} is feasible for IHI^{H} if i) for each doctor dd, we have |Yd|1|Y_{d}|\leq 1, ii) for each hospital hh, |Yh|qh|Y_{h}|\leq q_{h}, and iii) for any doctor dd and any hospital hh, we have dDh(Y)d\in D_{h}(Y) if and only if Hd(Y)={h}H_{d}(Y)=\{h\}.

Hospital-doctor matching with regional quotas

An instance IR\mathit{I^{R}} of hospital-doctor matching with regional quotas is a tuple (D,H,R,qH,δ¯,δ¯,𝒴,D,H,R)(D,H,R,q_{H},\overline{\delta},\underline{\delta},\mathcal{Y},\succ_{D},\succ_{H},\succ_{R}) with additional entries R,δ¯R,\overline{\delta}, δ¯\underline{\delta} and R\succ_{R}.

Let R={r1,,rj}R=\{r_{1},...,r_{j}\} denote a set of regions where each region riRr_{i}\in R is a subset of HH, i.e., riHr_{i}\subseteq H. A collection of regions PRP\subseteq R forms a partition of a subset of hospitals HHH^{\prime}\subseteq H, if rPr=H\bigcup_{r\in P}r=H^{\prime} and for any two different regions r,rPr,r^{\prime}\in P, we have rr=r\cap r^{\prime}=\emptyset. A collection of regions FRF\subseteq R forms a hierarchy of hospitals HHH^{\prime}\subseteq H, if rFr=H\bigcup_{r\in F}r=H^{\prime} and for any two regions r,r(r)Fr,r^{\prime}(\neq r)\in F, one of the three conditions holds: i) rr=r\cap r^{\prime}=\emptyset, ii) rrr\subseteq r^{\prime}, or iii) rrr^{\prime}\subseteq r.

Let δ¯=(δ¯r)rR\overline{\delta}=(\overline{\delta}_{r})_{r\in R} denote a vector consisting of each region’s maximum quota where δ¯r\overline{\delta}_{r} is region rr’s maximum quota. Similarly δ¯=(δ¯r)rR\underline{\delta}=(\underline{\delta}_{r})_{r\in R} is a vector of each region’s minimum quota.

For any Y𝒴,Y\subseteq\mathcal{Y}, let Yr=hrYhY_{r}=\bigcup_{h\in r}Y_{h} be the set of contracts involving region rr and let Dr(Y)={dD|(d,h)Yr}D_{r}(Y)=\{d\in D|(d,h)\in Y_{r}\} denote the set of doctors matched to region rr in the outcome YY.

The introduction of regional priorities was intended to resolve the conflicts when a region confronts more applicants of doctors than it could accommodate (Kamada and Kojima, 2017a, 2018).111In (Kamada and Kojima, 2017a, 2018), regional preferences are defined in a weaker way that regions only concern about how many doctors are matched rather than which doctors are matched. We followed this idea and assume that each region rr has a strict priority ordering over 𝒴r\mathcal{Y}_{r} and a contract (d,h)𝒴r(d,h)\in\mathcal{Y}_{r} is acceptable to region rr if (d,h)h(,h)(d,h)\succ_{h}(\emptyset,h) holds.222Similar ideas to regional preferences have already been considered. Kamada and Kojima (2015) studied the model where all regions form a partition of hospitals HH and they assume that each region specifies a precedence ordering over hospitals. Let R\succ_{R} denote the priority profile of all regions.

An outcome Y𝒴Y\subseteq\mathcal{Y} is feasible for IRI^{R} with regional quotas if YY is feasible for IHI^{H} and it respects regional quotas, i.e., for any region rr we have δ¯r|Dr(Y)|δ¯r\underline{\delta}_{r}\leq|D_{r}(Y)|\leq\overline{\delta}_{r}.

The following stability concept captures the idea that a blocking pair is not considered as legitimate if it does not take regional priorities into account.

Definition 2 (Stability with regional priorities)

Given a feasible outcome Y𝒴Y\subseteq\mathcal{Y} for instance IRI^{R} with regional quotas, a doctor dDd\in D and a hospital hHh\in H with (d,h)Y(d,h)\notin Y form a blocking pair with regional priorities if (d,h)(d,h) d\succ_{d} (d,Yd)(d,Y_{d}) and there exists a set of doctors DDh(Y)D^{\prime}\subseteq D_{h}(Y) such that i) for each doctor dDd^{\prime}\in D^{\prime}, we have (d,h)h(d,h)(d,h)\succ_{h}(d^{\prime},h), ii) for each doctor dDd^{\prime}\in D^{\prime} and for each region rr with hrh\in r, we have (d,h)r(d,h)(d,h)\succ_{r}(d^{\prime},h), and iii) the new outcome Y{(d,h)}(dD{(d,h)}Yd})Y\cup\{(d,h)\}\setminus(\bigcup_{d^{\prime}\in D^{\prime}}\{(d^{\prime},h)\}\cup Y_{d}\}) is feasible for instance IRI^{R}. A feasible outcome YY is stable with regional priorities if it is individually rational and it admits no blocking pair with regional priorities.

Definition 2 states that given a feasible outcome YY, doctor dd and hospital hh will form a blocking pair if doctor dd prefers hospital hh to his assigned hospital Hd(Y)H_{d}(Y) (which could be empty), and there exists a set of doctors DD^{\prime} (which could be empty) that are matched to hospital hh in the outcome XX such that i) each doctor dDd^{\prime}\in D^{\prime} has lower priority than doctor dd, ii) for each region rr that is associated with hospital hh, each doctor dDd^{\prime}\in D^{\prime} has lower regional priority than doctor dd and iii) hospital hh could admit doctor dd by (possibly) removing the set of doctors DD^{\prime}.

The difference between Definition 1 and Definition 2 is that a blocking pair for an instance of matching with regional quotas should respect the priorities of both hospitals and regions. When distribution constraints do not exist, both definitions collapse to the original stability concept (Roth, 1985).

We can decompose the blocking pair in the Definition 2 into two cases: when the set of doctors DD^{\prime} is non-empty, we say doctor dd has justified envy towards DD^{\prime} with regional priorities and an outcome is fair with regional priorities if it admits no justified envy. When the set of doctors DD^{\prime} is empty, we say the outcome is wasteful.

3 Transformation from Diversity Constraints to Regional Quotas

In this section, we explore the relation between (1) school choice with diversity constraints and (2) hospital-doctor matching with regional quotas in terms of feasibility and stability. We show how to convert an instance of (1) into an instance of (2) in polynomial time and how the feasibility and stability of corresponding matchings are preserved under the reduction.

In a recent work, Kamada and Kojima (2017b) illustrate how to associate one instance of (1) with another of (2) when each student belongs to exactly one type. The idea is straightforward: Each student corresponds to a doctor and each school corresponds to a region. For each region, create multiple hospitals such that each hospital is associated with one type and each doctor only considers the hospital of the same type acceptable. However, we cannot directly extend this idea to the general case allowing for overlapping types. The main issue is that a doctor should not be assigned to several hospitals corresponding to the types to which he belongs.

The crux of the transformation is how to eliminate overlapping types among students. We can just create a new type space 𝒯={t1,,t2|T|}\mathcal{T}=\{t_{1}^{\prime},...,t_{2^{|T|}}^{\prime}\} such that each unique type vector τs\tau_{s} corresponds to a new type tτst_{\tau_{s}}^{\prime}. It is not necessary to consider the whole type space 𝒯\mathcal{T} when 2|T|2^{|T|} is larger than |S||S|, because only the distinct type vectors that appear in type matrix τ\tau matter, of which the maximum number is no more than min(|S|,2|T|)\min(|S|,2^{|T|}), bounded by the number of students |S||S|. Let 𝒯\mathcal{T}^{*} be such a new type space induced from τ\tau with |𝒯|min(|S|,2|T|)|\mathcal{T}^{*}|\leq\min(|S|,2^{|T|}). Then we can assign a student ss with type vector τs\tau_{s} one new type tτs𝒯t_{\tau_{s}}^{\prime}\in\mathcal{T^{*}} and no two students have overlapping types.

Now we proceed to the polynomial-time reduction from an instance of diversity constraints ID=(S,C,qC,T,τ,η¯,η¯,𝒳,S,C)I^{D}=(S,C,q_{C},T,\tau,\overline{\eta},\underline{\eta},\mathcal{X},\succ_{S},\succ_{C}) to a corresponding instance of regional quotas IR=(D,H,R,qH,δ¯,δ¯I^{R}=(D,H,R,q_{H},\overline{\delta},\underline{\delta}, 𝒴,D,H,R)\mathcal{Y},\succ_{D},\succ_{H},\succ_{R}).

For each student sjSs_{j}\in S, create a corresponding doctor djd_{j}. Let D=sjSdjD=\bigcup_{s_{j}\in S}d_{j} denote the set of doctors. For each school ciCc_{i}\in C and each unique type vector τs\tau_{s} from type matrix τ\tau, create a hospital hciτsh^{\tau_{s}}_{c_{i}} with capacity qciq_{c_{i}}. Let Hci𝒯H^{\mathcal{T}^{*}}_{c_{i}} be the set of hospitals induced from school cic_{i}. Denote the set of hospitals as H=ciCHci𝒯H=\bigcup_{c_{i}\in C}H^{\mathcal{T}^{*}}_{c_{i}} with capacity vector qH=(qh)hHq_{H}=(q_{h})_{h\in H}.

For each school ciCc_{i}\in C, create a set of regions of size (|T||T| ++ 11) denoted by Rci={ri,ri1,,ri|T|}R_{c_{i}}=\{r_{i},r_{i}^{1},...,r_{i}^{|T|}\} where region rir_{i} corresponds to school cic_{i} and region rijr_{i}^{j} corresponds to type tjt_{j} at school cic_{i}. Region rir_{i} contains all hospitals induced from cic_{i}, i.e., Hri=Hci𝒯H_{r_{i}}=H^{\mathcal{T}^{*}}_{c_{i}} and each region rijr_{i}^{j} contains the hospitals induced from school cic_{i} that are associated with type tjt_{j}, i.e., Hrij={hiτsHi𝒯|τstj=1}H_{r_{i}^{j}}=\{h^{\tau_{s}}_{i}\in H^{\mathcal{T}^{*}}_{i}|\tau_{s}^{t_{j}}=1\}. The maximum and minimum regional quotas for each region are described in Table 1. Let R=ciCRciR=\bigcup_{c_{i}\in C}R_{c_{i}} denote the set of regions with δ¯=(δ¯r)rR\overline{\delta}=(\overline{\delta}_{r})_{r\in R} and δ¯=(δ¯r)rR\underline{\delta}=(\underline{\delta}_{r})_{r\in R}.

Table 1: Regions RciR_{c_{i}} induced from school cic_{i}
rir_{i} rijr_{i}^{j}
maximum quota δ¯ri=qci\overline{\delta}_{r_{i}}=q_{c_{i}} δ¯rij=η¯citj\overline{\delta}_{r_{i}^{j}}=\overline{\eta}_{c_{i}}^{t_{j}}
minimum quota δ¯ri=0\underline{\delta}_{r_{i}}=0 δ¯rij=η¯citj\underline{\delta}_{r_{i}^{j}}=\underline{\eta}_{c_{i}}^{t_{j}}
related hospitals ri=Hci𝒯r_{i}=H^{\mathcal{T}^{*}}_{c_{i}} rij={hciτsHci𝒯|τstj=1}r_{i}^{j}=\{h^{\tau_{s}}_{c_{i}}\in H^{\mathcal{T}^{*}}_{c_{i}}|\tau_{s}^{t_{j}}=1\}

For each contract x=(s,ci)𝒳x=(s,c_{i})\in\mathcal{X}, create a new contract y=(d,hciτs)y=(d,h^{\tau_{s}}_{c_{i}}) with doctor dd corresponding to student ss and hospital hciτsh^{\tau_{s}}_{c_{i}} corresponding to type vector τs\tau_{s} of student ss. Let 𝒴=x𝒳{y}\mathcal{Y}=\bigcup_{x\in\mathcal{X}}\{y\} be the set of available contracts. Each doctor dd’s preference ordering d\succ_{d} corresponds to s\succ_{s}. For each hospital hHci𝒯h\in H^{\mathcal{T}^{*}}_{c_{i}} and each region rRcir\in R_{c_{i}}, the preference orderings of h\succ_{h} and r\succ_{r} are consistent with ci\succ_{c_{i}} over corresponding contracts involving hospital hh and region rr. The preference profiles of doctors, hospitals and regions are denoted by D,H,R\succ_{D},\succ_{H},\succ_{R} respectively.

Proposition 3.1

The reduction takes time O(|C||S|2)O(|C|\cdot|S|^{2}).

Proof

The running time of the construction depends on the number of all induced elements. The number of doctors, hospitals and regions is O(|S|+|C||𝒯|+|C|(|T|+1))O(|S|+|C|\cdot|\mathcal{T}^{*}|+|C|\cdot(|T|+1)). The number of capacities and type-specific quotas is bounded by O(|H|+|R|)O(|H|+|R|). The number of contracts is at most |C||S||C|\cdot|S|. The number of preference orderings of doctors, hospitals and regions is O(|C||S|+|C||𝒯||S|+|C||S||T|)O(|C|\cdot|S|+|C|\cdot|\mathcal{T}^{*}|\cdot|S|+|C|\cdot|S|\cdot|T|) where |𝒯||S||\mathcal{T}^{*}|\leq|S| and |T||S||T|\leq|S|. Thus the running time of the reduction is O(|C||S|2)O(|C|\cdot|S|^{2}).

Example 1

We illustrate the reduction with the following example. Consider an instance IDI^{D} with diversity constraints:

S={s1,s2,s3,s4},C={c},qc=2,T={t1,t2},τs1=(0,0),τs2=(0,1),τs3=(1,0),τs4=(1,1),η¯c=(1,1),η¯c=(1,0),𝒳={(s1,c),(s2,c),(s3,c),(s4,c)},sS(s,c)s(s,),(s1,c)c(s2,c)c(s3,c)c(s4,c)c(,c).\begin{array}[]{l}S=\{s_{1},s_{2},s_{3},s_{4}\},C=\{c\},q_{c}=2,T=\{t_{1},t_{2}\},\tau_{s_{1}}=(0,0),\\ \tau_{s_{2}}=(0,1),\tau_{s_{3}}=(1,0),\tau_{s_{4}}=(1,1),\overline{\eta}_{c}=(1,1),\underline{\eta}_{c}=(1,0),\\ \mathcal{X}=\{(s_{1},c),(s_{2},c),(s_{3},c),(s_{4},c)\},\forall s\in S\ (s,c)\succ_{s}(s,\emptyset),\\ (s_{1},c)\succ_{c}(s_{2},c)\succ_{c}(s_{3},c)\succ_{c}(s_{4},c)\succ_{c}(\emptyset,c).\end{array}

Create a corresponding instance IRI^{R} with regional quotas as follows, where region rr corresponds to school cc, region r1r_{1} corresponds to type t1t_{1} and region r2r_{2} corresponds to type t2t_{2} at school cc.

D={d1,d2,d3,d4},H={h00,h01,h10,h11},hHqh=2,R={r,r1,r2},r=H,r1={h10,h11},r2={h01,h11},δr¯=2,δr1¯=δr2¯=1,δr¯=δr2¯=0,δr1¯=1,𝒴={(d1,h00),(d2,h01),(d3,h10),(d4,h11)},dD𝒴dd(d,),hH𝒴hh(,h),(d1,h00)r(d2,h01)r(d3,h11)r(d4,h11),(d3,h10)r1(d4,h11),(d2,h01)r2(d4,h11).\begin{array}[]{l}D=\{d_{1},d_{2},d_{3},d_{4}\},H=\{h_{00},h_{01},h_{10},h_{11}\},\forall h\in H\ q_{h}=2,\\ R=\{r,r_{1},r_{2}\},r=H,r_{1}=\{h_{10},h_{11}\},r_{2}=\{h_{01},h_{11}\},\\ \overline{\delta_{r}}=2,\overline{\delta_{r_{1}}}=\overline{\delta_{r_{2}}}=1,\underline{\delta_{r}}=\underline{\delta_{r_{2}}}=0,\underline{\delta_{r_{1}}}=1,\\ \mathcal{Y}=\{(d_{1},h_{00}),(d_{2},h_{01}),(d_{3},h_{10}),(d_{4},h_{11})\},\\ \forall d\in D\ \mathcal{Y}_{d}\succ_{d}(d,\emptyset),\forall h\in H\ \mathcal{Y}_{h}\succ_{h}(\emptyset,h),\\ (d_{1},h_{00})\succ_{r}(d_{2},h_{01})\succ_{r}(d_{3},h_{11})\succ_{r}(d_{4},h_{11}),\\ (d_{3},h_{10})\succ_{r_{1}}(d_{4},h_{11}),(d_{2},h_{01})\succ_{r_{2}}(d_{4},h_{11}).\end{array}

The relationship between two instances is shown in Figure 2. The index of each hospital is in binary corresponding to each distinct type vector. Note that four copy schools {c00,c01,c10,c11}\{c_{00},c_{01},c_{10},c_{11}\} are used for interpretation only which are not actually involved in the reduction.

ccc00c_{00}c10c_{10}c11c_{11}c01c_{01}t1t_{1}t2t_{2}rrh00h_{00}h10h_{10}h11h_{11}h01h_{01}r1r_{1}r2r_{2}
Figure 2: An example of reduction.

Next we show how the feasibility and stability of corresponding outcomes are preserved under the reduction. Given an outcome XX of instance IDI^{D} with diversity constraints, create an outcome YY of induced instance IRI^{R} with regional quotas by adding a corresponding contract yy to YY for each contract xXx\in X.

Proposition 3.2

The outcome XX is feasible for IDI^{D} with diversity constraints if and only if the induced outcome YY under the reduction is feasible for IRI^{R} with regional quotas.

Proof

If outcome XX is feasible for IDI^{D} with diversity constraints, then for each school cCc\in C, we have |Sc(X)|qc|S_{c}(X)|\leq q_{c} and η¯csSc(X)τsη¯c\underline{\eta}_{c}\leq\sum_{s\in S_{c}(X)}\tau_{s}\leq\overline{\eta}_{c}. Since each school cic_{i} corresponds to a region rir_{i}, then no more than qciq_{c_{i}} doctors are matched to region rir_{i} in outcome YY, which implies that each hospital hciτsrih^{\tau_{s}}_{c_{i}}\in r_{i} admits no more than qciq_{c_{i}} doctors. Since each type jj at school cic_{i} corresponds to a region rijr_{i}^{j}, then we have δ¯rij|Drij(Y)|δ¯rij\underline{\delta}_{r_{i}^{j}}\leq|D_{r_{i}^{j}}(Y)|\leq\overline{\delta}_{r_{i}^{j}}, i.e., the outcome YY respects the regional quotas of each region rijr_{i}^{j}.

If outcome YY is feasible for IRI^{R} with regional quotas, then for each region rr we have δ¯r|Dr(Y)|δ¯r\underline{\delta}_{r}\leq|D_{r}(Y)|\leq\overline{\delta}_{r}. Since each region rir_{i} corresponds to a school cic_{i}, then no more qciq_{c_{i}} students are matched to cic_{i} in outcome XX. Since each region rijr_{i}^{j} corresponds to one type jj at school cic_{i}, then all type-specific quotas of school cic_{i} are satisfied in outcome XX.

Proposition 3.3

The outcome XX is stable for instance IDI^{D} with diversity constraints if and only if the induced outcome YY under the reduction is stable with regional priorities for instance IRI^{R} with regional quotas.

Proof

If outcome XX is stable, for the sake of contradiction, suppose outcome YY admits a blocking pair (d,h)(d,h) with regional priorities induced from student ss and school cc respectively. Let DDh(Y)D^{\prime}\subseteq D_{h}(Y) be the set of doctors such that each dDd^{\prime}\in D^{\prime} has lower priority than dd at hospital hh as well as at each region rr with hrh\in r, and hospital hh can admit doctor dd by removing doctors DD^{\prime}. Let the set of students SS^{\prime} correspond to doctors DD^{\prime}. Then student ss and school cc could form a blocking pair, since school cc could admit student ss by removing SS^{\prime}, a contradiction.

If outcome YY is stable with regional priorities, suppose outcome XX admits a blocking pair (s,c)(s,c) and let SS^{\prime} be the set of students such that each sSs^{\prime}\in S^{\prime} has lower priority than ss at school cc and school cc could admit student ss by removing SS^{\prime}. Let doctor dd, region rr and a set of doctors DD^{\prime} correspond to student ss, school cc and set SS^{\prime} of students respectively. Then doctor dd and hospital hh could form a blocking pair with regional priorities through a set of doctors DD^{\prime}, since all induced hospitals HcH_{c} and regions RcR_{c} from school cc have the consistent priority orderings as school cc, a contradiction.

4 Transformation from regional minimum quotas to regional maximum quotas

In this section, we further show how to transform an instance of regional maximum and minimum quotas into a corresponding instance of regional maximum quotas only in terms of feasibility.

Goto et al. (2014, 2016) considered how to represent regional minimum quotas with regional maximum quotas in a restrictive setting where any doctor is acceptable to any hospital and vice versa. In addition, the total capacity of all hospitals is at least the number of doctors and no doctors are unmatched in any feasible outcome. Their idea works as follows: If region rr requires at least δ¯r\underline{\delta}_{r} doctors, then the number of doctors that can be assigned to other hospitals which do not belong to region rr cannot exceed |D|δ¯r|D|-\underline{\delta}_{r}. However, this does not hold in general if we relax these requirements.

Example 2

There are two doctors d1,d2d_{1},d_{2}, two hospitals h1,h2h_{1},h_{2} and two regions r1={h1},r2={h2}r_{1}=\{h_{1}\},r_{2}=\{h_{2}\} with δ¯r1=δ¯r2=δ¯r1=δ¯r2=1\underline{\delta}_{r_{1}}=\underline{\delta}_{r_{2}}=\overline{\delta}_{r_{1}}=\overline{\delta}_{r_{2}}=1. Following the reduction of (Goto et al., 2014, 2016), after removing regional minimum quotas, the regional quotas for the induced instance become δ¯r1=δ¯r2=1\overline{\delta}_{r_{1}}=\overline{\delta}_{r_{2}}=1. Then an empty outcome is feasible for the induced instance but not for the original one.

Next we generalize their idea to general setting without any assumption by adding an additional null hospital. Since we consider feasibility only, the preference and priority orderings of doctors, hospitals and regions are not necessary. Given a simplified instance of hospital-doctor matching with regional quotas IR=(D,H,R,qH,δ¯,δ¯,𝒴)I^{R}=(D,H,R,q_{H},\overline{\delta},\underline{\delta},\mathcal{Y}), construct an instance with regional maximum quotas only IR+=(D,H,R,qH,δ¯,𝒴)I^{R+}=(D^{\prime},H^{\prime},R^{\prime},q_{H^{\prime}},\overline{\delta^{\prime}},\mathcal{Y^{\prime}}) as follows:

The set of doctors remains the same and a null hospital h0h_{0} is added to HH, i.e., D=DD^{\prime}=D and H=H{h0}H^{\prime}=H\cup\{h_{0}\}. Let qH=(qh)hHq_{H^{\prime}}=(q_{h})_{h\in H^{\prime}} with qh0=|D|q_{h_{0}}=|D^{\prime}|. For each region rRr\in R, create a new region r^=H{r}\hat{r}=H^{\prime}\setminus\{r\} with δ¯r^=|D|δ¯r\overline{\delta}_{\hat{r}}=|D^{\prime}|-\underline{\delta}_{r}. The set of all regions is denoted as R=rR{r,r^}R^{\prime}=\bigcup_{r\in R}\{r,\hat{r}\}. For each doctor dDd\in D^{\prime}, add a new contract (d,h0)(d,h_{0}) to 𝒴\mathcal{Y}, i.e., 𝒴=dD{(d,h0)}𝒴\mathcal{Y}^{\prime}=\bigcup_{d\in D^{\prime}}\{(d,h_{0})\}\cup\mathcal{Y}.

Proposition 4.1

The reduction takes time O(|R|+|D||H|)O(|R|+|D|\cdot|H|).

Proof

The time of copying instance IRI^{R} is bounded by O(|D|+|H|+|R|+|D||H|)O(|D|+|H|+|R|+|D|\cdot|H|). In addition, we create one new hospital h0h_{0} and a set of new regions and contracts, whose total number is bounded by O(|R|+|D||H|)O(|R|+|D|\cdot|H^{\prime}|) where |H|=|H|+1|H^{\prime}|=|H|+1. Thus the construction of instance IR+I^{R+} takes time O(|R|+|D||H|)O(|R|+|D|\cdot|H|).

Proposition 4.2

An outcome YY is feasible for instance IRI^{R} with both regional minimum and maximum quotas if and only if the outcome YY is feasible for the induced instance IR+I^{R^{+}} with regional maximum quotas only.

Proof

If outcome YY respects regional quotas for IRI^{R}, then for each rRr\in R, the number of doctors matched to H{r}H\setminus\{r\} does not exceed δ¯r^=|D|δ¯r\overline{\delta}_{\hat{r}}=|D^{\prime}|-\underline{\delta}_{r}, which implies the region r^\hat{r} respects regional quotas. If YY respects regional quotas for IR+I^{R^{+}}, then for each r^R\hat{r}\in R^{\prime}, the number of doctors matched to r^\hat{r} does not exceed δ¯r^=|D|δ¯r\overline{\delta}_{\hat{r}}=|D^{\prime}|-\underline{\delta}_{r}, which implies at least δ¯r\underline{\delta}_{r} doctors are matched with region rr.

Our reduction reveals an important distinction between regional minimum quotas and regional maximum quotas: In an instance with regional maximum quotas only, any number of doctors can be placed at the null hospital without violating feasibility. But when we translate regional minimum quotas into regional maximum quotas, we limit the maximum number of doctors that can be matched to the null hospital by imposing regional caps to some regions that contain the null hospital.

5 Complexity results

This section is devoted to the complexity results on checking the existence of a feasible and stable outcome for both settings. We first prove NP-completeness for setting (1), then by reduction from setting (1) to setting (2), it implies that these NP-completeness results also hold for setting (2). In view of this, we help unify some complexity results that were already proved in previous literature, and we further show these NP-completeness results still hold under more restrictive settings.

Complexity of computing a feasible outcome

Next we provide a polynomial-time reduction from (3,3)-Set cover problem to school choice problem with diversity constraints. Gonzalez (1985) has proved that (3,3)-Set cover is NP-complete.

(3,3)-Set cover
Input: A collection FF of subsets of a finite set UU and a positive integer kk where each uUu\in U occurs in at most three subsets of FF and each fFf\in F contains at most three elements of UU. Question: Is there a subset FFF^{\prime}\subseteq F of size at most kk such that fFf=U\bigcup_{f\in F^{\prime}}f=U?

Proposition 5.1

It is NP-complete to check the existence of a feasible outcome for school choice problem with diversity constraints, even if there is only one school, each student belongs to at most three types, each type contains at most three students and there is no upper bound for any type.

Proof

Given an instance IDI^{D} with diversity constraints, to decide whether IDI^{D} admits a feasible outcome or not is in NP, since we can guess an outcome XX and check whether XX satisfies feasibility in polynomial-time.

Given an instance (F,U)(F,U) of (3,3)-Set Cover, create a corresponding instance IDI^{D} with diversity constraints as follows: For each element uiUu_{i}\in U, create a type tit_{i}. For each subset fjFf_{j}\in F, create a student sjs_{j}. A student sjs_{j} belongs to type tit_{i} if uifju_{i}\in f_{j}. Create one school cc with capacity qc=kq_{c}=k and minimum quota ηt¯=1\underline{\eta_{t}}=1 for each type tTt\in T. For each student sSs\in S, create a contract (s,c)(s,c) which is acceptable to both ss and cc. Create an arbitrary priority ordering c\succ_{c}.

If (F,U)(F,U) admits a Yes-instance FF^{\prime} of size at most kk, then let SS^{\prime} == fjF\bigcup_{f_{j}\in F^{\prime}} sjs_{j} denote the corresponding set of students. The outcome X=sjS(sj,c)X=\bigcup_{s_{j}\in S^{\prime}}(s_{j},c) is feasible for IDI^{D} with diversity constraints, since school cc admits at most kk students and each minimum type-specific quota is satisfied.

If IDI^{D} with diversity constraints admits a feasible outcome XX, let F=sjSc(X)fjF^{\prime}=\bigcup_{s_{j}\in S_{c}(X)}f_{j} denote the corresponding subsets of FF. Then we have a Yes-instance of (3,3)-Set Cover, since we have |F|k|F^{\prime}|\leq k and fFf=U\bigcup_{f\in F^{\prime}}f=U.

Although Goto et al. (2014) proved it is NP-complete to check the existence of feasible outcomes for setting (2), their original reduction requires both regional minimum quotas and regional maximum quotas (which needs to be equal for each region). Under the assumption that no doctor is unmatched in any feasible outcome (as discussed in last section), they infer that “checking whether a feasible matching exists or not is NP-complete where there are only regional minimum quotas or only regional maximum quotas”. However, if we relax the assumption, it can be done in polynomial-time to check whether a feasible matching exists when there are only regional maximum quotas, since an empty matching satisfies feasibility. With the help of our reduction, we further show this NP-completeness result even holds for a more restrictive setting when there are only minimum quotas in the following corollary:

Corollary 5.1

It is NP-complete to check the existence of a feasible outcome for hospital-doctor matching with regional quotas, even if each doctor is matched to at most three regions, each region admits at most three doctors and contains at most three hospitals, except for one region that contains all hospitals and is matched with all doctors.

In another recent work on public housing allocation with diversity constraints, Benabbou et al. (2018) showed that it is NP-complete to check whether there exists a feasible assignment with maximum social welfare, which is different from ours.

Complexity of computing a stable outcome

Now we move on to the complexity question of deciding the existence of a stable outcome. In previous work, Huang (2010) showed that it is NP-hard to compute a stable matching for school choice with diversity constraints when both minimum and maximum quotas exist. Next we show that if there are no minimum quotas, it is NP-complete to check whether a stable outcome exists under strict preferences. The following reduction is inspired by the work on hospital-doctor matching with couples by Ronn (1990) and McDermid and Manlove (2010).

Proposition 5.2

Given an instance of school choice with diversity constraints in which there are no minimum quotas, it is NP-complete to decide whether there exists a stable outcome under strict preferences, even if there are only two types, the capacity and type-specific maximum quotas for each school are at most 2 and the length of any preference / priority ordering is at most 4. 333In the previous version, we implicitly assume that there are no minimum quotas.

Proof

First we prove that if there are no minimum quotas, deciding whether a stable outcome exists is in NP. We can guess an outcome XX and check whether XX admits a blocking pair in polynomial time as follows: For each student ss and each school cc such that (s,c)sXs(s,c)\succ_{s}X_{s}, if school cc can admit student ss by removing all students who are matched to school cc with lower priority than student ss, then outcome XX is not stable. This is because removing a set of students with lower priority than student ss does not violate feasibility requirement.

Next we show it is NP-hard by reduction from a restricted version of 3-SAT where each literal appears exactly twice, which is NP-complete (Berman et al., 2003).

student type vector preference
s1is^{i}_{1} (1,1)(1,1) c1ict1ic^{i}_{1}\ c^{i}_{t_{1}}
s2is^{i}_{2} (1,1)(1,1) c2icf1ic^{i}_{2}\ c^{i}_{f_{1}}
s3is^{i}_{3} (1,0)(1,0) c1ic2ic^{i}_{1}\ c^{i}_{2}
s4is^{i}_{4} (0,1)(0,1) c2ic1ic^{i}_{2}\ c^{i}_{1}
s5is^{i}_{5} (0,0)(0,0) c1ict2ic^{i}_{1}\ c^{i}_{t_{2}}
s6is^{i}_{6} (0,0)(0,0) c2icf2ic^{i}_{2}\ c^{i}_{f_{2}}
t1it^{i}_{1} (1,0)(1,0) ct1io(t1i)β1,3ic^{i}_{t_{1}}\ o(t^{i}_{1})\ \beta^{i}_{1,3}
t2it^{i}_{2} (0,1)(0,1) ct2io(t2i)β2,3ic^{i}_{t_{2}}\ o(t^{i}_{2})\ \beta^{i}_{2,3}
f1if^{i}_{1} (1,0)(1,0) cf1io(f1i)β3,3ic^{i}_{f_{1}}\ o(f^{i}_{1})\ \beta^{i}_{3,3}
f2if^{i}_{2} (0,1)(0,1) cf2io(f2i)β4,3ic^{i}_{f_{2}}\ o(f^{i}_{2})\ \beta^{i}_{4,3}
αk,1i\alpha^{i}_{k,1} (0,1)(0,1) βk,2iβk,1i\beta^{i}_{k,2}\ \beta^{i}_{k,1}
αk,2i\alpha^{i}_{k,2} (1,0)(1,0) βk,1iβk,2i\beta^{i}_{k,1}\ \beta^{i}_{k,2}
αk,3i\alpha^{i}_{k,3} (1,1)(1,1) βk,3iβk,1i\beta^{i}_{k,3}\ \beta^{i}_{k,1}
Table 2: Students induced from variable uiu_{i} where k[1,4]k\in[1,4]
school capacity maximum quotas priority ordering
c1ic^{i}_{1} 2 (1,1)(1,1) s4is1is3is5is^{i}_{4}\ s^{i}_{1}\ s^{i}_{3}\ s^{i}_{5}
c2ic^{i}_{2} 2 (1,1)(1,1) s3is2is4is6is^{i}_{3}\ s^{i}_{2}\ s^{i}_{4}\ s^{i}_{6}
ct1ic^{i}_{t_{1}} 1 (1,1)(1,1) s1it1is^{i}_{1}\ t^{i}_{1}
ct2ic^{i}_{t_{2}} 1 (1,1)(1,1) s5it2is^{i}_{5}\ t^{i}_{2}
cf1ic^{i}_{f_{1}} 1 (1,1)(1,1) s2if1is^{i}_{2}\ f^{i}_{1}
cf2ic^{i}_{f_{2}} 1 (1,1)(1,1) s6if2is^{i}_{6}\ f^{i}_{2}
βk,1i\beta^{i}_{k,1} 2 (1,1)(1,1) αk,1iαk,3iαk,2i\alpha^{i}_{k,1}\ \alpha^{i}_{k,3}\ \alpha^{i}_{k,2}
βk,2i\beta^{i}_{k,2} 1 (1,1)(1,1) αk,2iαk,1i\alpha^{i}_{k,2}\ \alpha^{i}_{k,1}
βk,3i\beta^{i}_{k,3} 1 (1,1)(1,1) θαk,3i\theta\ \alpha^{i}_{k,3}
Table 3: Schools induced from variable uiu_{i}

Given an instance (U,W)(U,W) of 3-SAT in which each literal appears exactly twice, let U={u1,,uk}U=\{u_{1},\ldots,u_{k}\} denote a set of variables and W={w1,,wl}W=\{w_{1},\ldots,w_{l}\} be a set of clauses. Create an instance IDI^{D} with diversity constraints as follows.

For each variable uiUu_{i}\in U, create a gadget consisting of 22 students and 18 schools as shown in Table 3 and 3. Students t1it^{i}_{1} and t2it^{i}_{2} stand for the first and second occurrence of literal uiu_{i}. Students f1if^{i}_{1} and f2if^{i}_{2} stand for the first and second occurrence of literal u¯i\bar{u}_{i}.

For each clause wjWw_{j}\in W, create exactly one school ojo^{j} with capacity 22 and maximum quota 22 for two types. Let s(l1)ojs(l2)ojs(l3)s(l_{1})\succ_{o^{j}}s(l_{2})\succ_{o^{j}}s(l_{3}) denote the priority ordering of school ojo^{j} where s(lk)s(l_{k}) denotes the corresponding student of literal lkl_{k} that appears in clause wjw_{j}.

Note that o(t1i)o(t^{i}_{1}) and o(t2i)o(t^{i}_{2}) in the preference of student t1it^{i}_{1} and t2it^{i}_{2} stand for two schools induced by the clauses in which uiu_{i} appears for the first and second time. Similarly, o(f1i),o(f2i)o(f^{i}_{1}),o(f^{i}_{2}) correspond to two schools induced from the clauses in which u¯i\bar{u}_{i} appears for the first and second time. In the priority ordering of school βk,3i\beta^{i}_{k,3}, θ\theta stands for tkit^{i}_{k} when k[1,2]k\in[1,2] and stands for fk2if^{i}_{k-2} when k[3,4]k\in[3,4].

Lemma 5.1

If there exists a satisfying assignment α:U{false,true}\alpha:U\rightarrow\{false,true\} of instance (U,W)(U,W) of 3-SAT, then the induced instance IDI^{D} of school choice admits a stable outcome.

Proof

For each gadget induced from variable uiu_{i}, if the value of uiu_{i} is true in the assignment α\alpha, then select the outcome XTiX^{i}_{T}:

XTi={(s1i,c1i),(s2i,cf1i),(s3i,c2i),(s4i,c2i),(s5i,c1i),\displaystyle X^{i}_{T}=\{(s^{i}_{1},c^{i}_{1}),(s^{i}_{2},c^{i}_{f_{1}}),(s^{i}_{3},c^{i}_{2}),(s^{i}_{4},c^{i}_{2}),(s^{i}_{5},c^{i}_{1}),
(s6i,cf2i),(t1i,ct1i),(t2i,ct2i),(f1i,o(f1i)),(f2i,o(f2i)),\displaystyle(s^{i}_{6},c^{i}_{f_{2}}),(t^{i}_{1},c^{i}_{t_{1}}),(t^{i}_{2},c^{i}_{t_{2}}),(f^{i}_{1},o(f^{i}_{1})),(f^{i}_{2},o(f^{i}_{2})),
(αk,1i,βk,2i),(αk,2i,βk,1i),(αk,3i,βk,3i)};\displaystyle(\alpha^{i}_{k,1},\beta^{i}_{k,2}),(\alpha^{i}_{k,2},\beta^{i}_{k,1}),(\alpha^{i}_{k,3},\beta^{i}_{k,3})\};

otherwise select outcome XFiX^{i}_{F}:

XFi={(s1i,ct1i),(s2i,c2i),(s3i,c1i),(s4i,c1i),(s5i,ct2i),\displaystyle X^{i}_{F}=\{(s^{i}_{1},c^{i}_{t_{1}}),(s^{i}_{2},c^{i}_{2}),(s^{i}_{3},c^{i}_{1}),(s^{i}_{4},c^{i}_{1}),(s^{i}_{5},c^{i}_{t_{2}}),
(s6i,c2i),(t1i,o(t1i)),(t2i,o(t2i)),(f1i,cf1i),(f2i,cf2i)\displaystyle(s^{i}_{6},c^{i}_{2}),(t^{i}_{1},o(t^{i}_{1})),(t^{i}_{2},o(t^{i}_{2})),(f^{i}_{1},c^{i}_{f_{1}}),(f^{i}_{2},c^{i}_{f_{2}})
(αk,1i,βk,2i),(αk,2i,βk,1i),(αk,3i,βk,3i)}.\displaystyle(\alpha^{i}_{k,1},\beta^{i}_{k,2}),(\alpha^{i}_{k,2},\beta^{i}_{k,1}),(\alpha^{i}_{k,3},\beta^{i}_{k,3})\}.

For the school ojo^{j} induced from clause wj=(l1,l2,l3)w_{j}=(l_{1},l_{2},l_{3}), if the value of literal lkl_{k} is false in the assignment α\alpha, then match the student s(lk)s(l_{k}) corresponding to literal lkl_{k} to school ojo^{j}.

Next we show that none of induced schools would be part of any blocking pair. First consider any school ojo^{j} induced from clause wjw_{j}. Since the assignment α\alpha is satisfying, no more than two students will be matched to ojo^{j}, otherwise corresponding clause is false. School ojo^{j} would not be part of any blocking pair, since it can admit any two students without violating feasibility. Then consider the schools in the gadget induced by variable uiu_{i}. If we select outcome XTiX^{i}_{T}, then s1i,s4i,s5i,t1i,t2i,c2i,cf1i,cf2i,αk,1i,αk,2i,αk,3is^{i}_{1},s^{i}_{4},s^{i}_{5},t^{i}_{1},t^{i}_{2},c^{i}_{2},c^{i}_{f_{1}},c^{i}_{f_{2}},\alpha^{i}_{k,1},\alpha^{i}_{k,2},\alpha^{i}_{k,3} are matched with their top choices, which implies they cannot be part of any blocking pair. Then we can infer c1ic^{i}_{1} cannot form a blocking pair with s4is^{i}_{4}, ct1ic^{i}_{t_{1}} cannot form a blocking pair with s1is^{i}_{1} and ct2ic^{i}_{t_{2}} cannot form a blocking pair with s5is^{i}_{5}. If we select outcome XFiX^{i}_{F}, then s2i,s3i,s6i,f1i,f2i,c1i,ct1i,ct2i,αk,1i,αk,2i,αk,3is^{i}_{2},s^{i}_{3},s^{i}_{6},f^{i}_{1},f^{i}_{2},c^{i}_{1},c^{i}_{t_{1}},c^{i}_{t_{2}},\alpha^{i}_{k,1},\alpha^{i}_{k,2},\alpha^{i}_{k,3} are matched with their top choices. Then we can infer c2ic^{i}_{2} cannot form a blocking pair with s3is^{i}_{3}, cf1ic^{i}_{f_{1}} cannot form a blocking pair with s2is^{i}_{2} and cf2ic^{i}_{f_{2}} cannot form a blocking pair with s6is^{i}_{6}. Thus any induced school would not be part of any blocking pair.

Lemma 5.2

In any stable outcome XX for IDI^{D}, if (t1i,ct1i)(t^{i}_{1},c^{i}_{t_{1}}) X\in X, then (t2i,ct2i)(t^{i}_{2},c^{i}_{t_{2}}) X\in X; if (f1i,cf1i)(f^{i}_{1},c^{i}_{f_{1}}) X\in X, then (f2i,cf2i)(f^{i}_{2},c^{i}_{f_{2}}) X\in X.

Proof

If (t1i,ct1i)X(t^{i}_{1},c^{i}_{t_{1}})\in X, then we have (s1i,c1i)X(s^{i}_{1},c^{i}_{1})\in X, otherwise student s1is^{i}_{1} and school ct1ic^{i}_{t_{1}} will form a blocking pair. Then we can infer (s5i,c1i)X(s^{i}_{5},c^{i}_{1})\in X, otherwise student s5is^{i}_{5} and school ct1ic^{i}_{t_{1}} will form a blocking pair. Thus (t2i,ct2i)X(t^{i}_{2},c^{i}_{t_{2}})\in X holds, otherwise they will form a blocking pair. Similarly, if (f1i,cf1i)X(f^{i}_{1},c^{i}_{f_{1}})\in X, then we have (s2i,c2i)X(s^{i}_{2},c^{i}_{2})\in X, otherwise student s2s_{2} and school cf1ic^{i}_{f_{1}} will form a blocking pair. Then we can infer (s6i,c2i)X(s^{i}_{6},c^{i}_{2})\in X and (f2i,cf2i)X(f^{i}_{2},c^{i}_{f_{2}})\in X holds, otherwise they will form a blocking pair.

Lemma 5.3

For any stable outcome XX for induced instance IDI^{D}, students t1it^{i}_{1}, t2it^{i}_{2}, f1if^{i}_{1}, f2if^{i}_{2} must be matched with their first two choices.

Proof

Given any stable outcome XX, student t1it^{i}_{1} cannot be unmatched, otherwise t1it^{i}_{1} will form a blocking pair with school βk,3i\beta^{i}_{k,3}. Suppose student t1it^{i}_{1} is matched to βk,3i\beta^{i}_{k,3}. Then there will be two cases: either i) student αk,3i\alpha^{i}_{k,3} is matched to βk,1i\beta^{i}_{k,1} or ii) αk,3i\alpha^{i}_{k,3} is unmatched. For case i), we have student αk,1i\alpha^{i}_{k,1} is matched to βk,2i\beta^{i}_{k,2}, otherwise αk,1i\alpha^{i}_{k,1} and βk,1i\beta^{i}_{k,1} will form a blocking pair. However, αk,2i\alpha^{i}_{k,2} and βk,2i\beta^{i}_{k,2} will form a blocking pair, leading to a contradiction. For case ii), if both students αk,1i\alpha^{i}_{k,1} and αk,2i\alpha^{i}_{k,2} are matched to βk,1i\beta^{i}_{k,1}, then αk,1i\alpha^{i}_{k,1} and βk,2i\beta^{i}_{k,2} will form a blocking pair. If αk,1i\alpha^{i}_{k,1} is matched to βk,2i\beta^{i}_{k,2}, αk,3i\alpha^{i}_{k,3} and βk,1i\beta^{i}_{k,1} will form a blocking pair. Thus for any stable outcome XX, student t1it^{i}_{1} cannot be matched to βk,3i\beta^{i}_{k,3}, and student t1it^{i}_{1} must be matched to his first two choices. The similar arguments work for students t2it^{i}_{2}, f1if^{i}_{1}, f2if^{i}_{2}.

Lemma 5.4

For any stable outcome XX for induced instance IDI^{D}, either (t1i,ct1i)X(t^{i}_{1},c^{i}_{t_{1}})\in X, (f1i,cf1i)X(f^{i}_{1},c^{i}_{f_{1}})\notin X holds or (t1i,ct1i)X(t^{i}_{1},c^{i}_{t_{1}})\notin X, (f1i,cf1i)X(f^{i}_{1},c^{i}_{f_{1}})\in X holds.

Proof

For the sake of contradiction, suppose (t1i,ct1i)X(t^{i}_{1},c^{i}_{t_{1}})\notin X and (f1i,cf1i)X(f^{i}_{1},c^{i}_{f_{1}})\notin X first. Then we have (s1i,ct1i)(s^{i}_{1},c^{i}_{t_{1}}) X\in X and (s2i,cf1i)(s^{i}_{2},c^{i}_{f_{1}}) X\in X, otherwise students t1it^{i}_{1} and f1if^{i}_{1} will form blocking pairs with schools ct1ic^{i}_{t_{1}} and cf1ic^{i}_{f_{1}} respectively. Then we can infer that (s3i,c2i)(s^{i}_{3},c^{i}_{2}) X\in X, (s4i,c1i)(s^{i}_{4},c^{i}_{1}) X\in X, otherwise students s1is^{i}_{1} and s2is^{i}_{2} will form blocking pairs with c1ic^{i}_{1} and cf1ic^{i}_{f_{1}} respectively. However, XX is not stable, since s3is^{i}_{3} and s4is^{i}_{4} can form blocking pairs with c1ic^{i}_{1} and c2ic^{i}_{2} respectively.

Then suppose (t1i,ct1i)(t^{i}_{1},c^{i}_{t_{1}}) X\in X and (f1i,cf1i)(f^{i}_{1},c^{i}_{f_{1}}) X\in X. Then we have (s1i,c1i)(s^{i}_{1},c^{i}_{1}) X\in X and (s2i,c2i)(s^{i}_{2},c^{i}_{2}) X\in X, otherwise s1is^{i}_{1} and s2is^{i}_{2} will form blocking pairs with ct1ic^{i}_{t_{1}} and cf1ic^{i}_{f_{1}} respectively. However, outcome XX is not stable, since students s4is^{i}_{4} and s3is^{i}_{3} can form blocking pairs with schools c1ic^{i}_{1} and c2ic^{i}_{2} respectively.

Lemma 5.5

If there is a stable outcome XX for IDI^{D}, then there is a satisfying assignment α\alpha for the 3-SAT instance (U,W)(U,W).

Proof

For each variable uiu_{i}, if (t1i,ct1i)X(t^{i}_{1},c^{i}_{t_{1}})\in X, then set uiu_{i} to be true; if (f1i,cf1i)X(f^{i}_{1},c^{i}_{f_{1}})\in X, then set uiu_{i} to be false. By Lemmas 5.2, 5.3 and 5.4, this assignment is consistent. Suppose there is a clause wjw^{j} with all literals set to false. Then the corresponding school ojo^{j} must accommodate three students, exceeding the capacity of ojo^{j}, a contradiction.

This concludes the proof of the Proposition 5.2.

Corollary 5.2

It is NP-complete to check whether there exists a stable outcome with regional priorities for hospital-doctor matching without regional minimum quotas, even if each region contains at most 4 hospitals, the capacity of each hospital and the maximum quota of each region is at most 2, and the length of any preference / priority ordering is at most 4.

6 Algorithm design for Weaker Stability

In this section we discuss the second implication of the reduction: how positive results for setting (2) hospital-doctor matching with regional quotas could lead to corresponding results for setting (1) school choice with diversity constraints. In general, stability with regional preferences is too strong to guarantee the existence of stable outcomes and it is NP-complete to decide whether one exists even if there are only regional maximum quotas.

Suppose there exists an algorithm ϕ\phi that takes an instance IRI^{R} with regional quotas as input and returns a feasible and weaker stable outcome with respect to priorities of hospitals and regions. Given an instance IDI^{D} with diversity constraints, we can first convert IDI^{D} into a corresponding instance IRI^{R} with regional quotas in polynomial time. Then apply the algorithm ϕ\phi designed for matching with regional quotas to instance IRI^{R} to obtain some feasible and weaker stable outcome XRX^{R}. By the corresponding relation, we can restore the outcome XRX^{R} to an outcome XDX^{D} of instance IDI^{D} that also satisfies feasibility and some form of weaker stability.

Now the problem boils down to designing a weaker stable concept for setting (2) which should be weak enough to guarantee the existence of feasible and stable outcomes but still strong enough to lead to reasonable outcomes. One possible way, which has been considered in the mechanism design of matching markets Fragiadakis et al. (2016); Goto et al. (2016); Kurata et al. (2017), is to decompose stability into fairness and non-wastefulness, and then weaken one or both of them to obtain some weaker stable concept.

However, as far as we know, there is no convincing stability solution that works for a general instance of matching with regional quotas. Previous work mainly concentrates on special cases where regions form a partition or a hierarchy of hospitals, and they additionally assume there exists a strict master list over doctors that is used to determine which doctor should be matched when conflicts occur Goto et al. (2014, 2015, 2016). Note that a master list is equivalent to imposing unified regional priorities on all regions. Based on such a strict master list over doctors, we can easily design a strategy-proof algorithm that always returns a feasible and stable outcome with regional priorities when there are only regional maximum quotas. For example, we can just employ serial dictatorship, which lets each doctor choose their favorite hospital without violating hospital capacity and regional maximum quotas in the order of master list, to obtain a stable outcome with regional priorities where regional priorities are consistent with master list Goto et al. (2014); Fragiadakis et al. (2016).

The following fairness concept for school choice with diversity constraints is a weaker version of fairness by considering master list. Compared to original fairness definition, it additionally requires that student ss could have justified envy towards a non-empty set of students SS^{\prime} by master list if each student sSs^{\prime}\in S^{\prime} has lower master list priority than student ss.

Definition 3 (Fairness by Master List)

Given an instance IDI^{D} with diversity constraints, a feasible outcome XX for IDI^{D} and a strict master list ML\succ_{ML}, a student ss has justified envy toward a non-empty set of students SSc(X)S^{\prime}\subseteq S_{c}(X) by master list, if the following conditions hold: i) (s,c)s(s,Xs)(s,c)\succ_{s}(s,X_{s}), ii) for each sSs^{\prime}\in S^{\prime}, we have (s,c)c(s,c)(s,c)\succ_{c}(s^{\prime},c) and (s,c)ML(s,c)(s,c)\succ_{ML}(s^{\prime},c) and iii) X{(s,c)}(sS{(s,c)}Xs})X\cup\{(s,c)\}\setminus(\bigcup_{s^{\prime}\in S^{\prime}}\{(s^{\prime},c)\}\cup X_{s}\}) is a feasible outcome for IDI^{D}. A feasible outcome XX is fair by master list if XX does not admit justified envy by master list.

Consider an instance IDI^{D} of school choice without type-specific minimum quotas and a strict master list. If we apply serial dictatorship to its induced instance IRI^{R} with regional quotas, then we obtain a feasible and stable outcome XRX^{R} with regional priorities where each regional priority ordering is consistent with the master list. When we restore the outcome XRX^{R} to outcome XDX^{D} of the original instance IDI^{D}, we have a feasible, non-wasteful and fairness outcome by master list. This serves well for the illustration of how positive results for setting (2) could lead to corresponding results for setting (1). Designing an appropriate stability concept for matching with regional quotas requires more exploration. One key conclusion of our results is that further developments on axiomatic and algorithmic aspects of hospital-doctor matching with regional quotas will result in corresponding results in school choice with diversity constraints.

7 Summary and Discussion

In this paper we provide a formal connection between two important forms of distributional constraints via a polynomial-time reduction. Our reduction has two implications: First, if we have NP-completeness results in the model of school choice with diversity constraints, then these complexity results also carry over to the model with regional quotas. Second, positive results, such as polynomial-time algorithms that guarantee the existence of some weakly stable outcomes for the model with regional quotas, imply corresponding results for school choice with diversity constraints.

Note that our reduction can be generalized to new models appearing in recent literature on school choice with diversity constraints. Matching with slot-specific priorities was proposed in (Kominers and Sönmez, 2016) where each slot (an extension of type) could have different priority orderings. This requires that the induced hospitals HiH_{i} and regions RiR_{i} should have different priority orderings from school cic_{i}. To overcome non-existence of feasible outcomes and improve efficiency, diversity constraints may be regarded as soft bounds and schools could admit more students than the type-specific quotas allows if some seats are unoccupied (Kojima, 2012; Hafalir et al., 2013; Ehlers et al., 2014; Kurata et al., 2017). This implies that in the induced instance, each doctor could have contracts with multiple hospitals at the same region. We can still convert an instance with these complicated diversity constraints into a corresponding instance with regional quotas, since the mapping relationship between two models does not change. Further development on matching with regional quotas will shed light on the problem of school choice with diversity constraints. Finally, it will be interesting to explore similar connections with other matching models (see e.g., (Aziz et al., 2018; Jones and Teytelboym, 2016)).

Acknowledgments

Aziz is supported by the UNSW Scientia Fellowship. Gaspers is the recipient of an Australian Research Council (ARC) Future Fellowship (FT140100048) and he also acknowledges support under the ARC’s Discovery Projects funding scheme (DP150101134). Sun receives support through Australian Government Research Training Program Scholarship. Walsh is funded by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme via grant AMPLIFY 670077.

References

  • Abdulkadiroğlu [2005] A. Abdulkadiroğlu. College admissions with affirmative action. International Journal of Game Theory, 33(4):535–549, 2005.
  • Abdulkadiroğlu and Sönmez [2003] A. Abdulkadiroğlu and T. Sönmez. School choice: A mechanism design approach. American Economic Review, 93(3):729–747, 2003.
  • Aziz et al. [2018] H. Aziz, J. Chen, S. Gaspers, and Z. Sun. Stability and pareto optimality in refugee allocation matchings. In Proceedings of the 17th International Conference on Autonomous Agents and Multiagent Systems (AAMAS), pages 964–972, 2018. URL http://dl.acm.org/citation.cfm?id=3237842.
  • Benabbou et al. [2018] N. Benabbou, M. Chakraborty, X. Ho, J. Sliwinski, and Y. Zick. Diversity constraints in public housing allocation. In Proceedings of the 17th International Conference on Autonomous Agents and Multiagent Systems (AAMAS), pages 973–981, 2018.
  • Berman et al. [2003] P. Berman, M. Karpinski, and A. D. Scott. Approximation hardness of short symmetric instances of MAX-3SAT. Electronic Colloquium on Computational Complexity (ECCC), 049, 2003.
  • Biró et al. [2010] P. Biró, T. Fleiner, R. W. Irving, and D. F. Manlove. The college admissions problem with lower and common quotas. Theoretical Computer Science, 411(34):3136 – 3153, 2010. ISSN 0304-3975.
  • Ehlers et al. [2014] L. Ehlers, I. E. Hafalir, M. B. Yenmez, and M. A. Yildirim. School choice with controlled choice constraints: Hard bounds versus soft bounds. Journal of Economic Theory, 153:648–683, 2014.
  • Fragiadakis et al. [2016] D. Fragiadakis, A. Iwasaki, P. Troyan, S. Ueda, and M. Yokoo. Strategyproof matching with minimum quotas. ACM Transactions on Economics and Computation, 4(1):1–40, 2016.
  • Gonzalez [1985] T. F. Gonzalez. Clustering to minimize the maximum intercluster distance. Theoretical Computer Science, 38:293 – 306, 1985. ISSN 0304-3975.
  • Goto et al. [2014] M. Goto, N. Hashimoto, A. Iwasaki, Y. Kawasaki, S. Ueda, Y. Yasuda, and M. Yokoo. Strategy-proof matching with regional minimum quotas. In Proceedings of the 13th International Conference on Autonomous Agents and Multiagent Systems (AAMAS), pages 1225–1232, 2014.
  • Goto et al. [2015] M. Goto, R. Kurata, N. Hamada, A. Iwasaki, and M. Yokoo. Improving fairness in nonwasteful matching with hierarchical regional minimum quotas. In Proceedings of the 14th International Conference on Autonomous Agents and Multiagent Systems (AAMAS), pages 1887–1888, 2015.
  • Goto et al. [2016] M. Goto, A. Iwasaki, Y. Kawasaki, R. Kurata, Y. Yasuda, and M. Yokoo. Strategyproof matching with regional minimum and maximum quotas. Artificial intelligence, 235:40–57, 2016.
  • Hafalir et al. [2013] I. E. Hafalir, M. B. Yenmez, and M.A. Yildirim. Effective affirmative action in school choice. Theoretical Economics, 8(2):325–363, 2013.
  • Hatfield and Kojima [2008] J. W. Hatfield and F. Kojima. Matching with contracts: Comment. American Economic Review, 98(3):1189–94, 2008.
  • Hatfield and Milgrom [2005] J. W. Hatfield and P. R. Milgrom. Matching with contracts. American Economic Review, 95(4):913–935, 2005.
  • Huang [2010] C. Huang. Classified stable matching. In Proceedings of the 21th Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 1235–1253, 2010.
  • Jones and Teytelboym [2016] W. Jones and A. Teytelboym. Choices, preferences and priorities in a matching system for refugees. Forced Migration Review, 51:81–82, 2016.
  • Kamada and Kojima [2012] Y. Kamada and F. Kojima. Stability and strategy-proofness for matching with constraints: A problem in the japanese medical match and its solution. American Economic Review, 102(3):366–70, 2012.
  • Kamada and Kojima [2015] Y. Kamada and F. Kojima. Efficient matching under distributional constraints: Theory and applications. The American Economic Review, 105(1):67–99, 2015.
  • Kamada and Kojima [2017a] Y. Kamada and F. Kojima. Recent developments in matching with constraints. The American Economic Review, 107(5):200–204, 2017a.
  • Kamada and Kojima [2017b] Y. Kamada and F. Kojima. Stability concepts in matching under distributional constraints. Journal of Economic Theory, 168:107–142, 2017b.
  • Kamada and Kojima [2018] Y. Kamada and F. Kojima. Stability and strategy-proofness for matching with constraints: A necessary and sufficient condition. Theoretical Economics, 13(2):761–793, 2018.
  • Kojima [2012] F. Kojima. School choice: Impossibilities for affirmative action. Games and Economic Behavior, 75(2):685–693, 2012.
  • Kominers and Sönmez [2016] S. D. Kominers and T. Sönmez. Matching with slot-specific priorities: Theory. Theoretical Economics, 11(2):683–710, 2016.
  • Kurata et al. [2017] R. Kurata, N. Hamada, A. Iwasaki, and M. Yokoo. Controlled school choice with soft bounds and overlapping types. Journal of Artificial Intelligence Research, 58:153–184, 2017.
  • McDermid and Manlove [2010] E. J. McDermid and D. F. Manlove. Keeping partners together: algorithmic results for the hospitals/residents problem with couples. Journal of Combinatorial Optimization, 19(3):279–303, 2010.
  • Ronn [1990] E. Ronn. NP-complete stable matching problems. Journal of Algorithms, 11(2):285–304, 1990.
  • Roth [1985] A. E. Roth. The college admissions problem is not equivalent to the marriage problem. Journal of Economic Theory, 36(2):277–288, 1985.