arXiv is now an independent nonprofit! Learn more
License: arXiv.org perpetual non-exclusive license
arXiv:1906.00074v1 [cs.SI] 31 May 2019

Balancing spreads of influence in a social network

Ruben Becker    Federico Corò    Gianlorenzo D’Angelo    Hugo Gilbert
Gran Sasso Science Institute (GSSI), L’Aquila, Italy
Abstract

The personalization of our news consumption on social media has a tendency to reinforce our pre-existing beliefs instead of balancing our opinions. This finding is a concern for the health of our democracies which rely on an access to information providing diverse viewpoints. To tackle this issue from a computational perspective, Garimella et al. (NIPS’17) modeled the spread of these viewpoints, also called campaigns, using the well-known independent cascade model introduced by Kempe et al. (KDD’03) and studied an optimization problem that aims at balancing information exposure in a social network when two opposing campaigns propagate in the network. The objective in their 𝑁𝑃\mathit{NP}-hard optimization problem is to maximize the number of people that are exposed to either both or none of the viewpoints. For two different settings, one corresponding to a model where campaigns spread in a correlated manner, and a second one, where the two campaigns spread in a heterogeneous manner, Garimella et al. provide constant ratio approximation algorithms. In this paper, we investigate a more general formulation of this problem. That is, we assume that μ\mu different campaigns propagate in a social network and we aim to maximize the number of people that are exposed to either ν\nu or none of the campaigns, where μν2\mu\geq\nu\geq 2. We provide dedicated approximation algorithms for both the correlated and heterogeneous settings. Interestingly, while the problem can still be approximated within a constant factor in the correlated setting for any constant values of μ\mu and ν\nu, for the heterogeneous setting with ν3\nu\geq 3, we give reductions leading to several approximation hardness results. Maybe most importantly, we obtain that the problem cannot be approximated within a factor of ng(n)n^{-g(n)} for any g(n)=o(1)g(n)=o(1) assuming the Gap-ETH hypothesis, denoting with nn the number of nodes in the social network. For ν4\nu\geq 4, we furthermore show a stronger hardness of approximation bound under a different condition, that is, if a certain class of one-way functions exists, then there is no nϵn^{-\epsilon}-approximation algorithm where ϵ>0\epsilon>0 is a given constant which depends on ν\nu. This complements our finding of an approximation algorithm for the heterogeneous case that for arbitrary μ\mu and ν=3\nu=3 leads to an approximation ratio of order n1/2n^{-1/2}.

1 Introduction

One of the promises of a highly connected world is that of an impartial spread of opinions driven by free and unbiased sources of information. As a consequence, any opinion could have been equitably exposed to the wide public. On the contrary, the social network platforms that are currently governing news diffusion, while offering many seemingly-desirable features like searching, personalization, and recommendation, are reinforcing the centralization of information spreading and the creation of what is often termed echo chambers and filter bubbles [GDFMGM18]. Stated differently, algorithmic personalization of news diffusion are likely to create homogeneous polarized clusters where users get less exposure to conflicting viewpoints. A good illustration of this issue was given by Conover et al. [CRF+11] who studied the Twitter network during the 2010 US congressional midterm elections. The authors demonstrated that the retweet network had a highly segregated partisan structure with extremely limited connectivity between left-wing and right-wing users. Consequently, instead of giving users a diverse perspective and balancing users opinions by exposing them to challenging ideas, social media platforms are likely to make users more extreme by only exposing them to views that reinforce their pre-existing beliefs [CRF+11, DVBZ+16].

To address this issue from a computational perspective, Garimella et al. [GGPT17] introduced the problem of balancing information exposure in a social network. Following the influence maximization paradigm going back to the seminal work of Kempe et al. [KKT03, KKT15], their problem involves two opposing viewpoints or campaigns that propagate in a social network following the independent cascade model. Given initial seed sets for both campaigns, a centralized agent is then responsible for selecting a small number of additional seed users for each campaign in order to maximize the number of users that are reached by either both or none of the campaigns. The authors study this problem in two different settings, namely the heterogeneous and correlated settings. The heterogeneous setting corresponds to the general case in which there is no restriction on the probabilities with which the campaigns propagate. Contrarily, in the correlated setting, the probability distributions for different campaigns are identical and completely correlated. After proving that the optimization problem of balancing information exposure is 𝑁𝑃\mathit{NP}-hard, the authors designed efficient approximation algorithms with an approximation ratio of (11/eϵ)/2(1-1/e-\epsilon)/2 for both settings.

Our Contribution.

We address the main open problem in [GGPT17], that is we generalize their optimization problem to a setting with possibly more than two campaigns. More precisely, let μ\mu and ν\nu be fixed constants such that 2νμ2\leq\nu\leq\mu. In our general problem, there are μ\mu opposing campaigns and the task is to maximize the number of nodes in the network that are reached by at least ν\nu campaigns or remain oblivious to all of them. We term this problem the μ\mu-ν\nu-Balance problem. Interestingly, we obtain results that surprisingly differ from the ones of Garimella et al. [GGPT17]. Indeed, while we show that any μ\mu-ν\nu-Balance problem can be approximated within a constant factor in the correlated setting (Section 5), we obtain strong approximation hardness results in the heterogeneous setting. In particular, when ν3\nu\geq 3, we show that under the Gap Exponential Time Hypothesis [Man17], there is no ng(n)n^{-g(n)}-approximation algorithm with g(n)=o(1)g(n)=o(1) for the μ\mu-ν\nu-Balance problem where nn is the number of nodes. Moreover, when ν4\nu\geq 4, we show that if a certain class of one-way functions exists [App13], there is no nϵn^{-\epsilon}-approximation algorithm for the μ\mu-ν\nu-Balance problem where ϵ>0\epsilon>0 is a constant which depends on ν\nu (Section 3). We mitigate these hardness results by designing an algorithm with an approximation factor of Ω(n1/2)\Omega(n^{-1/2}) for the case where ν=3\nu=3 and μ\mu is an arbitrary constant (Section 4).

Related work.

There is a large literature on influence maximization, we refer the interested reader to [BBCL14, KKT15] and references therein. Here we focus on the literature about multiple campaigns running simultaneously on the same network. Budak et al. [BAEA11] studied the problem of limiting as much as possible the spread of a “bad” campaign by starting the spreading of another “good” campaign that blocks the first one. The two campaigns compete on the nodes that they reach: once a node becomes active in one campaign it cannot change campaign. They prove that the objective function is monotone and submodular and hence they obtain a constant approximation ratio. Similar concepts of competing cascades in which a node can only participate in one campaign have been studied in several works [AM11, BKS07, CNWVZ07, DGDM06, KOW08, LCL15, ML12]. Game theoretic aspects like the existence of Nash equilibria have been also investigated in this case [AFPT10, GHK14, TAM12]. Borodin et al. [BBLO17] consider the problem of controlling the spread of multiple campaigns by a centralized authority. Each campaign has its own objective function to maximize associated with its spread and the aim of a central authority is to maximize the social welfare defined as the sum of the selfish objective function of each campaign. They propose a truthful mechanism to achieve theoretical guarantees on the social welfare.

Two other works closely related to ours are the ones of Aslay et al. [AMGG18] and Matakos et al. [MG18]. The former work tackles an item-aware information propagation problem in which a centralized agent must recommend some articles to a small set of seed users such that the spread of these articles maximizes the expected diversity of exposure of the agents. The diversity exposure is measured by a sum of agent-dependent functions that takes into account user leanings. The authors show that the 𝑁𝑃\mathit{NP}-hard problem they define amounts to optimizing a monotone and submodular function under a matroid constraint and design a constant factor approximation algorithm. The latter paper models the problem of maximizing the diversity of exposure in a social network as a quadratic knapsack problem. Here also the problem amounts to recommending a set of articles to some users in order to maximize a diversity index taking into account users’ leanings and the strength of their connections in the social network. The authors show that the resulting diversity maximization problem is inapproximable and design a polynomial algorithm without an approximation guarantee.

2 Preliminaries

2.1 Independent Cascade model

We introduce the well-known Independent Cascade model. We mostly follow the terminology and notation from Kempe et al. [KKT15]. Given a directed graph G=(V,E)G=(V,E), probabilities p:E[0,1]p:E\rightarrow[0,1] and an initial node set AVA\subseteq V called a set of seed nodes. Define A0=AA_{0}=A. For t0t\geq 0, we call a node vAtv\in A_{t} active at time tt. If a node vv is active at time t0t\geq 0 but was not active at time t1t-1, i.e., vAtAt1v\in A_{t}\setminus A_{t-1} (formally let A1=A_{-1}=\emptyset), it tries to activate each neighbor ww, independently, with a probability of success equal to pvwp_{vw}. In case of success ww becomes active at step t+1t+1, i.e., wAt+1w\in A_{t+1}. If at some time t0t^{*}\geq 0, we have that At=At+1A_{t^{*}}=A_{t^{*}+1} we say that the process has quiesced and call tt^{*} the time of quiescence. For an initial set AA, we denote with σ(A)=E[|At|]\sigma(A)=\E[|A_{t^{*}}|] the expected number of nodes activated at the time of quiescence when running the process with seed nodes AA. Kempe et al. showed that this process is equivalent to what is referred to as the Triggering Model, see [KKT15, Proof of Theorem 4.5]. For a node vVv\in V, let NvN_{v} denote all in-neighbors of vv. Here, every node independently picks a triggering set TvNvT_{v}\subseteq N_{v} according to a distribution over subsets of its in-neighbors, namely Tv=ST_{v}=S with probability uSpuvuNvS(1puv)\prod_{u\in S}p_{uv}\cdot\prod_{u\in N_{v}\setminus S}(1-p_{uv}). For a possible outcome X=(Tv)vVX=(T_{v})_{v\in V} of triggering sets for the nodes VV, let ρX(A)\rho_{X}(A) be the set of nodes reachable from AA in the outcome XX. Note that after sampling XX, the quantity ρX(A)\rho_{X}(A) is deterministic. According to Kempe et al. [KKT15], this model is equivalent to the Independent Cascade model and it holds that σ(A)=EX[|ρX(A)|]\sigma(A)=\E_{X}[|\rho_{X}(A)|], where the expectation is over the outcome profile XX. While it is not feasible to compute ρX(A)\rho_{X}(A) for all outcome profiles XX, it is possible to obtain a (1±ϵ)(1\pm\epsilon)-approximation to σ(A)\sigma(A), with probability at least 1δ1-\delta, by sampling Ω(|V|2log(1/δ)/ϵ2)\Omega(|V|^{2}\log(1/\delta)/\epsilon^{2}) possible outcomes XX and computing the average over the corresponding values |ρX(A)||\rho_{X}(A)|, see [KKT15, Proposition 4.1].

2.2 The μ\mu-ν\nu-Balance problem

Inspired by the work of Garimella et al. [GGPT17], we consider several information spread processes, we also call them “campaigns”, unfolding in parallel, each following the Independent Cascade model described above. Formally, we are given a graph G=(V,E)G=(V,E) and μ\mu probability functions (pi)i[μ](p_{i})_{i\in[\mu]}, where each pip_{i} is a probability function as in the Independent Cascade model described above, i.e., pi:E[0,1]p_{i}:E\rightarrow[0,1].11 1 For nn\in\mathbb{N}, we use [n][n] to denote the set {1,,n}\{1,\ldots,n\}. For an index i[μ]i\in[\mu], let Xi=(Tv)vVX_{i}=(T_{v})_{v\in V} be a possible outcome sampled using probabilities pip_{i}. Then for a seed set AVA\subseteq V, we denote with ρXi(i)(A)\rho^{(i)}_{X_{i}}(A) the set of nodes reachable from AA in outcome XiX_{i}. For an arbitrary sequence =(Ri)i[μ]\mathcal{R}=(R_{i})_{i\in[\mu]} of subsets of VV, we define

NoSMμ,ν():=|(Vi[μ]Ri)M[μ]:|M|νiMRi|\textstyle\NoSM(\mathcal{R}):=\big|(V\setminus\bigcup_{i\in[\mu]}R_{i})\cup\bigcup_{M\subseteq[\mu]:|M|\geq\nu}\bigcap_{i\in M}R_{i}\big|

to be the number of nodes that are contained in none or in sufficiently many, i.e., in at least ν\nu, of the sets in \mathcal{R}. Let 𝒳=(Xi)i[μ]\mathcal{X}=(X_{i})_{i\in[\mu]} be an outcome profile by letting XiX_{i} be a possible outcome according to distribution pip_{i}. Then, for 𝒜=(Ai)i[μ]\mathcal{A}=(A_{i})_{i\in[\mu]} with AiVA_{i}\subseteq V, we denote with ρ𝒳(𝒜)=(ρXi(i)(Ai))i[μ]\rho_{\mathcal{X}}(\mathcal{A})=(\rho^{(i)}_{X_{i}}(A_{i}))_{i\in[\mu]} the set of reached nodes in the outcome 𝒳\mathcal{X} from seed sets 𝒜\mathcal{A}. For two sequences of sets 𝒜\mathcal{A}, 𝒜\mathcal{A}^{\prime}, and a set AA, we let 𝒜𝒜=(AiAi)i[μ]\mathcal{A}\cup\mathcal{A}^{\prime}=(A_{i}\cup A^{\prime}_{i})_{i\in[\mu]} be the element-wise union and 𝒜A=(AiA)i[μ]\mathcal{A}\cap A=(A_{i}\cap A)_{i\in[\mu]} be the element-wise intersection with the set AA.

For constant integers μν2\mu\geq\nu\geq 2, we consider the following optimization problem:

μ\mu-ν\nu-Balance Input: Graph G=(V,E)G=(V,E), probabilities 𝒫=(pi)i[μ]\mathcal{P}=(p_{i})_{i\in[\mu]}, seed sets =(Ii)i[μ]\mathcal{I}=(I_{i})_{i\in[\mu]}, and k2k\geq 2. Find: sets 𝒮=(Si)i[μ]\mathcal{S}=(S_{i})_{i\in[\mu]} with i[μ]|Si|k\sum_{i\in[\mu]}|S_{i}|\leq k, such that Φμ,ν(𝒮)\Phi^{\mathcal{I}}_{\mu,\nu}(\mathcal{S}) is maximum, where Φμ,ν(𝒮):=E𝒳[NoSMμ,ν(ρ𝒳(𝒮))].\Phi^{\mathcal{I}}_{\mu,\nu}(\mathcal{S}):=\E_{\mathcal{X}}[\NoSM(\rho_{\mathcal{X}}(\mathcal{I}\cup\mathcal{S}))].

We refer to the objective function simply by Φ(𝒮)\Phi(\mathcal{S}), in case \mathcal{I}, μ\mu, and ν\nu are clear from the context. We assume kν|V|k\leq\nu|V| as otherwise the problem becomes trivial by choosing Si=VS_{i}=V for every i[ν]i\in[\nu]. Moreover, we assume w.l.o.g. that |V|μ|V|\geq\mu and kνk\geq\nu, since |V||V| and kk are input parameters and μ\mu and ν\nu are constant numbers. Following Garimella et al. [GGPT17], we distinguish two settings. (1) The heterogeneous setting corresponds to the general case in which there is no restriction on 𝒫\mathcal{P}. (2) In the correlated setting, the distributions pip_{i} are identical and completely correlated for all i[μ]i\in[\mu]. That is, if an edge (u,v)(u,v) propagates a campaign to vv, it propagates all campaigns that reach uu to vv.

Decomposing the Objective Function.

In all of our algorithms, we use the approach of decomposing the objective function into summands and approximating the summands separately. For an outcome profile 𝒳\mathcal{X}, and seed sets =(Ii)i[μ]\mathcal{I}=(I_{i})_{i\in[\mu]}, we define V𝒳,VV^{\ell,\mathcal{I}}_{\mathcal{X}}\subseteq V, for =0,,μ\ell=0,\ldots,\mu, to be the set of nodes that are reached by exactly \ell campaigns from the seed sets \mathcal{I}. Formally, for any value [μ]\ell\in[\mu],

V𝒳,:=τ([μ])(iτρXi(i)(Ii)j[μ]τρXj(j)(Ij)),\textstyle V_{\mathcal{X}}^{\ell,\mathcal{I}}:=\bigcup_{\tau\in\binom{[\mu]}{\ell}}\Big(\bigcap_{i\in\tau}\rho_{X_{i}}^{(i)}(I_{i})\setminus\bigcup_{j\in[\mu]\setminus\tau}\rho_{X_{j}}^{(j)}(I_{j})\Big),

where ([μ])\binom{[\mu]}{\ell} denotes the set {τ[μ]:|τ|=}\{\tau\subseteq[\mu]~:~|\tau|=\ell\}. We write V𝒳V_{\mathcal{X}}^{\ell}, if the initial seed sets \mathcal{I} are clear from the context. In the above definition, by convention an empty union is the empty set, while an empty intersection is the whole universe, here VV. Accordingly, we define

Φ(𝒮):=E𝒳[NoSMμ,ν(ρ𝒳(𝒮)V𝒳,)].\Phi^{\ell}(\mathcal{S}):=\E_{\mathcal{X}}\Big[\NoSM(\rho_{\mathcal{X}}(\mathcal{I}\cup\mathcal{S})\cap V^{\ell,\mathcal{I}}_{\mathcal{X}})\Big].

Note that Φ(𝒮)\Phi^{\ell}(\mathcal{S}) is the expected number of nodes that are reached by 00 or at least ν\nu campaigns, resulting from nodes that have been reached by exactly \ell campaigns from \mathcal{I}. Now, the objective function decomposes as

