Equity in the Police Districting Problem: Balancing Territorial and Racial Fairness in Patrolling Operations

The Police Districting Problem concerns the definition of patrol districts that distribute police resources in a territory in such a way that high-risk areas receive more patrolling time than low-risk areas, according to a principle of territorial fairness. This results in patrolling configurations that are efficient and effective at controlling crime but that, at the same time, might exacerbate racial disparity in police stops and arrests. In this paper, an Equitable Police Districting Problem that combines crime-reduction effectiveness with racial fairness is proposed. The capability of this model in designing patrolling configurations that find a balance between territorial and racial fairness is assessed. Also, the trade-off between these two criteria is analyzed. The Equitable Police Districting Problem is defined as a mixed-integer program. The objective function is formulated using Compromise Programming and Goal Programming. The model is validated on a real-world case study on the Central District of Madrid, Spain, and its solutions are compared to standard patrolling configurations currently used by the police. A trade-off between racial fairness and crime control is detected. However, the experiments show that including the proposed racial criterion in the optimization of patrol districts greatly improves racial fairness with limited detriment to the policing effectiveness. Also, the model produces solutions that dominate the patrolling configurations currently in use by the police. The results show that the model successfully provides a quantitative evaluation of the trade-off between the criteria and is capable of defining patrolling configurations that are efficient in terms of both racial and territorial fairness.


Introduction
Predictive policing (Perry 2013) consists in the use of analytical quantitative techniques to identify potential criminal activity. Its most popular application is predictive hot spot policing, that is, the allocation of police resources to geographical points which are experiencing an increasing criminal trend, according to predictive models. Optimization models can be used in concert with predictive models to support decision-makers in distributing police resources based on risk maps. In particular, the Police Districting Problem (PDP) (Liberatore et al. 2020) defines efficient and effective patrol sectors according to the distribution of the crime risk in the territory considered. Specifically, PDP produces configurations where the total risk is evenly distributed between the sectors, thus, ensuring that each r-district 1 receives an amount of patrolling time that is proportional to its risk-level.
Targeting high crime areas has proven to be successful at reducing crime (Braga et al. 2014). Another consequence of having areas with a higher police presence is that more individuals will be stopped and, eventually, arrested in those areas. It is a well-known fact that high crime locations tend to be correlated with minority and disadvantaged populations (Sampson and Wilson 1995). Therefore, allocating more police presence at hot spots of crime or high crime risk areas has the likely consequence of contributing to racial disparity in police stops and arrests (Andrejevic and Gates 2014;Fagan 2017;Rinehart Kochel 2011). Also, the heightened exposition of minorities to police increases the feeling of pressure in these communities, subsequently producing negative repercussions for both the individuals and society as a whole. On the other hand, reducing police presence in the highest crime areas is not viable, as depriving minority communities of police protection would be discriminatory and could lead to further victimization (Mohler et al. 2015).
In this paper, a model that balances territorial and racial fairness in patrolling operations is proposed, building upon previous research from the authors (Camacho-Collados and Liberatore 2015; Liberatore and Camacho-Collados 2016). The model considers a territory divided into r-districts and a population segmented into groups (e.g., racial groups). Ideally, an efficient patrolling configuration assigns each r-district an amount of patrolling time that is proportional to its crime risk. At the same time, racial equity is achieved only if each population group has a police contact that is proportional to its size. Therefore, two criteria are maximized: the minimum proportion of service received across all r-districts and the minimum proportion of police exposure across all the population groups. Maximizing the minimum proportions increases the patrolling time received by each r-district and population group while reducing the variability as much as possible, thus ensuring efficiency and fairness. Given the distribution of the crime risk and the population on the territory, the model assigns r-districts to patrol sectors and optimizes both criteria according to a weight coefficient which expresses the decision-maker preference for one criterion over the other. By changing its value it is possible to generate the set of efficient solutions which comprise the Pareto frontier. 2 The proposed model presents some assumptions and limitations. First, policing involves a large number of different operations; this paper focuses on proactive patrolling operations in a territory. Second, the model does not explicitly incorporate the crimes' severity. On the other hand, the model makes use of the crime risk, which is a versatile measure that can be used to model the stakeholders' preferences as well as the priorities of the police department. For this reason, guidelines on how the crime risk is set are intentionally not provided. However, it could be easily defined to summarise the total level of crime in an area or to implement hot-spot policing. Finally, the proposed approach does not consider game-theoretical elements (Fu and Wolpin 2018;Galiani et al. 2018;Maheshri and Mastrobuoni 2019). As the aim of this paper is to explore the impact of including racial fairness in patrolling decisions, this is left as future research.
The contributions of this article are many-fold. A review of the latest papers on PDP is presented. This review expands and complements the literature review on the subject by Liberatore et al. (2020). The main contribution of this paper lies in the optimization model. To the best of the author's knowledge, this is the first PDP that explicitly considers racial fairness. The model is applied to a real-world case study which allows drawing insights on the usefulness and applicability of the model. The results show that in the dataset considered, a small loss in racial fairness leads to near-optimal solutions in terms of service level. More in general, the analysis of the trade-off between racial fairness and service level is necessary to inform decision-makers and to help them identify solutions that are socially equitable and highly performant.
The rest of the paper is organized as follows. In the next section (Sect. 2) the relevant literature is reviewed. The model is introduced in detail in Sect. 3. Next, in Sect. 4, the case study is introduced and the results of the model are discussed and analyzed. The paper concludes in Sect. 5 with a summary of the findings and a discussion of their implications.

