Eisenbrand et al [1] proposed an interesting technique called cycle connection game. In the setting, there is a set of nn core points connected by an undirected cycle. Each core point is associated with a client. There is a pair of opposite directed edges between each core point and its associated client. Each edge has a weight. We denote GG be the resulting graph.

An example is shown as follows:

Given an instance, we would like to sample a set of clients randomly (we call the sampled clients the "marked" ones). Then each client sends a unit capacity flow to its nearest marked client (could be itself). The cost of each flow is the sum of edges weights it travels. We would like to bound the expected cost of sending the flows.

The selection process requires a probability pp as input and consists of two steps:

  1. Pick a client uniformly at random. This step ensures that at least one client is selected.
  2. For the rest clients, pick each of them independently with some probability pp.

To analyze the expected cost, we divide the edges into three categories:

  1. Directed edges from the clients to the core points, denoted as Ein\mathcal{E}_{in};
  2. Directed edges from the core points to the clients, denoted as Eout\mathcal{E}_{out};
  3. Edges on the cycle, denoted as C\mathcal{C}.

The following theorem holds:

Theorem. The expected cost is bounded by

∑e∈Ein∪Eoutwe+12p∑e∈Cwe \sum_{e \in \mathcal{E}_{in} \cup \mathcal{E}_{out} } w_e + \frac{1}{2p} \sum_{e \in \mathcal{C} } w_e

Proof.
For each edge ee, define a random variable XeX_e to be the flow value of ee. For e∈Eine \in \mathcal{E}_{in}, it holds that Xe∈{0,1}X_e \in \{0, 1\}. However, for e∈Eout∪Ce \in \mathcal{E}_{out} \cup \mathcal{C}, it is possible that Xe>1X_e > 1. By linearity of expectation,