Φ(𝒮)=E𝒳[NoSMμ,ν(ρ𝒳(𝒮))]=E𝒳[[μ]NoSMμ,ν(ρ𝒳(𝒮)V𝒳)]=[μ]Φ(𝒮),\Phi(\mathcal{S})=\E_{\mathcal{X}}\Big[\NoSM(\rho_{\mathcal{X}}(\mathcal{I}\cup\mathcal{S}))\Big]=\E_{\mathcal{X}}\Big[\sum_{\ell\in[\mu]}\NoSM(\rho_{\mathcal{X}}(\mathcal{I}\cup\mathcal{S})\cap V^{\ell}_{\mathcal{X}})\Big]=\sum_{\ell\in[\mu]}\Phi^{\ell}(\mathcal{S}),

using linearity of expectation and that sets V𝒳V_{\mathcal{X}}^{\ell} are disjoint. Furthermore, we will denote by

Φ(𝒮):=i=μΦi(𝒮)=E𝒳[NoSMμ,ν(ρ𝒳(𝒮)(j=01V𝒳j))].\Phi^{\geq\ell}(\mathcal{S}):=\sum_{i=\ell}^{\mu}\Phi^{i}(\mathcal{S})=\E_{\mathcal{X}}\Big[\NoSM(\rho_{\mathcal{X}}(\mathcal{I}\cup\mathcal{S})\setminus(\cup_{j=0}^{\ell-1}V^{j}_{\mathcal{X}}))\Big].

Again, Φ(𝒮)\Phi^{\geq\ell}(\mathcal{S}) denotes the expected number of nodes that are reached by sufficiently many campaigns or none of them resulting from nodes that have previously been reached by at least \ell campaigns. Clearly, Φ(𝒮)=Φ0(𝒮)\Phi(\mathcal{S})=\Phi^{\geq 0}(\mathcal{S}). For convenience, in what follows, we will often refer to 𝒮\mathcal{S} as a set of pairs in V^:=V×[μ]\hat{V}:=V\times[\mu], where picking pair (v,i)(v,i) into 𝒮\mathcal{S} corresponds to picking vv into set SiS_{i}. We fix the following observations:

  • For =0\ell=0, Φ0(𝒮)\Phi^{0}(\mathcal{S}) is optimal when 𝒮=()i[μ]\mathcal{S}=(\emptyset)_{i\in[\mu]}. The achieved value is the expected size of V𝒳0V^{0}_{\mathcal{X}}: Φ0(𝒮)=E𝒳[NoSMμ,ν(ρ𝒳(()i[μ])V𝒳0)]=E𝒳[|V𝒳0|]\Phi^{0}(\mathcal{S})=\E_{\mathcal{X}}[\NoSM(\rho_{\mathcal{X}}(\mathcal{I}\cup(\emptyset)_{i\in[\mu]})\cap V^{0}_{\mathcal{X}})]=\E_{\mathcal{X}}[|V^{0}_{\mathcal{X}}|].

  • For =ν1\ell=\nu-1, the function Φν1(𝒮)=i=ν1μΦi(𝒮)\Phi^{\geq\nu-1}(\mathcal{S})=\sum_{i=\nu-1}^{\mu}\Phi^{i}(\mathcal{S}) is monotone and submodular.

A First Structural Lemma.

When applying the standard greedy hill climbing algorithm to finding a set of size kk maximizing a submodular set function the key property that is used in the analysis is the following. At any stage of the greedy algorithm there exists an element which leads to an improvement that is at least a fraction of kk of the difference of the optimal and the current solution, compare for example [Hoc97, Lemma 3.13]. Maybe the most important structural lemma underlying our algorithms is a very similar result for the functions Φ\Phi^{\geq\ell}.

Lemma 1.

Let [1,ν1]\ell\in[1,\nu-1] and 𝒮V^\mathcal{S}\subseteq\hat{V} with |𝒮|k(ν)|\mathcal{S}|\leq k-(\nu-\ell) and define U:={τV^,|τ|=ν}U:=\{\tau\subseteq\hat{V},|\tau|=\nu-\ell\}. Then, τ=𝚊𝚛𝚐max{Φ(Sτ):τU}\tau^{*}=\mathtt{arg}\max\{\Phi^{\geq\ell}(S\cup\tau):\tau\in U\} satisfies Φ(𝒮τ)Φ(𝒮)(Φ(𝒮)Φ(𝒮))/(kν),\Phi^{\geq\ell}(\mathcal{S}\cup\tau^{*})-\Phi^{\geq\ell}(\mathcal{S})\geq(\Phi^{\geq\ell}(\mathcal{S}^{*}_{\geq\ell})-\Phi^{\geq\ell}(\mathcal{S}))/\binom{k}{\nu-\ell}, where 𝒮\mathcal{S}^{*}_{\geq\ell} is an optimal solution of size kk to maximizing Φ\Phi^{\geq\ell}.

Proof.

Let 𝒳\mathcal{X} be an outcome profile and let vv be an arbitrary node in V:=Vj=01V𝒳jV^{\prime}:=V\setminus\bigcup_{j=0}^{\ell-1}V^{j}_{\mathcal{X}}. Let us denote by 𝟙𝒳𝒮(v)\mathds{1}^{\mathcal{S}}_{\mathcal{X}}(v) the indicator function that is one if vv is reached by at least ν\nu campaigns in outcome profile 𝒳\mathcal{X} from seed sets 𝒮\mathcal{I}\cup\mathcal{S} and zero otherwise. We note that Φ(𝒮)=E𝒳[vV𝟙𝒳𝒮(v)]\Phi^{\geq\ell}(\mathcal{S})=\E_{\mathcal{X}}\big[\sum_{v\in V^{\prime}}\mathds{1}_{\mathcal{X}}^{\mathcal{S}}(v)\big]. Now, define Y:={τ𝒮:|τ|=ν}Y:=\{\tau\subseteq\mathcal{S}_{\geq\ell}^{*}:|\tau|=\nu-\ell\}, i.e., YY are the sets of nodes in 𝒮\mathcal{S}^{*}_{\geq\ell} of size ν\nu-\ell. We now argue that the following inequality holds for vv and 𝒳\mathcal{X}:

𝟙𝒳𝒮(v)𝟙𝒳𝒮(v)τY(𝟙𝒳𝒮τ(v)𝟙𝒳𝒮(v)).\displaystyle\mathds{1}_{\mathcal{X}}^{\mathcal{S}^{*}_{\geq\ell}}(v)-\mathds{1}_{\mathcal{X}}^{\mathcal{S}}(v)\leq\sum\nolimits_{\tau\in Y}(\mathds{1}_{\mathcal{X}}^{\mathcal{S}\cup\tau}(v)-\mathds{1}_{\mathcal{X}}^{\mathcal{S}}(v)). (1)

If the left hand side is not positive, the inequality holds, since the right hand side cannot be negative by monotonicity. Hence, assume that the left hand side is positive. In that case it holds that 𝟙𝒳𝒮(v)=1\mathds{1}_{\mathcal{X}}^{\mathcal{S}^{*}_{\geq\ell}}(v)=1, but 𝟙𝒳𝒮(v)=0\mathds{1}_{\mathcal{X}}^{\mathcal{S}}(v)=0, i.e., in outcome profile 𝒳\mathcal{X}, vv is reached by at least ν\nu campaigns from seed sets 𝒮\mathcal{I}\cup\mathcal{S}^{*}_{\geq\ell} but not from seed sets 𝒮\mathcal{I}\cup\mathcal{S}. For such vv, there must be a set τY\tau\in Y such that adding τ\tau to 𝒮\mathcal{S} results in vv being reached by ν\nu campaigns (recall that vVv\in V^{\prime} and thus vv is already reached by at least \ell campaigns). Thus, there exists a set in YY that contributes a value of 1 on the right hand side and we may conclude that (1) holds. Now, using linearity of expectation and (1), we obtain

Φ(𝒮)Φ(𝒮)\displaystyle\Phi^{\geq\ell}(\mathcal{S}^{*}_{\geq\ell})-\Phi^{\geq\ell}(\mathcal{S}) =E𝒳[vV(𝟙𝒳𝒮(v)𝟙𝒳𝒮(v))]E𝒳[vVτY(𝟙𝒳𝒮τ(v)𝟙𝒳𝒮(v))].\displaystyle\!=\!\E_{\mathcal{X}}\Big[\sum_{v\in V^{\prime}}(\mathds{1}_{\mathcal{X}}^{\mathcal{S}^{*}_{\geq\ell}}(v)-\mathds{1}_{\mathcal{X}}^{\mathcal{S}}(v))\Big]\!\leq\!\E_{\mathcal{X}}\Big[\sum_{v\in V^{\prime}}\sum_{\tau\in Y}(\mathds{1}_{\mathcal{X}}^{\mathcal{S}\cup\tau}(v)-\mathds{1}_{\mathcal{X}}^{\mathcal{S}}(v))\Big].

Using linearity of expectation again, we obtain that the right hand side above is equal to τY(Φ(𝒮τ)Φ(𝒮))\sum_{\tau\in Y}(\Phi^{\geq\ell}(\mathcal{S}\cup\tau)-\Phi^{\geq\ell}(\mathcal{S})). Then, the statement follows by the maximality of τ\tau^{*} and the fact that |Y|(kν)|Y|\leq\binom{k}{\nu-\ell}. ∎

The Correlated Case.

For the correlated setting, where probability functions are identical for all campaigns and the cascade processes are completely correlated, we introduce an additional function called Ψ\Psi. First note that in the correlated setting, the outcome profile 𝒳\mathcal{X} in the definition of Φ(𝒮)\Phi(\mathcal{S}) satisfies X1==XμX_{1}=\ldots=X_{\mu}. In order to define Ψ\Psi, we introduce an additional fictitious campaign, call it campaign 00, that spreads with the same probability p0=p1==pμp_{0}=p_{1}=\ldots=p_{\mu} as the other μ\mu campaigns. We extend the outcome 𝒳=(Xi)i[μ]\mathcal{X}=(X_{i})_{i\in[\mu]} with X1==XμX_{1}=\ldots=X_{\mu} to contain also an identical copy X0X_{0} and define Ψ:2V×{0}[n]\Psi:2^{V\times\{0\}}\rightarrow[n] by

Ψ(𝒯):=E𝒳[|(ρX0(0)(𝒯)j=1ν1V𝒳j)j=νμV𝒳j|].\Psi(\mathcal{T}):=\E_{\mathcal{X}}\big[\big|\big(\rho_{X_{0}}^{(0)}(\mathcal{T})\cap\bigcup_{j=1}^{\nu-1}V_{\mathcal{X}}^{j}\big)\cup\bigcup_{j=\nu}^{\mu}V_{\mathcal{X}}^{j}\big|\big].

Observe that Ψ(𝒯)\Psi(\mathcal{T}) measures the expected number of nodes that are either (1) reached by more than ν\nu campaigns from \mathcal{I} or (2) are reached by at least one campaign from \mathcal{I} and are reached by the fictitious campaign 00 from 𝒯\mathcal{T}. Note that nodes from (1) are already reached by sufficiently many campaigns while nodes from (2) have been reached by some campaign from \mathcal{I} and, as witnessed by Ψ\Psi, can be reached from the nodes in 𝒯\mathcal{T}. Note that Ψ\Psi is monotone and submodular in 𝒮\mathcal{S} which follows directly from σ\sigma having these properties.

Approximating Ψ\Psi and Φ\Phi^{\geq\ell}.

As mentioned above, already in the standard independent cascade process, it is not feasible to evaluate the function σ\sigma exactly. However, σ\sigma can be approximated to within a factor of (1±ϵ)(1\pm\epsilon) by sampling a polynomial number of times. A very similar approach works for approximating the functions Ψ\Psi and Φ\Phi^{\geq\ell} for [0,ν]\ell\in[0,\nu]. That is, there is an algorithm approx(f,𝒮,,ν,ϵ,δ)(f,\mathcal{S},\mathcal{I},\nu,\epsilon,\delta) that, for f{Ψ,Φ0,,Φν}f\in\{\Psi,\Phi^{\geq 0},\ldots,\Phi^{\geq\nu}\}, sets 𝒮\mathcal{S} and \mathcal{I}, and parameters ν,ϵ,δ\nu,\epsilon,\delta returns a (1±ϵ)(1\pm\epsilon)-approximation of f(S)f(S) with probability 1δ1-\delta. We prove this fact in Appendix A in Lemma 14. The proof relies on a Chernoff bound and is very similar to the original proof of Proposition 4.1 in [KKT15] for the σ\sigma-function.

All of our algorithms are of a greedy flavor, that is, we greedily choose sets in order to build the output set 𝒮\mathcal{S}. We investigate the impact of the approximation on this approach in the following lemma. To this end, let ff be a function from {Ψ,Φ1,,Φν}\{\Psi,\Phi^{\geq 1},\ldots,\Phi^{\geq\nu}\} and, for some 0<ϵ10<\epsilon\leq 1, let f~\tilde{f} be a (1±ϵ)(1\pm\epsilon^{\prime})-approximation of ff with ϵ:=ϵ/(e(kλ(f)))\epsilon^{\prime}:=\epsilon/(e\cdot\binom{k}{\lambda(f)}), where λ(f)\lambda(f) depends on ff, namely λ(f):=ν\lambda(f):=\nu-\ell for f=Φf=\Phi^{\geq\ell} and λ(f):=1\lambda(f):=1 for f=Ψf=\Psi. We denote with DfD_{f} the universe over which ff is defined, i.e., Df:=V^D_{f}:=\hat{V} for f=Φf=\Phi^{\geq\ell}, while Df:=V×{0}D_{f}:=V\times\{0\} for f=Ψf=\Psi.

Lemma 2.

Let ff and f~\tilde{f} be as above for some 0<ϵ10<\epsilon\leq 1. Let U:={τDf,|τ|=λ(f)}U:=\{\tau\subseteq D_{f},|\tau|=\lambda(f)\}, 𝒮Df\mathcal{S}\subseteq D_{f} with |𝒮|kλ(f)|\mathcal{S}|\leq k-\lambda(f), and let 𝒮\mathcal{S}^{*} denote a set maximizing ff of size kk. Then, either

f(𝒮)(11e)f(𝒮) or f(𝒮τ~)f(𝒮)(1ϵ)(f(𝒮τ)f(𝒮)),f(\mathcal{S})\geq\big(1-\frac{1}{e}\big)\cdot f(\mathcal{S}^{*})\quad\text{ or }\quad f(\mathcal{S}\cup\tilde{\tau})-f(\mathcal{S})\geq(1-\epsilon)\cdot(f(\mathcal{S}\cup\tau^{*})-f(\mathcal{S})),

where τ:=𝚊𝚛𝚐max{f(Sτ):τU}\tau^{*}:=\mathtt{arg}\max\{f(S\cup\tau):\tau\in U\}, and τ~:=𝚊𝚛𝚐max{f~(Sτ):τU}\tilde{\tau}:=\mathtt{arg}\max\{\tilde{f}(S\cup\tau):\tau\in U\}.

We defer the proof to Appendix A. In summary: either 𝒮\mathcal{S} already yields a (11/e)(1-1/e)-approximation of the optimum of ff or a set τ\tau of size λ(f)\lambda(f) maximizing an approximation f~\tilde{f} of ff can lead to a progress of at least an (1ϵ)(1-\epsilon)-fraction of the maximum progress possible.

Maximizing Φν1\Phi^{\geq\nu-1} and Ψ\Psi.

Here, we fix the result that the standard greedy hill climbing algorithm, we refer to it as Greedy(f,ϵ,δ,,ν,k)(f,\epsilon,\delta,\mathcal{I},\nu,k), can be applied in order to approximate both f{Φν1,Ψ}f\in\{\Phi^{\geq\nu-1},\Psi\} to within a factor of 11/eϵ1-1/e-\epsilon for any 0<ϵ<10<\epsilon<1 with probability at least 1δ1-\delta for any 0<δ1/20<\delta\leq 1/2. This is based on the fact that these functions are submodular and monotone set functions. See Appendix A for a pseudo-code implementation and a proof of the submodularity property. Since we can only evaluate Φν1\Phi^{\geq\nu-1} and Ψ\Psi approximately, we obtain the additive ϵ\epsilon-term.

Lemma 3.

Let f{Φν1,Ψ}f\in\{\Phi^{\geq\nu-1},\Psi\} and let 0<ϵ<10<\epsilon<1 and 0<δ1/20<\delta\leq 1/2. With probability at least 1δ1-\delta, Greedy(f,ϵ,δ,,ν,k)(f,\epsilon,\delta,\mathcal{I},\nu,k) returns 𝒮\mathcal{S} satisfying f(𝒮)(11/eϵ)f(𝒮),f(\mathcal{S})\geq(1-1/e-\epsilon)\cdot f(\mathcal{S}^{*}), where 𝒮\mathcal{S}^{*} is an optimal solution of size kk to maximizing ff.

3 Hardness of Approximation for the Heterogeneous Case

We now let d2d\geq 2 be a constant. In this section, we show that in the heterogeneous setting for νd+1\nu\geq d+1, the μ\mu-ν\nu-Balance problem is as hard to approximate as the Densest-kk-Sub-dd-hypergraph problem [CDK+18]. Notably, this result has the following consequences: if d=2d=2 there is no ng(n)n^{-g(n)}-approximation algorithm with g(n)=o(1)g(n)=o(1) for μ\mu-ν\nu-Balance under the Gap Exponential Time Hypothesis (Gap-ETH). For general d3d\geq 3, we get that there is no nϵn^{-\epsilon}-approximation algorithm for a given constant ϵ>0\epsilon>0 which depends on dd under the assumption that a particular class of one way functions exists [App13]. We recall the definition of the Densest-kk-Sub-dd-hypergraph problem.