Literature Review
Territory districting (or design) problems are a subfield of discrete optimization related to partitioning decisions (Ríos-Mercado 2020). In its most generic definition, a set of atomic geographic units must be divided into districts or territories, according to specific requirements which depend on the application context. This family of problems has been applied to a wide number of fields, including politics (Kim and Kim 2020;Ricca and Scozzari 2020), health-care (Yanık and Bozkaya 2020; Enayati et al. 2020), sales (Moya-García and Salazar-Aguilar 2020) and, of course, policing.
Since the seminal papers by Mitchell (1972) and Bodily (1978), several police districting models have been presented in the literature. The interested reader is referred to 1 3 the extensive literature review and annotated bibliography on the subject by Liberatore et al. (2020), which is complemented in the rest of this section by presenting alternative approaches to the problem of police patrolling. Leigh et al. (2019) combine predictive policing and optimization models to target highcrime areas. Hot spots are identified using kernel density estimation, while police agents' positions are determined by a maximum coverage location problems. Another line of research concerns the definition of patrol routes for police agents on a street map. For a review on the subject, the interested reader is referred to Dewinter et al. (2020). Chen et al. (2017) propose a heuristic online Bayesian ant-colony optimization algorithm which generates routes that are efficient, flexible, unpredictable, scalable, and robust. In a subsequent paper (Chen et al. 2019), the authors present a street network police districting problem where the territorial atomic units are street segments rather than r-districts. The objective is to generate patrol sectors that are balanced in terms of workload, which is defined as a combination of multiple factors, i.e., risk, area, and diameter, following previous contributions in the literature (Liberatore et al. 2020). The problem is formulated as a Mixed Integer Programming (MIP) problem and a tabu search meta-heuristic algorithm (Glover 1989(Glover , 1990 is proposed for its solution. In their paper, Chen et al. recognize that the street map approach, although it leads to a more straightforward definition of the police patrol routes, is an obstacle to the incorporation of census data in the model, as street segments are incompatible with census districts. A possible solution to this issue is given by Kim (2018) that proposes and tests methods for re-distribute census data to street segments.
To the best of the authors' knowledge, the only optimization model in the context of police patrolling that considers issues of racial inequalities is proposed by Wheeler (2019a). The author tackles the problem of police resources allocation in a fixed districting configuration. The problem is formulated as a Linear Programming model which includes a constraint that imposes an upper bound to the expected police contact of minorities. In another paper, Wheeler (2019b) formulates a police districting model based on the p-median problem. Additional constraints have been added by the author to ensure that the clusters are connected. However, these constraints are incorrect as they still allow a sector to have nonconnected subclusters, i.e., the model allows a sector to be formed by geographically disconnected areas. On the other hand, as the focus of the model is on answering calls rather than on proactive patrolling, this does not hinder its usefulness and applicability.

