Distribution equality as an optimal epidemic mitigation strategy ================================================================ * Adar Hacohen * Reuven Cohen * Sol Efroni * Ido Bachelet * Baruch Barzel ## Abstract Upon the development of a drug or vaccine, a successful response to a global pandemic, such as COVID-19, requires the capacity for efficient distribution at a global scale. Considering constraints on production and shipping, most existing strategies seek to maximize the outflow of therapeutics, hence optimizing for rapid dissemination. Surprisingly, we find that this intuitive approach is counterproductive. The reason is that focusing strictly on the quantity of disseminated therapeutics, such strategies disregard their specific spreading patterns, most crucially – they overlook the interplay of these spreading patterns with those of the pathogens. This results in a discrepancy between supply and demand, that prohibits efficient mitigation even under optimal conditions of superfluous drug/vaccine flow. Therefore, here, we design a dissemination strategy that naturally follows the predicted spreading patterns of the epidemic, optimizing not just for supply volume, but also for its congruency with the anticipated demand. Specifically, we show that epidemics spread relatively uniformly across all destinations, and hence we introduce an equality constraint into our dissemination that prioritizes supply homogeneity. This strategy may, at times, slow down the supply rate in certain locations, however, thanks to its egalitarian nature, which mimics the flow of the viral spread, it provides a dramatic leap in overall mitigation efficiency, saving more lives with orders of magnitude less resources. As international mobility continues to grow, so does the threat of global pandemics1–6 as clearly demonstrated in the worldwide inpact of COVID-197–12. Originating in a random outbreak, the pathogens spread internationally though air-travel, then locally, at each destination via contagion dynamics13,14. When such scenario transpires, the challenge is to (i) develop *ad hoc* a drug or vaccine; then (ii) design an efficient strategy for its dissemination. The point is that even if a therapeutic is avalibale, eliminating challenge (i), we must still ship it to multiple destinations worldwide – potentially stretching our transportation resources to their maximal capacity. Hence, to outrun the epidemic, we crucially need optimal dissemination schemes, that achieve the most efficient mitigation, given the constraints imposed by our production/shipping limitations15–19. The common approach to this challenge treats the dissemination as a commodity flow problem20–23. The therapeutics (commodity) are shipped from one or few sources, propagating along transportation routes – each with a given shipping capacity – aiming to fill the demand of all network nodes. In this framework one seeks the shipping sequence that optimally utilizes the network – namely, allows for the maximal *volume* of commodity to flow through the network per unit time24,25 **(Fig. 1a)**. This well-established approach, however, overlooks the interplay between the supply, *i.e*. the therapeutics, and the demand, here generated by the spreading pathogens. *Indeed, our true goal is not just to generate the maximal outflow of commodities, but rather to obtain the optimal mitigation of the disease*. This does not depend only on the *volume* of commodity flow, but also on the manner in which this volume is distributed across all potential destinations, *i.e*. its *spreading pattern*26,27. Consequently, the most efficient shipping sequence might be one that compromises maximal flow, but generates the optimal spreading patterns, compatible with those of the pathogens13,28–32. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/18/2020.09.15.20194506/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/09/18/2020.09.15.20194506/F1) Figure 1. Optmizing therapeutic dissemination for efficient mitigation We treat the theraputic spread as a commodity flow problem. (a) Commodities flow from the source *s* to fill the demands *Dn* **(***t***)** at all destinations *a* through *e*. The transportation network imposes constraints though the restricted shipping capacity of all links. For example, the route *s → a* can transport at most 5 units per day (edge labels). The shipping sequence determines the quotas shipped via each link *Lnm*, and hence the rates by which all nodes are supplied. Demands are updated as therapeutics are supplied – for instance, at *t* = 1 *b*’s demand *D*2(*t*) is updated to zero, as *b*’s demand is filled. Under *max-flow* the goal is to achieve the highest net outflow of commodities from *s, i.e*. maximize *Ls→*(*t*). Here, nodes were supplies within two days, a net flow L*s→*(1) *=* 9 in day 1 and *Ls→*(*t*) *=* 8, in day 2. (b) During an epidemic, however, the goal is to achieve optimal *mitigation*, not just optimal commodity flow. This requires us to consider the interplay between our dissemination and the spread of the epidemic. Each node’s demand is determined by its infection levels (*jn*(*t*), red), leading us to seek the shipping sequence which best adresses the anticiated spread of infections – aiming for supply (*qn* (*t*), blue) that is most compatible with demand. (c) We find that the optimal strategy is to design an *egalitarian* flow, in which supply is homogeneously spread across all destinations. Using the same network as in (a), our egalitarian strategy yields an even shipping sequence, in which nodes are supplied concurrently at roughly equal rates. For example, instead of fully supplying node *b* in *t* = 1 and only then shipping to *d* at *t* = 2, as done under max-flow in panel (a), our egalitarian algorithm favors the shipping sequence where *b* and *d* are simulataneasly supplied at equal rates. Under global demand, we show, such egalitarian supply is orders of magnitude more effective in terms of mitigation and resource disposal. With this in mind, we take here a network dynamics approach, and analyse the spatiotemporal propagation patterns of both the commodities and the pathogens **(Fig. 1b)**. We find that while the disease impacts the majority of nodes uniformly, maximum commodity flow optimization naturally yields a highly uneven spread, conflicting a concurrent and homogeneous global demand with an extremely heterogeneous supply pattern. As a result, not only does maximum flow not guarantee optimal mitigation, it is, in fact, counterproductive, indeed, generating a desirable *volume* of available therapeutics, but at the same time a highly undesirable *spreading pattern*. We, therefore, introduce an *equality* criterion to the optimization, that prioritizes dissemination sequences with homogeneous spreading patterns **(Fig. 1c)**. This ensures that the therapeutic flow is not just rapid, but also properly disseminated, generating supply patterns that are compatible with demand. The resulting dissemination, we find, is orders of magnitude more efficient, achieving higher mitigation levels even under extremely prohibitive shipping constraints. ## MITIGATION VIA MAXIMUM-FLOW To examine our response to a global epidemic, we consider different epidemiological scenarios, from mildly contagious to extremely virulent, in which a disease spreads globally via air-travel, under the susceptible-infected-recovered (SIR) epidemic model13,14,33 **(Box I)**. We used empirical data on human aviation to evaluate the flux of passengers between *N* = 1,292 local populations (nodes), each with *Mn* individuals (*n* = 1*, …,N*), and quantified the impact of the epidemic through its global *coverage* ![Formula][1] where *rn*(*t*) is the fraction of recovered individuals in *n* and ![Graphic][2] represents the global population. Hence, *R*(*t*) in (1) captures the fraction of impacted individuals worldwide. In **Fig. 2a** we show the unmitigated *R*(*t*) under three different scenarios: extremely contagious, where *R*(*t* → *∞*) → 1 (dark red), medium (*R*(*t* → *∞*) ≈ 0.5, red) and mild (*R*(*t* → *∞*) ≈ 0.3, light red). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/18/2020.09.15.20194506/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/09/18/2020.09.15.20194506/F2) Figure 2. Outrunning a contageuous epidemic using max-flow vs. egalitarian dissemination (a) The global coverage *R* (*t*) vs. *t* following an outbreak at Burundi (BJM) under three scenarios: severe (dark red ℛ *=* 10), intermediate (red, ℛ *=* 2) and mild (light red, ℛ *=* 1.5). (b) The state of the epidemic at *t* = 56 days, directly before drug dissemination begins. The local coverage *rn*(*t*) in each node and the flux of infected individuals along each link (air-route) are represented by their red color depth. (c) *Max-flow dissemination:* starting at *t*R *=* 56 we begin drug dissemination from Osaka (ITM) via max-flow; Eqs. (2)-(5). Drug fluxes *Lnm*(*t*) (links) and drug availability *qn*(*t*) (nodes circumference) are represented by blue color depth. We observe a *race* between the therapeutic and the disease, both spreading along similar routes, ending in a significant fraction of infected individuals, as indicated by the prevalence of red nodes at *t* = 100. (d) *Egalitarian dissemination:* under the same scenario, the egalitarian strategy of Eqs. (11)-(13) achieves a much more efficient mitigation, where by *t* = 100, the epidemic is almost eliminated. (e)-(g) *R*(*t*) at the three timepoints: unmitigated (red), under max-flow (yellow) and under egalitarian (green). By *t* = 100 egalitarian has reduced the global coverage by a factor of ~ 20, from *R*(*t*) *=* 0.65 to *R*(*t*) *=* 0.04, while max-flow has achieved a mere ~ 1.5 factor reduction. Three crucial parameters may impact our mitigation: the global transportation capacity *C*, our response time *t*R and the severity of the contagion ℛ. (h) Mitigation efficiency *E* vs. *C* under egalitarian (green) and max-flow (yellow). The former achieves higher efficiency with orders of magnitude less resources. For example, egalitarian provides *E* = 0.8 under *C ~* 10−3, as compared to *C ~* 10−1, a two order of magnitude gap, required under max-flow for similar efficiency. (i) *E* vs. *tR*. Again, we find that egalitarian shows a significantly higher performance in the face of a late response. (j) *E* vs. ℛ, confirming, again, the consistent advantage of egalitarian mitigation. In each panel we vary one parameter and list the set values of the other two. Following the initial outbreak at *t* = 0, we define the response time *t*R as the time required to begin the distribution of a therapy. This therapy can be in the form of a vaccine, designed to prevent the infection of *susceptible* individuals, or a drug, facilitating the recovery of *infected* individuals. As we consider a scenario in which the disease is rapidly spreading, and hence many individuals may already be infected, we focus below on the dissemination of drugs, which remain relevant at all stages of the spread34,35. Examining the efficiency of *dissemination*, rather than production, we assume that the therapeutic is already stockpiled in sufficient quantities at a specific source node *s*. Therefore, the main challenge is to optimally distribute it via the air-transportation network to all destinations. The network is characterized by the *carrying capacities Bnm* that quantify the daily volume of drugs that can be shipped through each air-route *m → n, i.e. Bnm = Bn←m*. The challenge is, therefore, to obtain the optimal shipping sequence from s to all other nodes, under the constraints imposed by *Bnm*. Max-flow mitigation (Fig. 2a-c). In the maximum flow dissemination strategy we seek to maximize the daily volume of drug doses spreading from s to the rest of the network. Denoting the number of doses shipped from *m* to *n* on day *t* by *Lnm*(*t*), our optimization translates to ![Formula][3] where ![Formula][4] captures the net daily outflow of drugs from *n*. Hence, in Eq. (2), by maximizing *Ls→*(*t*)*, i.e*. the outflow from the source, we seek to enhance the total volume of drugs introduced into the network in each day. This maximization is subject to two constraints: ![Formula][5] ![Formula][6] The first constraint (4) ensures that the daily flux along each route *m → n* is within the bounds of the route carrying capacity *Bnm*. The second constraint (5) restricts nodes from accumulating doses in excess of their remaining demand at time *t*. Hence *n’s* net *incoming* flux *Ln←*(*t*) is bounded by *n’s* current demand, *Dn*(*t*)*;* note that *Ln←*(*t*) *= −Ln→*(*t*). We consider three strategies for setting the demands: (i) *Population based*. At *t* = 0 we set *Dn*(0) to be proportional to each node’s population *Mn*, then update these demands daily according to *n’s* supply, as *Dn*(*t* + l) = *Dn*(*t*) *− Ln←*(*t*) • (ii) *Impact based*. Setting *Dn*(0) according to the projected number of infected individuals at *n* • (iii) *Urgency based*. Updating *Dn*(*t +* l) dynamically based on *n’s* supply gap at *t*, such that highly impacted (or undersupplied) nodes are prioritized. While (i) is simplest, it is also inefficient, tending to over-estimate the actual demand at *n;* (iii), on the other hand, is highly complex, relying on real-time tracking and updating of the demands, but, in some cases, benefits, by design, from a high level of congruency with the spread of the pathogens. An expanded discussion of these three strategies appears in **Supplementary Section 2.4**. **Network capacity**. Applying (2) – (5), the maximal dissemination rate one can achieve is restricted by the out-degree of the source node, ![Graphic][7], representing an upper-bound on the volume of drugs that can be introduced into the network per day. We, therefore, define the network’s normalized *capacity* with respect to each source *s* as ![Formula][8] normalized to *Cs* = 1 in case *s* can disseminate sufficient doses to meet the entire global demand (Ω) in one day. The network capacity can be controlled by rescaling all individual carrying capacities *Bnm*, to describe affluent vs. restrictive dissemination scenarios **(Supplementary Section 2.3)**. In **Fig. 2b,c** we present the evolution of the epidemic at four selected time-points. At *t* = 0 we simulate an outbreak (red) at Burundi (BJM), emulating the 2013 Ebola, which originated in Africa36,37, then track its spread through air-travel, setting ℛ *=* 2. The node infection levels and the epidemic fluxes, *i.e*. the daily volume of infected passengers on each route, are represented by red color depth. Drug dissemination via max-flow optimization (blue) begins at *t*R *=* 57 days in Osaka (ITM), using blue color depth to signify the availability/flux of drugs in each node/route. We set the network capacity in Eq. (6) to *C* = 10*−*2, a dissemination capability of 1% of the global demand per day, and assign demands *Dn*(*t*) using the impact based strategy. We find, through the long-term prevalence of infections (**Fig. 2c**, red), that under these conditions, mitigation falls short. Indeed, in **Fig. 2e-g** we observe that the coverage *R*(*t*) in (1) is only slightly affected by our max-flow mitigation (red vs. yellow), illustrating the failure to effectively suppress the epidemic. For a more systematic assessment of our dissemination strategy, we track the mitigation efficiency via15 ![Formula][9] where *R∞*(*C*) is the long term coverage, *i.e. R*(*t* → *∞*), of the epidemic under mitigation with network capacity *Cs = C*. In case mitigation is effective we expect *R∞*(*C*) ≪ *R∞*(0), a significant reduction in the disease coverage, which translates to *E →* 1. Conversely, a failed mitigation leaves infection levels almost unchanged, leading in (7) to *E →* 0. Testing *E* vs. *C* under maximum flow dissemination, we find that for a broad range of *C* levels the epidemic is almost unaffected, a consistently inefficient mitigation in which *E* ≪ 1 (**Fig. 1h**, yellow). Effective mitigation is only achieved around *C* ≳ 0.1, a limit in which s is capable of shipping doses in the order of the entire global demand in few days. Such optimal conditions are not only unlikely, but mainly, they indicate the inefficiency of this dissemination strategy, requiring an extreme volume of therapeutics shipped in a highly constrained timeframe in order to achieve a measurable impact on the epidemic. It seems, therefore, that the max-flow optimization strategy is inadequate for the containment of a globally spreading epidemic. In **Supplementary Section 2** we analyze the population and urgency based strategies, observing similar challenges. ### The roots of the max-flow inefficiency Our mitigation is, in its essence, a spreading competition between the therapeutics and the pathogens, both progressing along the same underlying network, *i.e*. air-transportation19. It seems, therefore, that winning this competition is a matter of shipping capacity: we must generate sufficient therapeutic fluxes to outrun the spread of the disease, namely we must increase *Cs*. However, the analysis above indicates, that there is an intrinsic deficiency in the spread of therapeutics, that cannot be easily compensated by simply increasing shipping rates. Next, we show that the two competing processes – viruses vs. therapeutics – indeed, lead to fundamentally different spreading patterns, in which the viruses benefit from an intrinsic advantage. **Propagation of pathogens (Fig. 3a-c)**. Viruses spread via diffusion coupled with local SIR dynamics, as captured by Eqs. (15) in **Box I**. In this process, upon penetration of node *n*, the viruses reproduce locally through SIR, until reaching peak infection at ![Graphic][10]. In a network environment, since the majority of nodes are at the mean distance from the initial outbreak, we find that after a limited propagation time, most nodes reach peak infection approximately simultaneously. We track this via the infection rate ![Formula][11] capturing the average *speed* by which the pathogens propagate to *n*. In **Fig. 3c** we show ![Graphic][12], the probability density for a randomly selected node *n* to have ![Graphic][13]. Indeed, we find that ![Graphic][14] follows a bounded distribution, where the majority of nodes are infected at a similar rate, and hence their peak infection occurs at roughly at the same time. Such uniform propagation patterns lead to a simultaneous global peak infection, and, as a consequence, to a concurrent global demand for therapeutics. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/18/2020.09.15.20194506/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/09/18/2020.09.15.20194506/F3) Figure 3. Why max-flow provides sub-optimal mitigation (a) Pathogens spread via contagion dynamics, reproducing in each desination independently via the SIR process. The result is a relatively homogeneous spread, in which the majority of destinations reach peak infection (red) at approximately the same time. Here this is observed by the almost uniform red shade and size of all nodes, capturing the instantaneuous infection level *jn*(*t*) during the time of global peak infection. (b) To quantify this we measure the infection rate ![Graphic][15] of all nodes as obtained from Eq. (8), and plot its distribution ![Graphic][16]. The bounded form of P(I) indicates the homogeneous nature of the pathogen spread. Indeed, we find that 80% of the nodes are impacted wthin a small margin around the average (shaded area). (c) Using max-flow we track the availability of the therapeutics *qn*(*t*) at all destinations (node size/blue color depth). In contrast to the pathogens, the therapeutics naturally spread extremely unevenly, with a small minority of early supplied nodes (large, dark blue), coexisting alogside a majority of delayed destinations (small, light blue). (d) The supply rate distribution ![Graphic][17] as obtained via max-flow (circles). As predicted ![Graphic][18] can be approximated by a power-law ![Graphic][19], with *ν* = 2 (solid line). This captures an extremely heterogeneous dissemination pattern, in which 80% of the nodes are supplied at a below average rate (shaded area). Hence, max-flow confronts a homogeneous demand (*jn*(*t*), red) with an extremely uneven supply (*qn*(*t*), blue). (e) The local efficiency *En* vs. supply rate ![Graphic][20] as obtained via max-flow dissemination. As expected we find that ![Graphic][21] crucially impacts the mitgation effectiveness at *n*. Therefore, the fat-tailed nature of ![Graphic][22], in which the majority of nodes have a small ![Graphic][23] (80%, shaded) creates an intrinsic mitigation defficiency, in which most destinations experience insufficient supply. The challenge, we emphasize, is not in the production/shipping capacities, as indeed we have ![Graphic][24], namely, that on average, supply rate can meet the evolving demand. Rather, is is rooted in the patterns of spread of the existing supply, which are fundamentally incongruent with those of the pathogens. **Dissemination of therapeutics (Fig. 3d-f, Supplementary Section 4)**. The dynamics of commodity flow, as provided by Eqs. (2) – (5) are fundamentally different. The therapeutics flow from a single source node *s*, and undergo dilution as they disperse across the exponentially growing number of pathways38,39. Such flow patterns naturally lead to a highly heterogeneous supply pattern across all nodes. To observe this we consider the therapeutic supply time ![Graphic][25] for node *n* to fill its therapeutic demand, and its consequent supply *rate* ![Formula][26] where ![Graphic][27], is the elapsed time from initial dissemination (*t*R) to *n*’s final supply ![Graphic][28]. This rate ![Graphic][29] quantifies how rapidly *n* receives treatment. In **Supplementary Section 4** we show that under max-flow mitigation the probability density ![Graphic][30] follows a power law of the form ![Graphic][31], with *ν* = 2 **(Fig. 3f)**. Hence, in contrast to the homogeneity of the epidemic spread, *supply* is extremely heterogeneous, with a vast majority of undersupplied nodes (small ![Graphic][32], shaded 80%), and a selected privileged minority of well-treated destinations (large ![Graphic][33]). Together, this combination of *homogeneous* demand and *heterogeneous* supply, a consequence of the naturally occurring spreading patterns of pathogens versus commodities, creates a crucial gap in our ability to effectively mitigate a global pandemic. To understand this, consider ideal flow conditions, where ![Graphic][34], namely that, on average, dissemination is more rapid than the viral propagation. Still, the fat-tailed nature of ![Graphic][35] vs. the bounded form of ![Graphic][36] indicates that while most nodes have ![Graphic][37], when it comes to supply, the typical node has ![Graphic][38]. Indeed, in a scale-free distribution the majority of entries are below the mean. This is illustrated in **Fig. 3**, where we show that 80% of nodes are within a 20*%* margin of ![Graphic][39], *i.e*. they are all impacted within a narrow time window around the mean (**Fig. 3b**, shaded). Yet, in contrast, a similar 80% have a therapeutic rate below ![Graphic][40], *i.e*. they witness supply deficiency (**Fig. 3d**, shaded). Hence we are confronted with a reality in which almost all nodes require treatment, and yet only a small minority receives sufficient, and at times, even superfluous, supply. This uneven distribution ![Graphic][41] directly impacts our mitigation efficiency. To observe this, consider the local efficiency ![Formula][42] where *rn,∞*(*C*) *= rn*(*t* → ∞) under mitigation with capacity *Cs = C*. Hence, while *E* in (7) captures the global mitigation efficiency, *En* focuses specifically on the effect observed at *n*. In **Fig. 3e** we show *En* vs. ![Graphic][43]. As predicted, we observe that ![Graphic][44] is a crucial determinant of mitigation efficiency. Indeed, nodes with low ![Graphic][45] (shaded), by far the majority of nodes according to the power-law structure of ![Graphic][46], exhibit an almost vanishing *En*. This illustrates, once again, the flaw of max-flow: most nodes, the deprived 80%, have *En →* 0, while a small minority benefit from *En →* 1. Taken together, our analysis shows that the common approach of maximizing flow is insufficient, as it treats the *mean* flow, but not its *distribution*. Most essentially, it disregards the interplay with the competing flow of the pathogens. When confronting these two spreading processes, we observe that the main challenge is rooted in equality ![Graphic][47] rather than in quantity ![Graphic][48], namely not the volume of the outflux, *Ls→*(*t*), but rather the way this volume distributes across the network. Hence, below we complement our optimization for max-flow with an additional requirement, seeking supply homogeneity. ## OPTIMIZING FOR HOMOGENEITY **Egalitarian mitigation (Fig. 1c)**. The optimization of Eqs. (2) – (5) is designed to increase the volume of drugs introduced into the network, but as we have shown, it leads to an extremely sub-optimal dispersion pattern. To remedy this, we introduce a homogeneity criterion, replacing the maximization in (2) by (**Supplementary Section 2**) ![Formula][49] where the homogeneity coefficient *H* determines the relative balance between our optimization for max-flow (*Ls→*(*t*)) vs. our demand for egalitarian spread (*λ*(*t*)). In addition to constrains (4) and (5) above, we now also require ![Formula][50] and ![Formula][51] forcing, at each step, to supply, at the least, a fraction *λ*(*t*) of every node’s current demand. This strategy seeks to maximize flow, but at the same time favors solutions that supply all nodes simultaneously. Indeed, constraint (12) prohibits shipping sequences that at any step supply exclusively a subset of the nodes. This additional constraint may, in general, lead to a diminished net flow, however, thanks to its egalitarian nature, we find that it dramatically enhances mitigation efficiency. To observe this, we revisit the epidemic spread of **Fig. 2**, leaving all conditions unchanged, except that now, instead of the max-flow optimization of (2) – (5) we mitigate the disease via the egalitarian strategy of (11) – (13). The results are striking: egalitarian drug dissemination practically eliminates the epidemic, as observed by the dominance of blue nodes/links at *t > t*R (**Fig. 2d**). Indeed, compared to the ~ 30*%* reduction in *R*(*t* → ∞) afforded by max-flow, egalitarian achieves, under the same capacity *Cs*, a ~ 95*°%* reduction, representing a practically perfect mitigation **(Fig. 2g)**. A systematic comparison of the two mitigation strategies consistently supports the crucial role of dissemination homogeneity. In **Fig. 2h** we show the efficiency *E* in (7) vs. *C* under egalitarian dissemination (green). We find that *E* approaches 80% already at *C ~* 10−3, representing an extreme case of supply deficiency. A similar *E* under max-flow was only observed for *Cs ~* 10−1 (yellow), a two order of magnitude advantage for egalitarian. A crucial factor impacting our mitigation efficiency, is the response time *t*R, required to identify the threat and initiate a response. To observe this, in **Fig. 2i** we present *E* vs. *t*R for both max-flow (yellow) and egalitarian (green) mitigation. As expected, we find that *E* declines with *t*R, however, for the entire range of response times egalitarian consistently outperforms max-flow. In **Fig. 2j** we further show that the egalitarian advantage is consistently maintained for a range of contagion levels ℛ. Hence, examined against a series of potential challenges, from the disease parameters (ℛ) to the performance level of our response (tR, *Cs*), we find that egalitarian allows enhanced mitigation, in some cases, by a dramatic margin. This, again, is thanks to egalitarian’s homogeneous dissemination patterns, which, under the price of a potentially diminished supply *rate*, afford us a more desirable supply *distribution*, and consequently, an optimized mitigation. ### The roots of the egalitarian advantage The efficiency of egalitarian mitigation may seem, at first glance, implausibly high. For example, **Fig. 2h** indicates that effective mitigation is already achieved with capacity as low as *C* ~ 10−3 day-1. This represents an extreme scenario, in which we require an order of 103 days to supply the global demand, a time scale that by far exceeds the final spread of the epidemic. Such extraordinary efficiency is rooted in the fact that even a tiny influx of drugs introduced at an early enough stage of the epidemic can dramatically impact its future development. To understand this, consider a node *n* at the early stage of the spread, when *jn*(*t*) ≪ 1. As time advances, this small seed of infections reproduces until *n* reaches its final infection levels *rn*(*t → ∞*). Suppressing this seed when it is small is, therefore, the optimal strategy, allowing us, with a tiny supply of therapeutics, to eliminate not just the current *jn*(*t*), but also terminate its potential future reproduction. Hence, maintaining a continuous influx of therapeutics, even at a very low rate, keeps *jn*(*t*) subdued, never allowing the local infections to reach their potential coverage. Egalitarian mitigation is precisely designed to provide such supply patterns: instead of rapidly supplying a select group of nodes at a time, it favors a slow influx that is spread evenly across all nodes, starting from the *get-go* at *t = t*R. In contrast, max-flow generates a similar (or even higher) total flux, but spreads it out sequentially – few nodes receive all their supply early on, and are hence successfully treated, while the majority fill their demand at a much later time, when *jn*(*t*) is already harder to contain. To observe this, we revisit the simulated mitigation of **Fig. 2**, this time focusing on a pair of specific nodes **(Fig. 4a)**: SDJ which is among the few early supplied nodes (![Graphic][52], red), and CGQ, a typical node, which represents the majority destinations, which receive their supply later on (![Graphic][53], Color). In **Fig. 4b,c** we show the drug availability *qn*(*t*) in both nodes (blue), together with their unmitigated (red) vs. mitigated (yellow) disease coverage. Indeed, we find that early supply is crucial, with SDJ benefiting from an almost perfect mitigation, and CGQ being almost unaffected by the late arriving drugs. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/18/2020.09.15.20194506/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/09/18/2020.09.15.20194506/F4) Figure 4. The roots of the egalitarian advantage (a) We focus on two specific nodes, SDJ (right), who, under max-flow, is among the early supplied nodes ![Graphic][54], and on CGQ (left), representing the majority of undersupplied destinations ![Graphic][55]. (b) Therapeutic availability *qn*(*t*) vs. *t* at CGQ (blue) as obtained under max-flow dissemination. Due to the sequential nature of this dissemination scheme, the therapeutic reaches CGQ too late, and consequently its infection level *jn*(*t*) remains high, both without (red) and under the effect of max-flow mitigation (yellow). (c) In contrast, SDJ is among the first to fill its demand, hence the sharp rise in *qn*(*t*) at an early stage (blue). As a result this node benefits from highly efficienct mitigation, with *rn*(*t → ∞*) reduced dramatically from 80% (unmitigated, red) to a mere 15% (mitigated, yellow). The challenge is that while CGQ is representative of the majority of destinations, the highly benefited SDJ belongs to a select minoroty under the max-flow framework. (d)-(e) Egalitarian generates a homogeneous supply pattern, hence now both CGQ and SDJ share a similar (low) supply rate of ~ 10−3. Despite that they both benefit from a near perfect mitigation (unmitigated – red; mitigated – green). This is thanks to the fact that despite their *equal but slow* supply, they now both begin to receive the therapeutic at an early stage (*qn*(*t*), blue). Indeed, even SDJ, who is now, under egalitarian, experiencing a dramtic drop in its supply rate, from 0.2 to 10−3, continues to benefit from highly efficient mitigation. Hence, the simulataneous, albeit uniformly slow, treatment afforded via egalitarian, not only benefits the max-flow deprived majority, but also continues to effectively treat its superflous minority. (f) The global consumption *Q*(*t*) vs. *t*, approaching unity, *i.e*. all demands are filled, when *t → ∞*. (g) Global mitigation efficiency *E* vs. the volume of disseminated therapeutics *Q*. Under max-flow (yellow) efficient mitigation *E →* 1 requires *Q →* 1 (yellow). In contrast, egalitarian has *E ~* 1 already at *Q ~* 0.1, achieving, thanks to its concurrent supply patterns, an almost perfect mitigation with only 10% of the resources. Under egalitarian mitigation the supply pattern is fundamentally different: indeed, as per the equality constraint, both nodes settle for the slower rate ![Graphic][56]. Hence, egalitarian provides no overall rate advantage, and, in fact, disadvantages the high rate SDJ. However, now, instead of *qn*(*t*) remaining idle and then, in due time, receiving a sharp supply boost, it is continuously supplied from the outset at *tR* **(Fig. 4 d,e)**, a direct consequence of the homogeneity constraint of (11) – (13). The advantage is that now both SDJ and CGQ receive treatment early on, mitigating the disease while it is still small, and thus suppressing its potential to grow. Therefore, while egalitarian impedes the supply of the high rate minority, it still allows us to successfully treat both that minority and the low rate majority. This is because its simultaneous, even if slow, supply, affords early treatment for *all* nodes, overcoming *jn*(*t*) while it is still at its embryonic stage. ### Mitigation under restricted supply A crucial implication of the above analysis touches on the global resources required for mitigation40,41. Clearly, running the max-flow/egalitarian algorithm until all demands are filled will inevitably cost a total of ~ Ω drug units, supplied at a time-scale of ~ Ω*/Cs* days. However, since egalitarian begins supplying all nodes while *jn*(*t*) is still small, it has the potential to overcome the pandemic even before all initial demands *Dn*(0) are completely satiated. To observe this, we run both strategies – max-flow and egalitarian – for a limited time, from *t* = 0 to *T*. We then measured the efficiency *E* (7), in function of the total resource consumption ![Formula][57] Equation (14) quantifies the total outflux of therapeutics from *s* as a fraction of the global demand Ω, throughout the mitigation period, *t* = 0 to *T*. The greater is *Q*(*T*), the more therapeutic resources that are consumed. As predicted we find in **Fig. 4g** that egalitarian, not only provides a higher mitigation efficiency, but achieves this in a fraction of the resources. For max-flow (green) we observe a linear gain with *Q*, reaching peak efficiency only around *Q* → 1*, i.e*. when supply fills *all* demands. Egalitarian, in contrast, thanks to its simultaneous treatment of all nodes, suppresses the epidemic at an early stage, when it is still manageable, with limited supply, therefore enabling mitigation of above 90*%* (grey dashed line) with total consumption as low as *Q ~* 0.1, *i.e*. a mere 10*%* of the global demand. ## DISCUSSION Commodity flow problems are at the heart of many crucial applications, from communications to supply chains, seeking to optimize a network’s capacity to distribute information or goods20-25,42-44. Most often the target optimization function, *e.g*., max-flow, arises naturally from the purpose of the distribution. However, in the context of disease mitigation, the optimization is a consequence not just of the dissemination scheme, but also of its interplay with the viral spread. Indeed, our true goal is not just distribution efficiency in and of itself, but rather *mitigation* efficiency. Hence, to assess a dissemination protocol we must couple it with the spread of the epidemic, and observe its effectiveness in terms of the actual observed reduction in infection levels. Here we have shown that the naive approach of maximizing the commodity flow fails this test, indeed, providing rapid supply, but at the same time a highly inefficient mitigation. On the other hand – optimizing for homogeneity dramatically improves mitigation efficiency. Our findings are relevant in case the epidemic spreads globally – a scenario averted in many recent pandemics, but experienced quite severely during COVID-19. In such scenario, the peak infection occurs approximately simultaneously at logarithmic time-scales, due to the short network paths27,45. Our egalitarian protocol, by spreading simultaneously, albeit slowly, to all nodes outlines the most effective strategy, catching the disease at its early stage in all locations. To fully benefit from this advantage, it is crucial to activate the egalitarian dissemination as early as possible, to truly attack the epidemic at its seed. Otherwise, social distancing measures must be employed to delay the viral spread and keep it at manageable levels until egalitarian dissemination can be instigated46–52. In this regard, the recent successes in epidemic forecasting53 are key: first, predict whether the pathogens *will* spread globally, then instigate egalitarian dissemination. Of course, a false alarm might result in a waste of resources, however our analysis indicates that early intervention allows mitigation with relatively small amount of doses, hence we believe that the potential gain by far outweighs the risk. Battling a global spread, it is crucial to appropriately set the demands *Dn*(*t*), which determine how much supply will be shipped to each node. On the side of caution, the natural tendency is to overestimate these demands, *i.e*. generate redundancy. However, in the face of a virulent epidemic, where the distribution network is stretched to its maximal capacity, such strategy may over-allocate resources to some destinations, while depriving others. On the other hand, fine-tuning the demands to meet the exact supply gap of each node, as we do, for example in our *urgency-based* strategy, is highly complex, not always practical, and, if inaccurate, may lead to undersupply in some destinations. Our egalitarian optimization addresses this trade-off quite naturally, by simultaneously supplying all nodes. To understand this consider two nodes *n* and *m*. Under max-flow, supply will be typically sequential, *e.g., n* and then *m*. Under these conditions, overestimating *Dn*(*t*) will lead to the shipment of more drugs to *n*, thus wrongly delaying the treatment of *m*. In contrast, such discrepancy under egalitarian, will have no bearing on *m*, as both nodes are simultaneously supplied from the start. This makes egalitarian highly robust against the specific assignment of *Dn*(*t*), allowing us to design a simple dissemination protocol, free of the need to accurately tailor and fine-tune all demands. Our analysis focuses solely on the mitigation efficiency, yet its conclusions touch upon the ethics of resource allocation – here, prioritizing the even distribution of drugs in a global emergency. While often one is confronted with a clash of values in such cases, here we find it encouraging that efficiency and equality are, in fact, compatible – representing a case in which scientific results go hand in hand with moral directives. ## Data Availability Upon publication, the authors plan to make all datasets and numerical codes publicly available. ## Author contribution All authors designed and conducted the research. AH analyzed the data and conducted the numerical simulations, AH, BB and RC lead the analytical derivations, and IB and SE advised on the biological applicability. BB was the lead writer. ## ADDITIONAL INFORMATION **Competing interests**. Authors declare no competing interests. **Data availability**. Upon publication, the authors plan to make all datasets and numerical codes publicly available. #### Box I. Modeling the mitigation of a global epidemic In a network of *N* coupled nodes *n* = 1*, …,N*, each with a population of *Mn* individuals, we use the SIR model to track the fraction of *Mn* who are susceptible (*sn*), infected (*jn*) or removed (*rn*), where the infected are divided as ![Graphic][58], among the treated (T) individuals, who have been provided a therapeutic and the untreated (U) individuals, who have not yet gained access to it. The epidemic dynamics is driven by (**Supp. Sec. 1**) ![Formula][59] where *α* is the infection rate, *β* is the mortality/recovery rate of the untreated individuals and ζ > *β* is the recovery rate under treatment. The epidemic reproduction number ℛ = *α/β* captures the level of contagion of the (untreated) disease. The therapeutic consumption rate *ρ*(*qn*) depends on the availability of the therapeutic *qn*(*t*) in *n*, as ![Formula][60] Hence *ρ*(*qn*) increases linearly with *qn*(*t*) as long as the demand (denominator) exceeds the supply, and saturates to unity when *n* has excess quantities of the therapeutic, avoiding over consumption (**Supp. Sec. 1.2**). The availability *qn*(*t*) is determined by our dissemination scheme, following either max-flow, as in Eq. (3), or egalitarian flow, as in Eq. (11). In (15) we introduce an invasion threshold *ε* through the sigmoidal function ![Formula][61] which activates the local SIR dynamics only when the local infection levels ![Graphic][62] exceed *ε*. The diffusion of individuals between nodes is mediated by *Anm*, derived from the empirical international air-travel network (**Supp. Sec. 2.1**). * Received September 15, 2020. * Revision received September 15, 2020. * Accepted September 18, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## REFERENCES 1. 1.Pastor-Satorras, R. & Vespignani, A. Epidemic Spreading in Scale-Free Networks. Phys. Rev. Lett. 86, 3200–3203 (2001). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/PhysRevLett.86.3200&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11290142&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000167866300072&link_type=ISI) 2. 2. R. Pastor-Satorras, C. Castellano, P. Van Mieghem and A. Vespignani. Epidemic processes in complex networks. Rev. Mod. Phys. 87, 925–958 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/RevModPhys.87.925&link_type=DOI) 3. 3.Barzel, B. & Barabási, A.-L. Universality in network dynamics. Nat. Phys. 9, 673–681 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nphys2741&link_type=DOI) 4. 4.Colizza, V. & Vespignani, A. Epidemic modeling in metapopulation systems with heterogeneous coupling pattern: Theory and simulations. J. Theor. Biol. 251, 450–467 (2008). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2007.11.028&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18222487&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000255082000006&link_type=ISI) 5. 5.Apolloni, A., Poletto, C., Ramasco, J. J., Jensen, P. & Colizza, V. Metapopulation epidemic models with heterogeneous mixing and travel behaviour. Theor. Biol. Med. Model. 11, 3 (2014). 6. 6.Tizzoni, M. et al. On the Use of Human Mobility Proxies for Modeling Epidemics. PLoS Comput. Biol. 10, 1003716 (2014). 7. 7.Rossman, H. et al. A framework for identifying regional outbreak and spread of COVID-19 from one-minute population-wide surveys. Nature Medicine 26, 634–638 (2020). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) 8. 8.Wu, F. et al. A new coronavirus associated with human respiratory disease in China. Nature 579, 265–269 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2008-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) 9. 9.Chinazzi, M. et al. The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak. Science 368, 395–400 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ4OS8zOTUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8xOC8yMDIwLjA5LjE1LjIwMTk0NTA2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 10. 10.Xu, B. et al. Epidemiological data from the COVID-19 outbreak, real-time case information. Sci. Data 7, 1–6 (2020). 11. 11.Zhang, J. et al. Evolving epidemiology and transmission dynamics of coronavirus disease 2019 outside Hubei province, China: a descriptive and modelling study. Lancet Infect. Dis. 20, 793–802 (2020). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) 12. 12.Dong, E., Du, H. & Gardner, L. An interactive web-based dashboard to track COVID-19 in real time. The Lancet Infectious Diseases 20, 533–534 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30120-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) 13. 13.Brockmann, D. & Helbing, D. The Hidden Geometry of Complex, Network-Driven Contagion Phenomena. Science 342, 1337–1342 (2013). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNDIvNjE2NC8xMzM3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMTgvMjAyMC4wOS4xNS4yMDE5NDUwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 14. 14.Sattenspiel, L. & Dietz, K. A structured epidemic model incorporating geographic-mobility among regions. Math. Biosci. 128, 71–91 (1995). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0025-5564(94)00068-B&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7606146&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995RC24500007&link_type=ISI) 15. 15.Hacohen, A., Cohen, R., Efroni, S., Barzel, B. & Bachelet, I. Digitizable therapeutics for decentralized mitigation of global pandemics. Sci. Rep. 9, 14345 (2019). 16. 16.Pastor-Satorras, R. & Vespignani, A. Immunization of complex networks. Phys. Rev. E 65, 036104 (2002). 17. 17.Oles, K., Gudowska-Nowak, E. & Kleczkowski, A. Understanding disease control: Influence of epidemiological and economic factors. PLoS One 7, 1–9 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0035495&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22506074&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) 18. 18.Arinaminpathy, N. & McLean, A. . Antiviral treatment for the control of pandemic influenza: some logistical constraints. J. R. Soc. Interface 5, 545–553 (2008). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rsif.2007.1152&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17725972&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) 19. 19.Woolley-Meza, O. et al. Complexity in human transportation networks: A comparative analysis of worldwide air transportation and global cargo-ship movements. European Physical Journal B 84, 589–600 (Springer, 2011). 20. 20.Ahuja, R. K., Magnanti, T. L. & Orlin, J. B. Network flows: theory, algorithms, and applications. (Prentice Hall, 1993). 21. 21.Dai, W., Zhang, J. & Sun, X. On solving multi-commodity flow problems: An experimental evaluation. Chinese J. Aeronaut. 30, 1481–1492 (2017). 22. 22.Barnhart, C., Krishnan, N. & Vance, P. H. Multicommodity Flow Problems. in Encyclopedia of Optimization 1583–1591 (Springer US, 2001). 23. 23.Fleischer, L. K. Faster algorithms for the quickest transshipment problem. SIAM J. Optim. 12, 18–35 (2001). 24. 24.Ford, L. R. & Fulkerson, D. R. Maximal flow through a network. in Classic Papers in Combinatorics 243–248 (Birkhauser Boston, 2009). 25. 25.Goldberg, A. V. & Tarjan, R. E. Efficient maximum flow algorithms. Commun. ACM 57, 82–89 (2014). 26. 26.Harush, U. & Barzel, B. Dynamic patterns of information flow in complex networks. Nat. Commun. 8, 2181 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-01916-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29259160&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) 27. 27.Hens, C., Harush, U., Haber, S., Cohen, R. & Barzel, B. Spatiotemporal signal propagation in complex networks. Nat. Phys. 15, 403–412 (2019). 28. 28.Brockmann, D., Hufnagel, L. & Geisel, T. The scaling laws of human travel. Nature 439, 462–465 (2006). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature04292&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16437114&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234859200044&link_type=ISI) 29. 29.Brockmann, D. & Helbing, D. The hidden geometry of complex, network-driven contagion phenomena. Science 342, 1337–1342 (2013). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNDIvNjE2NC8xMzM3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMTgvMjAyMC4wOS4xNS4yMDE5NDUwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 30. 30.Zhang, X., Ruan, Z., Zheng, M., Barzel, B. & Boccaletti, S. Epidemic spreading under infection-reduced-recovery. Chaos, Solitons and Fractals 140, 110130 (2020). 31. 31.Stegmaier, T., Oellingrath, E., Himmel, M. & Fraas, S. Differences in epidemic spread patterns of norovirus and influenza seasons of Germany: an application of optical flow analysis in epidemiology. Sci. Rep. 10, 14125 (2020). 32. 32.Buckee, C. O. et al. Aggregated mobility data could help fight COVID-19. Science 368, 145–146 (2020). [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE0OiIzNjgvNjQ4Ny8xNDUtYSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzE4LzIwMjAuMDkuMTUuMjAxOTQ1MDYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 33. 33.Barrat, A., Barthelemy, M. & Vespignani, A. Dynamical Processes on Complex Networks. (Cambridge University Press, 2008). 34. 34.Cohen, J. Designer antibodies could battle COVID-19 before vaccines arrive. Science. (2020). 35. 35.Régnier, S. A. & Huels, J. Drug versus vaccine investment: A modelled comparison of economic incentives. Cost Eff. Resour. Alloc. 11, 23 (2013). 36. 36.Cenciarelli, O. et al. Ebola virus disease 2013-2014 outbreak in west Africa: an analysis of the epidemic spread and response. Int. J. Microbiol. 2015, 769121 (2015). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25852754&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) 37. 37.WHO | Statement on the 1st meeting of the IHR Emergency Committee on the 2014 Ebola outbreak in West Africa. WHO (2014). Available at: [http://www.who.int/mediacentre/news/statements/2014/ebola-20140808/en/](http://www.who.int/mediacentre/news/statements/2014/ebola-20140808/en/). 38. 38. B. Barzel and A.-L. Barabási. Network link prediction by global silencing of indirect correlations. Nat. Biotechnol. 31, 720–725 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.2601&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23851447&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) 39. 39. B. Barzel and O. Biham. Quantifying the connectivity of a network: The network correlation function method. Phys. Rev. E 80, 46104–46115 (2009). 40. 40.Emanuel, E. J. et al. Fair allocation of scarce medical resources in the time of Covid-19. N. Engl. J. Med. 382, 2049–2055 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMsb2005114&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32202722&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) 41. 41.Laventhal, N. et al. The ethics of creating a resource allocation strategy during the COVID-19 Pandemic. Pediatrics 146, (2020). 42. 42.Ford, L. R. & Fulkerson, D. R. Constructing maximal dynamic flows from static flows. Oper. Res. 6, 419–433 (1958). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1287/opre.6.3.419&link_type=DOI) 43. 43.Fleischer, L. & Skutella, M. Quickest flows over time. SIAM J. Comput. 36, 1600–1630 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1137/S0097539703427215&link_type=DOI) 44. 44.Aronson, J. E. A survey of dynamic network flows. Ann. Oper. Res. 20, 1–66 (1989). 45. 45. R. Cohen and S. Havlin. Scale-free networks are ultrasmall. Phys. Rev. Lett. 90, 58701–58704 (2003). 46. 46.Steinegger, B., Arenas, A., Gomez-Gardenes, J. & Granell, C. Pulsating campaigns of human prophylaxis driven by risk perception palliate oscillations of direct contact transmitted diseases. Phys. Rev. Res. 2, 23181 (2020). 47. 47.Soriano-Paños, D., Ghoshal, G., Arenas, A. & Gómez-Gardeñes, J. Impact of temporal scales and recurrent mobility patterns on the unfolding of epidemics. J. Stat. Mech. Theory Exp. 2020, (2019). 48. 48.Mancastroppa, M., Burioni, R., Colizza, V. & Vezzani, A. Active and inactive quarantine in epidemic spreading on adaptive activity-driven networks. Phys. Rev. E 102, 020301 (2020). 49. 49.Abouk, R. & Heydari, B. The Immediate Effect of COVID-19 Policies on Social Distancing Behavior in the United States. SSRN Electron. J. (2020). 50. 50.Rader, B. et al. Crowding and the epidemic intensity of COVID-19 transmission. *medRxiv* 2020.04.15.20064980 (2020). 51. 51.Bianconi, G., Sun, H., Rapisardi, G. & Arenas, A. A message-passing approach to epidemic tracing and mitigation with apps. *arXiv:2007.05277* (2020). 52. 52.Meidan, D. et al. Alternating quarantine for sustainable mitigation of COVID-19. *arXiv2004.01453* (2020). 53. 53.Tizzoni, M. et al. Real-time numerical forecast of global epidemic spreading: case study of 2009 A/H1N1pdm. BMC Med. 10, 165 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1741-7015-10-165&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23237460&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20194506.atom) [1]: /embed/graphic-2.gif [2]: /embed/inline-graphic-1.gif [3]: /embed/graphic-4.gif [4]: /embed/graphic-5.gif [5]: /embed/graphic-6.gif [6]: /embed/graphic-7.gif [7]: /embed/inline-graphic-2.gif [8]: /embed/graphic-8.gif [9]: /embed/graphic-9.gif [10]: /embed/inline-graphic-3.gif [11]: /embed/graphic-10.gif [12]: /embed/inline-graphic-4.gif [13]: /embed/inline-graphic-5.gif [14]: /embed/inline-graphic-6.gif [15]: F3/embed/inline-graphic-7.gif [16]: F3/embed/inline-graphic-8.gif [17]: F3/embed/inline-graphic-9.gif [18]: F3/embed/inline-graphic-10.gif [19]: F3/embed/inline-graphic-11.gif [20]: F3/embed/inline-graphic-12.gif [21]: F3/embed/inline-graphic-13.gif [22]: F3/embed/inline-graphic-14.gif [23]: F3/embed/inline-graphic-15.gif [24]: F3/embed/inline-graphic-16.gif [25]: /embed/inline-graphic-17.gif [26]: /embed/graphic-12.gif [27]: /embed/inline-graphic-18.gif [28]: /embed/inline-graphic-19.gif [29]: /embed/inline-graphic-20.gif [30]: /embed/inline-graphic-21.gif [31]: /embed/inline-graphic-22.gif [32]: /embed/inline-graphic-23.gif [33]: /embed/inline-graphic-24.gif [34]: /embed/inline-graphic-25.gif [35]: /embed/inline-graphic-26.gif [36]: /embed/inline-graphic-27.gif [37]: /embed/inline-graphic-28.gif [38]: /embed/inline-graphic-29.gif [39]: /embed/inline-graphic-30.gif [40]: /embed/inline-graphic-31.gif [41]: /embed/inline-graphic-32.gif [42]: /embed/graphic-13.gif [43]: /embed/inline-graphic-33.gif [44]: /embed/inline-graphic-34.gif [45]: /embed/inline-graphic-35.gif [46]: /embed/inline-graphic-36.gif [47]: /embed/inline-graphic-37.gif [48]: /embed/inline-graphic-38.gif [49]: /embed/graphic-14.gif [50]: /embed/graphic-15.gif [51]: /embed/graphic-16.gif [52]: /embed/inline-graphic-39.gif [53]: /embed/inline-graphic-40.gif [54]: F4/embed/inline-graphic-41.gif [55]: F4/embed/inline-graphic-42.gif [56]: /embed/inline-graphic-43.gif [57]: /embed/graphic-18.gif [58]: /embed/inline-graphic-44.gif [59]: /embed/graphic-19.gif [60]: /embed/graphic-20.gif [61]: /embed/graphic-21.gif [62]: /embed/inline-graphic-45.gif