Densest-kk-Sub-dd-hypergraph Input: dd-Regular Hypergraph G=(V,E)G=(V,E), integer kdk\geq d. Find: set SVS\subseteq V with |S|k|S|\leq k, s.t. |E(S)||E(S)| is maximum, where E(S):={eE:eS}.E(S):=\{e\in E:e\subseteq S\}.

A dd-regular hypergraph is a hypergraph in which all hyperedges are composed of exactly dd vertices, where dd is a constant. When d=2d=2, Densest-kk-Sub-dd-hypergraph is known as the Densest-kk-Subgraph problem. For the hardness of approximation proof, we consider the following transform τ\tau of an instance (G=(V,E),k)(G=(V,E),k) of the Densest-kk-Sub-dd-hypergraph problem into an instance τ(G,k)=(G¯=(V¯,A¯),𝒫,,k¯)\tau(G,k)=(\overline{G}=(\overline{V},\overline{A}),\mathcal{P},\mathcal{I},\overline{k}) of the μ\mu-ν\nu-Balance problem.

  • Define V¯:=VV\overline{V}:=V_{\boxempty}\cup V_{\ocircle}, where V:=VV_{\boxempty}:=V, i.e., for each node vVv\in V, we get a node vv in V¯\overline{V}. Moreover, let J:=([μν+d]d)J:=\binom{[\mu-\nu+d]}{d}, and SdS_{d} be the set of permutations of [d][d]; we then define VV_{\ocircle} as V:={eι,πt:eE,ιJ,πSd,t[l]},V_{\ocircle}:=\{e^{t}_{\iota,\pi}:e\in E,\iota\in J,\pi\in S_{d},t\in[l]\}, i.e., for each edge eEe\in E, we create λl\lambda l nodes, where l:=|V|+1l:=|V|+1 and λ:=|Sd||J|=d!(μν+dd)\lambda:=|S_{d}|\cdot|J|=d!\binom{\mu-\nu+d}{d}. That is, each set ι\iota of dd campaigns in JJ, induces ll nodes eι,πt,t[l]e^{t}_{\iota,\pi},t\in[l] for each π\pi in SdS_{d}.

  • The arc set A¯\overline{A} and the probabilities are defined as shown in Figure 1 illustrating the case of d=3d=3 (a more detailed illustration is provided in Appendix B in Fig. 2). We get this scheme in G¯\overline{G} for every edge e={v1,,vd}Ee=\{v_{1},\ldots,v_{d}\}\in E, for each permutation π\pi in SdS_{d}, and for each set in JJ of dd campaigns.

  • The initial seed sets \mathcal{I} are defined as I1=I2==Iμν+d=I_{1}\!=\!I_{2}\!=\!\ldots\!=\!I_{\mu-\nu+d}\!=\!\emptyset, Iμν+d+1==Iμ=V¯I_{\mu-\nu+d+1}\!=\!\ldots\!=\!I_{\mu}\!=\!\overline{V}.

  • The budget is the same as in the Densest-kk-Sub-dd-hypergraph problem, i.e., k¯=k\overline{k}=k.

uuvvwweι,π1e^{1}_{\iota,\pi}eι,π2e^{2}_{\iota,\pi}eι,π3e^{3}_{\iota,\pi}\ldotseι,πl1e^{l-1}_{\iota,\pi}eι,πle^{l}_{\iota,\pi}pπ(i)=1p_{\pi(i)}=1pπ(j)=1p_{\pi(j)}=1pπ(k)=1p_{\pi(k)}=1pπ(i)=pπ(j)p_{\pi(i)}=p_{\pi(j)}=pπ(k)=1=p_{\pi(k)}=1pπ(i)=pπ(j)p_{\pi(i)}=p_{\pi(j)}=pπ(k)=1=p_{\pi(k)}=1pπ(i)=pπ(j)p_{\pi(i)}=p_{\pi(j)}=pπ(k)=1=p_{\pi(k)}=1
Figure 1: This figure illustrates the case d=3d=3. For an hyperedge e={u,v,w}e=\{u,v,w\} in GG, we get d!(μν+dd)d!\binom{\mu-\nu+d}{d} schemes of the above type, one for each set ι={i,j,k}J\iota=\{i,j,k\}\in J and for each way of ordering them given by a permutation πSd\pi\in S_{d}. Probabilities that are not given are equal to 0.

Note that each node in G¯\overline{G} is already covered by νd\nu-d campaigns and that the instance generated is deterministic, in the sense that probability values are either 0 or 1.

Let us now fix a μ\mu-ν\nu-Balance instance P=(G¯=(V¯,A¯),𝒫,,k¯)P=(\overline{G}=(\overline{V},\overline{A}),\mathcal{P},\mathcal{I},\overline{k}) resulting from the transform τ\tau as image of a Densest-kk-Sub-dd-hypergraph instance Q=(G=(V,E),k)Q=(G=(V,E),k). Clearly, V¯\overline{V} is of cardinality |V|+λl|E||V|+\lambda l|E| and A¯\overline{A} is of cardinality λ(l+d1)|E|\lambda(l+d-1)|E|. Let us denote by Σ\Sigma the set of feasible solutions for PP. For each 𝒮Σ\mathcal{S}\in\Sigma, it holds that the objective function Φ(𝒮)\Phi(\mathcal{S}) can be decomposed as Φ(𝒮)=Φ(𝒮)+Φ(𝒮),\Phi(\mathcal{S})=\Phi_{\boxempty}(\mathcal{S})+\Phi_{\ocircle}(\mathcal{S}), where

Φ(𝒮):=NoSMμ,ν(ρ𝒳(𝒮)V) and Φ(𝒮):=NoSMμ,ν(ρ𝒳(𝒮)V),\Phi_{\boxempty}(\mathcal{S}):=\NoSM(\rho_{\mathcal{X}}(\mathcal{I}\cup\mathcal{S})\cap V_{\boxempty})\quad\text{ and }\quad\Phi_{\ocircle}(\mathcal{S}):=\NoSM(\rho_{\mathcal{X}}(\mathcal{I}\cup\mathcal{S})\cap V_{\ocircle}),

for 𝒳\mathcal{X} being the only possible (deterministic) outcome profile. Now, let 𝒮\mathcal{S}^{*}, 𝒮\mathcal{S}^{*}_{\boxempty}, and 𝒮\mathcal{S}^{*}_{\ocircle} denote optimal solutions to the problem of maximizing Φ\Phi, Φ\Phi_{\boxempty}, and Φ\Phi_{\ocircle}, respectively, over Σ\Sigma. The following lemma whose proof can be found in Appendix B collects three statements. The first statement says that an optimal solution to Φ\Phi also maximizes Φ\Phi_{\ocircle}. The second statement says that there exists a feasible solution to PP which achieves at least a multiple of lpl\cdot p of the objective value in Densest-kk-Sub-dd-hypergraph with p=d!/ddp=d!/d^{d}. In the third statement, we observe that from a feasible solution to PP, we can construct a feasible solution to QQ while loosing only a factor of λl\lambda l in objective value.

Lemma 4.
  1. (1)

    An optimal solution to Φ\Phi also maximizes Φ\Phi_{\ocircle}, i.e., Φ(𝒮)=Φ(𝒮)\Phi_{\ocircle}(\mathcal{S}_{\ocircle}^{*})=\Phi_{\ocircle}(\mathcal{S}^{*}).

  2. (2)

    It holds that Φ(𝒮)lpDKSHd\Phi_{\ocircle}(\mathcal{S}^{*}_{\ocircle})\geq l\cdot p\cdot\dkshop^{*}_{d}, where DKSHd\dkshop^{*}_{d} is the optimal value of Densest-kk-Sub-dd-hypergraph in QQ and p=d!/ddp=d!/d^{d}.

  3. (3)

    Given 𝒮Σ\mathcal{S}\in\Sigma, we can, in polynomial time, build a feasible solution SS of QQ such that |E(S)|Φ(𝒮)/(λl)|E(S)|\geq\Phi_{\ocircle}(\mathcal{S})/(\lambda l).

We are now ready to show the following relations between the complexity of the two problems. Note that the assumption that gg is non-increasing is w.l.o.g.

Theorem 5.

Let d2d\geq 2, νd+1\nu\geq d+1, and p=d!/ddp=d!/d^{d}, then we have the following two cases:

Case d=2d=2: Let α(n)=ng(n)\alpha(n)=n^{-g(n)} with gg being non-increasing, g(n)=o(1)g(n)=o(1) and α(n)(0,1]\alpha(n)\in(0,1] and β(n)=pn6g(n)2λ\beta(n)=\frac{p\cdot n^{-6g(n)}}{2\lambda}.

Case d3d\geq 3: Let α(n)=nϵ(d)\alpha(n)=n^{-\epsilon(d)} where ϵ(d)>0\epsilon(d)>0 is a constant which depends on dd, α(n)(0,1]\alpha(n)\in(0,1] and β(n)=pnϵ(d)2λ\beta(n)=\frac{p\cdot n^{-\epsilon^{\prime}(d)}}{2\lambda}, with ϵ(d)=(d+4)ϵ(d)\epsilon^{\prime}(d)=(d+4)\cdot\epsilon(d).

In both cases the following statement holds: If there is an α(|V¯|)\alpha(|\overline{V}|)-approximate algorithm for the deterministic μ\mu-ν\nu-Balance problem, then there is a β(|V|)\beta(|V|)-approximate algorithm for Densest-kk-Sub-dd-hypergraph. Here |V¯||\overline{V}| and |V||V| denote the number of vertices in the μ\mu-ν\nu-Balance and the Densest-kk-Sub-dd-hypergraph problems, respectively and λ=d!|J|\lambda=d!|J|.

Proof.

Let Q=(G,k)Q=(G,k) be an instance of the Densest-kk-Sub-dd-hypergraph problem and let P:=(G¯=(V¯,E¯),𝒫,,k¯)=τ(G,k)P:=(\overline{G}=(\overline{V},\overline{E}),\mathcal{P},\mathcal{I},\overline{k})=\tau(G,k) be the instance of the μ\mu-ν\nu-Balance problem obtained by the transform τ\tau. For brevity, let n:=|V|n:=|V| and n¯:=|V¯|\overline{n}:=|\overline{V}|. Moreover, let 𝒮\mathcal{S} be an α(|V¯|)\alpha(|\overline{V}|)-approximate solution to PP, that is Φ(𝒮)α(|V¯|)Φ(𝒮)\Phi(\mathcal{S})\geq\alpha(|\overline{V}|)\Phi(\mathcal{S}^{*}). We show how to construct a β(n)\beta(n)-approximate solution SS to QQ.

Using Lemma 4, (3), we obtain a feasible solution SS to QQ with |E(S)|Φ(𝒮)/(λl)|E(S)|\geq\Phi_{\ocircle}(\mathcal{S})/(\lambda l). We proceed by lower-bounding Φ(𝒮)\Phi_{\ocircle}(\mathcal{S}). We can w.l.o.g. assume that 𝒮V=\mathcal{S}\cap V_{\ocircle}=\emptyset and that Φ(𝒮)l\Phi_{\ocircle}(\mathcal{S})\geq l. Indeed, if Φ(𝒮)<l\Phi_{\ocircle}(\mathcal{S})<l then Φ(𝒮)=0\Phi_{\ocircle}(\mathcal{S})=0 and we can build in polynomial-time a better solution by identifying one edge (v1,,vd)(v_{1},\ldots,v_{d}) and propagating campaign ii in viv_{i}. This further implies that Φ(𝒮)Φ(𝒮)\Phi_{\ocircle}(\mathcal{S})\geq\Phi_{\boxempty}(\mathcal{S}) as l>nΦ(𝒮)l>n\geq\Phi_{\boxempty}(\mathcal{S}). We obtain

Φ(𝒮)Φ(𝒮)2α(n¯)2Φ(𝒮)α(n¯)2Φ(𝒮)=α(n¯)2Φ(𝒮)α(n¯)lp2DKSHd,\Phi_{\ocircle}(\mathcal{S})\geq\frac{\Phi(\mathcal{S})}{2}\geq\frac{\alpha(\overline{n})}{2}\cdot\Phi(\mathcal{S}^{*})\geq\frac{\alpha(\overline{n})}{2}\cdot\Phi_{\ocircle}(\mathcal{S}^{*})=\frac{\alpha(\overline{n})}{2}\cdot\Phi_{\ocircle}(\mathcal{S}^{*}_{\ocircle})\geq\frac{\alpha(\overline{n})\cdot l\cdot p}{2}\cdot\dkshop^{*}_{d},

using Lemma 4, (1) and (2) in the last two steps. In summary, we have |E(S)|α(n¯)p2λDKSHd|E(S)|\geq\frac{\alpha(\overline{n})\cdot p}{2\lambda}\dkshop^{*}_{d}. Note that 2λ/p2\lambda/p is a constant.

Case d=2d=2: Since gg is non-increasing, we get α(n¯)=n¯g(n¯)=2g(n¯)log(n¯)2g(n)log(2λn3)26g(n)log(n)=n6g(n),\alpha(\overline{n})={\overline{n}}^{-g(\overline{n})}=2^{-g(\overline{n})\log(\overline{n})}\geq 2^{-g(n)\log(2\lambda n^{3})}\geq 2^{-6g(n)\log(n)}=n^{-6g(n)}, where we used 2nn¯2λn32\leq n\leq\overline{n}\leq 2\lambda n^{3} and λμ2n2\lambda\leq\mu^{2}\leq n^{2} (as d=2d=2). This completes this case.

Case d3d\geq 3: In this case α(n¯)=n¯ϵ(d)=2ϵ(d)log(n¯)2ϵ(d)log(2λn3)2(d+4)ϵ(d)log(n)=n(d+4)ϵ(d),\alpha(\overline{n})={\overline{n}}^{-\epsilon(d)}=2^{-\epsilon(d)\log(\overline{n})}\geq 2^{-\epsilon(d)\log(2\lambda n^{3})}\geq 2^{-(d+4)\epsilon(d)\log(n)}=n^{-(d+4)\epsilon(d)}, where we used 2nn¯2λn32\leq n\leq\overline{n}\leq 2\lambda n^{3} and λμdnd\lambda\leq\mu^{d}\leq n^{d}. This completes this case.∎

To sum up, our reduction shows that: (1) as Densest-kk-Sub-dd-hypergraph cannot be approximated within 1/nϵ1/n^{\epsilon} for some constant ϵ>0\epsilon>0 which depends on dd, if a particular class of one way functions exists [App13], we have shown that the same hardness result holds for any μ\mu-ν\nu-Balance problem with νd+14\nu\geq d+1\geq 4; (2) moreover as Densest-kk-Subgraph cannot be approximated within 1/no(1)1/n^{o(1)}, if the Gap-ETH holds [Man17], we have shown that the same hardness result holds for any μ\mu-ν\nu-Balance problem with ν3\nu\geq 3.

Other approximation hardness results exist for Densest-kk-Subgraph. We review them here, highlighting the hardness results that our reduction implies in each case.

  • Densest-kk-Subgraph cannot be approximated within any constant, if the Unique Games with Small Set Expansion (UGSSE) conjecture holds [RS10]. Therefore, under the UGSSE conjecture it is easy to prove that the reduction given above shows that any μ\mu-ν\nu-Balance problem with ν3\nu\geq 3 cannot be approximated within any constant.

  • Densest-kk-Subgraph cannot be approximated within n(loglogn)cn^{-(\log\log n)^{-c}}, for some constant cc if the exponential time hypothesis holds [Man17]. Under the same conjecture, our reduction implies the same hardness result for any μ\mu-ν\nu-Balance problem with ν3\nu\geq 3.

4 Approximation Algorithm for the Heterogeneous Case

Our approach for maximizing Φ(𝒮)\Phi(\mathcal{S}) decomposes it as Φ(𝒮)=Φ0(𝒮)+Φ1(𝒮)\Phi(\mathcal{S})=\Phi^{0}(\mathcal{S})+\Phi^{\geq 1}(\mathcal{S}) and works on each summand separately. In the following two subsections, we give two different algorithms for maximizing Φ1(𝒮)\Phi^{\geq 1}(\mathcal{S}). At the end of the section, we show how to combine them.

Greedily Picking Tuples.

In this paragraph, we present GreedyTuple(ϵ,δ,,,ν,k)(\epsilon,\delta,\ell,\mathcal{I},\nu,k) that, for given \ell, computes a solution to maximizing Φ\Phi^{\geq\ell}. For =ν1\ell=\nu-1 the algorithm is identical to the standard greedy hill climbing algorithm. For the general case of ν1\ell\leq\nu-1, we will show the following theorem. The algorithm is inspired by a greedy algorithm, called Greedy1, due to [DOS18] for solving the so-called maximum coverage with pairs problem.

Theorem 6.

Let ϵ(0,1)\epsilon\in(0,1), δ1/2\delta\leq 1/2, and [1,ν1]\ell\in[1,\nu-1]. If k2ν/ϵk\geq 2\nu/\epsilon, with probability at least 1δ1-\delta, the algorithm GreedyTuple(ϵ,δ,,,ν,k)(\epsilon,\delta,\ell,\mathcal{I},\nu,k) returns a solution 𝒮\mathcal{S} satisfying Φ(𝒮)(11eϵ)/(k1ν1)Φ(𝒮),\Phi^{\geq\ell}(\mathcal{S})\geq(1-\frac{1}{e}-\epsilon)/\binom{k-1}{\nu-\ell-1}\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}), where 𝒮\mathcal{S}_{\geq\ell}^{*} is an optimal solution to Φ\Phi^{\geq\ell}.