Model
The problem considers a territory partitioned into non-overlapping areas, such as report or census districts. The population in the territory is divided into population groups, that could be defined according to a characteristic of interest, e.g., race, ethnicity or origin. The crime distribution and the population groups distribution are known and, therefore, each area is characterized by a crime level and by the size of each population group. The main decision involves assigning the areas to the patrol districts (multiple assignments are allowed). Then, the time capacity available at each patrol district must be distributed between the tasks of travelling between the areas and patrolling the areas in the district. The patrolling time spent in the areas is used to compute two scores: the service score is the proportion of demand received by an area as patrolling time, while the contact score is the ratio between the police contact time of a population group and its goal. Territorial and racial fairness are achieved by maximizing the minimum service score and the minimum contact score, respectively. In fact, Maximizing the minimum score produces solutions with a high average value and low variance.
The rest of the section provides technical details of the methodology. Readers only interested in the practical outcomes of this work may omit this section. A summary of the notation used in the model is provided in Table 1. The parameters of the model are presented in the following. Let G = (N, E) be an undirected graph where N is the set of areas and E the set of edges connecting them, and let K be the number of districts (indexed by k) and P the set of population groups (indexed by p). Each district has a time capacity, cap, and each population group has a normalized police contact goal (e.g., proportional to the size of the group), goal p , such that ∑ p∈P goal p = K ⋅ cap . Each area i ∈ N is characterized by: -a population distribution d ip representing the percentage of the population of area i that belongs to group p, being ∑ p∈P d ip = 1, ∀i ∈ N. -a normalized service demand (e.g., proportional to the expected crime risk in the area), dem i , such that ∑ i∈N dem i = K ⋅ cap. Finally, the edges have an associated distance, dist ij , which represents the time required to go from area i to j.
In the problem considered, the main decisions are represented in the model by the following variables: Each district k defines a subgraph G k ⊂ G generated by the areas in the district and their incident edges. This is illustrated through an example in Fig. 1a. In the figure, the circled portion of the graph corresponds to a district; the areas that comprise it and the edges that connect them (represented in red) identify a subset, G k , of the original graph G. It is important to notice that exactly defining the route that a patrol would use is out of scope and, also, it is not recommendable as it would make the patrols predictable. Therefore, inside of each district k, the patrol route is approximated by a Spanning Tree (ST) of the graph G k . ST is a subgraph which connects all of the vertices of the original graph at minimum cost. That is, ST identifies the set of edges having minimum total distance that connects all the areas in the district (see Fig. 1b for an example). Implementing ST has the additional advantage of ensuring that the districts are connected, i.e., it is possible to connect every pair of areas in the district using only edges that are part of the district. The interested reader can find the complete formulation of ST in Appendix A.
Let ST k ≥ 0 be the length of ST corresponding to district k. It is assumed that a complete patrol route would traverse ST twice as to complete a round trip (i.e., a loop that visits each area at least once and ends where the patrol started) each edge needs to be travelled The number next to each edge represent its distance and the edges comprising the ST are drawn in red. The total distance (i.e., the sum of the distances of the edges in ST) is minimum. Also, it is possible to connect any pair of nodes in G k using only edges in ST (Color figure online) two times. Therefore, the capacity of the district is reduced by twice the length of the corresponding ST (i.e., ST k ) and the remaining time can be allocated to its areas. In the following, the constraints and objectives of the model are presented according to their scope.

Districting Constraints
These constraints concern the definition of feasible districts and the allocation of the service time to the areas. The first two sets of constraints (1 and 2) state that no district can be empty and all areas must be assigned to at least one district, respectively. The allocation constraints (3) state that an area can receive patrol time only from districts to which it is assigned. Constraints (4) specify that, for each district, the sum of the patrol time allocated to the areas cannot exceed the district capacity reduced by twice the length of the corresponding ST. The last two sets of constraints present the existence conditions for variables assign ik and time ik . Below, the criteria that allow to explore the trade-off between territorial and racial fairness are presented.

Optimization Criteria
Let us consider an area i ∈ N and its demand, dem i . The service level received by an area can be computed as the percentage of demand satisfied, i.e., the ratio between the patrolling time assigned to the area and its demand: In terms of patrolling efficiency, it is desirable that all the areas receive the highest level of service with little variability among them, to patrol the areas as homogeneously as possible according to their demands and the available capacity. Thus, the territorial fairness criterion is formulated as: which is linearized by introducing an additional continuous variable, W s , representing the lowest service level across all the areas: On the other hand, the police contact of a population group p ∈ P can be computed as the ratio between the patrolling time received by the population type and the contact goal: The equation assumes that all the citizens in an area receive an identical amount of patrolling time, i.e., police contact is assumed to follow a uniform probability distribution. Therefore, the expected police contact of a population group is the sum across the areas of the proportion of time received by the area relative to the size of the group. Police contact should be as homogeneous as possible among the groups to ensure racial fairness. At the same time, it is still important to provide each group with the highest possible level of police contact, to avoid social abandonment and victimization from crime. Therefore, the racial fairness criterion can be formulated as: which is linearized by introducing an additional continuous variable, W c , representing the lowest exposition level across all the groups:

