key: cord-0839470-dlam5xyo authors: Santini, Alberto title: Optimising the assignment of swabs and reagent for PCR testing during a viral epidemic date: 2020-09-22 journal: Omega DOI: 10.1016/j.omega.2020.102341 sha: 2789c618ad4ac0e086079a1c978ef4e789f6fd79 doc_id: 839470 cord_uid: dlam5xyo Early large-scale swab testing is a fundamental tool for health authorities to assess the prevalence of a virus and enact appropriate mitigation measures during an epidemic. The COVID-19 pandemic has shown that the availability of chemical reagent required to carry out the tests is often a bottleneck in increasing a country’s testing capacity. Further, demand is unevenly spread between more affected regions (which require more tests they can perform) and less affected ones (which have spare capacity). These issues hint at the opportunity of increasing test capacity via the optimal allocation of swabs and reagent to laboratories. We prove that this is the case, proposing an Integer Programming formulation to maximise the number of tests a country can perform and validating our approach on both real-life data from Italy and synthetic instances. Our results show that increased inter-regional collaboration and a steadier supply of reagent (i.e., coming from local production sites rather than international shipments) can dramatically increase testing capacity. Accordingly, we propose short-term and long-term recommendations for policy makers and health authorities. Viruses are pathogens that replicate after penetrating the living cells of other organisms. Outside cells, a virus is nothing more than genetic material (DNA or RNA) surrounded by protective layers of proteins and lipids. Once inside a host cell, the virus uses the cell's structures to replicate its genetic material and assembly its protective layers, thus creating a copy of itself. At the end, it kills the cell to release both the original and the copy. Viruses cause a host of human diseases, ranging from the common cold to AIDS. Of particular interest for this work is the recent Coronavirus disease 2019 (COVID- 19) pandemic. Under particular circumstances, a virus (and the disease it causes) can spread to a large part of a population in a short time, leading to an epidemic. When the epidemic spreads across national borders and infects people worldwide, it is termed a pandemic. SARS-CoV-2, the virus spreading COVID-19, for example, possessed the right characteristics to turn its associated disease into a pandemic: people carrying the virus do not show symptoms for an average of ve days, during which they can spread the virus to others [44] ; its basic reproduction number has been estimated between 1.4 and 5.7, i.e., each infected person in turns infects an average of up to 5.7 people [33] ; up to 44% of patients show no symptoms at all during the infection period, making their diagnosis di cult [29] . Massive testing of the population is a fundamental part of the response against epidemics and pandemics. When carriers of the virus remain asymptomatic while infecting others, as for COVID-19, testing is one of the main tools health authorities can deploy to contain the spread of the disease [37] . Viruses such as SARS-CoV-2 are RNA viruses, i.e., they contain a single strand of nucleotides rather than the "double helix" typical of DNA. For this brief introduction, we call the four nucleotides composing the RNA chain with their initials: A, C, G, U. The quickest available test for RNA viruses uses real-time reverse transcription polymerase chain reaction (rRT-PCR). A popular test consists in collecting a sample of nasal secretions on a swab [7] , looking for an RNA subsequence which is unique to the virus. To identify this subsequence, the tester uses inverse RNA: a sequence of nucleotides which pairs those they are looking for (remember that A pairs with U and C with G). If the sample contains viral RNA, it will pair with the chain injected by the tester. Once the nucleotides are paired, they e ectively transform the single-strand RNA into double-strand DNA. To detect a signi cant presence of the viral DNA, the genetic material goes through an ampli cation process, the polymerase chain reaction [6] . Both transforming RNA into DNA, and amplifying it, cannot happen without certain enzymes, commonly called reagents. While this simpli ed description is enough to introduce our problem, we refer the reader to the reviews of Bustin [4] and Freeman, Walker, and Vrana [16] for a more accurate description of the rRT-PCR technique. Figure 1 is an example of a protocol for rRT-PCR tests used by the United Kingdom's National Health Service. During the COVID-19 pandemic, laboratory capacities and reagent availability have become the bottleneck for increasing the number of tests in much of Europe, North America and Asia [13, 1, 43, 27] . Therefore, health authorities need to optimise the allocation of resources to their test laboratories, starting from the distribution of reagent and the assignments of testing tasks. With this work, we use tools from Operational Research to help decision makers maximise the number of tests they can conduct while limited by reagent availability and logistic constraints. Our main contributions are the following. • We introduce the problem of maximising a country's test capacity for rRT-PCR tests, during the outbreak of an epidemic. In such a situation, (i) reagent scarcity is often a bottleneck for testing, and (ii) some regions are foci of the disease, requiring more tests than others. These two conditions provide an opportunity to optimise the allocation of reagent and swabs, with the aim of maximising the number of swabs tested and reducing the wait time for results. • We provide an Integer Programming formulation for this problem, extend it to consider the real-world case in which the health-care system is decentralised at the regional level, and consider a hierarchical multi-objective version of the problem. Because commercial software such as Gurobi solves the formulation in a few minutes, this model is suitable for direct use by public health decision makers and to run what-if analyses. • We conduct two case studies to validate our approach. In the rst one, we use real-life data referring to the Italian situation during the COVID-19 pandemic. Our analysis con rms that fostering inter-regional collaboration can increase by up to +40% the number of swabs tested during a 13-day planning horizon. In the second one, we use synthetic data to determine the most important factors which impact testing capacity. We conclude that a strong industrial network able to provide a steady supply of reagent provides the greatest advantage, vs. procuring reagent via sporadic shipments, even when the total purchased quantity stays the same. Further, countries with limited lab infrastructure should prioritise acquiring more machines and training more personnel, while increasing inter-regional coordination helps those countries with a developed laboratory network. • We make available both our data and the code we used, under an open-source license [35] . In our problem, a set of laboratories L = {1, . . . , |L|} has to perform rRT-PCR tests on swabs during a time horizon T = {1, . . . , |T |} (each time unit corresponds to a day). We assume, without loss of generality, that for each swab a lab needs one unit of reagent. This allows to consider di erent testing routines, such as group testing, by rescaling the amount of reagent needed to test one swab. Each lab l ∈ L starts with a reserve ρ l 0 of reagent at the beginning of the time horizon. The labs can receive further units of reagent from factories R = {1, . . . , |R|}, limited by their production capacity. Each factory r ∈ R produces f r t units of reagent on day t ∈ T . Factories also store an initial amount ρ r 0 of reagent at the beginning of the planning horizon. A factory can also model other types of facilities; for example, a warehouse receiving new reagent once per week, or even a one-time shipment from a foreign country. Each laboratory l has a capacity Q l of swabs it can test during one day. This capacity applies even if the lab has a larger amount of reagent, due to limitations on available machinery and workforce. The e ective testing capacity of l, then, is the minimum between Q l and the amount of reagent available at l. A lab l ∈ L is tasked with testing m l t swabs on day t, according to a prede ned schedule to meet epidemiological needs. If a laboratory does not have enough reagent to test all the swabs, the decision-maker has three options: (i) moving some reagent from a factory to the lab, (ii) moving some swabs to another lab, or (iii) storing the swabs and schedule their testing for another day. Remark. We assume the planner has already chosen the laboratory to which he assigns the swabs. More generally, though, a planner might need that swabs be tested in a given region, without any constraint on which speci c lab performs the test, as long as they are geographically close to the point where the swabs were collected. In Section 5.1, we will extend the model to take into account this general case, making the number of swabs assigned to each lab a decision variable. We also refer the reader to Section 5.2, in which we introduce a fourth possibility for the decision maker: to move some reagent from a lab which has a surplus to one which is lacking. The amount of reagent and swabs that a planner can move between locations each day is bounded by q rg and q sw , respectively. Remark. We assume that quantities q rg and q sw are global capacity limits on the total amount of reagent and swab movements. Such an assumption rests on the observation that, during an epidemic, a central decision-maker manages logistic resources and can reallocate them from one area to another. If this were not the case, it would not be hard to enforce local limits, e.g., on the number of swabs that the planner can move out of a single laboratory. In a model extensions presented in Section 5.1, locations are partitioned in regions and we impose region-speci c capacity limits. Given this input data, a planner must determine, for each day of the planning horizon, (i) how many units of reagent to move from the factories to the laboratories, and (ii) how many swabs to move between laboratories, with the aim of maximising the number of tests carried out. Remark. While the main objective of the planner is to maximise the number of swabs tested, it is possible to account for secondary objectives. For example, because e ective testing must be both large-scale and quick, one might want to minimise the average time swabs spend waiting at lab facilities. In Section 5.3 we propose a multi-objective extension of the problem in which, among all solutions which maximise the number of tests carried out, we select the one with lowest average waiting time. Our problem, described in Section 2, resembles super cially problems in the area of production planning, material requirement planning and supply chain optimisation (for a complete overview on these topics see, e.g., Pochet and Wolsey [30] ). In our case, carrying out tests corresponds to manufacturing a product using two raw materials (swabs and reagent); for us, though, the nished product demand (tested swabs) always matches one of the raw materials' availability (swabs), while the other (reagent) acts as a bottleneck for production. As in production planning, we deal with the ow of goods over time, including shipments of raw material and inventory carry-over. We also allow the creation of backlogs, i.e., the storage of swabs which cannot be tested on the day of collection. Other characteristics of our problem, however, are not frequently considered in production planning optimisation: the geographic position of suppliers and production facilities (in our case, laboratories) and their regional subdivision; the possibility of reassigning orders between facilities, subject to logistics constraints; the uneven demand between regions. Most importantly, unlike problems in production planning, we disregard costs and focus exclusively on maximising production (in our case, number of tests). To this end, we do not consider transport, production, storage, raw material, order placement, personnel, nor opportunity costs, which instead have a prominent role when optimising commercial supply chains and production processes. In the standard case, in fact, the aim of the planner is pure cost minimisation [30] and satisfying the demand in full could even be sub-optimal, if the cost of producing one more unit is larger than the opportunity cost of missing a sale. Other approaches commonly used in production planning include multi-objective optimisation in which cost minimisation is tackled together with secondary objectives, such as demand satisfaction [26] , steadiness of production rates [40] or environmental performance [2] . Even in the healthcare sector, if considering times when authorities are not responding to an emergency, costs are taken into account. For example, Zhao et al. [48] consider a drugs production planning problem in which they minimise both costs (avoiding over-production of drugs) and the duration of the trial. Other works (e.g., [10, 14] ) have considered the production and supply optimisation problems in the clinical area, but have focused on production and distribution of drugs in "normal" times, rather than during an emergency response to a sudden epidemic event. Literature focusing speci cally on allocation of medical resources during emergencies mainly concentrates on the distribution of (life-saving) medical devices. As such, most models optimise the assignment of existing material to hospitals or individual patients, rather than the production and allocation of new material. Xiang and Zhuang [46] , for example, consider the problem of assigning medical resources to patients, taking into account that the conditions of patients left untreated will deteriorate over time. Cao and Huang [5] use a simulation approach in which they consider di erent strategies to allocate scarce medical resources (e.g., milder patients rst, critical patients rst, randomly) and complement an analysis of the utility of each strategy with ethical considerations on prioritising some patients over others. In the wake of the COVID-19 epidemics, the optimisation community has rushed to provide decision support tools for decision-makers involved in the ght against the virus. Similar to our work, researchers at MIT's Operations Research Center [9, 8] have shown the advantages of inter-regional collaboration in sharing automatic ventilators across U.S. states. Mehrotra et al. [23] consider the same problem and use a stochastic model to show how sharing ventilators across states, even when each state is averse to the risk of not having enough machines for its own patients, can reduce the probability of shortages. The problem of allocating ventilator during in uenza epidemics was already considered in the literature, e.g., by Huang et al. [19] , Meltzer et al. [24] , and Zaza et al. [47] . Lampariello and Sagratella [22] study a problem that, like ours, concerns diagnostic tests during the COVID-19 pandemics. The authors consider the problem of deciding how many swabs each region should collect to provide an accurate estimate of the prevalence of the disease among its population, while restricted by a nation-wide capacity limit on the number of tests performed. The problem is motivated by the practice employed in Italy during the early period of the epidemic, of double-checking at one single central laboratory (located in Rome) the results of each positive test conducted by local labs. The capacity of the central lab limits, in this case, the number of total expected positive results coming from all other labs in the country. Note how a decision maker working under the assumptions of [22] can use the model therein to generate regional demands m it , which they can then feed to our model. Finally, Seccia [36] considers the problem of devising balanced rosters for nurses in hospitals while, due to the virus outbreak, there is both a peak in demand for nurses' services and fewer available health professionals, because many are COVID-19-positive and con ned at home. We propose an Integer Programming model using the following sets of variables. The number of units of reagent moved from factory r ∈ R to lab l ∈ L on day t ∈ T is denoted by variable x r l t ∈ Z + . Variables l 1 l 2 t ∈ Z + represents the number of swabs going from laboratory l 1 ∈ L to l 2 ∈ L on day t ∈ T . The quantities of reagent stored, respectively, at factory r ∈ R and lab l ∈ L on day t ∈ T are denoted by variables ρ r t ∈ Z + and ρ l t ∈ Z + . Variables z l t ∈ Z + represent the number of swabs stored at lab l ∈ L during day t ∈ T . Laboratories store swabs if they cannot carry out the tests on day t because of a lack of reagent necessary for the rRT-PCR procedure. We assume without loss of generality that, at the beginning of the time horizon, labs do not have any stored swab (z l 0 = 0 for all l ∈ L). Maximising the number of tests corresponds to minimising the number of swabs stored at labs at the end of the time horizon. The objective function, then, is: The bounds on the quantities of reagent and swabs moved each day translate into the following two constraints: The next constraints limit the number of tests carried out at each lab on a given day, based on the quantity of reagent available and lab capacities. To this end, it is convenient to introduce an auxiliary variable w l t ∈ Z + representing the number of tests carried out at lab l ∈ L on day t ∈ T . Quantities m l t (number of tests requested each day), l ∈L, l l ll t (swabs moving from l to other labs), l ∈L, l l l l t (swabs moving from other labs to l), z l ,t −1 (swabs backlog from the previous day) and z l t (swabs stored at the end of the day) uniquely determine the value of w l t , according to the linking relation Figure 2 shows this relation visually. In the gure, each node represents lab l on a di erent day and arrows display the ow of swabs in and out of the laboratory. With linking variables w l t , we can constrain the number of tests performed with the following inequalities: Next, we consider ow balance equations analogous to eq. (4), relative to the number of reagent stored at factories and laboratories: Constraint (7) states that the quantity of reagent in the factory's inventory at the beginning of the day, plus the quantity produced, is equal to the quantity in inventory at the end of the day, plus the quantity shipped to the labs. Constraint (8) equates, on the left-hand side, the quantity of reagent in the lab's inventory at the start of the day and the total amount received with, on the right-hand side, the inventory at the end of the day, plus the quantity of reagent used (or, which is the same, the number of tests performed). Because all variables are non-negative, constraint (8) makes constraint (5) redundant. Model (1)-(8), together with the following variable domain de nitions, is the Base Model for our problem. x r l t ∈ Z + ∀r ∈ R, ∀l ∈ L, ∀t ∈ T (9) A disadvantage of the Base Model (1)- (12) is that it su ers from symmetry. For example, consider any solution with slack capacity to transport swabs between two labs l 1 , l 2 ∈ L during a day t ∈ T . One can get another solution of the same cost in which, on day t, l 1 sends one more swab to l 2 and l 2 sends one more swab to l 1 . To reduce symmetry, we impose two restrictions on the transport of swabs between labs: (i) a lab l can ship out swabs only if it is already working at full capacity, i.e., it is testing Q l swabs or it ran out of reagent; (ii) a lab can only send or receive swabs during a given day, but not both. To tackle the rst restriction, we introduce indicator variables γ l t ∈ {0, 1}, which take value 1 i lab l ∈ L is running at full capacity during day t ∈ T . These variables link to variables w and ρ via the following indicator constraints (see, e.g., [3] ): The restriction is, then, enforced with another indicator constraint: Constraint (16) guarantees that laboratory l cannot ship out swabs unless it is at full capacity. To address the second restriction, we can add two more sets of indicator variables and constraints, keeping track of whether a lab sends out or receives swabs. Note that the black-box commercial solver we use (Gurobi) converts all indicator constraints presented in this section into linear constraints. The new variables are γ + l t , γ − l t ∈ {0, 1}, which assume value 1 i , respectively, lab l sends or receives swabs on day t. The necessary constraints to link these variables and enforce the restriction are: Constraint (19) enforces that each lab, on each day, cannot both send and receive swabs. This restriction also avoids that a laboratory works as a relay, receiving and instantly re-shipping (i.e., "teleporting") swabs between otherwise distant labs. The solutions produced by the Base Model can have another property which is undesirable in real-life applications: consider two laboratories l 1 , l 2 ∈ L and two reagent factories r 1 , r 2 ∈ R such that l 1 is closer to r 1 than to r 2 and l 2 is closer to r 2 than to r 1 . The model allows l 1 to buy reagent from the more distant factory r 2 while l 2 buys reagent from r 1 . To avoid this counter-intuitive situation which represents a sub-optimal use of logistic resources, we introduce parameter d lr ∈ R + representing the distance between lab l ∈ L and factory r ∈ R and add the following constraint to the model: While (20) refers to reagent, we add an analogous constraint to prevent moving swabs ine ciently: where d ll ∈ R + is the distance between labs l ∈ L and l ∈ L. In this section we propose extensions to the model presented above, to account for more realistic scenarios. We present an Extended Model, taking into account these considerations, in Appendix A. In Section 2 we mentioned two limitations of the Base Model which might not hold true in real-life applications. The rst limitation is that there is a xed quantity of swabs assigned to each laboratory. In most countries, central or regional health authorities oversee and plan testing. Swabs collected in a region, for example, might be sent for testing to any laboratory in the same region. The second limitation is that tests (and reagent) can move overnight from any location to any other location subject to global capacity limits. When applying the model to large countries, though, logistic constraints might impose that movement of material only happen between geographically close locations. To address these aspects, we extend the initial problem formulation as follows. The set of laboratories is partitioned as L = L 1 ∪ . . . ∪ L n (with L i ∩ L j = ∅ for any two i, j ∈ {1, . . . , n}, i j). We call each set composing the partition a region. The planner receives as input the number of swabs to test each day in each region L i , denoted as m it , but must determine how to split the swabs between the region's labs. A parameter δ r l ∈ {0, 1} determines whether a factory r ∈ R can ship reagent overnight to laboratory l ∈ L. Analogously, parameter µ l 1 l 2 ∈ {0, 1} determines whether two labs l 1 , l 2 ∈ L can send each other swabs overnight. In our model, logistic capacities refer to each region. We consider, respectively, quantities q rg i and q sw i of reagent and swabs that can be shipped to region i ∈ {1, . . . , n} in a day. Depending on the logistic infrastructure, a planner could specify bounds at di erent levels of aggregation and for both inbound and outbound movements. For example, one might consider a maximum quantity of reagent shipped out of a factory, or a group of factories. Although we present here regional-level inbound capacities, our formulation is exible enough to allow a wide array of modelling options. A limitation of the Base Model is that it allows to move swabs but not reagent between laboratories. This re ects the conservative stance of laboratories, which are often unwilling to part with a scarce resource even if they have it in surplus, out of fear of running out later and not seeing their "favour" returned. Without capacities, moving swabs and reagent would be equivalent. For example, a lab l 1 could send a swab it cannot test to another lab l 2 , because l 1 does not have enough reagent to test it; or, which is the same, l 2 could send one unit of reagent to l 1 , which could test their swab by themselves. In the presence of capacities determined by equipment and personnel availability, though, the two options are not interchangeable. In the example above, l 2 might not have the capacity to test one more swab, while l 1 could have it if provided with the necessary reagent. In light of the above considerations, we propose an extension of the model in which reagent can be forwarded between labs: a practice which goes by the name of transshipment. A planner can also use transshipments as a way to extend the radius at which factories send their products, by rst storing them at a lab to which the factory ships and, when the need arises, forwarding them to other labs -even if these labs were not originally reachable from the factory with an overnight shipment. Because reagent take one day to move between labs, we are indirectly introducing distance-based lead times in the reagent supply chain. To account for transshipment of reagent, we extend the model introducing a new set of variables x l 1 l 2 t ∈ Z + indicating the amount of reagent moved from lab l 1 ∈ L to lab l 2 ∈ L on day t ∈ T . We assume that transshipments happen between laboratories which are close enough (µ l 1 l 2 = 1) and that a unit sent on day t is available on day t + 1. Figure 3 : Flow of reagent for a lab l ∈ L at time intervals from t −1 to t +1, when accounting for transshipments. Figure 3 shows the ow of reagent for a lab l ∈ L when accounting for transshipments. Each day, the lab has in its availability the reagent already stored there, plus that received from factories and from other laboratories. At the same time, it uses some units of reagent to perform tests and eventually sends some to other labs; the remaining units go to the inventory for next day. In the model, we replace constraint (8) with the following: where x l ,l ,t −1 = 0 if t = 0, and constraint (2) with: Lastly, we add a constraint analogous to (20) and (21) to ensure that reagent move e ciently between labs: x l 2 l 1 t = 0 ∨ x l 1 l 2 t = 0 ∀t ∈ T , ∀l 1 , l 1 , l 2 , l 2 ∈ L : d l 1 l 1 < d l 2 l 1 and d l 2 l 2 < d l 1 l 2 . As discussed in Section 2, the number of tests carried out is not the only important parameter for an e ective testing campaign. Having fast results also helps assessing the disease's spread and enact appropriate and well-timed con nement and mitigation measures. We propose, then, a Multi-Objective Optimisation (MOO) approach to the problem. For a recent review of MOO applications in optimisation, we refer the reader to Gunantara [17] and Kalyanmoy [20] . In the following, we consider the hierarchical objective approach, in which the planner rst optimises for their main objective (number of tests performed) and then, among all the solutions which give an optimal objective value, looks for the one optimising a secondary objective (mean test waiting time). Let Z * be the optimal objective value of model (27)-(53). We can obtain the solution to the hierarchical problem by solving the following Hierarchical Model: where eq. (25) minimises the sum of all amounts of swabs stored at labs waiting to be tested and eq. (26) imposes that the considered solutions match the optimal value with respect to objective function (27) . Objective function (25) is equivalent to minimising the average swab waiting time because (26) xes the total number of swabs tested. When minimising the primary objective, two untested swabs collected, respectively, on the rst and last day of the time horizon would contribute the same penalty. Adding the secondary objective avoids solution where old swabs remain untested for a long time, while newer ones get tested quickly. A planner could also add further objectives. For example, if he wants to guarantee a more equitable distribution of swabs to laboratories, he might want to minimise their avereage workload. To this end, we introduce a new variable W ∈ R + such that W ≥ 1 Q l t ∈T w l t for all labs l ∈ L. Minmising W corresponds to minimising the highest average workload across all labs. Alternatively, if W ≥ w l t Q l for all l ∈ L and t ∈ T , then minimising W amounts to minimising the peak workload over all labs and all days of the planning horizon. We want to assess the impact of optimising the allocation of reagent and swabs to laboratories on the testing capacity of a country. We present two case studies. In the rst, we perform an analysis on a realistic dataset relative to the Italian response to the COVID-19 pandemic, during the period from April 1 st to April 13 th , 2020. In the second, we use synthetic data to do a sensitivity analysis highlighting the decisions a planner can take to most increase testing capacity and lower waiting times. The datasets and the code used are available at Santini [35] . We use the Hierarchical Model, which we solve using commercial solver Gurobi [18] with a time limit of 15 minutes for the primary and 5 minutes for the secondary objective, on a laptop equipped with a 4-core Intel i7 processor running at 1.6GHz. This time limit proved enough to solve to optimality the majority of instances, while those that remain open all show gaps smaller than 5% and, often, smaller than 1%. We create a set of instances based on data sources relative to the COVID-19 pandemic in Italy [32, 25, 28, 38, 12] . We use data for the period 01-13 April, 2020 for the entire national territory. Figure 4a shows the distribution of laboratories and reagent factories in Italy [25] : dots represents labs which are o cially authorised to execute COVID-19 swab tests by the Italian Health Ministry, while stars represent factories whose testing kits (including reagent, extraction kits, probes, negative controls) have been certi ed for use in COVID-19 swab testing. Figure 4b reports the number of swabs analysed in Italy during the period 01-13 April, 2020, on a regional basis. We note that there is no publicly available data on the average wait times between swab collection and analysis. (a) Map of rRT-PCR laboratories (dots) and reagent factories (stars). (b) Number of swabs analysed per region, over the period from April, 1 to April, 13, 2020. We analyse the number of swabs tested, day by day and region by region, using ve di erent scenarios. The rst uses the o cial numbers released by Italy's "Dipartimento della Protezione Civile", which we refer to as Real Data. Next, we consider the output obtained by the Hierarchical Model, under the following hypotheses. Swabs collected in a region cannot be sent to other regions, but can be allocated to any laboratory within the collection region (for this reason we label this data as Model -Regional). In modelling terms, this corresponds to setting µ l 1 l 2 = 1 if and only if, for two labs l 1 , l 2 ∈ L, there is a region L k such that l 1 , l 2 ∈ L k . We use this hypothesis to model how testing happens in Italy, with local laboratories and regional "reference laboratories" collecting swabs from all-over the region and providing double checks before non-urgent hospital admissions (emergency patients are admitted even before test results are ready) and discharges. Also remind that the Italian healthcare system is managed on a regional basis, with ample autonomy given to the local health authorities but little inter-regional interaction [15] . We assume that labs can procure reagent from their geographically closest factory, and that procurement is optimised by a planner deciding the quantities of reagent sent to each laboratory, to maximise the number of tested swabs. This assumption re ects how Italian regions were procuring supplies autonomously, before the central government created and rendered operational a COVID-19 Commissariat which also acted as a national procurement centre [11, Art. 122 ]. Because news sources and press releases by regional health authorities (we found statements by authorities of Basilicata, Campania, Emilia-Romagna, Marche, Lombardia, Piemonte, Puglia, Sardegna, Sicilia, Toscana, Umbria and Veneto, i.e., 12 out of the 20 Italian regions) hint that reagent was a major bottleneck in expanding test capabilities, we expect this model to be able to increase the number of swabs tested even if it does not allow inter-regional swab reassignment, just by optimising reagent assignment. Figure 5a shows how laboratories cluster within each region: an edge between two laboratories means that they belong to the same region. Finally, we consider the output obtained from the Hierarchical Model, but allowing for reallocation of swabs between laboratories in di erent regions. This assumption e ectively means that regional daily demand can be shared across regions. Because swabs would need to move within a short time (e.g., overnight) we put limitations on the inter-laboratory distance that allows swab transfer. We consider three thresholds of 100Km, 200Km and 400Km and we denote the corresponding data as Model -100Km, Model -200Km and Model -400Km. (We make an exception for the island region of Sardinia, for which moving swabs outside of the region would be unfeasible even if there are other laboratories within 400Km.) Figure 5b shows how laboratories can transfer swabs when using the 200Km threshold. Compare Figures 5a and 5b to note how inter-regional collaboration gives a central planner more opportunities to move swabs between regions in case there should be a day with a demand peak in a particular area. We refer the reader to [35] for the full description of the instance generation process and to [34] for an interactive dashboard presenting the results of the analysis. Figure 6 shows the number of swabs tested, day-by-day, when using our proposed models Regional, 100Km, 200Km, and 400Km. Note how Model -400Km fully uses the system's capacity during the last period of the time horizon, where the number of tested swabs becomes a at line. During some days, a model with fewer lab transfer capabilities gives more tested swabs than a model in which more transfers are allowed. For example, on April 6, using Model -400Km results in fewer tests than Model -200Km and Model -100Km. This is not surprising if we consider that the models maximise the number of swabs tested over the whole planning horizon, even if this results in fewer tests on any particular day. Table 1 reports detailed results of the analysis. Each row corresponds to one day, and the columns list the real data and the results obtained using the model. Columns #tests report the number of swabs tested during that day. Columns %gain are the percentage gains obtained using the models, compared to the number of swabs actually tested in Italy in the described period. The last row gives aggregate results for the whole country. Note how using optimisation techniques results in a testing capacity increase between 14% and 39%, under the hypotheses explained above. In particular, assuming that tests in Italy were largely limited by the unavailability of reagent due to a poor allocation of supplies and low inter-regional collaboration, optimising these two aspects can allow to perform up to 250 000 more tests (the increase between Real Data and Model -400Km) during 13 days. Figure 7 shows example results from the Regional and 400Km models, for four representative Italian regions from the north (Piemonte and Emilia-Romagna), the centre (Lazio) and the south (Sicilia) of the country. Each Model with intra-regional (orange column) and inter-regional (green, pink and blue lines) swab transfers. line inter-regional transfer radii considered are 100, 200 and 400Km. chart reports the number of swabs analysed per day. Note how, for three of the four regions, the results of the Regional model agree with the real data, validating our modelling assumptions and hinting that regions managed resources locally e ciently. In all four regions, during the last period of the planning horizon, the models use the full laboratory capacity and their curves appear at, i.e., the limiting factor is not reagent availability anymore, but structural limitations such as testing machines or lab personnel. Note how, in Sicily during the middle period, the 400Km model provides a solution with fewer tested swabs than the Regional model and the Real Data, for two days. However, there was a catch-up on the tests during the last part of the planning horizon, resulting in more tests carried out in total. By conducting a scenario analysis on synthetic data, we want to learn how some aspects of the healthcare system of the area hit by an epidemic can a ect its ability to perform large-scale testing of the population. The factors we analyse are: 1. The e ects of having a healthcare system with a large level of independent decision-making at the sub-national level (as is the case, e.g., of Italian regions, Spanish autonomous communities, German lander or the United States) and the bene ts of increased cross-regional collaboration, common planning, and resource sharing. 2. The capillarity and production capacity of bio-chemical and pharmaceutical industries able to prepare the reagent needed to carry out rRT-PCR tests, especially in light of the increasing trend to o -shore such manufacturing capabilities to developing countries (see, e.g., [21, 42] ). 3. Reagent-independent laboratory capacity, i.e., the test capacity deriving by the number of available machines and personnel. This corresponds to the full laboratory capacity in the ideal case in which any amount of reagent were immediately available. 4. The logistic capabilities of the country, i.e., how labs can procure reagent from manufacturing facilities and transfer swabs away from saturated labs. We list in the following the characteristics of the instances we generate. We consider |L| = 100 laboratories, which we distribute in the square S = [0, 100] 2 in the euclidean plane as explained below. First, we de ne parameter λ L×R , which is the expected size of regions, i.e., the expected number of labs in each region. Next, we select n = |L|/λ L×R random points in [10, 90] 2 which we use as initial centroids for the regions (note how we excluded the external boundary of S). We place the labs in S by rst choosing a random centroid and then a random point in S which lies at distance of at most 20 from the centroid, and repeating this procedure |L| times. In this way, we aim at placing clusters in loose regional clusters. To complete the assignment of labs to regions, we assign each lab to its closest centroid. (a) Spatial set-up of the instance. Dots represent laboratories, with their size proportional to the average daily demand in their region. Stars represent factories and lines are regional boundaries. (b) Daily reagent production. Each smaller bar represents one factory and the larger bars represent the total daily production. Another parameter, λ F×R , determines how distributed is the country's reagent production capacity; this parameters represents the expected number of reagent factories per region. We place |R| = n · λ F×R factories at random in S. Figure 8a shows an example distribution of labs and factories. Dots represent laboratories, stars denote factories, and lines are the Voronoi boundaries (see, e.g., [31, Ch. 5] ) associated with regional centroids (i.e., the lines delimit the regions). To determine the demand of each region, we assume that 40% of the regions will have a criticality, such as a focus of the illness, and will thus require more tests. where ε is a noise term distributed according to a truncated normal distribution between 0.95 and 1.05. The de nition of m it is justi ed by the assumption that the number of swabs to test depends on regional population, which correlates positively with the number of laboratories in the region. The dot size in Figure 8a is directly proportional to the lab's region average demand in the example instance; note how two bordering regions are "critical". Lab capacities are then de ned as a fraction of the ideal number of tests assigned to each lab: Q l =mλ LC , where λ LC is a parameter. Values of the parameter larger than 1 correspond to a situation in which reagent availability is the main bottleneck for testing, while smaller values indicate that laboratories have structural de ciencies which limit the number of swabs they can test, even when reagent is abundant. Each lab can procure reagent from its λ LF closest factories, and can send swabs to other laboratories in the same region. A parameter λ LL gives the maximum distance at which two labs in di erent regions can exchange swabs. In particular, when λ LL = 0 we forbid cross-regional collaboration and sharing of testing capacity. Regional inbound capacities for reagent (q rg i ) and for swabs (q sw i ) are determined as: for 80% of the regions m · |L i | · 0.75 for 20% "impervious" regions q sw i = m · |L i | · 0.25 for 80% of the regions m · |L i | · 0.10 for 20% "impervious" regions, We assume here that the "impervious" regions have tighter logistical constraints; for example, islands, mountainous areas, or regions with poorer infrastructure. We set the initial amount of reagent stored at labs and factories to a minimal amount, as follows: where ε is distributed according to a truncated normal distribution between 0 and 0.25. Such a choice corresponds to labs and factories having a small fraction of their daily, respectively needed and produced, reagent in the initial stock -a situation that can commonly occur during an epidemic, when safety stocks are quickly depleted. Next, we describe how we determine factories' daily production. We propose two production patterns (selected via parameter λ PP ): steady and bumpy. Under the steady pattern, each factory r produces f r t =m · λ L×R λ F×R · λ FP · ε units of reagent each day, where ε follows a truncated normal distribution between 0.95 and 1.05 and λ FP is a parameter. Note how λ FP < 1 drives the total factory production capacity to be smaller than the total demand, indicating a situation in which there is less reagent than the nation needs. Values of λ FP ≥ 1 drive the total reagent production towards meeting the demand; in this case, tests not performed because of lack of reagent are more likely due to poor resource allocation over the national territory, rather than an absolute scarcity of reagent. Under the bumpy pattern, for each factory we rst choose two random days of the week τ 1 r , τ 2 r ∈ {0, . . . , 6}, which we call "release days". The factory will produce reagent at the same rate as under Table 3 : Results of a linear regression analysis using the parameters as independent variables and the primary (number of swabs tested, expressed as a percentage of the total demand) and secondary (average waiting time of a swab before it is tested) objectives as dependent variables. Columns coef report the linear regression coe cients, after centring and standardising the inputs. All independent variables resulted signi cant, with p-values always smaller than 0.001. the steady pattern, but will not make them available until the current day t corresponds to one of the two τ 1 r , τ 2 r (i.e., t = 0 mod τ 1 r or t = 0 mod τ 2 r ). At that point, it is ready to ship all the reagent produced since the last release day. This second pattern simulates a situation in which it is infeasible, due to logistic constraints, to have daily shipments out of the factories. Figure 8b shows an example of bumpy production. The coloured bars correspond to the quantities released by the factories on their release days, while grey background bars count the total amount of reagent released on that day, by all factories. The gure shows that, even if factories have, in principle, a steady production rate (only uctuating between −5% and +5% due to the noise term ε), making new reagent available twice per week produces an uneven total daily distribution. Last, we introduce transshipment parameter λ TS , which takes values yes or no. This parameter does not a ect instance generation; rather, it enables or disables the transshipment of reagent between laboratories in the model. If transshipments are not allowed, it forces all variables x ll t to take value zero. Table 2 summarises all instance generation parameters we used. Because in our scenario analysis we vary each parameter independently, we ran the Hierarchical Model over the full grid of 60 480 instances. We ran the model twice for each instance, once enabling and once disabling transshipments. We report the results of our scenario analysis, highlighting how the parameters a ect the testing capacity and commenting on which real-life decisions a planner can take to increase the number of swabs tested. A better testing capacity depends on both structural factors (such as a strong reagent industrial base) which require long-term decisions, and organisational factors (such as more inter-regional collaboration) which are in uenced by operational decisions. Figure 9 shows the percentage of swabs tested within the end of the time horizon, over all the swabs collected and available for testing. In other words, if the planner manages to allocate swabs to laboratories in such a way that all swabs are tested within the considered time horizon, then this percentage would be 100%. Each chart corresponds to one parameter; we use line charts for the seven numeric parameters, and box plots for λ PP and λ TS , which are categorical. In each chart, the x axis reports the value taken by the parameter. The axis reports the percentage of swabs tested when solving instances generated with the corresponding parameter. For the line charts, the solid line denotes the average percentage of swabs tested over all instances generated with the given parameter value, while the shaded area is the 95% con dence interval obtained via bootstrapping (using 1000 iterations). For the box plots, each box denotes the two central quartiles, the horizontal line is the median, and whiskers extend to the rest of the distribution. Figure 9 : Variation in the percentage of swabs tested ( axis), when varying each parameter (x axis). In each chart, the line shows the average over all instances generated with the given parameter value and the shaded area around the line is the 95% con dence interval. For parameters λ PP and λ TS , as they take categorical values, we show a box plot instead of a line plot. The central line indicates the median, the box spans the two central quartiles, while whiskers extend to the rest of the distribution (excluding outliers). Two linear regression analyses provide additional ndings. They use the instance parameters as independent variables. We centre and standardise numeric parameters, and transform λ PP and λ TS into numeric parameters, respectively, with value 0 if the production was steady and 1 if bumpy, and 0 when transshipments are not allowed and 1 when they are. The percentage of swabs tested and the average time a swab spent in a lab waiting to be tested are, respectively, the dependent variables in the two analyses. Table 3 reports the results, while in the following we comment on the insights obtained. Steadiness of reagent supply. Parameter λ PP is the one with the largest impact in our analysis: a steady availability of reagent allows the laboratories to operate without interruptions, testing more swabs and with shorter delays. The average percentage of swabs tested when λ PP = steady is double than when λ PP = bumpy. Recall that a factory in our model can represent a real production site, but also a warehouse or a shipment from abroad. A physical factory will be able to maintain a steadier production rate, while relying on shipments will resemble more the "bumpy" scenario. The analysis, thus, hints to the fact that a solid industrial sector and a well-function supply chain can have the greatest impact on testing capacity; not necessarily because they allow to procure more reagent in total, but because they provide a steadier ow of it. A country that cannot count on such production capacity can mitigate the adverse e ects of irregular supply with larger orders (see the positive impact of parameter λ FP ), more inventory capacity and a quick dispatch of reagent to labs. Capillarity of the reagent industry. Having more small reagent production sites works better than having few, larger ones (parameter λ F×R ). When λ F×R moves from 0.1 to 1.0, the average percentage of swabs tested climbs from 40.83% to 51.46%. An e cient supply chain which is able to move material over longer distances (parameter λ LF ) also has a positive impact on testing capacity. When each lab can procure from its closest factory (λ LF = 1), the analysis shows that 41.36% of swabs are tested. This number increases to 47.86% and then to 49.61%, when λ LF increases to, respectively, 2 and 3. Figure 10 : Variation in the percentage of swabs tested ( axis), when varying parameter λ LL (x axis). The chart on the left refers to instances with λ LC < 1, while that on the right, to instances with λ LC ≥ 1. Lab capacity and inter-regional capacity sharing. Other two factors which impact the amount of swabs tested are lab capacity and the degree to which regions can share this capacity by transferring swabs between each other. Parameter λ LC clearly in uences testing capabilities: for λ LC < 1 largely less than half of the swabs are tested within the end of the time horizon. When λ LC ≥ 1, i.e., when the total lab capacity is at least as large as the average demand, gains are marginal, as exhibited by the curve attening in Figure 9 . It is at this point that we would expect capacity sharing via increased inter-regional collaboration to have the greatest e ect. The extent at which capacity is shared depends on parameter λ LL . It might appear, looking at Figure 9 , that the impact of λ LL is small. However, if we divide the analysis of its impact between instances with λ LC < 1 and instances with λ LC ≥ 1, the situation is di erent, as reported in Figure 10 . Note how the curve referring to instances with λ LC ≥ 1 is much steeper, corresponding to a larger percentage increase in swabs tested. We can conclude that countries with insu cient lab infrastructures should primarily focus on procuring test machinery and training personnel to carry out rRT-PCR tests. Countries which already have a developed test lab network, on the other hand, bene t more from inter-regional collaboration and centralised planning than they do from increasing lab capacity. Other factors. Allowing transshipments produces an increase in swabs tested, from 43.30% to 45.85% on average (an increase of +5.03%), placing parameter λ TS 's impact behind several others. We suspect that the reason is that there are few cases in which moving reagent in one direction gives di erent results compared to moving swabs in the opposite direction, as mentioned in Section 5.2. Transshipments' e ect is larger if we focus on instances with bumpy reagent production. In this case they increase the number of swabs tested by +11.14%. This result suggests that transshipments can smooth irregular production patterns: a country which cannot procure reagent steadily should increase internal collaboration by moving both swabs and reagent. The percentage of swabs tested decreases when there are fewer, larger regions (parameter λ L×R ). This is because, with fewer regions, the impact of one region being critical, i.e., with larger demands m it , is proportionally larger. The e ect of region size is rather modest, as indicated by a small regression coe cient associated with λ L×R . Regarding the planning horizon length |T |, its impact is unclear. On one hand, a longer time horizon allows more opportunities for optimisation. On the other hand, there is a reduction in tested swabs when moving from |T | = 5 to |T | = 7, which is hard to explain. In practice, a planner would want to balance the advantages usually conferred by longer-term planning with the quality of the demand predictions that he can obtain from epidemiological models. Finally, while the signs of the regression coe cients in the two regression analyses tend to coincide (meaning that parameters which have a positive impact on the number of swabs tested also have a positive impact on shortening wait times), the signs disagree for parameters λ L×R and |T |, making it di cult to derive actionable recommendations for decision makers based on these two parameters. Table 4 : Relation between the percentage of swabs moved between laboratories (column "% swabs moved") over all swabs collected, and the success of the testing campaign as indicated by the percentage of swabs tested over all swabs collected (column "avg % swabs tested"). The second column contains averages over all values of the rst column belonging to each interval. To summarise the impact of resource sharing, Table 4 shows how the quantity of swabs tested changes depending on the quantity of swabs moved. Indicator "% swabs moved" is the percentage of swabs sent from one lab to another, over the total number of swabs collected. The table shows the correlation between swabs sharing and the e ectiveness of the test campaign, independently from how swab sharing was achieved. In our simulations swab sharing depends on instance generation parameters, but in real life it can depend on actions taken by the planner. We have proposed a model for the problem of assigning reagent and swabs to laboratories in a country which wants to maximise the number of rRT-PCR tests during a viral epidemic. After describing an Integer Programming model for a basic version of the problem, we have strengthened and extended the formulation to include real-life aspects of the organisation of healthcare systems. The complete model takes into account that many countries manage healthcare at a regional level and that labs can function as reagent transshipment points, and allows to optimise two hierarchical objectives: primarily, maximising the number of tests performed and, secondarily, minimising the average time a swab waits before it is tested. We performed two computational studies: the rst based on real-life data from Italy and the second on synthetic instances. Under the assumptions we used, the main one being that reagent availability is often a bottleneck for testing more swabs (see, e.g., [13, 1, 43] ), the analysis has revealed two main ndings: • Increased inter-regional coordination has the potential of increasing the testing capabilities of the country as a whole, as shown by the analysis of Italian data. In particular, moving swabs from critical regions to regions with spare capacity produced an increase of up to 40% in the number of swabs tested in the Italian scenario. The study performed on synthetic data con rmed this nding, but found a more modest e ect. The recommendation, in fact, is that countries which do not have adequate lab infrastructure (i.e., which would have insu cient test capacity even if fully stocked with reagent) rst focus on expanding their lab network, e.g., by purchasing new machines and training more lab personnel. Countries with an advanced laboratory infrastructure, instead, will bene t more from inter-regional collaboration. • Steady availability of reagent has a large impact on the number of swabs tested and their waiting time. A solid chemical industrial base, with many (although possibly small) factories producing reagent daily gives an advantage compared with having to rely on a discontinuous supply, e.g., when purchasing from international markets. If the supply is not steady, then we recommend to allow labs to work as small warehouses by letting them forward reagent to other labs, i.e., enabling transshipments. On the short term, increasing inter-regional collaboration requires political will and an e cient supply chain which can be used to move swabs and reagent. On the other hand, increasing reagent production capacity involves longer-term investments and balancing the economic interests of pharmaceutical companies with national interests. Future research on optimising the delivery of mass-scale test programmes can develop integrating other test types, such as immunological tests, i.e., tests which detect the presence of antibodies which a person has developed in response to the viral infection. These tests are too slow for early detection of positive patients; e.g., T and B cells response to COVID-19 were detected around one week after the onset of symptoms [39] , while swab tests can identify the virus even before symptoms appear. However, because antibodies persist in the blood even after a patient does not host the virus any longer, health authorities use immunological tests for mass-screening of populations (including asymptomatic patients) with the aim of assessing the prevalence of the virus, the level of population immunity and other epidemiological dynamics. It would be interesting to study not only how to increase capacity for immunological tests, but also to study the problem of selecting the patients to test with higher priority, in order to create an accurate epidemiological picture of the population while using a restricted supply of antibody tests such as ELISA (enzyme-linked immunosorbent assay) [41] or neutralization assays [45] . Why Widespread Coronavirus Testing Isn't Coming Anytime Soon Assessing alternative production options for eco-e cient food supply chains using multiobjective optimization On mathematical programming with indicator constraints Quanti cation of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems Principles of scarce medical resource allocation in natural disaster relief: a simulation approach Revision: 03. U.S. Department of Health and Human Services Processing of Sputum Specimens for Nucleic Acid Extraction Optimization can solve the ventilator shortage Ventilator Pooling: Formulation and Data Sources Simulation-optimization approach to clinical trial supply chain management with demand scenario forecast Misure di potenziamento del Servizio sanitario nazionale e di sostegno economico per famiglie, lavoratori e imprese connesse all'emergenza epidemiologica da COVID-19. Gazzetta U ciale della Repubblica Italiana Open Data Coronavirus European Centre for Disease Prevention and Control. Coronavirus disease 2019 (COVID-19) pandemic: increased transmission in the EU/EEA and the UK. Tech. rep. Seventh update Positioning Inventory in Clinical Trial Supply Chains The Italian health-care system Quantitative RT-PCR: pitfalls and potential A review of multi-objective optimization: Methods and its applications Gurobi Optimization LLC. Gurobi Optimizer Reference Manual Stockpiling ventilators for in uenza pandemics Multi-Objective Optimization The Geography of Prescription Pharmaceuticals Supplied to the US: Levels, Trends and Implications E ectively managing diagnostic tests to monitor the COVID-19 outbreak in Italy A Model of Supply-Chain Decisions for Resource Sharing withan Application to Ventilator Allocation to Combat COVID-19 Estimates of the demand for mechanical ventilation in the United States during an in uenza pandemic Pandemia di COVID-19 -Aggiornamento delle indicazioni sui test diagnostici e sui criteri da adottare nella determinazione delle priorità. Aggiornamento delle indicazioni relative alla diagnosi di laboratorio A multi-objective robust optimization model for multi-product multi-site aggregate production planning in a supply chain under uncertainty UNICEF provides COVID-19 testing reagents to support national pandemic response in Thailand Open Data Covid19 Asymptomatic cases in a family cluster with SARS-CoV-2 infection Production Planning by Mixed Integer Programming Computational Geometry: an Introduction COVID-19 Italia: Monitoraggio della Situazione High Contagiousness and Rapid Spread of Severe Acute Respiratory Syndrome Coronavirus 2 Increase COVID-19 Swab Testing Capacity via Optimisation Optimisation tools for the COVID-19 pandemic. Github repository The Nurse Rostering Problem in COVID-19 emergency scenario Recent advances in diagnostic testing for viral infections Istituto Superiore di Sanità. Epidemia COVID-19, Aggiornamento Nazionale The trinity of COVID-19: immunity, in ammation and intervention Integrated production planning and control: a multi-objective optimization model Unique epidemiological and clinical features of the emerging 2019 novel coronavirus pneumonia (COVID-19) implicate special control measures Safeguarding Pharmaceutical Supply Chains in a Global Economy Advice on the use of point-of-care immunodiagnostic tests for COVID-19 The SARS-CoV-2 outbreak: what we know Neutralizing antibody responses to SARS-CoV-2 in a COVID-19 recovered patient cohort and their implications A medical resource allocation model for serving emergency victims with deteriorating health conditions A conceptual framework for allocation of federally stockpiled ventilators during large-scale public health emergencies A multi-objective production planning problem with the consideration of time and cost in clinical trials The author would like to thank Michele Iacono of Roche Diagnostics S.p.A. and Xavier Jiménez Fàbrega of the Health Department of Catalonia for comments and fruitful discussions on this work. The author would also like to thank two anonymous reviewers for suggesting important improvements to this paper.Funding: the author was partially supported by grant "RTI2018-095197-B-I00" from the Spanish Ministry of Economy and Competitiveness.Declarations of interest: none. The extended model, which incorporates the considerations made in Section 5, uses a new variable u l t ∈ Z + , which represents the number of new swabs assigned to lab l ∈ L on day t ∈ T . The model reads as follows:x r l t = 0 ∀r ∈ R, ∀l ∈ L : δ r l = 0, ∀t ∈ T (30)x r l t ∀l ∈ L, ∀t ∈ T (35)x r l t ∀r ∈ R, ∀t ∈ T (37)x ll t ∀l ∈ L, ∀t ∈ T (38)x r 2 l 1 t = 0 ∨ x r 1 l 2 t = 0 ∀t ∈ T , ∀l 1 , l 2 ∈ L, ∀r 1 , r 2 ∈ R : (46)ρ l t , z l t , w l t , u l t ∈ Z + ∀l ∈ L, ∀t ∈ T (52)Without loss of generality, one can linearise indicator constraints of the type = 0 → ax ≤ b, for a binary variable and a non-negative integer variable x, introducing a non-negative slack variable s and transforming inequality ax ≤ b into ax + s ≤ b. Then the "big-M" inequality s ≤ M enforces the indicator constraint.Similarly, one can linearise constraints of type (x = 0)∨( = 0), where both x and are integer an non-negative. Consider two binary variables s x , s which will take value 1 if, respectively, x and are strictly positive. Then the following three inequalities linearise the original constraint: x ≤ Ms x , ≤ Ms , s x + s ≤ 1. Alberto Santini: everything. Other authors: no other authors. 26