t:=kν(|V^|ν)t:=\lceil\frac{k}{\nu-\ell}\rceil\binom{|\hat{V}|}{\nu-\ell}, δδ/t\delta^{\prime}\leftarrow\delta/t, ϵϵ/(2e(kν))\epsilon^{\prime}\leftarrow\epsilon/(2e\cdot\binom{k}{\nu-\ell}), 𝒮\mathcal{S}\leftarrow\emptyset 1
while |𝒮|k(ν)|\mathcal{S}|\leq k-(\nu-\ell) do 2
Compute τ𝚊𝚛𝚐maxτV^,|τ|=ν{approx(Φ,𝒮τ,,ν,ϵ,δ)}\tau\leftarrow\mathtt{arg}\max_{\tau\subseteq\hat{V},|\tau|={\nu-\ell}}\{\texttt{approx}(\Phi^{\geq\ell},\mathcal{S}\cup\tau,\mathcal{I},\nu,\epsilon^{\prime},\delta^{\prime})\}, set 𝒮𝒮τ\mathcal{S}\leftarrow\mathcal{S}\cup\tau return 𝒮\mathcal{S}
Algorithm 1 GreedyTuple(ϵ,δ,,,ν,k)(\epsilon,\delta,\ell,\mathcal{I},\nu,k)

We let 𝒮i\mathcal{S}^{i} denote the set 𝒮\mathcal{S} at the end of iteration ii of the algorithm. The main idea underlying the analysis of GreedyTuple is very much related to the analysis of the standard greedy algorithm. That is (ignoring the approximation issue), every step of the algorithm incurs a factor of (1(11/(kν))CLOSE(1-(1-1/\binom{k}{\nu-\ell}). For =ν1\ell=\nu-1, this coincides with the standard case.

Lemma 7.

Let 0<ϵ<10<\epsilon<1, δ1/2\delta\leq 1/2, and [1,ν1]\ell\in[1,\nu-1]. With probability at least 1δ1-\delta, after each iteration ii of Algorithm 1, it either holds that

Φ(𝒮i)(1(11ϵ2(kν))i)Φ(𝒮) or Φ(𝒮i)(11e)Φ(𝒮).\Phi^{\geq\ell}(\mathcal{S}^{i})\geq\Big(1-\Big(1-\frac{1-\frac{\epsilon}{2}}{\binom{k}{\nu-\ell}}\Big)^{i}\Big)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*})\quad\text{ or }\quad\Phi^{\geq\ell}(\mathcal{S}^{i})\geq\big(1-\frac{1}{e}\big)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}).

The proof of this lemma can be found in Appendix C.1, it uses Lemmata 1 and 2. We are now ready to give the proof of Theorem 6.

Proof of Theorem 6.

Let 𝒮\mathcal{S} denote the set returned by the algorithm. Clearly, Φ(𝒮)Φ(𝒮ι)\Phi^{\geq\ell}(\mathcal{S})\geq\Phi^{\geq\ell}(\mathcal{S}^{\iota}), where ι\iota denotes the number of iterations of the while loop in the algorithm. By assumption k2ν/ϵk\geq 2\nu/\epsilon and thus ι=kνkν1(1ϵ2)kν\iota=\lfloor\frac{k}{\nu-\ell}\rfloor\geq\frac{k}{\nu-\ell}-1\geq(1-\frac{\epsilon}{2})\cdot\frac{k}{\nu-\ell}. Using Lemma 7 for 𝒮ι\mathcal{S}^{\iota} yields that either Φ(𝒮ι)(11/e)Φ(𝒮)\Phi^{\geq\ell}(\mathcal{S}^{\iota})\geq(1-1/e)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}) or

Φ(𝒮ι)(1(11ϵ2(kν))ι)Φ(𝒮)(1(11ϵ2(kν))(1ϵ2)kν)Φ(𝒮).\displaystyle\Phi^{\geq\ell}(\mathcal{S}^{\iota})\geq\Big(1-\Big(1-\frac{1-\frac{\epsilon}{2}}{\binom{k}{\nu-\ell}}\Big)^{\iota}\Big)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*})\geq\Big(1-\Big(1-\frac{1-\frac{\epsilon}{2}}{\binom{k}{\nu-\ell}}\Big)^{(1-\frac{\epsilon}{2})\frac{k}{\nu-\ell}}\Big)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}).

For the former case, note that 11/e1-1/e is greater than the approximation factor required by the theorem. For the latter case note that, as 1xexp(x)1-x\leq\exp(-x) for any real xx, we have

1(11ϵ2(kν))(1ϵ2)kν1exp((1ϵ2)2(k1ν1))11eϵ(k1ν1),1-\Big(1-\frac{1-\frac{\epsilon}{2}}{\binom{k}{\nu-\ell}}\Big)^{(1-\frac{\epsilon}{2})\frac{k}{\nu-\ell}}\geq 1-\exp\Big(\frac{-(1-\frac{\epsilon}{2})^{2}}{\binom{k-1}{\nu-\ell-1}}\Big)\geq\frac{1-\frac{1}{e}-\epsilon}{\binom{k-1}{\nu-\ell-1}},

where the last inequality uses that 1exp(x)x(1exp(1))1-\exp(-x)\leq x\cdot(1-\exp(-1)) and (11/e)(1x)11/ex(1-1/e)(1-x)\leq 1-1/e-x for any x0x\geq 0. This completes the proof. ∎

Being Iteratively Greedy.

Recall that, at the beginning of this section, we have defined Φ(𝒮):=E𝒳[NoSMμ,ν(ρ𝒳(𝒮)(j=01V𝒳j,))]\Phi^{\geq\ell}(\mathcal{S}):=\E_{\mathcal{X}}[\NoSM(\rho_{\mathcal{X}}(\mathcal{I}\cup\mathcal{S})\setminus(\cup_{j=0}^{\ell-1}V^{j,\mathcal{I}}_{\mathcal{X}}))]. We now extend this notation by letting

Φβ(,𝒮):=E𝒳[NoSMμ,β(ρ𝒳(𝒮)j=01V𝒳j,)]\textstyle\Phi^{\geq\ell}_{\beta}(\mathcal{R},\mathcal{S}):=\E_{\mathcal{X}}[\NoSMs_{\mu,\beta}(\rho_{\mathcal{X}}(\mathcal{R}\cup\mathcal{S})\setminus\bigcup_{j=0}^{\ell-1}V^{j,\mathcal{R}}_{\mathcal{X}})]

where [ν1]\ell\in[\nu-1] and β[ν]\beta\in[\nu]; we will mainly be working with the case β=+1\beta=\ell+1. The function measures the expected number of nodes that are reached by at least β\beta campaigns from 𝒮\mathcal{R}\cup\mathcal{S} within the set of nodes that have originally been reached by at least \ell campaigns from \mathcal{R}. Our goal now is to maximize Φ()\Phi(\cdot) through the following iterative scheme: for \ell from 11 to ν1\nu-1, we find sets 𝒮[]\mathcal{S}^{[\ell]} of size k/(ν1)\lfloor k/(\nu-1)\rfloor maximizing Φ+1([],)\Phi^{\geq\ell}_{\ell+1}(\mathcal{R}^{[\ell]},\cdot), where []:=j=11[j]\mathcal{R}^{[\ell]}:=\mathcal{I}\cup\bigcup_{j=1}^{\ell-1}\mathcal{R}^{[j]}. That is, in the th\ell^{th} iteration, we maximize the number of nodes reached by +1\ell+1 campaigns that have previously been reached by at least \ell campaigns. The approach is motivated by the observation that, for any [ν1]\ell\in[\nu-1] and initial sets \mathcal{R}, the function Φ+1(,𝒮)\Phi^{\geq\ell}_{\ell+1}(\mathcal{R},\mathcal{S}) is monotone and submodular in 𝒮\mathcal{S}, compare with Section 2.2 where we used this fact for =ν1\ell=\nu-1. Using Lemma 3 applied to Φ+1(,)\Phi^{\geq\ell}_{\ell+1}(\mathcal{R},\cdot) with ν=+1\nu=\ell+1 we get that the standard greedy algorithm can be used in order to obtain a (11/eϵ)(1-1/e-\epsilon)-approximate solution when maximizing Φ+1(,)\Phi^{\geq\ell}_{\ell+1}(\mathcal{R},\cdot). Note that our algorithm, called GreedyIter is inspired by a similar greedy algorithm called Greedy2 from [DOS18] that is used there for the maximum coverage with pairs problem. We will prove the following theorem in this section.

Theorem 8.

Let 0<ϵ<10<\epsilon<1 and δ1/2\delta\leq 1/2. With probability 1δ1-\delta, GreedyIter(ϵ,δ,,ν,k)(\epsilon,\delta,\mathcal{I},\nu,k) returns 𝒮\mathcal{S} satisfying Φ(𝒮)(11eϵ)ν1ν2ν3(k2|V|)ν2Φ1(,𝒮1),\Phi(\mathcal{S})\geq\frac{(1-\frac{1}{e}-\epsilon)^{\nu-1}}{\nu^{2\nu-3}}(\frac{k}{2|V|})^{\nu-2}\cdot\Phi^{\geq 1}(\mathcal{I},\mathcal{S}^{*}_{\geq 1}), where 𝒮1\mathcal{S}_{\geq 1}^{*} is a set of cardinality kk maximizing Φ1(,)\Phi^{\geq 1}(\mathcal{I},\cdot).

δδ/ν\delta^{\prime}\leftarrow\delta/\nu, ϵϵ/2\epsilon^{\prime}\leftarrow\epsilon/2, [1]\mathcal{R}^{[1]}\leftarrow\mathcal{I} 1
for =1,,ν1\ell=1,\ldots,\nu-1 do 2
𝒮[]Greedy(Φ+1([],),ϵ,δ,[],+1,k/(ν1))\mathcal{S}^{[\ell]}\leftarrow\textsc{Greedy}(\Phi^{\geq\ell}_{\ell+1}(\mathcal{R}^{[\ell]},\cdot),\epsilon^{\prime},\delta^{\prime},\mathcal{R}^{[\ell]},\ell+1,\lfloor k/(\nu-1)\rfloor), [+1][]𝒮[]\mathcal{R}^{[\ell+1]}\leftarrow\mathcal{R}^{[\ell]}\cup\mathcal{S}^{[\ell]} 3
 return i=1ν1𝒮[i]\bigcup_{i=1}^{\nu-1}\mathcal{S}^{[i]}
Algorithm 2 GreedyIter(ϵ,δ,,ν,k)(\epsilon,\delta,\mathcal{I},\nu,k)

The proof of Theorem 8 relies on the following two lemmata whose proofs are given in the appendix, see Section C.2. In a sense the first lemma quantifies the loss in approximation of the first iteration of GreedyIter, while the second lemma quantifies the loss of the later iterations. Both proofs rely on the submodularity of Φ+1(,)\Phi^{\geq\ell}_{\ell+1}(\mathcal{R},\cdot).

Lemma 9.

Let ϵ>0\epsilon>0 and assume that k2(ν1)/ϵk\geq 2(\nu-1)/\epsilon. If 𝒮[1]V^\mathcal{S}^{[1]}\subseteq\hat{V} is the set of cardinality k/(ν1)\lfloor k/(\nu-1)\rfloor selected in the first iteration of GreedyIter(ϵ,δ,,ν,k)(\epsilon,\delta,\mathcal{I},\nu,k), then, with probability at least 1δ/ν1-\delta/\nu, it holds that Φ21(,𝒮[1])11eϵνΦ1(,𝒮1),\Phi^{\geq 1}_{2}(\mathcal{I},\mathcal{S}^{[1]})\geq\frac{1-\frac{1}{e}-\epsilon}{\nu}\cdot\Phi^{\geq 1}(\mathcal{I},\mathcal{S}^{*}_{\geq 1}), where 𝒮1\mathcal{S}^{*}_{\geq 1} is a set of cardinality kk maximizing Φ1(,)\Phi^{\geq 1}(\mathcal{I},\cdot).

Lemma 10.

Let ϵ>0\epsilon>0, 2\ell\geq 2 and assume that k2(ν1)/ϵk\geq 2(\nu-1)/\epsilon. If 𝒮[]V^\mathcal{S}^{[\ell]}\subseteq\hat{V} is the set of cardinality k/(ν1)\lfloor k/(\nu-1)\rfloor selected in the \ell’th iteration of GreedyIter(ϵ,δ,,ν,k)(\epsilon,\delta,\mathcal{I},\nu,k), then, with probability at least 1δ/ν1-\delta/\nu, it holds that Φ+1([],𝒮[])(11eϵ)k2(+1)(ν1)|V|Φ1([1],𝒮[1]).\Phi^{\geq\ell}_{\ell+1}(\mathcal{R}^{[\ell]},\mathcal{S}^{[\ell]})\!\geq\!\frac{(1-\frac{1}{e}-\epsilon)k}{2(\ell+1)(\nu-1)|V|}\cdot\Phi^{\geq\ell-1}_{\ell}(\mathcal{R}^{[\ell-1]},\mathcal{S}^{[\ell-1]}).

Proof of Theorem 8.

Since 𝒮=i=1ν1𝒮[i]\mathcal{S}=\bigcup_{i=1}^{\nu-1}\mathcal{S}^{[i]}, we obtain Φ(𝒮)Φνν1([ν1],𝒮ν1)\Phi(\mathcal{S})\geq\Phi^{\geq\nu-1}_{\nu}(\mathcal{R}^{[\nu-1]},\mathcal{S}^{\nu-1}). Using the union bound, ν2\nu-2 times Lemma 10 and then Lemma 9 yield that, with probability at least 1δ1-\delta, it holds that

Φ(𝒮)\displaystyle\Phi(\mathcal{S}) ((11eϵ)k2ν(ν1)|V|)ν2Φ21([1],𝒮[1])(11eϵ)ν1ν2ν3(k2|V|)ν2Φν1(,𝒮1).\displaystyle\geq\Big(\frac{(1-\frac{1}{e}-\epsilon)k}{2\nu(\nu-1)|V|}\Big)^{\nu-2}\Phi^{\geq 1}_{2}(\mathcal{R}^{[1]},\mathcal{S}^{[1]})\geq\frac{(1-\frac{1}{e}-\epsilon)^{\nu-1}}{\nu^{2\nu-3}}\Big(\frac{k}{2|V|}\Big)^{\nu-2}\Phi^{\geq 1}_{\nu}(\mathcal{I},\mathcal{S}^{*}_{\geq 1}).\qed
Algorithm for General Heterogeneous μ\mu-ν\nu-Balance problem.

Our approach to solving the general μ\mu-ν\nu-Balance problem is now to use both algorithms presented above. According to Theorem 6, using GreedyTuple(ϵ,δ/2,1,,ν,k)(\epsilon,\delta/2,1,\mathcal{I},\nu,k), we obtain a set 𝒮1\mathcal{S}^{1} that with probability 1δ/21-\delta/2 satisfies Φ1(𝒮1)α1Φ1(𝒮1)\Phi^{\geq 1}(\mathcal{S}^{1})\geq\alpha_{1}\cdot\Phi^{\geq 1}(\mathcal{S}_{\geq 1}^{*}), where 𝒮1\mathcal{S}_{\geq 1}^{*} denotes an optimal solution of size kk to maximizing Φ1\Phi^{\geq 1} and α1=(11eϵ)/(k1ν2)\alpha_{1}=(1-\frac{1}{e}-\epsilon)/\binom{k-1}{\nu-2}. According to Theorem 8, using GreedyIter(ϵ,δ/2,,ν,k)(\epsilon,\delta/2,\mathcal{I},\nu,k), we obtain a set 𝒮2\mathcal{S}^{2} that with probability 1δ/21-\delta/2 satisfies Φ(𝒮2)α2Φ1(𝒮1)\Phi(\mathcal{S}^{2})\geq\alpha_{2}\cdot\Phi^{\geq 1}(\mathcal{S}^{*}_{\geq 1}), where 𝒮1\mathcal{S}_{\geq 1}^{*} is as above and α2=(11eϵ)ν1(k(2|V|CLOSE)ν2/ν2ν3\alpha_{2}=(1-\frac{1}{e}-\epsilon)^{\nu-1}(\frac{k}{(2|V|})^{\nu-2}/\nu^{2\nu-3}. Now, we define 𝒮\mathcal{S}^{\prime} to be the solution that achieves the maximum max{Φ(𝒮1),Φ(𝒮2)}\max\{\Phi(\mathcal{S}^{1}),\Phi(\mathcal{S}^{2})\} and 𝒮\mathcal{S} to be the solution that achieves the maximum max{Φ(),Φ(𝒮)}\max\{\Phi(\emptyset),\Phi(\mathcal{S}^{\prime})\}. We obtain 2Φ(𝒮)Φ()+Φ(𝒮)Φ0()+α1α2Φ1(𝒮1),2\cdot\Phi(\mathcal{S})\geq\Phi(\emptyset)+\Phi(\mathcal{S}^{\prime})\geq\Phi^{0}(\emptyset)+\sqrt{\alpha_{1}\cdot\alpha_{2}}\cdot\Phi^{\geq 1}(\mathcal{S}_{\geq 1}^{*}), using that the maximum Φ(𝒮)\Phi(\mathcal{S}^{\prime}) is lower bounded by the geometric mean of Φ(𝒮1)\Phi(\mathcal{S}^{1}) and Φ(𝒮2)\Phi(\mathcal{S}^{2}), which are in turn lower bounded by α1Φ1(𝒮1)\alpha_{1}\cdot\Phi^{\geq 1}(\mathcal{S}_{\geq 1}^{*}) and α2Φ1(𝒮1)\alpha_{2}\cdot\Phi^{\geq 1}(\mathcal{S}^{*}_{\geq 1}), respectively. Now, let 𝒮\mathcal{S}^{*} be an optimal solution of size kk to maximizing Φ\Phi. Using that the empty set maximizes Φ0\Phi^{0}, we have Φ0()Φ0(𝒮)\Phi^{0}(\emptyset)\geq\Phi^{0}(\mathcal{S}^{*}). Furthermore Φ1(𝒮1)Φ1(𝒮)\Phi^{\geq 1}(\mathcal{S}_{\geq 1}^{*})\geq\Phi^{\geq 1}(\mathcal{S}^{*}), thus

OPENΦ(𝒮)α1α22(Φ0(𝒮)+Φ1(𝒮)))=α1α22Φ(𝒮).\Phi(\mathcal{S})\geq\frac{\sqrt{\alpha_{1}\alpha_{2}}}{2}\cdot(\Phi^{0}(\mathcal{S}^{*})+\Phi^{\geq 1}(\mathcal{S}^{*})))=\frac{\sqrt{\alpha_{1}\alpha_{2}}}{2}\cdot\Phi(\mathcal{S}^{*}).