Payoff Matrix and Compromise Programming Model
To verify the trade-off between the criteria, the payoff matrix is computed: The matrix presents for each criterion (i.e., W s and W c ) the ideal and the anti-ideal values, identified by a '+' and a '−' superscript, respectively. The ideal and the anti-ideal values are the best and worst possible values obtainable for each criterion, respectively. The payoff matrix is computed following Algorithm 1.
Each step of the algorithm computes one of the values in the payoff matrix (15). First, the ideal territorial fairness ( W + s ) is calculated by maximizing the minimum service level across the districts. Then, the anti-ideal racial fairness ( W − c ) is identified by maximizing the minimum police contact while imposing that the territorial fairness must be at least as high as its ideal. This completes the first row of the payoff matrix. Next, the ideal racial fairness ( W + c ) is determined by maximizing the minimum police contact and, finally, the anti-ideal territorial fairness ( W − s ) is assessed by maximizing the minimum service level across the districts while imposing that the racial fairness must be at least as high as its ideal. After this step, the full payoff matrix is achieved.
To obtain a solution to this multi-criteria problem, Compromise Programming (CP) (Ringuest 1992) is used. CP is a decision making technique that finds the best solution by minimizing the linear combination of the normalized distance of the criteria values to the ideal point. In this model the distance value is calculated as: where ∈ [0, 1] is a parameter that specifies the preference of the decision-maker on the territorial fairness criterion over the racial fairness criterion. The complete model, called the Equitable Police Districting Problem (EquPDP), is given below. [EquPDP] min Z Compromise programming distance: (16) Optimization criteria: (7), (10), (11), (14) Partitioning constraints: (1)-(6) Spanning tree constraints: (23)- (35) Let K be the number of districts to be defined, |N| the number of areas, |P| the number of population groups, and |E| the number of edges. Then, the model includes: (2K|N| + K|E|) binary variables, (3 + K + |P| + |N| + 2K|N| + K|E|) real variables, and (1 + 4K + 2|P| + 3|N| + 9K|N| + 5K|E|) constraints (including the existence conditions for the variables). The solutions provided by EquPDP are efficient and by varying the value of it is possible to obtain the Pareto Frontier.

Case Study
EquPDP is applied to a real case study on the Central District of Madrid, Spain. The Central District of Madrid is approximately 5.23 km 2 (2.02 sq mi) in size, it has a population of 149,718 people and a population density of 28,587/km 2 (74,040/sq mi). Administratively, the district is divided into six wards (also called barrios): Cortes, Embajadores, Justicia, Universidad, Palacio, and Sol. The crime risk definition given below has been determined according to the objectives of the Spanish National Police Corps (SNPC) and the Central District of Madrid PD in particular. According to the OECD Better Life Index, Spain actually has a lower incidence of crime than in many other OECD countries (OECD 2020), e.g., Spain has a homicide rate of 0.6 murders per 100,000 inhabitants, compared to the OECD average of 3.6 per 100,000. In particular, the capital city of Madrid and its surrounding region is the safest of the large population centers; in fact, the region of Madrid contains 13.75% of Spanish residents but it is responsible for just 10.41% of all crime. Theft is the most frequent type of crime committed in Spain and one of the main priorities for the SNPC is its reduction. For the above reasons, in this case study the crime risk of an area is the number of thefts reported.