E[∑e∈Ein∪Eoutwe⋅Xe+∑e∈Cwe⋅Xe]=∑e∈Ein∪Eoutwe⋅E[Xe]+∑e∈Cwe⋅E[⋅Xe]. \mathbb{E} \left[ \sum_{e \in \mathcal{E}_{in} \cup \mathcal{E}_{out} } w_e \cdot X_e + \sum_{e \in \mathcal{C} } w_e \cdot X_e \right] = \sum_{e \in \mathcal{E}_{in} \cup \mathcal{E}_{out} } w_e \cdot \mathbb{E} \left[ X_e \right] + \sum_{e \in \mathcal{C} } w_e \cdot \mathbb{E} \left[ \cdot X_e \right].
  1. For e∈Eine \in \mathcal{E}_{in}, there is a unit flow on ee if only the associated client is not sampled. This happens with probability (1−1n)(1−p)(1 - \frac{1}{n})(1 - p), where nn is the number of clients. Therefore, E[Xe⋅we]=(1−1n)(1−p)⋅we\mathbb{E} [X_e \cdot w_e] = (1 - \frac{1}{n})(1 - p) \cdot w_e and

    E[∑e∈EinXe⋅we]=(1−1n)(1−p)⋅∑e∈Einwe. \mathbb{E} \left[ \sum_{e \in \mathcal{E}_{in} } X_e \cdot w_e \right] = (1 - \frac{1}{n})(1 - p) \cdot \sum_{ e \in \mathcal{E}_{in} } w_e.
  2. For e∈Eoute \in \mathcal{E}_{out}, the analysis is a little tricky. For each particular sample instance, it holds that ∑e∈EoutXe=∑e∈EinXe\sum_{e \in \mathcal{E}_{out}} X_e = \sum_{e \in \mathcal{E}_{in}} X_e, by flow conservation. Therefore,

    E[∑e∈EoutXe]=E[∑e∈EinXe]=(1−1n)(1−p)n. \mathbb{E} \left[ \sum_{e \in \mathcal{E}_{out}} X_e \right] = \mathbb{E} \left[ \sum_{e \in \mathcal{E}_{in}} X_e \right] = (1 - \frac{1}{n})(1 - p) n.

    By symmetry, {Xe}e∈Eout\{ X_e \}_{ e \in \mathcal{E}_{out} } are identical distributed random variables. Therefore E[Xe]=(1−1n)(1−p)\mathbb{E} [X_e] = (1 - \frac{1}{n})(1 - p). It concludes that

    E[∑e∈EoutXe⋅we]=(1−1n)(1−p)⋅∑e∈Eoutwe. \mathbb{E} \left[ \sum_{e \in \mathcal{E}_{out} } X_e \cdot w_e \right] = (1 - \frac{1}{n})(1 - p) \cdot \sum_{ e \in \mathcal{E}_{out} } w_e.
  3. For e∈Ce \in \mathcal{C}, the analysis of E[Xe]\mathbb{E} [X_e] is the hardest part. We look at the suboptimal flow routing scheme such that each un-marked client sends a unit flow to a nearest marked client, with respect to unit cycle edge weight. Let YeY_e be the number of flows on edge ee. Clearly,

    E[∑e∈CXe⋅we]≤E[∑e∈CYe⋅we]. \mathbb{E} \left[ \sum_{e \in \mathcal{C} } X_e \cdot w_e \right] \le \mathbb{E} \left[ \sum_{e \in \mathcal{C} } Y_e \cdot w_e \right].

    Further, for each client jj, let the random variable YjY_j be the number of cycle edges travelled by its flow, then

    ∑e∈CYe=∑jYj. \sum_{e \in \mathcal{C} } Y_e = \sum_j Y_j.

    For a fixed jj and some integer 0≤k≤⌊(n−2)/2⌋0 \le k \le \lfloor (n - 2) / 2 \rfloor, the event Yj>kY_j > k happens if and only if jj is not marked and none of the first kk client neighbors to the left and right of jj is marked, the probability of which is given by

    Pr⁡[Yj>k]=(1−2k+1n)q2k+1. \Pr[Y_j > k] = \left( 1 - \frac{2k + 1}{n} \right) q^{2k + 1}.

    where q=(1−p)q = (1 - p). Now the expectation of YjY_j is given by

    E[Yj]=∑k=0⌊(n−2)/2⌋(1−2k+1n)q2k+1≤e−1nq∑k=0⌊(n−2)/2⌋e−k2nqk2≤e−1nq∑k=0⌊(n−2)/2⌋(e−2nq2)k=e−1nq(1−(e−2nq2)⌊(n−2)/2⌋)1−(e−2nq2)≤e−1nq1−(e−2nq2)≤12(1−e−1nq). \begin{aligned} \mathbb{E} [Y_j] &= \sum_{k = 0}^{\lfloor (n - 2) / 2 \rfloor } \left( 1 - \frac{2k + 1}{n} \right) q^{2k + 1} \\ &\le e^{ - \frac{ 1 }{n} } q \sum_{k = 0}^{\lfloor (n - 2) / 2 \rfloor } e^{ - k \frac{ 2 }{n} } q^{k2} \\ &\le e^{ - \frac{ 1 }{n} } q \sum_{k = 0}^{\lfloor (n - 2) / 2 \rfloor } (e^{ - \frac{ 2 }{n} } q^{2} )^k \\ &= e^{ - \frac{ 1 }{n} } q \frac{ (1 - (e^{ - \frac{ 2 }{n} } q^{2} )^{\lfloor (n - 2) / 2 \rfloor } ) }{1 - (e^{ - \frac{ 2 }{n} } q^{2} ) } \\ &\le \frac{ e^{ - \frac{ 1 }{n} } q }{ 1 - (e^{ - \frac{ 2 }{n} } q^{2} ) } \\ &\le \frac{1}{ 2 (1 - e^{ - \frac{ 1 }{n} } q) }. \end{aligned}

    The final inequality holds as 2x−2x2≤1−x22x - 2x^2 \le 1 - x^2. Replacing qq with 1−p1 - p, we have

    E[Yj]≤12(1−e−1n(1−p))≤12(1−e−1n+e−1np)≤12p. \mathbb{E} [Y_j] \le \frac{1}{ 2 (1 - e^{ - \frac{ 1 }{n} } (1 - p)) } \le \frac{1}{ 2 (1 - e^{ - \frac{ 1 }{n} } + e^{ - \frac{ 1 }{n} } p) } \le \frac{1}{2p}.

    Remark: a more coarse calculation which leads to the same result would be

    E[Yj]≤∑k=0∞(1−p)2k+1≤1−p1−(1−p)2≤12p. \mathbb{E} [Y_j] \le \sum_{k = 0}^\infty (1 - p)^{2k + 1} \le \frac{1 - p}{1 - (1 - p)^2 } \le \frac{1}{2p}.

    By symmetry, the YjY_j's are identical distributed random variables. So are the variable YeY_e's for e∈Ce \in \mathcal{C}. It follows

    E[Ye]=E[Yj]≤12p. \begin{aligned} \mathbb{E} [Y_e] = \mathbb{E} [Y_j] &\le \frac{1}{2p}. \end{aligned}

    and

    E[∑e∈CYe⋅we]≤12p∑e∈Cwe. \begin{aligned} \mathbb{E} \left[ \sum_{e \in \mathcal{C}} Y_e \cdot w_e \right] &\le \frac{1}{2p} \sum_{e \in \mathcal{C}} w_e. \end{aligned}

■\blacksquare

Reference:

  1. Eisenbrand, Friedrich, Fabrizio Grandoni, Thomas Rothvoß, and Guido Schäfer. "Approximating connected facility location problems via random facility sampling and core detouring." In Proceeding of Nineteenth annual ACM-SIAM Symposium (SODA'08), no. CONF, pp. 1174-1183. 2008.

  2. Eisenbrand, F., Grandoni, F., Rothvoß, T., Schäfer, G., 2010. Connected facility location via random facility sampling and core detouring. Journal of Computer and System Sciences 76, 709–726.