Plugging in α1\alpha_{1} and α2\alpha_{2} and using kν2(k1ν2)k^{\nu-2}\geq\binom{k-1}{\nu-2}, we get the following theorem.

Theorem 11.

Let 0<ϵ<10<\epsilon<1 and δ1/2\delta\leq 1/2. There is an algorithm that, with probability 1δ1-\delta, outputs a solution 𝒮\mathcal{S} that satisfies Φ(𝒮)(11eϵ)ν2(12|V|)ν22ν2ν32Φ(𝒮),\Phi(\mathcal{S})\geq(1-\frac{1}{e}-\epsilon)^{\frac{\nu}{2}}(\frac{1}{2|V|})^{\frac{\nu-2}{2}}\nu^{-\frac{2\nu-3}{2}}\cdot\Phi(\mathcal{S}^{*}), where 𝒮\mathcal{S}^{*} denotes an optimal solution of size kk to maximizing Φ()\Phi(\cdot).

Note that for ν=3\nu=3, we obtain an algorithm with an approximation ratio of order n1/2n^{-1/2}.

5 Approximation Algorithm for the Correlated Case

We now turn to the correlated case. Recall that here the probability functions are identical for all campaigns, i.e., p1(e)==pμ(e)p_{1}(e)=\ldots=p_{\mu}(e) for every edge eEe\in E. Moreover, the cascade processes are completely correlated, that is, for any edge (u,v)(u,v), if node uu propagates campaign ii to vv, then node uu also propagates all other campaigns that reach it to vv.

We will consider the same decomposition of the objective function as in the heterogeneous case, i.e., Φ(𝒮)\Phi(\mathcal{S}) as Φ(𝒮)=Φ1(𝒮)+Φ0(𝒮)\Phi(\mathcal{S})=\Phi^{\geq 1}(\mathcal{S})+\Phi^{0}(\mathcal{S}) for a solution 𝒮\mathcal{S}. Recall that Φ1(𝒮)\Phi^{\geq 1}(\mathcal{S}) counts the number of nodes that are reached by sufficiently many, i.e. ν\nu, campaigns from 𝒮\mathcal{I}\cup\mathcal{S} and have been reached by at least one campaign from \mathcal{I}. Similarly, Φ0(𝒮)\Phi^{0}(\mathcal{S}) counts nodes that are reached by sufficiently many campaigns or none and have previously been reached by no campaign from \mathcal{I}. Clearly, as in the heterogeneous case, Φ0(𝒮)\Phi^{0}(\mathcal{S}) is optimal when 𝒮=\mathcal{S}=\emptyset. Differently from the heterogeneous case however, we will see that in the correlated setting, there is an approximation algorithm for Φ1\Phi^{\geq 1} that achieves a constant factor, namely, (11/eϵ)/(ν+1)(1-1/e-\epsilon)/(\nu+1). The idea is to pick ν\nu campaigns and propagate them in the same k/ν\lfloor k/\nu\rfloor nodes, exploiting that all campaigns spread in an identical manner.

To that end, we consider the problem of maximizing influence spread with one fictitious campaign, say campaign 00, spreading with the same probabilities as the others. We will consider the nodes reached by campaign 00 among the nodes that were (a) reached by at least one campaign from \mathcal{I} and were (b) reached by no more than ν\nu campaigns from \mathcal{I}. For this purpose, we had defined the function Ψ\Psi in Section 2. Recall

Ψ(𝒯):=E𝒳[|(ρX0(0)(𝒯)j=1ν1V𝒳j)j=νμV𝒳j|]\Psi(\mathcal{T}):=\E_{\mathcal{X}}\Big[\big|(\rho_{X_{0}}^{(0)}(\mathcal{T})\cap\bigcup_{j=1}^{\nu-1}V_{\mathcal{X}}^{j})\cup\bigcup_{j=\nu}^{\mu}V_{\mathcal{X}}^{j}\big|\Big]

and observe that Ψ(𝒯)\Psi(\mathcal{T}) measures the expected number of nodes that are either (1) reached by more than ν\nu campaigns from \mathcal{I} or (2) are reached by campaign 00 from 𝒯\mathcal{T} and were reached by at least one campaign from \mathcal{I}. Recall that we had shown that Ψ\Psi is submodular and that the greedy hill-climbing algorithm leads to an approximation factor of at least 11/eϵ1-1/e-\epsilon for any ϵ>0\epsilon>0 when applied to maximizing Ψ\Psi. The following lemma whose proof can be found in Appendix D collects three statements. The first statement relates the optimum of Ψ\Psi to the optimum of Φ1\Phi^{\geq 1}. The second statement says that we loose a factor of roughly ν\nu when choosing a set of size k/ν\lfloor k/\nu\rfloor instead of kk when maximizing Ψ\Psi (this is due to submodularity). The last statement shows that a certain solution 𝒮\mathcal{S}^{\prime} for Φ1\Phi^{\geq 1} constructed from a solution 𝒯\mathcal{T} to Ψ\Psi achieves the same value.

Lemma 12.
  1. (1)

    If 𝒯kV×{0}\mathcal{T}^{k}\subseteq V\times\{0\} is a solution of size kk maximizing Ψ\Psi and 𝒮kV^\mathcal{S}^{k}\subseteq\hat{V} is a solution of size kk maximizing Φ1\Phi^{\geq 1}, then Ψ(𝒯k)Φ1(𝒮k)\Psi(\mathcal{T}^{k})\geq\Phi^{\geq 1}(\mathcal{S}^{k}).

  2. (2)

    Let ϵ>0\epsilon>0 and kν/ϵk\geq\nu/\epsilon. If 𝒯kV×{0}\mathcal{T}^{k}\subseteq V\times\{0\} is a solution of size kk maximizing Ψ\Psi and 𝒯k/νV×{0}\mathcal{T}^{\lfloor k/\nu\rfloor}\subseteq V\times\{0\} is a solution of size k/ν\lfloor k/\nu\rfloor maximizing Ψ\Psi, then Ψ(𝒯k/ν)1ϵν+1Ψ(𝒯k)\Psi(\mathcal{T}^{\lfloor k/\nu\rfloor})\geq\frac{1-\epsilon}{\nu+1}\cdot\Psi(\mathcal{T}^{k}).

  3. (3)

    Let 𝒯V×{0}\mathcal{T}\subseteq V\times\{0\} be of size k/ν\lfloor k/\nu\rfloor. Then 𝒮:={(v,j)|(v,0)𝒯,j[ν]}V^\mathcal{S}^{\prime}:=\{(v,j)|(v,0)\in\mathcal{T},j\in[\nu]\}\subseteq\hat{V} is a set of size at most kk such that Φ1(𝒮)=Ψ(𝒯)\Phi^{\geq 1}(\mathcal{S}^{\prime})=\Psi(\mathcal{T}).

Now let 0<ϵ<10<\epsilon<1, 0<δ1/20<\delta\leq 1/2, 𝒯:=\mathcal{T}:= Greedy(Ψ,ϵ/2,δ,,ν,k/ν)(\Psi,\epsilon/2,\delta,\mathcal{I},\nu,\lfloor k/\nu\rfloor), and assume that k2ν/ϵk\geq 2\nu/\epsilon. Furthermore, let 𝒮:={(v,j)|(v,0)𝒯,j[ν]}V^\mathcal{S}^{\prime}:=\{(v,j)|(v,0)\in\mathcal{T},j\in[\nu]\}\subseteq\hat{V} be as in Lemma 12, (3). Then, according to Lemma 3, it holds that Φ1(𝒮)=Ψ(𝒯)αΨ(𝒯k/ν)\Phi^{\geq 1}(\mathcal{S}^{\prime})=\Psi(\mathcal{T})\geq\alpha^{\prime}\cdot\Psi(\mathcal{T}^{\lfloor k/\nu\rfloor}) with α:=11/eϵ/2\alpha^{\prime}:=1-1/e-\epsilon/2. Using Lemma 12, (2) and (1), we obtain

Φ1(𝒮)α1ϵ2ν+1Ψ(𝒯k)αΦ1(𝒮1),\Phi^{\geq 1}(\mathcal{S}^{\prime})\geq\alpha^{\prime}\cdot\frac{1-\frac{\epsilon}{2}}{\nu+1}\cdot\Psi(\mathcal{T}^{k})\geq\alpha\cdot\Phi^{\geq 1}(\mathcal{S}^{*}_{\geq 1}),

with α:=(11/eϵ)/(ν+1)\alpha:=(1-1/e-\epsilon)/(\nu+1) and 𝒮1\mathcal{S}^{*}_{\geq 1} being an optimal solution to Φ1\Phi^{\geq 1} of size kk. Now let 𝒮\mathcal{S} be the set among 𝒮\mathcal{S}^{\prime} and \emptyset that achieves the maximum out of Φ(𝒮)\Phi(\mathcal{S}^{\prime}) and Φ()\Phi(\emptyset). Then 𝒮\mathcal{S} satisfies

2Φ(𝒮)Φ()+Φ(𝒮)Φ0(𝒮)+αΦ1(𝒮)αΦ(𝒮),2\cdot\Phi(\mathcal{S})\geq\Phi(\emptyset)+\Phi(\mathcal{S}^{\prime})\geq\Phi^{0}(\mathcal{S}^{*})+\alpha\cdot\Phi^{\geq 1}(\mathcal{S}^{*})\geq\alpha\cdot\Phi(\mathcal{S}^{*}),

where 𝒮\mathcal{S}^{*} is an optimal solution of size kk to maximizing Φ\Phi. Thus we get the following theorem.

Theorem 13.

Let 0<ϵ<10<\epsilon<1 and δ1/2\delta\leq 1/2. In the correlated setting, there is an algorithm that, with probability 1δ1-\delta, outputs a solution 𝒮\mathcal{S} that satisfies Φ(𝒮)11/eϵ2(ν+1)Φ(𝒮),\Phi(\mathcal{S})\geq\frac{1-1/e-\epsilon}{2(\nu+1)}\cdot\Phi(\mathcal{S}^{*}), where 𝒮\mathcal{S}^{*} denotes an optimal solution of size kk to maximizing Φ\Phi.

6 Conclusion and Future Works

In this paper, we introduced the μ\mu-ν\nu-Balance problem which is the generalization of the problem of balancing information exposure in a social network defined by [GGPT17]. We studied two settings called the correlated and the heterogeneous setting. While we designed an approximation algorithm with a constant approximation factor in the correlated setting, we obtained an approximation hardness result in the heterogeneous setting stating that it is unlikely to find an ng(n)n^{-g(n)}-approximation algorithm with g(n)=o(1)g(n)=o(1) if ν3\nu\geq 3 or even a nϵn^{-\epsilon}-approximation algorithm where ϵ\epsilon is a constant depending on ν\nu if ν4\nu\geq 4. In this setting, we designed an approximation algorithm with approximation ratio Ω(n1/2)\Omega(n^{-1/2}) for the case when ν=3\nu=3.

Several directions of future work are conceivable. First, it is interesting to improve the approximation guarantee for the μ\mu-ν\nu-Balance problem in both settings, most importantly for the heterogeneous case with ν>3\nu>3. Second, since the ν\nu parameter in the problem is of a threshold flavor, it would be interesting to investigate a smoother objective function by considering various ν\nu values, with different weights, such that a node reached by ν1\nu_{1} campaigns contributes more to the objective function than a node reached by ν2<ν1\nu_{2}<\nu_{1} campaigns, etc.

References

  • [AFPT10] Noga Alon, Michal Feldman, Ariel D Procaccia, and Moshe Tennenholtz. A note on competitive diffusion through social networks. Information Processing Letters, 110(6):221–225, 2010.
  • [AM11] Krzysztof R Apt and Evangelos Markakis. Diffusion in social networks with competing products. In International Symposium on Algorithmic Game Theory, pages 212–223. Springer, 2011.
  • [AMGG18] Cigdem Aslay, Antonis Matakos, Esther Galbrun, and Aristides Gionis. Maximizing the diversity of exposure in a social network. In 2018 IEEE International Conference on Data Mining (ICDM), pages 863–868. IEEE, 2018.
  • [App13] Benny Applebaum. Pseudorandom generators with long stretch and low locality from random local one-way functions. SIAM Journal on Computing, 42(5):2008–2037, 2013.
  • [BAEA11] Ceren Budak, Divyakant Agrawal, and Amr El Abbadi. Limiting the spread of misinformation in social networks. In Proceedings of the 20th international conference on World wide web, pages 665–674. ACM, 2011.
  • [BBCL14] Christian Borgs, Michael Brautbar, Jennifer Chayes, and Brendan Lucier. Maximizing social influence in nearly optimal time. In Proceedings of the Twenty-fifth Annual ACM-SIAM Symposium on Discrete Algorithms, SODA ’14, pages 946–957, 2014.
  • [BBLO17] Allan Borodin, Mark Braverman, Brendan Lucier, and Joel Oren. Strategyproof mechanisms for competitive influence in networks. Algorithmica, 78(2):425–452, 2017.
  • [BKS07] Shishir Bharathi, David Kempe, and Mahyar Salek. Competitive influence maximization in social networks. In International Workshop on Web and Internet Economics, pages 306–311. Springer, 2007.
  • [CDK+18] Eden Chlamtác, Michael Dinitz, Christian Konrad, Guy Kortsarz, and George Rabanca. The densest k-subhypergraph problem. SIAM J. Discrete Math., 32(2):1458–1477, 2018.
  • [CNWVZ07] Tim Carnes, Chandrashekhar Nagarajan, Stefan M Wild, and Anke Van Zuylen. Maximizing influence in a competitive social network: a follower’s perspective. In Proceedings of the ninth international conference on Electronic commerce, pages 351–360. ACM, 2007.
  • [CRF+11] Michael D. Conover, Jacob Ratkiewicz, Matthew R. Francisco, Bruno Gonçalves, Filippo Menczer, and Alessandro Flammini. Political polarization on twitter. In Proceedings of the Fifth International Conference on Weblogs and Social Media, Barcelona, Catalonia, Spain, July 17-21, 2011, 2011.
  • [DGDM06] Pradeep Dubey, Rahul Garg, and Bernard De Meyer. Competing for customers in a social network: The quasi-linear case. In Proceedings of 2nd International Workshop on Internet and Network Economics, pages 162–173, 2006.
  • [DOS18] Gianlorenzo D’Angelo, Martin Olsen, and Lorenzo Severini. Coverage centrality maximization in undirected networks. CoRR, abs/1811.04331, 2018. To Appear in AAAI 2019.
  • [DVBZ+16] Michela Del Vicario, Alessandro Bessi, Fabiana Zollo, Fabio Petroni, Antonio Scala, Guido Caldarelli, H. Eugene Stanley, and Walter Quattrociocchi. The spreading of misinformation online. Proceedings of the National Academy of Sciences, 113(3):554–559, 2016.
  • [GDFMGM18] Kiran Garimella, Gianmarco De Francisci Morales, Aristides Gionis, and Michael Mathioudakis. Political discourse on social media: Echo chambers, gatekeepers, and the price of bipartisanship. In Proceedings of the 2018 World Wide Web Conference, WWW ’18, pages 913–922, 2018.
  • [GGPT17] Kiran Garimella, Aristides Gionis, Nikos Parotsidis, and Nikolaj Tatti. Balancing information exposure in social networks. In Isabelle Guyon, Ulrike von Luxburg, Samy Bengio, Hanna M. Wallach, Rob Fergus, S. V. N. Vishwanathan, and Roman Garnett, editors, Advances in Neural Information Processing Systems 30: Annual Conference on Neural Information Processing Systems 2017, 4-9 December 2017, Long Beach, CA, USA, pages 4666–4674, 2017.
  • [GHK14] Sanjeev Goyal, Hoda Heidari, and Michael Kearns. Competitive contagion in networks. Games and Economic Behavior, 2014.
  • [Hoc97] Dorit S. Hochbaum, editor. Approximation Algorithms for NP-hard Problems. PWS Publishing Co., Boston, MA, USA, 1997.
  • [KKT03] David Kempe, Jon Kleinberg, and Éva Tardos. Maximizing the spread of influence through a social network. In Proceedings of the Ninth ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, KDD ’03, pages 137–146. ACM, 2003.
  • [KKT15] David Kempe, Jon M. Kleinberg, and Éva Tardos. Maximizing the spread of influence through a social network. Theory of Computing, 11:105–147, 2015.
  • [KOW08] Jan Kostka, Yvonne Anne Oswald, and Roger Wattenhofer. Word of mouth: Rumor dissemination in social networks. In International Colloquium on Structural Information and Communication Complexity, pages 185–196. Springer, 2008.
  • [LCL15] Wei Lu, Wei Chen, and Laks VS Lakshmanan. From competition to complementarity: comparative influence diffusion and maximization. Proceedings of the VLDB Endowment, 9(2):60–71, 2015.
  • [Man17] Pasin Manurangsi. Almost-polynomial ratio eth-hardness of approximating densest k-subgraph. In Proceedings of the 49th Annual ACM SIGACT Symposium on Theory of Computing, pages 954–961. ACM, 2017.
  • [McD98] Colin McDiarmid. Concentration, pages 195–248. Springer Berlin Heidelberg, 1998.
  • [MG18] Antonis Matakos and Aristides Gionis. Tell me something my friends do not know: Diversity maximization in social networks. In 2018 IEEE International Conference on Data Mining (ICDM), pages 327–336. IEEE, 2018.
  • [ML12] Seth A Myers and Jure Leskovec. Clash of the contagions: Cooperation and competition in information diffusion. In Data Mining (ICDM), 2012 IEEE 12th International Conference on, pages 539–548. IEEE, 2012.
  • [RS10] Prasad Raghavendra and David Steurer. Graph expansion and the unique games conjecture. In Proceedings of the forty-second ACM symposium on Theory of computing, pages 755–764. ACM, 2010.
  • [TAM12] Vasileios Tzoumas, Christos Amanatidis, and Evangelos Markakis. A game-theoretic analysis of a competitive diffusion process over social networks. In International Workshop on Internet and Network Economics, pages 1–14. Springer, 2012.