Dataset
The topological data has been obtained from official sources (NOMECALLES 2020). Census districts have been chosen as the atomic territorial units. Each census district is represented by its centroid, which translates into an area in the graph on which the model is formulated. Two areas are connected by an edge if the corresponding census districts share part of their perimeters. The territory and the resulting graph, comprised of 111 areas and 308 edges, are shown in Fig. 2. For each edge (i, j) ∈ E , let length ij be the great-circle distance 3 between the areas in meters. Then, the edges' distance is calculated as follows: which corresponds to the time in minutes that it takes to walk the edge assuming a walking speed of 100 m/min (equivalent to 6 km/h).
The crime data has been provided by the SNPC. Specifically, the dataset includes all the reported thefts occurred in each area during the following shifts: The heat maps for each shift are represented in Fig. 3. Let risk i be the number of thefts in area i. Then, the areas' demand is set to where cap is the duration of the shift in minutes, i.e., 600 for SATT3 and 420 for SUNT1 and MONT2. Concerning the number of districts, the values considered are K = {2, 6} . In fact, the standard patrolling configurations for the Central District of Madrid adopted by SNPC (represented in Fig. 4) partition the territory into either two districts (i.e., north/ south of the main artery, the Gran Via) or six (i.e., according to the wards), which are each assigned to multiple officers. The population data has been obtained from the 2011 Spanish Census (INE 2011) which provides data on a census district level and segments the population according to the following geographical regions of birth: (a) Spain; (b) other EU country; (c) European non-EU country; (d) Africa; (e) Caribbean, Central and South America; (f) North America; (g) Asia; (h) Oceania. Following this categorization, the population groups set P has eight elements, one for each region of birth. The distribution of each population group is illustrated in Fig. 5.
Let pop ip be the number of people living in i that were born in region p. Then, the population distribution over the territory and the patrolling goals for the population groups are calculated as follows:  (Gurobi Optimization 2020). All the experiments have been run on a Dell Precision 5540, equipped with Intel Ⓡ Core TM i9-9880H CPU @ 2.30GHz × 16 and 16GB RAM. The standard configuration of Gurobi has been used, which applies multithreading. A CPU time limit of 3600s has been set on all the optimization processes. Overall, six instances are considered, obtained by combining three shifts (i.e., SATT3, SUNT1, and MONT2) with two values of K (i.e., K = 2, 6 ). For each instance the payoff matrix is computed and the Pareto frontier is sampled at the following values of the parameter = {0.0, 0.25, 0.5, 0.75, 1.0} . This entails solving nine optimization models per instance, i.e., four to compute the payoff matrix and one for each value of . Due to the size of the graph, none of the optimization processes could complete within the time limit. Therefore, the results of EquPDP presented below are sub-optimal. On the other hand, the evaluation of the standard configurations by SNPC involves solving EquPDP with fixed assignment variables (i.e., assign ik ) and takes less than one second to complete to optimality. Table 2 illustrates the payoff matrices for the instances considered for K = 2 and K = 6 . It is possible to identify a trade-off between the criteria. More in detail, when K = 6 it can be appreciated that the territorial fairness criterion presents a larger gap between the ideal and the anti-ideal value than the racial fairness criterion. Also, when optimizing with respect to one criterion, the model achieves a high level of efficiency, higher than 0.85 for both criteria. On the other hand, when K = 2 , the criteria present similar values and gaps between the ideal and the anti-ideal points.

Racial Fairness Criterion Analysis
To verify the ability of the racial fairness criterion to generate solutions with homogeneous police contact levels across all the populations groups, the solutions obtained by the steps 2 and 4 of Algorithm 1 are considered and compared. In particular, these steps are equivalent to solving EquPDP with the following objective functions: and where "lex" stands for "lexicographic". 4 Equation (21) prioritizes W s over W c , yielding solutions that are more efficient yet less equitable in terms of population contact with the police, while Equation (22) does the opposite. Table 3 illustrates the average value and the standard deviations for the population groups' contact levels with the police on all the problem instances considered. It is evident  (21) and (22) The columns present the following information: shift, number of districts (K), average and standard deviation of contact level (variables contact p ) for objective function (21)  that prioritizing the racial fairness criteria results in solutions that have almost no variability and high contact level on average. On the other hand, optimizing with respect to territorial fairness yields solutions that are highly unequal in terms of police contact, with standard deviations ranging from 0.12 to 0.24, approx. Therefore, it can be concluded that the racial fairness criterion indeed produces the desired results. Figure 6 illustrates the sampled Pareto frontiers for the considered instances and the corresponding SNPC patrolling configurations (see Fig. 4). The figure is comprised of six plots, one per problem instance (i.e., SATT3, SUNT1, and MONT2) and value of K = 2, 6 . Each plot is a scatterplot of the solution values found by EquPDP and of the SNPC solution at different values of the preference coefficient = 0, 0.25, 0.5, 0.75, 1 . The axis correspond to the fairness criteria values (i.e., the territorial fairness, W s , and the racial fairness, W c , criteria). In each scatterplot, the solution points are drawn differently depending of whether they are generated by EquPDP (dot) or are the SNPC solution (triangle). Also, the points are connected by lines (continuous and dotted, respectively) to show the corresponding frontier. It can be seen that, in every plot, the frontier corresponding to EquPDP is always on top of the SNPC's frontier; this means that the solutions defined by EquPDP always dominate the SNPC's standard configuration. This means that EquPDP's solutions are always better than the SNPC's standard configuration and, therefore, the proposed model clearly improves on the quality of the patrolling operations currently in use. On a different note, when considering only the solutions found by EquPDP it can be seen that the frontiers present a very clear elbow. More importantly, the curves show that non-extreme solutions (i.e., the solutions obtained when 0.25 ≤ ≤ 0.75 ) provide a good compromise between the two criteria. Given the proximity of the non-extreme solution points in the plots, the value of the parameter does not seem to have a major impact on the criteria values, as long as it is 0.25 ≤ ≤ 0.75.

