Motivation
When social-distancing rules cap how many members of an organization may be physically present, a school, office or plant of 2K members can host only K at a time. A common response is rotation: split the members into two fixed groups and alternate which group is on site. How often to alternate is a design choice, and it matters because an infection that enters one group spreads only while that group is active, and is detected only after a delay, when symptoms appear or tests are run.
Ely, Galeotti and Steiner (2021) model this choice and show that the optimal rotation is extreme: either the groups alternate as often as possible, or one group is kept active throughout. Which extreme is optimal is decided by a single comparison between the organization's reaction time and the time a logistic epidemic needs to infect half a group. The model is a small continuous-time control problem with a switching activity coefficient, a logistic state equation and an expectation over the infection date; its main theorem is an exact optimality statement rather than an approximation.
Setting
Two groups i∈{0,1} each have measure K>1. A rotation schedule is a pair of measurable functions r(0,⋅),r(1,⋅) with values in [0,1] and r(0,t)+r(1,t)=1 for every time t; r(i,t) is the share of time group i is active around t. The frequent rotation is r(i,t)=1/2 for both groups and all t (the limit of alternating very often); activating group i throughout is r(i,t)=1, r(1−i,t)=0.
At a random time t∗≥0 with an arbitrary law μ, one member of a group ω becomes infected; ω is 0 or 1 with probability one half each, independently of t∗. The prevalence pt of the infected group starts at pt∗=1/K and follows the logistic equation
p˙t=γr(ω,t)pt(1−pt),t≥t∗,(2)
with transmission rate γ>0; the other group stays uninfected. The organization reacts at time t∗+R, where R≥0 is the reaction delay, and the rotation problem is
rmin Ept∗+Rω(ω),(P1)
the expectation being over t∗ and ω.
Two derived objects organize the analysis. The exposure time τ=∫t∗t∗+Rr(ω,t)dt is how long the infected group was active before the reaction. The loss and the half-infection time are
ℓ(t)=1+eγ(t−b)eγ(t−b),b=γ1ln(K−1);
ℓ(t) is the prevalence after t units of uninterrupted spread, and b is the time at which it reaches one half.
Formalization targets
Goal: Proposition 1 (p. 10)
2R<b⟹r(0,t)=r(1,t)=21 solves (P1),2R>b⟹r(i,t)=1 for all t solves (P1), i∈{0,1}.
Optimality is over all rotation schedules and holds for every law of t∗ on [0,∞). The boundary case R/2=b is not part of the statement.
Milestones, in the paper's order
- ℓ solves ℓ˙=γℓ(1−ℓ) with ℓ(0)=1/K (§2.1, p. 7).
- ℓ(b)=1/2; ℓ is convex below b and concave above (§2.1, p. 7).
- Lemma 1 (p. 7): the prevalence at the reaction time is pt∗+Rω(ω)=ℓ(τ).
- Lemma 2 (p. 8, in corrected form): Eτ=R/2 and Pr(τ≤T)=Pr(τ≥R−T).
- The exposure times of the two groups sum to R, and Eℓ(τ)=21E[ℓ(τ)+ℓ(R−τ)] (p. 8).
- Lemma 3 (pp. 9–10): ℓ(τ)+ℓ(R−τ) is U-shaped around R/2 if R/2<b, hump-shaped if R/2>b.
- The relaxed problem (P2), minFEF[ℓ(τ)+ℓ(R−τ)] subject to EFτ=R/2 over laws on [0,R], is solved by the point mass at R/2 when R/2<b and by 21δ0+21δR when R/2>b (proof of Proposition 1, pp. 10–11).
- The frequent rotation gives expected terminal prevalence ℓ(R/2), and activating one group gives (ℓ(0)+ℓ(R))/2 (p. 11).
Significance
The result turns a control problem over all measurable switching schedules into a comparison of two numbers, R/2 and b. It says that intermediate rotation frequencies are never strictly better than both extremes, and it identifies the curvature of the logistic curve over the reaction window as the quantity that decides between them: frequent rotation equalizes the two groups' exposure and so insures against exposure-time risk, which is worth doing only while prevalence is still in its convex phase. The paper's Corollary 1 (comparative statics in R, γ and K) follows from it.
The proposition is proved in the paper, in a few lines that rely on informal steps: solving (2) by a change of variable, the symmetry of the exposure-time distribution, and the "standard concavification argument" for (P2). To our knowledge none of it has been machine-checked. A formalization makes each step precise for measurable schedules, for which (2) holds only almost everywhere, and records one correction: the symmetry identity of Lemma 2 as printed, F(τ)=1−F(R−τ), is false at atoms of the distribution, and the mission states the form the proof actually establishes.
Difficulty
The obvious argument treats (2) as a classical differential equation and substitutes t^=∫r. For a merely measurable schedule this needs absolute continuity, a chain rule for a Lipschitz time change, and uniqueness of solutions of a Carathéodory equation; Lemma 1 has to supply all three, including that every solution stays the composite ℓ(∫r). The second difficulty is that (P1) optimizes over schedules, not over distributions: the reduction to (P2) uses only the mean and the symmetry of τ, and the optimum of (P2) must then be shown to be attained by an actual schedule. Lemma 3 holds on all of R, but the solution of (P2) in the hump-shaped case uses that exposure times lie in [0,R]: off that interval ℓ(τ)+ℓ(R−τ) falls below its endpoint values.
Formalization scope
All quantities are real numbers. The parameters γ>0, K>1, R≥0 are bundled as standing assumptions; R≥0 is added (a delay). Groups are Fin 2. Schedules are measurable functions on all of R satisfying r(i,t)∈[0,1] and r(0,t)+r(1,t)=1 at every t; the page's domain is R+, and the two readings agree for (P1) because t∗≥0. The law μ of t∗ is any probability measure on R with μ((−∞,0))=0, and the objective is 21∑ω∫pt∗+Rω,t∗dμ(t∗).
Equation (2) is used in integral form: a prevalence path is continuous and [0,1]-valued on [t∗,∞) and satisfies pt=1/K+∫t∗tγr(ω,s)ps(1−ps)ds; this admits the paper's periodic schedules, which have no classical derivative. A schedule is optimal when it has a prevalence family (one path for every ω and every t∗≥0) and, for every such family, its objective is at most that of any schedule with any family; the existence clause keeps optimality from holding vacuously. Terminal prevalence is the value of the solution of (2), never defined as ℓ(τ): that identification is Lemma 1 and is part of the work. Lemma 2 is stated as Pr(τ≤T)=Pr(τ≥R−T) for all real T, not as the printed distribution-function identity.
A complete development needs: the logistic function and its derivative; FTC and change of variables for absolutely continuous functions; a Grönwall-type uniqueness argument for the integral equation; and an elementary concavification bound for a function that is monotone on each side of R/2. The logistic and uniqueness parts are reusable beyond this mission. Proofs of any milestone, and alternative routes to the goal, are welcome.
Not formalized here: Corollary 1 (p. 11), the numerical calibration remark of p. 11, and §4 on mixing between groups (Lemma 4, Propositions 2 and 3), which is a different model.
Selected references