Appendix A Deferred Proofs for Section 2

Approximating Φ(𝒮)\Phi^{\geq\ell}(\mathcal{S}) and Ψ\Psi.

We use the following algorithm for approximating f{Ψ,Φ0,,Φν}f\in\{\Psi,\Phi^{\geq 0},\ldots,\Phi^{\geq\nu}\}.

// Note that, if f=Ψf=\Psi, then 𝒮V×{0}\mathcal{S}\subseteq V\times\{0\}, otherwise 𝒮V^\mathcal{S}\subseteq\hat{V}.
T|V|2ln(1/δ)/ϵ2T\leftarrow|V|^{2}\ln(1/\delta)/\epsilon^{2} 1
for t=1,,Tt=1,\ldots,T do 2
Sample outcome profile 𝒳\mathcal{X} 3
 if f=Ψf=\Psi then 4
   Compute RρX0(0)(𝒮)R\leftarrow\rho_{X_{0}}^{(0)}(\mathcal{S}) and nt|(Rj=1ν1V𝒳j)j=νμV𝒳j|n_{t}\leftarrow|(R\cap\bigcup_{j=1}^{\nu-1}V_{\mathcal{X}}^{j})\cup\bigcup_{j=\nu}^{\mu}V_{\mathcal{X}}^{j}| 5
    else 6
   Compute (ρ𝒳i(i)(IiSi))i[μ]\mathcal{R}\leftarrow(\rho_{\mathcal{X}_{i}}^{(i)}(I_{i}\cup S_{i}))_{i\in[\mu]} and ntNoSMμ,ν((j=01V𝒳j))n_{t}\leftarrow\NoSM(\mathcal{R}\setminus(\cup_{j=0}^{\ell-1}V^{j}_{\mathcal{X}})) 7
    return 1Tt=1Tnt\frac{1}{T}\sum_{t=1}^{T}n_{t}
Algorithm 3 approx(f,𝒮,,ν,ϵ,δ)(f,\mathcal{S},\mathcal{I},\nu,\epsilon,\delta)

We show the following lemma. As a condition for the lemma, we have the requirement that f(𝒮)1f(\mathcal{S})\geq 1 for the set 𝒮\mathcal{S} that we evaluate ff on. We argue at the end of this section, see Lemma 15 that we can assume Φ(𝒮)1\Phi^{\geq\ell}(\mathcal{S})\geq 1 for any [0,ν]\ell\in[0,\nu] and 𝒮\mathcal{S} as well as Ψ(𝒯)1\Psi(\mathcal{T})\geq 1 for any 𝒯\mathcal{T} at the cost of an arbitrarily small ϵ\epsilon in the approximation guarantee.

Lemma 14.

Let f{Ψ,Φ0,,Φν}f\in\{\Psi,\Phi^{\geq 0},\ldots,\Phi^{\geq\nu}\} and let 𝒮\mathcal{S} be such that f(𝒮)1f(\mathcal{S})\geq 1. Let f~(𝒮):=approx(f,𝒮,,ν,ϵ,δ)\tilde{f}(\mathcal{S}):={\texttt{approx}}(f,\mathcal{S},\mathcal{I},\nu,\epsilon,\delta) for some 0<δ1/20<\delta\leq 1/2 and 0<ϵ<10<\epsilon<1, then f~(𝒮)\tilde{f}(\mathcal{S}) is a (1±ϵ)(1\pm\epsilon)-approximation of f(𝒮)f(\mathcal{S}) with probability at least 1δ1-\delta.

Proof.

The proof is very similar to the proof of Proposition 4.1 in [KKT15], it is a straightforward application of a Chernoff bound, we use Theorem 2.3 from [McD98] here. Let us define TT random variables, one for each iteration of the algorithm, Y1,,YTY_{1},\ldots,Y_{T} by Yt:=nt/|V|Y_{t}:=n_{t}/|V|. Note that the YtY_{t} are independent and Yt[0,1]Y_{t}\in[0,1]. Let ST:=t=1TYtS_{T}:=\sum_{t=1}^{T}Y_{t} and μ:=E[ST]\mu:=\E[S_{T}], then ST=Tf~(𝒮)/|V|S_{T}=T\cdot\tilde{f}(\mathcal{S})/|V| and μ=Tf(𝒮)/|V|\mu=T\cdot f(\mathcal{S})/|V|. Thus, setting γ:=ϵf(𝒮)/|V|\gamma:=\epsilon f(\mathcal{S})/|V|, the Chernoff bound yields

Pr[|f(𝒮)f~(𝒮)|ϵf(𝒮)]=Pr[|STμ|Tγ]2e2Tγ2=2e2Tϵ2f(𝒮)2|V|2δ,\Pr[|f(\mathcal{S})-\tilde{f}(\mathcal{S})|\geq\epsilon f(\mathcal{S})]=\Pr[|S_{T}-\mu|\geq T\gamma]\leq 2e^{-2T\gamma^{2}}=2e^{-\frac{2T\epsilon^{2}f(\mathcal{S})^{2}}{|V|^{2}}}\leq\delta,

since T=|V|2ln(1/δ)/ϵ2T=|V|^{2}\ln(1/\delta)/\epsilon^{2} and f(𝒮)1f(\mathcal{S})\geq 1. ∎

Motivated by Lemma 1, we now investigate how the error in approximating f{Ψ,Φ1,f\in\{\Psi,\Phi^{\geq 1}, ,Φν}\ldots,\Phi^{\geq\nu}\} affects the error of the difference f(Sv)f(S)f(S\cup v)-f(S). In other words, we quantify how much we loose while maximizing ff by picking an element τ\tau with respect to an approximation of ff only: For some 0<ϵ10<\epsilon\leq 1, let f~\tilde{f} be a (1±ϵ)(1\pm\epsilon^{\prime})-approximation of ff with ϵ:=ϵ/(e(kλ(f)))\epsilon^{\prime}:=\epsilon/(e\cdot\binom{k}{\lambda(f)}), where λ(f)\lambda(f) is a constant that depends on ff, namely λ(f):=ν\lambda(f):=\nu-\ell for f=Φf=\Phi^{\geq\ell} and λ(f):=1\lambda(f):=1 for f=Ψf=\Psi. We get the following lemma.

See 2

Proof.

We distinguish two cases. First, assume that f(𝒮τ)f(𝒮)f(𝒮)/(e(kν))f(\mathcal{S}\cup\tau^{*})-f(\mathcal{S})\leq f(\mathcal{S}^{*})/(e\cdot\binom{k}{\nu-\ell}). If f=Φf=\Phi^{\ell} for some \ell, Lemma 1 yields that τ\tau^{*} satisfies

f(𝒮τ)f(𝒮)1(kλ(f))(f(𝒮)f(𝒮)).\displaystyle f(\mathcal{S}\cup\tau^{*})-f(\mathcal{S})\geq\frac{1}{\binom{k}{\lambda(f)}}\cdot(f(\mathcal{S}^{*})-f(\mathcal{S})). (2)

For f=Ψf=\Psi, we note that λ(f)=1\lambda(f)=1 and thus (kλ(f))=k\binom{k}{\lambda(f)}=k, so we get inequality (2) by the submodularity of Ψ\Psi. Thus, in both case by combining the two inequalities, we get f(𝒮)(11/e)f(𝒮)f(\mathcal{S})\geq(1-1/e)\cdot f(\mathcal{S}^{*}), which concludes this case.

On the other hand, assume f(𝒮τ)f(𝒮)>f(𝒮)/(e(kλ(f)))f(\mathcal{S}\cup\tau^{*})-f(\mathcal{S})>f(\mathcal{S}^{*})/(e\cdot\binom{k}{\lambda(f)}). Using the approximation guarantee of f~\tilde{f}, the definition of τ~\tilde{\tau}, and again the approximation guarantee, we get

f(𝒮CLOSE\displaystyle f(\mathcal{S}\cup OPENτ~)f(𝒮)1ϵ1+ϵf(𝒮τ)f(𝒮)\displaystyle\tilde{\tau})-f(\mathcal{S})\geq\frac{1-\epsilon^{\prime}}{1+\epsilon^{\prime}}f(\mathcal{S}\cup\tau^{*})-f(\mathcal{S})
=f(𝒮τ)f(𝒮)2ϵ1+ϵf(𝒮τ)\displaystyle=f(\mathcal{S}\cup\tau^{*})-f(\mathcal{S})-\frac{2\epsilon^{\prime}}{1+\epsilon^{\prime}}f(\mathcal{S}\cup\tau^{*})
=(1ϵ)(f(𝒮τ)f(𝒮))+ϵ(1+ϵ)2ϵ1+ϵf(𝒮τ)ϵf(𝒮).\displaystyle=(1-\epsilon)\cdot(f(\mathcal{S}\cup\tau^{*})-f(\mathcal{S}))+\frac{\epsilon(1+\epsilon^{\prime})-2\epsilon^{\prime}}{1+\epsilon^{\prime}}\cdot f(\mathcal{S}\cup\tau^{*})-\epsilon\cdot f(\mathcal{S}).

Thus, it remains to argue that the latter two summands are non-negative. From the case assumption and the optimality of 𝒮\mathcal{S}^{*}, we have f(𝒮)f(𝒮τ)/(1+1/(e(kλ(f))))f(\mathcal{S})\leq f(\mathcal{S}\cup\tau^{*})/(1+1/(e\cdot\binom{k}{\lambda(f)})) and thus the above latter two summands can be lower bounded by