Pareto Frontiers and Comparison with SNPC Configurations
It is important to notice that, as mentioned above, the solutions obtained in this experiments by the EquPDP are sub-optimal, which means that their value (variable Z, Eq. 16) is an upper bound to the optimum ( Z ⋆ ). In a plot where the axes are W s and W c (such as in Fig. 6), the point corresponding to a sub-optimal solution would be closer to the origin than the point corresponding to the optimal one. Therefore, the Pareto frontier identified by optimal solutions would be farther away from the origin than the current one and, as a consequence, the dominance of the solutions found by EquPDP over the solution of the SNPC would be even greater. This is a further proof of the superiority of EquPDP over the standard configurations currently in use. Table 4 compares the criteria values for all instances at = 0.0 and = 0.25 . The results displayed in the table confirm that, by sacrificing approximately 0.01 in racial fairness, it is possible to improve the service level by more than 0.50 for K = 2 , and 0.75 for K = 6 . This outcome is specific to the dataset considered; in fact, the gains in terms of racial fairness depend on the level of segregation among population groups. In particular, it is expected for these gains to be lower for higher levels of segregation. Regardless of that, what these results manifest is that the trade-off between racial and service equity needs to be assessed to make informed decisions, as high level of racial fairness could still lead to nearly-optimal solutions in terms of service level. Figures 7 and 8 shows the patrolling configurations corresponding to the solutions produced by EquPDP for K = 2, 6 and = 0.25 . Despite EquPDP allows to assign an area to multiple district (see Eq. 2) only two configurations have overlapping patrol districts. The solutions for K = 2 (Fig. 7) are quite similar, with small variations to adapt to the idiosyncrasies of each shift. On the other hand, for K = 6 ( Fig. 8) very different configurations are obtained: -In shift SATT3, most crimes are concentrated in the middle section of the district and in the north-eastern corner (see Fig. 4a), which represent approximately 50% of the ter- ritory and correspond to busy nightlife hot spots. In the solution identified by EquPDP this area is covered by as many as five different sectors. -Concerning shift SUNT1, most crimes are located in the center and in the south. The census district with the highest crime rate corresponds to a flea market that takes place every Sunday morning. In the solution, the yellow sector covers the center while the dark blue sector focuses on the area of the flea market. Three sectors (green, light blue, and red) provide support to cover the south, while the pink sector is dedicated exclusively to the north. -During shift MONT2, crime is concentrated almost exclusively in the middle section of the district, corresponding to about one third of the territory. EquPDP produces a solution where this area is covered by five different sectors, while the red sector focuses on the north, which is almost entirely free of crime.

Conclusions
In this paper, a police districting problem that balances between territorial and racial fairness is proposed. The problem is formulated as a MIP model that defines patrolling configurations that are operationally efficient and, at the same time, equitable in terms of the contact between the police force and the different population groups. The model is tested on a real-world case study on the Central District of Madrid, Spain, which allows drawing some insights regarding the model that are presented in the following.