f(Sτ)(ϵ(1+ϵ)2ϵ1+ϵϵ(1+1/(e(kλ(f)))CLOSE).\displaystyle f(S\cup\tau^{*})\cdot\Big(\frac{\epsilon(1+\epsilon^{\prime})-2\epsilon^{\prime}}{1+\epsilon^{\prime}}-\frac{\epsilon}{(1+1/(e\cdot\binom{k}{\lambda(f)})}\Big).

The latter is non-negative, since

(ϵ(1+ϵ)2ϵ)(1+1e(kλ(f)))ϵ(1+ϵ)=ϵ(1+ϵ)2ϵe(kλ(f))0(\epsilon(1+\epsilon^{\prime})-2\epsilon^{\prime})\Big(1+\frac{1}{e\cdot\binom{k}{\lambda(f)}}\Big)-\epsilon(1+\epsilon^{\prime})=\epsilon^{\prime}(1+\epsilon^{\prime})-\frac{2\epsilon^{\prime}}{e\cdot\binom{k}{\lambda(f)}}\geq 0

by the choice of ϵ:=ϵ/(e(kλ(f)))\epsilon^{\prime}:=\epsilon/(e\cdot\binom{k}{\lambda(f)}) and 0<ϵ10<\epsilon\leq 1. This concludes the proof. ∎

Lower Bound on Φ\Phi^{\geq\ell} and Ψ\Psi.

Our goal in this section is to argue that there is a transform τ\tau that takes an instance (G,𝒫,,k)(G,\mathcal{P},\mathcal{I},k) of the μ\mu-ν\nu-Balance problem and outputs a (slightly modified) instance (G,𝒫,,k):=τ(G,𝒫,,k)(G^{\prime},\mathcal{P},\mathcal{I}^{\prime},k):=\tau(G,\mathcal{P},\mathcal{I},k) such that the function f(𝒮)f^{\prime}(\mathcal{S}) is at least 11 for every argument 𝒮\mathcal{S} in the new instance for any f{Φ0,,Φν,Ψ}f\in\{\Phi^{\geq 0},\ldots,\Phi^{\geq\nu},\Psi\}. Moreover, given an approximation algorithm for Φ\Phi with approximation ratio α\alpha, we will show that applying this algorithm on the transformed instance τ(G,𝒫,,k)\tau(G,\mathcal{P},\mathcal{I},k) leads to a solution of approximation ratio at least αϵ\alpha-\epsilon for the original instance, for any ϵ>0\epsilon>0.

The transform τ\tau is defined as follows. Obtain GG^{\prime} by adding an isolated node vv to GG and extend \mathcal{I} to \mathcal{I}^{\prime} by adding vv to IiI_{i} for every i[ν]i\in[\nu]. Now clearly, for every solution 𝒮\mathcal{S}, it holds that f(𝒮)=f(𝒮)+11f^{\prime}(\mathcal{S})=f(\mathcal{S})+1\geq 1, where the 1 originates from the additional node vv that is initially covered by ν\nu\geq\ell campaigns. Moreover, we get the following lemma.

Lemma 15.

Let ϵ>0\epsilon>0. Then, for instances (G,𝒫,,k)(G,\mathcal{P},\mathcal{I},k) with k2ν/ϵk\geq 2\nu/\epsilon, the following holds: Let 𝒮\mathcal{S}^{\prime} be a solution in (G,𝒫,,k):=τ(G,𝒫,,k)(G^{\prime},\mathcal{P},\mathcal{I}^{\prime},k):=\tau(G,\mathcal{P},\mathcal{I},k) such that Φ(𝒮)αΦ(𝒮)\Phi^{\prime}(\mathcal{S}^{\prime})\geq\alpha\cdot\Phi^{\prime}(\mathcal{S}^{\prime*}), where 𝒮\mathcal{S}^{\prime*} denotes an optimal solution for maximizing Φ\Phi^{\prime} in the new instance (G,𝒫,,k)(G^{\prime},\mathcal{P},\mathcal{I}^{\prime},k). Then 𝒮:=𝒮{v}\mathcal{S}:=\mathcal{S}^{\prime}\setminus\{v\} satisfies Φ(𝒮)(αϵ)Φ(𝒮)\Phi(\mathcal{S})\geq(\alpha-\epsilon)\cdot\Phi(\mathcal{S}^{*}), where 𝒮\mathcal{S}^{*} denotes an optimal solution for maximizing Φ\Phi in the original instance (G,𝒫,,k)(G,\mathcal{P},\mathcal{I},k).

Proof.

First note that Φ(𝒮)k/νk/ν11/ϵ\Phi(\mathcal{S}^{*})\geq\lfloor k/\nu\rfloor\geq k/\nu-1\geq 1/\epsilon or equivalently 1ϵΦ(𝒮)1\leq\epsilon\Phi(\mathcal{S}^{*}). This yields the claim, since

Φ(𝒮)=Φ(𝒮)1αΦ(𝒮)1αΦ(𝒮)1(αϵ)Φ(𝒮).\Phi(\mathcal{S})=\Phi^{\prime}(\mathcal{S})-1\geq\alpha\cdot\Phi^{\prime}(\mathcal{S}^{\prime*})-1\geq\alpha\cdot\Phi(\mathcal{S}^{*})-1\geq(\alpha-\epsilon)\cdot\Phi^{\geq\ell}(\mathcal{S}^{*}).\qed
Maximizing Φν1\Phi^{\geq\nu-1} and Ψ\Psi.

Our goal here is to show that the standard greedy hill climbing algorithm, we refer to it as Greedy(f,ϵ,δ,,ν,k)(f,\epsilon,\delta,\mathcal{I},\nu,k), can be applied in order to approximate both Φν1\Phi^{\geq\nu-1} and Ψ\Psi to within a factor of 11/eϵ1-1/e-\epsilon for any 0<ϵ<10<\epsilon<1 with probability at least 1δ1-\delta for any 0<δ1/20<\delta\leq 1/2. We first formally prove that these functions are submodular.

Lemma 16.

The functions Ψ\Psi and Φν1\Phi^{\geq\nu-1} are monotone and submodular.

Proof.

The monotonicity of Ψ\Psi and Φν1\Phi^{\geq\nu-1} is straightforward. We argue the submodularity of Ψ\Psi (Φν1\Phi^{\geq\nu-1}) in a similar way as we argued in the proof of Lemma 1. To this end, let D(Ψ)=V×{0}D(\Psi)=V\times\{0\} and D(Φν1)=V^=V×[μ]D(\Phi^{\geq\nu-1})=\hat{V}=V\times[\mu] denote the domain of Ψ\Psi and Φν1\Phi^{\geq\nu-1}, respectively, and let 𝒮\mathcal{S} and 𝒮\mathcal{S}^{\prime} be subsets of D(Ψ)D(\Psi) (D(Φν1)D(\Phi^{\geq\nu-1})) such that 𝒮𝒮\mathcal{S}\subseteq\mathcal{S}^{\prime}, and let τ\tau be an element of the domain D(Ψ)D(\Psi) (D(Φν1)D(\Phi^{\geq\nu-1})). Furthermore, let 𝒳\mathcal{X} be an outcome profile w.r.t. the correlated (heterogeneous) probability distributions. Lastly, let vVV𝒳0v\in V\setminus V_{\mathcal{X}}^{0} (vVj=0ν2V𝒳jv\in V\setminus\bigcup_{j=0}^{\nu-2}V_{\mathcal{X}}^{j}) be a node that can contribute to the value of Ψ\Psi (Φν1\Phi^{\geq\nu-1}).22 2 Recall that V𝒳jV_{\mathcal{X}}^{j} is the set of nodes that was reached by jj campaigns from seed sets \mathcal{I}. We denote by 𝟙𝒳𝒮(v)\mathds{1}_{\mathcal{X}}^{\mathcal{S}}(v) the indicator function that is 1 if vv contributes to Ψ\Psi (Φν1\Phi^{\geq\nu-1}) in outcome profile 𝒳\mathcal{X} with initial seed sets \mathcal{I} and additional seed sets 𝒮\mathcal{S} and 0 otherwise. We now argue that the following inequality holds:.

𝟙𝒳𝒮τ(v)𝟙𝒳𝒮(v)𝟙𝒳𝒮τ(v)𝟙𝒳𝒮(v)\displaystyle\mathds{1}_{\mathcal{X}}^{\mathcal{S}^{\prime}\cup\tau}(v)-\mathds{1}_{\mathcal{X}}^{\mathcal{S}^{\prime}}(v)\leq\mathds{1}_{\mathcal{X}}^{\mathcal{S}\cup\tau}(v)-\mathds{1}_{\mathcal{X}}^{\mathcal{S}}(v) (3)

Note that the right-hand side cannot be negative by monotonicity and that, if the left-hand side is positive for Ψ\Psi (Φν1\Phi^{\geq\nu-1}), then it must hold that the node vv is reached by a subset M[μ]M\subseteq[\mu] of campaigns from \mathcal{I} with |M|[1,ν1]|M|\in[1,\nu-1] (|M|=ν1|M|=\nu-1). Furthermore, the node vv is not reached by campaign 00 (is not reached by a campaign j[μ]Mj\in[\mu]\setminus M) from 𝒮\mathcal{S}^{\prime}, but it is reached by campaign 00 (campaign jj) from τ\tau. Now, observe that 𝒮𝒮\mathcal{S}\subseteq\mathcal{S}^{\prime} and thus the node vv is not reached by campaign 0 (by campaign jj) from 𝒮\mathcal{S} neither. Hence it follows that the right-hand side is also positive. Taking the expected value on both sides of (3) yields that Ψ(𝒮τ)Ψ(𝒮)Ψ(𝒮τ)Ψ(𝒮)\Psi(\mathcal{S}^{\prime}\cup\tau)-\Psi(\mathcal{S}^{\prime})\leq\Psi(\mathcal{S}\cup\tau)-\Psi(\mathcal{S}) (OPENΦν1(𝒮τ)Φν1(𝒮)Φν1(𝒮τ)Φν1(𝒮))\Phi^{\geq\nu-1}(\mathcal{S}^{\prime}\cup\tau)-\Phi^{\geq\nu-1}(\mathcal{S}^{\prime})\leq\Phi^{\geq\nu-1}(\mathcal{S}\cup\tau)-\Phi^{\geq\nu-1}(\mathcal{S})) due to linearity of expectation. This establishes submodularity and concludes the proof. ∎

We now recall the following classical result concerning the greedy algorithm for maximizing a submodular function:

Lemma 17 (Theorem 3.9 in [Hoc97]).

The greedy hill-climbing algorithm, that at each step picks an element that leads to an increment being within factor β\beta of the optimal increment possible, achieves an approximation ratio of at least 1(1β/k)k>11/eβ1-(1-\beta/k)^{k}>1-1/e^{\beta}.

We have seen that both Φν1\Phi^{\geq\nu-1} and Ψ\Psi can be approximated within a (1±ϵ)(1\pm\epsilon)-factor using the approx-routine. In Lemma 2 we argued that using the approximations we can find an element vv (or a set τ\tau of cardinality λ(f)=1\lambda(f)=1) that when added to 𝒮\mathcal{S} leads to a progress of at least a factor of (1ϵ)(1-\epsilon) of the maximal progress possible. We prove Lemma 3.

δδ/(k|V^|)\delta^{\prime}\leftarrow\delta/(k|\hat{V}|), ϵϵ/(ek)\epsilon^{\prime}\leftarrow\epsilon/(ek), 𝒮\mathcal{S}\leftarrow\emptyset 1
while |𝒮|k|\mathcal{S}|\leq k do 2
Compute v𝚊𝚛𝚐max{approx(f,𝒮{v},,ν,ϵ,δ):vDf}v\leftarrow\mathtt{arg}\max\{\texttt{approx}(f,\mathcal{S}\cup\{v\},\mathcal{I},\nu,\epsilon^{\prime},\delta^{\prime}):v\in D_{f}\}, set 𝒮𝒮{v}\mathcal{S}\leftarrow\mathcal{S}\cup\{v\} 3
 return 𝒮\mathcal{S}
Algorithm 4 Greedy(f,ϵ,δ,,ν,k)(f,\epsilon,\delta,\mathcal{I},\nu,k)

See 3

Proof.

The union bound over all at most k|V^|k|\hat{V}| calls to approx, yields that, with probability at least 1δ1-\delta, each call resulted in a 1±ϵ1\pm\epsilon^{\prime}-approximation. Then Lemma 2 applied to ff guarantees that after each iteration ii either an element vv is picked such that the increment using vv is at least a (1ϵ)(1-\epsilon)-fraction of the optimal increment possible in this iteration or the current set 𝒮i\mathcal{S}^{i} is already a (11/e)(1-1/e)-approximation of the optimum set 𝒮ν1\mathcal{S}^{*}_{\nu-1}. In the latter case the lemma is fulfilled by the monotonicity of ff. In the former case we get an 𝒮\mathcal{S} having an approximation ratio of at least 1(1(1ϵ)/k)k11/e1ϵ1-(1-(1-\epsilon)/k)^{k}\geq 1-1/e^{1-\epsilon} according to Lemma 17. Since 11e1ϵ(1ϵ)(11e)11eϵ,1-\frac{1}{e^{1-\epsilon}}\geq(1-\epsilon)\cdot\big(1-\frac{1}{e}\big)\geq 1-\frac{1}{e}-\epsilon, this concludes the proof. ∎

Appendix B Deferred Proofs for Section 3

Further illustration of the reduction described in Section 3.

P[3],(2,3,1)P_{[3],(2,3,1)}P[3],(3,1,2)P_{[3],(3,1,2)}P[3],(3,2,1)P_{[3],(3,2,1)}P[3],(1,2,3)P_{[3],(1,2,3)}P[3],(1,3,2)P_{[3],(1,3,2)}P[3],(2,1,3)P_{[3],(2,1,3)}uuvvwwp2=1p_{2}=1p3=1p_{3}=1p1=1p_{1}=1p3=1p_{3}=1p1=1p_{1}=1p2=1p_{2}=1p3=1p_{3}=1p2=1p_{2}=1p1=1p_{1}=1p1=1p_{1}=1p2=1p_{2}=1p3=1p_{3}=1p1=1p_{1}=1p3=1p_{3}=1p2=1p_{2}=1p2=1p_{2}=1p1=1p_{1}=1p3=1p_{3}=1
Figure 2: For a set ι={i,j,k}J\iota=\{i,j,k\}\in J of dd campaigns and a permutation πSd\pi\in S_{d}, let Pι,πP_{\iota,\pi} stand for the path in Figure 1 of nodes eι,π1,,eι,πle_{\iota,\pi}^{1},\ldots,e_{\iota,\pi}^{l} connected by arcs (eι,πt,eι,πt+1)(e_{\iota,\pi}^{t},e_{\iota,\pi}^{t+1}) for t=1,,l1t=1,\ldots,l-1 with probabilities on these edges being one for pπ(i)p_{\pi(i)}, pπ(j)p_{\pi(j)}, and pπ(k)p_{\pi(k)} and zero for all other indices. The figure illustrates the case d=3d=3 and μ=ν=4\mu=\nu=4 and the portion of the network that is generated in the transform τ\tau of Section3 for one hyper-edge e={u,v,w}e=\{u,v,w\} and the only set ι=[3]J=(33)\iota=[3]\in J=\binom{3}{3}. Probabilities that are not given are equal to zero.

Figure 2 illustrates the scheme induced by an hyperedge e=(u,v,w)e=(u,v,w) when d=3d=3 and μ=ν=4\mu=\nu=4. In this case, J=([μν+d]d)J=\binom{[\mu-\nu+d]}{d} is only composed of set [3][3] and S3S_{3} is composed of 6 permutations. We use the standard tuple notation for permutations.

Deferred Proofs for Reducing Densest-kk-Sub-dd-hypergraph to μ\mu-ν\nu-Balance.

We start by defining the Multicld-Edge Densest-Sub-dd-hypergraph problem which is closely related to the Densest-kk-Sub-dd-hypergraph problem.

Multicld-Edge Densest-Sub-dd-hypergraph Input: dd-Regular Hypergraph G=(V,E)G=(V,E), integer kk Find: set SVS\subseteq V with |S|k|S|\leq k and a coloring function φ:S[d]\varphi:S\rightarrow[d], s.t. |Eφ(S)||E_{\varphi}(S)| is maximal, where Eφ(S):={eE:e=(v1,,vd)Sφ(vi)φ(vj),ij}.E_{\varphi}(S):=\{e\in E:e=(v_{1},\ldots,v_{d})\subseteq S~\wedge~\varphi(v_{i})\neq\varphi(v_{j}),\forall i\neq j\}.

Problem Multicld-Edge Densest-Sub-dd-hypergraph will be of interest to us due to the following results. We first prove a lemma showing the existence of an assignment φ\varphi^{\prime} such that at least a fraction pp of the hyperedges in the induced sub-hypergraph of a set SS have differently colored endpoints.

Lemma 18.

Let G=(V,E)G=(V,E) be a dd-regular hypergraph. For any set SVS\subseteq V, there is an assignment φ:S[d]\varphi^{\prime}:S\rightarrow[d] s.t. |Eφ(S)|p|E(S)||E_{\varphi^{\prime}}(S)|\geq p\cdot|E(S)| where p=d!/ddp=d!/d^{d}.

Proof.

Let SVS\subseteq V. Consider the probabilistic procedure in which, for each node, we assign a color from [d][d] uniformly at random and independently of the other nodes. This procedures yields a coloring φ\varphi. For any e=(v1,,vd)e=(v_{1},\ldots,v_{d}) in SS, the probability that φ(vi)φ(vj)\varphi(v_{i})\neq\varphi(v_{j}) for all iji\neq j is pp. This property is guaranteed if and only if (φ(v1),,φ(vd))(\varphi(v_{1}),\ldots,\varphi(v_{d})) corresponds to one of the d!d! permutation of [d][d]. In total, there are ddd^{d} ways of coloring ee. Hence, the expected value of |Eφ(S)||E_{\varphi}(S)| is p|E(S)|p\cdot|E(S)|. Consequently, the function φ\varphi^{\prime} that maximizes |Eφ(S)||E_{\varphi^{\prime}}(S)| satisfies |Eφ(S)|p|E(S)||E_{\varphi^{\prime}}(S)|\geq p\cdot|E(S)|. ∎

This leads to the following corollary.

Corollary 19.

Denoting with DKSHd(G,k)\dkshop^{*}_{d}(G,k) and MCDd(G,k)\mcdshop^{*}_{d}(G,k) the value of the optimal solution for Densest-kk-Sub-dd-hypergraph and Multicld-Edge Densest-Sub-dd-hypergraph on (G,k)(G,k), respectively, we have that DKSHd(G,k)MCDp(G,k)/p\dkshop^{*}_{d}(G,k)\leq\mcdshop^{*}_{p}(G,k)/p, where p=d!/ddp=d!/d^{d}.

Proof.

Let SS be a set that achieves DKSHd(G,k)=|E(S)|\dkshop^{*}_{d}(G,k)=|E(S)|, then DKSHd(G,k)=|E(S)||Eφ(S)|/pMCDd(G,k)/p\dkshop^{*}_{d}(G,k)=|E(S)|\leq|E_{\varphi^{\prime}}(S)|/p\leq\mcdshop^{*}_{d}(G,k)/p, where φ\varphi^{\prime} is as in Lemma 18. ∎

Recall that we fixed a μ\mu-ν\nu-Balance instance P=(G¯=(V¯,A¯),𝒫,,k¯)P=(\overline{G}=(\overline{V},\overline{A}),\mathcal{P},\mathcal{I},\overline{k}) resulting from the transform τ\tau as image of an Densest-kk-Sub-dd-hypergraph instance Q=(G=(V,E),k)Q=(G=(V,E),k). In what follows nodes in VV_{\boxempty} (resp. VV_{\ocircle}) are called rectangle-nodes (resp. circle-nodes). See 4

Proof.
  1. 1.

    We can w.l.o.g. assume that 𝒮V=\mathcal{S}^{*}\cap V_{\ocircle}=\emptyset and 𝒮V=\mathcal{S}_{\ocircle}^{*}\cap V_{\ocircle}=\emptyset. Then, it follows that both Φ(𝒮)\Phi_{\ocircle}(\mathcal{S}^{*}) and Φ(𝒮)\Phi_{\ocircle}(\mathcal{S}_{\ocircle}^{*}) are multiples of ll. Now, assume for the purpose of contradiction that Φ(𝒮)>Φ(𝒮)\Phi_{\ocircle}(\mathcal{S}_{\ocircle}^{*})>\Phi_{\ocircle}(\mathcal{S}^{*}). Then, Φ(𝒮)Φ(𝒮)+l\Phi_{\ocircle}(\mathcal{S}^{*}_{\ocircle})\geq\Phi_{\ocircle}(\mathcal{S}^{*})+l which leads to

    Φ(𝒮)=Φ(𝒮)+Φ(𝒮)Φ(𝒮)l+|V|<Φ(𝒮),\Phi(\mathcal{S}^{*})=\Phi_{\ocircle}(\mathcal{S}^{*})+\Phi_{\boxempty}(\mathcal{S}^{*})\leq\Phi_{\ocircle}(\mathcal{S}^{*}_{\ocircle})-l+|V|<\Phi(\mathcal{S}^{*}_{\ocircle}),

    using that l>|V|l>|V|. This is a contradiction to 𝒮\mathcal{S}^{*} being optimal.

  2. 2.

    Let (S,φ)(S^{*},\varphi^{*}) be an optimal solution to the Multicld-Edge Densest-Sub-dd-hypergraph problem induced by QQ. Construct a solution 𝒮\mathcal{S} for μ\mu-ν\nu-Balance by letting Si:={vV:φ(v)=i},i[d]S_{i}:=\{v\in V:\varphi(v)=i\},\forall i\in[d]. Clearly Φ(𝒮)=l|Eφ(S)|\Phi_{\ocircle}(\mathcal{S})=l|E_{\varphi^{*}}(S^{*})|. Thus, using Corollary 19: Φ(𝒮)lMCDdlpDKSHd.\Phi_{\ocircle}(\mathcal{S}^{*}_{\ocircle})\geq l\cdot\mcdshop^{*}_{d}\geq l\cdot p\cdot\dkshop^{*}_{d}.

  3. 3.

    Let S={vV:vSi for some i[μν+d]}V=VS=\{v\in V_{\boxempty}:v\in S_{i}\text{ for some }i\in[\mu-\nu+d]\}\subseteq V_{\boxempty}=V be the set of rectangle-nodes where 𝒮\mathcal{S} propagates at least one campaign in [μν+d][\mu-\nu+d]. Clearly, |S|k|S|\leq k. Let q=|E(S)|q=|E(S)| be the number of edges in the sub-graph of GG induced by SS. Then, Φ(𝒮)λlq\Phi_{\ocircle}(\mathcal{S})\leq\lambda lq, since each edge in GG can count for λl\lambda l circle-nodes if the dd corresponding rectangle-nodes propagate all campaigns in [μν+d][\mu-\nu+d]. It follows that |E(S)|Φ(𝒮)/(λl)|E(S)|\geq\Phi_{\ocircle}(\mathcal{S})/(\lambda l).∎

Appendix C Deferred Proofs for Section 4

C.1 Deferred Proofs for the Analysis of Algorithm GreedyTuple

The aim of this section is to prove the following Lemma. See 7

For this purpose, we will first prove the following lemma.

Lemma 20.

Let 0<ϵ<10<\epsilon<1, δ1/2\delta\leq 1/2, and [1,ν1]\ell\in[1,\nu-1]. With probability at least 1δ1-\delta, after each iteration ii of GreedyTuple(ϵ,δ,,,ν,k)(\epsilon,\delta,\ell,\mathcal{I},\nu,k), it either holds that

Φ(𝒮i)Φ(𝒮i1)1ϵ2(kν)(Φ(𝒮)Φ(𝒮i1)). or Φ(𝒮i)(11e)Φ(𝒮).\Phi^{\geq\ell}(\mathcal{S}^{i})-\Phi^{\geq\ell}(\mathcal{S}^{i-1})\geq\frac{1-\frac{\epsilon}{2}}{\binom{k}{\nu-\ell}}\cdot(\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*})-\Phi^{\geq\ell}(\mathcal{S}^{i-1})).\quad\text{ or }\quad\Phi^{\geq\ell}(\mathcal{S}^{i})\geq\big(1-\frac{1}{e}\big)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}).
Proof.

Algorithm GreedyTuple(ϵ,δ,,,ν,k)(\epsilon,\delta,\ell,\mathcal{I},\nu,k) calls algorithm approx at most tt times. Let us call EiE_{i} the event that the ii’th call to approx “succeeds”, i.e., that the call results in 1±ϵ1\pm\epsilon^{\prime}-approximation Φ~(𝒯)\tilde{\Phi}^{\geq\ell}(\mathcal{T}). That is, it holds that (1ϵ)Φ(𝒯)Φ~(𝒯)(1+ϵ)Φ1(𝒯)(1-\epsilon^{\prime})\Phi^{\geq\ell}(\mathcal{T})\leq\tilde{\Phi}^{\geq\ell}(\mathcal{T})\leq(1+\epsilon^{\prime})\Phi^{\geq 1}(\mathcal{T}). This event happens with probability at least 1δ=1δ/t1-\delta^{\prime}=1-\delta/t. Since there are at most tt many evaluations, using the union bound, we obtain that the probability that all evaluations succeed is at least 1δ1-\delta. Now the statement follows with Lemma 2. It states that either Φ(𝒮i)(11e)Φ(𝒮)\Phi^{\geq\ell}(\mathcal{S}^{i})\geq\big(1-\frac{1}{e}\big)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}) or, for the element τ\tau picked by the algorithm, it holds that Φ(𝒮i1τ)Φ(𝒮i1)(1ϵ2)/(kν)(Φ(𝒮)Φ(𝒮i1))\Phi^{\geq\ell}(\mathcal{S}^{i-1}\cup\tau)-\Phi^{\geq\ell}(\mathcal{S}^{i-1})\geq(1-\frac{\epsilon}{2})/\binom{k}{\nu-\ell}\cdot(\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*})-\Phi^{\geq\ell}(\mathcal{S}^{i-1})) using Lemma 1. ∎

We can now prove Lemma 7.

Proof of Lemma 7.

We show the statement by induction. For i=1i=1, we note that by Lemma 20, we either have Φ(𝒮1)Φ(𝒮0)(1ϵ2)/(kν)(Φ(𝒮)Φ(𝒮0))\Phi^{\geq\ell}(\mathcal{S}^{1})-\Phi^{\geq\ell}(\mathcal{S}^{0})\geq(1-\frac{\epsilon}{2})/\binom{k}{\nu-\ell}\cdot(\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*})-\Phi^{\geq\ell}(\mathcal{S}^{0})) or Φ(𝒮i)(11e)Φ(𝒮)\Phi^{\geq\ell}(\mathcal{S}^{i})\geq\big(1-\frac{1}{e}\big)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}). In the latter case the statement holds, in the former case, we get Φ(𝒮1)(1ϵ2)/(kν)Φ(𝒮)\Phi^{\geq\ell}(\mathcal{S}^{1})\geq(1-\frac{\epsilon}{2})/\binom{k}{\nu-\ell}\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}) and thus the statement follows in both cases. For i>1i>1, let us assume that the statement holds after iteration i1i-1. If Φ(𝒮i1)(11/e)Φ(𝒮)\Phi^{\geq\ell}(\mathcal{S}^{i-1})\geq(1-1/e)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}), we have Φ(𝒮i)(11/e)Φ(𝒮)\Phi^{\geq\ell}(\mathcal{S}^{i})\geq(1-1/e)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}) by monotonicity. In the other case, we have that

Φ(𝒮i1)(1(11ϵ2(kν))i1)Φ(𝒮).\displaystyle\Phi^{\geq\ell}(\mathcal{S}^{i-1})\geq\Big(1-\Big(1-\frac{1-\frac{\epsilon}{2}}{\binom{k}{\nu-\ell}}\Big)^{i-1}\Big)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}). (4)

Applying Lemma 20 yields that either Φ(𝒮i)(11/e)Φ(𝒮)\Phi^{\geq\ell}(\mathcal{S}^{i})\geq(1-1/e)\cdot\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}), in which case the statement holds, or we obtain

Φ(𝒮i)\displaystyle\!\!\Phi^{\geq\ell}(\mathcal{S}^{i}) =Φ(𝒮i1)+(Φ(𝒮i)Φ(𝒮i1))(ϵ2(kν1))Φ(𝒮i1)+ϵ2(kν1)Φ(𝒮).\displaystyle\!=\!\Phi^{\geq\ell}(\mathcal{S}^{i-1})\!+\!(\Phi^{\geq\ell}(\mathcal{S}^{i})\!-\!\Phi^{\geq\ell}(\mathcal{S}^{i-1}))\!\geq\!\Big(1\!-\!\frac{1\!-\!\frac{\epsilon}{2}}{\binom{k}{\nu-1}}\Big)\Phi^{\geq\ell}(\mathcal{S}^{i-1})\!+\!\frac{1\!-\!\frac{\epsilon}{2}}{\binom{k}{\nu-1}}\Phi^{\geq\ell}(\mathcal{S}_{\geq\ell}^{*}).

Applying (4) yields the claim. ∎

C.2 Deferred Proofs for the Analysis of Algorithm GreedyIter

In this section we prove lemmata 9 and 10 which are paramount in proving the approximation ratio of Algorithm GreedyIter. See 9

Proof.

Let 𝒮k/(ν1)\mathcal{S}_{\lfloor k/(\nu-1)\rfloor}^{*} and 𝒮k\mathcal{S}_{k}^{*} be sets of cardinality k/(ν1)\lfloor k/(\nu-1)\rfloor and kk, respectively, maximizing Φ21(,)\Phi^{\geq 1}_{2}(\mathcal{I},\cdot). Furthermore, let 𝒯\mathcal{T} be a subset of 𝒮k\mathcal{S}_{k}^{*} of cardinality k/(ν1)\lfloor k/(\nu-1)\rfloor that maximizes Φ21(,)\Phi^{\geq 1}_{2}(\mathcal{I},\cdot). Lemma 3 yields that for ϵ=ϵ/2\epsilon^{\prime}=\epsilon/2, with probability at least 1δ/ν1-\delta/\nu, we have that

Φ21(,𝒮[1])(11eϵ)Φ21(,𝒮k/(ν1))(11eϵ)Φ21(,𝒯).\displaystyle\Phi^{\geq 1}_{2}(\mathcal{I},\mathcal{S}^{[1]})\geq\big(1-\frac{1}{e}-\epsilon^{\prime}\big)\cdot\Phi^{\geq 1}_{2}(\mathcal{I},\mathcal{S}_{\lfloor k/(\nu-1)\rfloor}^{*})\geq\big(1-\frac{1}{e}-\epsilon^{\prime}\big)\cdot\Phi^{\geq 1}_{2}(\mathcal{I},\mathcal{T}). (5)

Using the submodularity and monotonicity of Φ21(,)\Phi^{\geq 1}_{2}(\mathcal{I},\cdot) and the maximum choice of 𝒯\mathcal{T} yields

Φ21(,𝒮k)kk/(ν1)Φ21(,𝒯)ν1ϵΦ21(,𝒯),\displaystyle\Phi^{\geq 1}_{2}(\mathcal{I},\mathcal{S}_{k}^{*})\leq\left\lceil\frac{k}{\lfloor k/(\nu-1)\rfloor}\right\rceil\cdot\Phi^{\geq 1}_{2}(\mathcal{I},\mathcal{T})\leq\frac{\nu}{1-\epsilon^{\prime}}\cdot\Phi^{\geq 1}_{2}(\mathcal{I},\mathcal{T}), (6)

as k(ν1)/ϵk\geq(\nu-1)/\epsilon^{\prime} implies kν+1(1ϵ)kk-\nu+1\geq(1-\epsilon^{\prime})k and thus kk/(ν1)k(ν1)kν+1ν11ϵν1ϵ\lceil\frac{k}{\lfloor k/(\nu-1)\rfloor}\rceil\leq\lceil\frac{k(\nu-1)}{k-\nu+1}\rceil\leq\lceil\frac{\nu-1}{1-\epsilon^{\prime}}\rceil\leq\frac{\nu}{1-\epsilon^{\prime}}. By combining the estimates from (5) and (6), we obtain

Φ21(,𝒮[1])(11eϵ)(1ϵ)νΦ21(,𝒮k)(11eϵ)νΦν1(,𝒮1),\displaystyle\Phi^{\geq 1}_{2}(\mathcal{I},\mathcal{S}^{[1]})\geq\frac{(1-\frac{1}{e}-\epsilon^{\prime})(1-\epsilon^{\prime})}{\nu}\cdot\Phi^{\geq 1}_{2}(\mathcal{I},\mathcal{S}_{k}^{*})\geq\frac{(1-\frac{1}{e}-\epsilon)}{\nu}\cdot\Phi^{\geq 1}_{\nu}(\mathcal{I},\mathcal{S}_{\geq 1}^{*}),

where the last step uses that x(1ϵ)xϵx(1-\epsilon^{\prime})\geq x-\epsilon^{\prime} for any x1x\leq 1, the definition of ϵ=ϵ/2\epsilon^{\prime}=\epsilon/2, and the fact that 𝒮k\mathcal{S}_{k}^{*} and 𝒮1\mathcal{S}_{\geq 1}^{*} are both of size kk and thus Φ21(,𝒮k)Φ21(,𝒮1)Φν1(,𝒮1)\Phi^{\geq 1}_{2}(\mathcal{I},\mathcal{S}_{k}^{*})\!\geq\!\Phi^{\geq 1}_{2}(\mathcal{I},\mathcal{S}_{\geq 1}^{*})\!\geq\!\Phi^{\geq 1}_{\nu}(\mathcal{I},\mathcal{S}_{\geq 1}^{*}). ∎

See 10

Proof.

We use the shorthand Φ[]():=Φ+1([],)\Phi^{[\ell]}(\cdot)\!:=\!\Phi^{\geq\ell}_{\ell+1}(\mathcal{R}^{[\ell]},\cdot) and similar Φ[1](𝒮):=Φ1([1],)\Phi^{[\ell-1]}(\mathcal{S})\!:=\!\Phi^{\geq\ell-1}_{\ell}(\mathcal{R}^{[\ell-1]},\cdot). We define U:=V×[+1]U:=V\times[\ell+1] and partition it into sets of cardinality k/(ν1)\lfloor k/(\nu-1)\rfloor plus a possible set of smaller size. The number of sets in the partition is t:=(+1)|V|k/(ν1)t:=\lceil\frac{(\ell+1)|V|}{\lfloor k/(\nu-1)\rfloor}\rceil. Denote these sets by U1,,UtU_{1},\ldots,U_{t}. Now, let 𝒯\mathcal{T} be any set of cardinality k/(ν1)\lfloor k/(\nu-1)\rfloor that maximizes Φ[]()\Phi^{[\ell]}(\cdot) and assume for the purpose of contradiction that

Φ[](𝒯)Φ[]()<1t(Φ[1](𝒮[1])Φ[]()).\displaystyle\Phi^{[\ell]}(\mathcal{T})-\Phi^{[\ell]}(\emptyset)<\frac{1}{t}\cdot(\Phi^{[\ell-1]}(\mathcal{S}^{[\ell-1]})-\Phi^{[\ell]}(\emptyset)). (7)

By definition of 𝒯\mathcal{T}, we have Φ[](𝒯)Φ[](Ui)\Phi^{[\ell]}(\mathcal{T})\geq\Phi^{[\ell]}(U_{i}) for i[t]i\in[t]. Hence, by submodularity we get

Φ[](U)Φ[]()i=1t(Φ[](Ui)Φ[]())t(Φ[](𝒯)Φ[]())<Φ[1](𝒮[1])Φ[]().\displaystyle\Phi^{[\ell]}(U)\!-\!\Phi^{[\ell]}(\emptyset)\!\leq\!\sum_{i=1}^{t}\!\big(\Phi^{[\ell]}(U_{i})\!-\!\Phi^{[\ell]}(\emptyset)\big)\!\leq\!t\big(\Phi^{[\ell]}(\mathcal{T})\!-\!\Phi^{[\ell]}(\emptyset)\big)\!<\!\Phi^{[\ell-1]}(\mathcal{S}^{[\ell-1]})\!-\!\Phi^{[\ell]}(\emptyset).

Since the maximum possible number of nodes, say NN are guaranteed to be reached by +1\ell+1 campaigns from sets UU, we have however that Φ[](V+1)=N\Phi^{[\ell]}(V^{\ell+1})=N. On the other hand we have Φ[1](𝒮[1])N\Phi^{[\ell-1]}(\mathcal{S}^{[\ell-1]})\leq N, which leads to a contradiction. From Lemma 3 we know that, with probability at least 1δ/ν1-\delta/\nu, it holds that Φ[](𝒮[])(11/eϵ)Φ[](𝒯)\Phi^{[\ell]}(\mathcal{S}^{[\ell]})\geq(1-1/e-\epsilon^{\prime})\cdot\Phi^{[\ell]}(\mathcal{T}) with ϵ=ϵ/2\epsilon^{\prime}=\epsilon/2. Thus, together with the converse of (7), we get

Φ[](𝒮[])\displaystyle\Phi^{[\ell]}(\mathcal{S}^{[\ell]}) 11eϵt(Φ[1](𝒮[1])Φ[]())+Φ[]()11eϵtΦ[1](𝒮[1]).\displaystyle\geq\frac{1-\frac{1}{e}-\epsilon^{\prime}}{t}\cdot(\Phi^{[\ell-1]}(\mathcal{S}^{[\ell-1]})-\Phi^{[\ell]}(\emptyset))+\Phi^{[\ell]}(\emptyset)\geq\frac{1-\frac{1}{e}-\epsilon^{\prime}}{t}\cdot\Phi^{[\ell-1]}(\mathcal{S}^{[\ell-1]}).

It remains to observe that k(ν1)/ϵk\geq(\nu-1)/\epsilon^{\prime} implies kν+1(1ϵ)kk-\nu+1\geq(1-\epsilon^{\prime})k and thus

t=(+1)|V|k/(ν1)(+1)(ν1)|V|kν+12(+1)(ν1)|V|(1ϵ)k\displaystyle t=\left\lceil\frac{(\ell+1)|V|}{\lfloor k/(\nu-1)\rfloor}\right\rceil\leq\left\lceil\frac{(\ell+1)(\nu-1)|V|}{k-\nu+1}\right\rceil\leq\frac{2(\ell+1)(\nu-1)|V|}{(1-\epsilon^{\prime})k}

where the last inequality follows since kν|V|k\leq\nu\cdot|V| and ν2\nu\geq 2 yield that the argument of the ceil-function is at least 1, and thus the error due to rounding is upper bounded by a factor of 2. The choice of ϵ=ϵ/2\epsilon^{\prime}=\epsilon/2 leads the result. ∎

Appendix D Deferred Proofs for Section 5

See 12

Proof.
  1. 1.

    Define 𝒯:={(v,0)|(v,i)𝒮k}\mathcal{T}:=\{(v,0)|(v,i)\in\mathcal{S}^{k}\} and observe that |𝒯|k|\mathcal{T}|\leq k. For a given outcome 𝒳\mathcal{X} a node that contributes to Φ1(𝒮k)\Phi^{\geq 1}(\mathcal{S}^{k}) is either reached by at least ν\nu campaigns in \mathcal{I} or has to be reached by a node in 𝒮k\mathcal{S}^{k}. In this case, for the same 𝒳\mathcal{X} this node will also contribute to Ψ(𝒯)\Psi(\mathcal{T}). Hence, we have Ψ(𝒯)Φ1(𝒮k)\Psi(\mathcal{T})\geq\Phi^{\geq 1}(\mathcal{S}^{k}). The optimality of 𝒯k\mathcal{T}^{k} concludes the proof.

  2. 2.

    First observe that kk/ν<kk/ν1+1ν+11ϵ\lceil\frac{k}{\lfloor k/\nu\rfloor}\rceil<\frac{k}{k/\nu-1}+1\leq\frac{\nu+1}{1-\epsilon} by the assumption on kk. Now, let 𝒯\mathcal{T} be a subset of 𝒯k\mathcal{T}^{k} of size k/ν\lfloor k/\nu\rfloor maximizing Ψ\Psi. By submodularity of Ψ\Psi, we have Ψ(𝒯k)kk/νΨ(𝒯)ν+11ϵΨ(𝒯)\Psi(\mathcal{T}^{k})\leq\lceil\frac{k}{\lfloor k/\nu\rfloor}\rceil\Psi(\mathcal{T})\leq\frac{\nu+1}{1-\epsilon}\Psi(\mathcal{T}). Using the optimality of 𝒯k/ν\mathcal{T}^{\lfloor k/\nu\rfloor} concludes the proof.

  3. 3.

    Since the cascade processes are completely correlated, given an outcome 𝒳\mathcal{X}, assume that a node contributes to Ψ(𝒯)\Psi(\mathcal{T}), then either it is reached by ν\nu campaigns from \mathcal{I} or it is reached by 𝒯\mathcal{T}. In the former case, the same node also contributes to Φ1(𝒮)\Phi^{\geq 1}(\mathcal{S}^{\prime}) as it is reached by ν\nu campaigns from \mathcal{I}. In the later case, it will be reached by all campaigns in [ν][\nu] by 𝒮\mathcal{S}^{\prime} and will therefore also contribute to Φ1(𝒮)\Phi^{\geq 1}(\mathcal{S}^{\prime}).∎