Results Summary
The following highlights the major findings obtained from the experiments.
-There is a clear trade-off between territorial and racial fairness.
-For K = 6 , the gap between the ideal and the anti-ideal values is larger for the territorial fairness criterion. Also, solutions that prioritize racial fairness result in a service level close to zero. On the other hand, prioritizing territorial fairness produces solutions with a racial fairness close to 0.50. For K = 2 these differences are not observed. -The racial fairness criterion effectively produces solutions with high average police contact level and extremely low variability (i.e., standard deviation close to zero). -The solutions identified by EquPDP dominate the standard patrolling configurations adopted by SNPC. -The Pareto frontier presents a well-defined elbow. Furthermore, non-extreme solutions are clustered in proximity of this elbow and are close to the ideal point. This implies that the value of parameter does not have a major impact on the criteria values and that, as long as it is strictly larger than zero and smaller than one, the solution found by EquPDP will be "good." -In terms of compromise between criteria, a very small decrease in racial fairness (0.01 approx.) results in a great improvement in territorial fairness (higher than 0.70). -Thus, it is extremely beneficial to optimize considering both criteria at the same time, as it results in solutions that are both very effective at controlling crime and equitable.
These insights are specific to the case study considered; however, they illustrate the applicability and usefulness of the methodology proposed. In particular, the results show that the Pareto frontier has a very sharp elbow. This indicates that most compromise solutions are clustered in the proximity of the ideal point and, therefore, choosing any among them does not have a major impact on performances. Optimizing exclusively for territorial fairness results in patrolling configurations that are highly imbalanced in terms of racial fairness and, vice-versa, prioritizing racial fairness produces very inefficient solutions. On the other hand, by sacrificing approximately 0.01 in racial fairness, it is possible to significantly improve the service level. This demonstrates the importance for the decision-makers of studying the trade-off between racial and territorial fairness.

Discussion
The experimental results verify that, as postulated by Kleinberg et al. (2018), the explicit incorporation of racial information into police decision making can lead to more fair and efficient outcomes. On the other hand, disregarding this information might lead to unfair results which have the potential to exacerbate racial disparity, due to the structural differences between population groups and the correlation between crime level and minorities presence.
The approach presented in this paper hinges on improving racial fairness at the expenses of effectiveness. Certainly, voluntarily depriving high-crime areas of much-needed security gives rise to an ethical dilemma. Although the methodology introduced allows quantifying the trade-off between effectiveness and racial fairness, it does not provide a solution to such predicament. What is usually recommended in these cases is that the decision-makers (e.g., law enforcement agencies and/or the public) should determine the degree of preference of racial fairness over effectiveness (Cohen 2017), that is represented in the model by the coefficient . However, it is the opinion of the authors that the ethical dilemma described above arises only when looking at the crime-reduction benefits of police operations in the short and medium term. Recent events (Wikipedia contributors 2020) have clearly shown the consequences that racial inequality in police interventions can have on society in the long term. These consequences include protest marches, social unrest and riots, which might lead to a surge in violence and crime. Factoring these long-term effects provides sufficient ground to estimate with objectivity and precision the most beneficial level of preference between the two criteria. This would require longitudinal and correlation study, and is left as future research.
The methodology presented in this paper does not solve the problem of racial inequality in police interventions. However, it can be used in conjunction with other strategies intended to mitigate long term harm to particular communities, such as problem-oriented policing (Scott and Clarke 2020), proactive policing (National Academies of Sciences 2018), and increasing the perceived legitimacy of police interactions with the public (Braga et al. 2019).
The authors hope that this work will be a useful source of ideas for future research on PDP and will contribute further in the development and solution of more complex and more realistic models in the context of policing, public security and safety.

Appendix A Spanning Tree Formulation
In the following, the mathematical programming model for the ST is provided. ST is formulated using a flow model for each district. One of the areas in the district is selected as the "root" that emits supply while the other areas must receive one unit of supply each. A flow of supply travels from area to area on the edges; all the non-root areas must satisfy the condition that the amount of entering flow must be equal to the amount of flow that leaves it plus one extra unit of supply that each area must receive. This ensures that the corresponding district is connected, as all the areas in ST must be reachable from the root of ST using only edges in the district (Shirabe 2005; Duque et al. 2011). A summary of the notation is given in Table 5. The variables required to model ST are: -ST k ≥ 0 , length of ST corresponding to district k.
root ik = 1 if area i is the root of ST of district k 0 otherwise cross ijk = 1 if edge (i, j) is in ST of district k 0 otherwise 1 3 -supply ik ≥ 0 , supply of area i in the flow model of district k. This variable takes value greater than zero only for the root area of ST, i.e., root ik = 1 . In this case, its value is the number of areas in k. -flow ijk ≥ 0 , flow on the edge (i, j) in district k. This variable can take value greater than zero only for the edges in ST, i.e., cross ijk = 1.