The Densest Subgraph Problem with a Convex/Concave Size Function

TL;DR

This paper generalizes the densest subgraph problem by introducing an $f$-density with convex/concave functions, providing efficient algorithms with approximation guarantees for different cases.

Contribution

It proposes new algorithms for the $f$-densest subgraph problem that handle size constraints via convex or concave functions, including approximation and exact solutions.

Findings

Developed a nearly-linear-time approximation algorithm for convex $f$-DS.

Designed an LP-based exact algorithm for concave $f$-DS.

Created a flow-based $O(|V|^3)$-time exact algorithm for unweighted graphs.

Abstract

In the densest subgraph problem, given an edge-weighted undirected graph $G=(V,E,w)$, we are asked to find $S\subseteq V$ that maximizes the density, i.e., $w(S)/|S|$, where $w(S)$ is the sum of weights of the edges in the subgraph induced by $S$. This problem has often been employed in a wide variety of graph mining applications. However, the problem has a drawback; it may happen that the obtained subset is too large or too small in comparison with the size desired in the application at hand. In this study, we address the size issue of the densest subgraph problem by generalizing the density of $S\subseteq V$. Specifically, we introduce the $f$-density of $S\subseteq V$, which is defined as $w(S)/f(|S|)$, where $f:\mathbb{Z}_{\geq 0}\rightarrow \mathbb{R}_{\geq 0}$ is a monotonically non-decreasing function. In the $f$-densest subgraph problem ($f$-DS), we aim to find $S\subseteq V$ that maximizes the $f$-density $w(S)/f(|S|)$. Although $f$-DS does not explicitly specify the size of the output subset of vertices, we can handle the above size issue using a convex/concave size function $f$ appropriately. For $f$-DS with convex function $f$, we propose a nearly-linear-time algorithm with a provable approximation guarantee. On the other hand, for $f$-DS with concave function $f$, we propose an LP-based exact algorithm, a flow-based $O(|V|^3)$-time exact algorithm for unweighted graphs, and a nearly-linear-time approximation algorithm.

Figures

Equations

$$\frac{f(|S^{\ast }|)}{|S^{\ast }|}=\frac{w(S^{\ast })/|S^{\ast }|}{w(S^{\ast })/f(|S^{\ast }|)}\le \frac{w(S_{\text{DS}}^{\ast })/|S_{\text{DS}}^{\ast }|}{w(S_{\text{DS}}^{\ast })/f(|S_{\text{DS}}^{\ast }|)}=\frac{f(|S_{\text{DS}}^{\ast }|)}{|S_{\text{DS}}^{\ast }|}.$$

$$f(x)=\{\begin{array}{ll}x& \text{if }x\le k,\\ \frac{w(V)}{w(e)/2}(x-k)+k& \text{otherwise}.\end{array}$$

$$\frac{w(S)}{f(|S|)}=\frac{w(S)}{\frac{w(V)}{w(e)/2}(|S|-k)+k}<\frac{w(S)}{\frac{w(V)}{w(e)/2}}\le \frac{w(e)}{2},$$

$$\min \left\{\frac{f(2)/2}{f(|S^{\ast }|)/|S^{\ast }{|}^2},\frac{2f(n)/n}{f(|S^{\ast }|)-f(|S^{\ast }|-1)}\right\}.$$

$$\frac{w(S_i^{\ast })}{\left(\genfrac{}{}{0ex}{}{i}{2}\right)}\ge \frac{w(S_{i+1}^{\ast }\setminus \{u\})}{\left(\genfrac{}{}{0ex}{}{i}{2}\right)}=\frac{w(S_{i+1}^{\ast })-d_{S_{i+1}^{\ast }}(u)}{\left(\genfrac{}{}{0ex}{}{i}{2}\right)}\ge \frac{(1-\frac{2}{i+1})\cdot w(S_{i+1}^{\ast })}{\left(\genfrac{}{}{0ex}{}{i}{2}\right)}=\frac{w(S_{i+1}^{\ast })}{\left(\genfrac{}{}{0ex}{}{i+1}{2}\right)},$$

$$\frac{w(S^{\ast })}{f(|S^{\ast }|)}\le \frac{2\cdot f(k)/k^2}{f(|S^{\ast }|)/|S^{\ast }{|}^2}\cdot \frac{w(S_k^{\ast })}{f(k)}.$$

$$\frac{w(S^{\ast })}{f(|S^{\ast }|)}$$

$$=\frac{1-1/|S^{\ast }|}{1-1/k}\cdot \frac{f(k)/k^2}{f(|S^{\ast }|)/|S^{\ast }{|}^2}\cdot \frac{w(S_k^{\ast })}{f(k)}\le \frac{2\cdot f(k)/k^2}{f(|S^{\ast }|)/|S^{\ast }{|}^2}\cdot \frac{w(S_k^{\ast })}{f(k)},$$

$$\frac{w(S^{\ast })}{f(|S^{\ast }|)}\le \frac{2f(n)/n}{f(|S^{\ast }|)-f(|S^{\ast }|-1)}\cdot \frac{w(S)}{f(|S|)}.$$

$$\frac{w(S^{\ast })}{f(|S^{\ast }|)}\ge \frac{w(S^{\ast }\setminus \{v\})}{f(|S^{\ast }|-1)}.$$

$$d_{S^{\ast }}(v)\ge (f(|S^{\ast }|)-f(|S^{\ast }|-1))\cdot \frac{w(S^{\ast })}{f(|S^{\ast }|)}.$$

$$\frac{w(S_l)}{f(l)}$$

$$\ge \frac{f(|S^{\ast }|)-f(|S^{\ast }|-1)}{2f(l)/l}\cdot \frac{w(S^{\ast })}{f(|S^{\ast }|)}\ge \frac{f(|S^{\ast }|)-f(|S^{\ast }|-1)}{2f(n)/n}\cdot \frac{w(S^{\ast })}{f(|S^{\ast }|)},$$

$$\min \left\{\frac{f(2)/2}{f(s)/s^2},\frac{2f(n)/n}{f(s)-f(s-1)}\right\}$$

$$\le \min \left\{2\cdot s^{2-\alpha },\frac{2n^{\alpha -1}}{s^{\alpha -1}}\right\}\le 2\cdot n^{(\alpha -1)(2-\alpha )}.$$

$$\frac{2\cdot f(k)/k^2}{f(s)/s^2}=2\cdot \frac{\lambda /k+(1-\lambda )}{\lambda /s+(1-\lambda )}\le 2\cdot \frac{\lambda /k+(1-\lambda )}{1-\lambda }=2+\frac{2\lambda }{(1-\lambda )k}.$$

$$\min \left\{\frac{f(2)/2}{f(s)/s^2},\frac{2f(n)/n}{f(s)-f(s-1)}\right\}$$

$$\le \min \left\{\frac{2(\lambda s+(1-\lambda ))}{1+\lambda },\frac{\frac{2n}{\lambda n+(1-\lambda )}}{\frac{s}{\lambda s+(1-\lambda )}}\right\}$$

$$\le \frac{2}{1+\lambda }\left(\lambda \cdot \frac{(1+\lambda )n}{\lambda n+(1-\lambda )}+(1-\lambda )\right)$$

$$\le \frac{2}{1+\lambda }\left(\lambda \cdot \frac{1+\lambda }{\lambda }+(1-\lambda )\right)=\frac{4}{1+\lambda },$$

$$\hat{\lambda }(|S|-k)+w(S_k^{\ast })\ge \frac{w(S_k^{\ast })}{f(k)}\cdot f(|S|)\ge \frac{w(S)}{f(|S|)}\cdot f(|S|)=w(S)$$

$${\mathrm{LP}}_k:$$

$$\sum _{v\in V}{y}_v=k,$$

$$x_e\le y_u,x_e\le y_v$$

$$x_e,y_v\in [0,1]$$

$$x_e=\{\begin{array}{ll}1& \text{if }e\in E(S),\\ 0& \text{otherwise},\end{array}\text{and}{y}_v=\{\begin{array}{ll}1& \text{if }v\in S,\\ 0& \text{otherwise}.\end{array}$$

$${\int }_0^{y_1^{\ast }}w(S^{k^{\ast }}(r))dr$$

$$=\sum _{e=\{u,v\}\in E}{\int }_0^{y_1^{\ast }}\left(w(e)\cdot [y_u^{\ast }\ge r\text{and}{y}_v^{\ast }\ge r]\right)dr$$

$$=\sum _{e=\{u,v\}\in E}w(e)\cdot \min \{y_u^{\ast },y_v^{\ast }\}=\sum _{e\in E}w(e)\cdot x_e^{\ast }\ge w(S^{\ast }),$$

$${\int }_0^{y_1^{\ast }}f(|S^{k^{\ast }}(r)|)dr$$

$$\le \sum _{h=1}^{k^{\ast }}(f(h)-f(h-1))=f(k^{\ast })-f(0)=f(k^{\ast }),$$

$$\frac{w(S^{\ast })}{f(k^{\ast })}$$

$$\le \frac{{\int }_0^{y_1^{\ast }}\left(\frac{w(S^{k^{\ast }}(r^{\ast }))}{f(|S^{k^{\ast }}(r^{\ast })|)}\cdot f(|S^{k^{\ast }}(r)|)\right)dr}{{\int }_0^{y_1^{\ast }}f(|S^{k^{\ast }}(r)|)dr}=\frac{w(S^{k^{\ast }}(r^{\ast }))}{f(|S^{k^{\ast }}(r^{\ast })|)}.$$

$$A_s$$

$$A_1$$

$$w_{\beta }(e)=\{\begin{array}{ll}d(v)/2& (e=(s,v)\in A_s),\\ \beta \cdot k\cdot a_k& (e=(p_k,t)\in A_t),\\ 1/2(=w(\{u,v\})/2)& (e=(u,v)\in A_1),\\ \beta \cdot a_k& (e=(v,p_k)\in A_2),\end{array}$$

$$a_k=\{\begin{array}{ll}2f(k)-f(k+1)-f(k-1)& (k=1,\dots ,n-1),\\ f(n)-f(n-1)& (k=n).\end{array}$$

Full text

1 Introduction

Finding dense components in a graph is an active research topic in graph mining. Techniques for identifying dense subgraphs have been used in various applications. For example, in Web graph analysis, they are used for detecting communities (i.e., sets of web pages dealing with the same or similar topics) [9] and spam link farms [13]. As another example, in bioinformatics, they are used for finding molecular complexes in protein–protein interaction networks [4] and identifying regulatory motifs in DNA [11]. Furthermore, they are also used for expert team formation [6, 20] and real-time story identification in micro-blogging streams [2].

To date, various optimization problems have been considered to find dense components in a graph. The densest subgraph problem is one of the most well-studied optimization problems. Let $G=(V,E,w)$ be an edge-weighted undirected graph consisting of $n=|V|$ vertices, $m=|E|$ edges, and a weight function $w:E\rightarrow\mathbb{Q}_{>0}$, where $\mathbb{Q}_{>0}$ is the set of positive rational numbers. For a subset of vertices $S\subseteq V$, let $G[S]$ be the subgraph induced by $S$, i.e., $G[S]=(S,E(S))$, where $E(S)=\{\{i,j\}\in E\mid i,j\in S\}$. The density of $S\subseteq V$ is defined as $w(S)/|S|$, where $w(S)=\sum_{e\in E(S)}w(e)$. In the (weighted) densest subgraph problem, given an (edge-weighted) undirected graph $G=(V,E,w)$, we are asked to find $S\subseteq V$ that maximizes the density $w(S)/|S|$.

The densest subgraph problem has received significant attention recently because it can be solved exactly in polynomial time and approximately in nearly linear time. In fact, there exist a flow-based exact algorithm [14] and a linear-programming-based (LP-based) exact algorithm [7]. Charikar [7] demonstrated that the greedy algorithm designed by Asahiro et al. [3], which is called the greedy peeling, obtains a $2$-approximate solution222 A feasible solution is said to be $\gamma$-approximate if its objective value times $\gamma$ is greater than or equal to the optimal value. An algorithm is called a $\gamma$-approximation algorithm if it runs in polynomial time and returns a $\gamma$-approximate solution for any instance. For a $\gamma$-approximation algorithm, $\gamma$ is referred to as an approximation ratio of the algorithm. for any instance. This algorithm runs in $O(m+n\log n)$ time for weighted graphs and $O(m+n)$ time for unweighted graphs.

However, the densest subgraph problem has a drawback; it may happen that the obtained subset is too large or too small in comparison with the size desired in the application at hand. To overcome this issue, some variants of the problem have often been employed. The densest $k$-subgraph problem (D$k$S) is a straightforward size-restricted variant of the densest subgraph problem [10]. In this problem, given an additional input $k$ being a positive integer, we are asked to find $S\subseteq V$ of size $k$ that maximizes the density $w(S)/|S|$. Note that in this problem, the objective function can be replaced by $w(S)$ since $|S|$ is fixed to $k$. Unfortunately, it is known that this size restriction makes the problem much harder to solve. In fact, Khot [16] proved that D$k$S has no PTAS under some reasonable computational complexity assumption. The current best approximation algorithm has an approximation ratio of $O(n^{1/4+\epsilon})$ for any $\epsilon>0$ [5].

Furthermore, Andersen and Chellapilla [1] introduced two relaxed versions of D$k$S. The first problem, the densest at-least-$k$-subgraph problem (Dal$k$S), asks for $S\subseteq V$ that maximizes the density $w(S)/|S|$ under the size constraint $|S|\geq k$. For this problem, Andersen and Chellapilla [1] adopted the greedy peeling, and demonstrated that the algorithm yields a $3$-approximate solution for any instance. Later, Khuller and Saha [17] investigated the problem more deeply. They proved that Dal$k$S is NP-hard, and designed a flow-based algorithm and an LP-based algorithm. These algorithms have an approximation ratio of $2$, which improves the above approximation ratio of $3$. The second problem is called the densest at-most-$k$-subgraph problem (Dam$k$S), which asks for $S\subseteq V$ that maximizes the density $w(S)/|S|$ under the size constraint $|S|\leq k$. The NP-hardness is immediate since finding a maximum clique can be reduced to it. Khuller and Saha [17] proved that approximating Dam$k$S is as hard as approximating D$k$S, within a constant factor.

1.1 Our Contribution

In this study, we address the size issue of the densest subgraph problem by generalizing the density of $S\subseteq V$. Specifically, we introduce the $f$-density of $S\subseteq V$, which is defined as $w(S)/f(|S|)$, where $f:\mathbb{Z}_{\geq 0}\rightarrow\mathbb{R}_{\geq 0}$ is a monotonically non-decreasing function with $f(0)=0$.333 To handle various types of functions (e.g., $f(x)=x^{\alpha}$ for $\alpha>0$), we set the codomain of the function $f$ to be the set of nonnegative real numbers. We assume that we can compare $p\cdot f(i)$ and $q\cdot f(j)$ in constant time for any $p,q\in\mathbb{Q}$ and $i,j\in\mathbb{Z}_{\geq 0}$.

Note that $\mathbb{Z}_{\geq 0}$ and $\mathbb{R}_{\geq 0}$ are the sets of nonnegative integers and nonnegative real numbers, respectively. In the $f$-densest subgraph problem ($f$-DS), we aim to find $S\subseteq V$ that maximizes the $f$-density $w(S)/f(|S|)$. For simplicity, we assume that $E\neq\emptyset$. Hence, any optimal solution $S^{*}\subseteq V$ satisfies $|S^{*}|\geq 2$. Although $f$-DS does not explicitly specify the size of the output subset of vertices, we can handle the above size issue using a convex size function $f$ or a concave size function $f$ appropriately. In fact, we can show that any optimal solution to $f$-DS with convex (resp. concave) function $f$ has a size smaller (resp. larger) than or equal to that of any densest subgraph (i.e., any optimal solution to the densest subgraph problem). For details, see Sections 2 and 3.

Here we mention the relationship between our problem and D$k$S. Any optimal solution $S^{*}\subseteq V$ to $f$-DS is a maximum weight subset of size $|S^{*}|$, i.e., $S^{*}\in\mathop{\rm argmax}\{w(S)\mid S\subseteq V,~{}|S|=|S^{*}|\}$, which implies that $S^{*}$ is also optimal to D$k$S with $k=|S^{*}|$. Furthermore, the iterative use of a $\gamma$-approximation algorithm for D$k$S leads to a $\gamma$-approximation algorithm for $f$-DS. Using the above $O(n^{1/4+\epsilon})$-approximation algorithm for D$k$S [5], we can obtain an $O(n^{1/4+\epsilon})$-approximation algorithm for $f$-DS.

In what follows, we summarize our results for both the cases where $f$ is convex and where $f$ is concave.

The case where $f$ is convex.

We first describe our results for the case where $f$ is convex. A function $f:\mathbb{Z}_{\geq 0}\rightarrow\mathbb{R}_{\geq 0}$ is said to be convex if $f(x)-2f(x+1)+f(x+2)\geq 0$ holds for any $x\in\mathbb{Z}_{\geq 0}$. We first prove the NP-hardness of $f$-DS with a certain convex function $f$ by constructing a reduction from Dam$k$S. Thus, for $f$-DS with convex function $f$, one of the best possible ways is to design an algorithm with a provable approximation guarantee.

To this end, we propose a $\min\left\{\frac{f(2)/2}{f(|S^{*}|)/|S^{*}|^{2}},~{}\frac{2f(n)/n}{f(|S^{*}|)-f(|S^{*}|-1)}\right\}$-approximation algorithm, where $S^{*}\subseteq V$ is an optimal solution to $f$-DS with convex function $f$. Our algorithm consists of the following two procedures, and outputs the better solution found by them. The first one is based on the brute-force search, which obtains an $\frac{f(2)/2}{f(|S^{*}|)/|S^{*}|^{2}}$-approximate solution in $O(m+n)$ time. The second one adopts the greedy peeling, which obtains a $\frac{2f(n)/n}{f(|S^{*}|)-f(|S^{*}|-1)}$-approximate solution in $O(m+n\log n)$ time. Thus, the total running time of our algorithm is $O(m+n\log n)$. Our analysis on the approximation ratio of the second procedure extends the analysis by Charikar [7] for the densest subgraph problem.

At the end of our analysis, we observe the behavior of the approximation ratio of our algorithm for three concrete size functions. We consider size functions between linear and quadratic because, as we will see later, $f$-DS with any super-quadratic size function is a trivial problem; in fact, it only produces constant-size optimal solutions. The first example is $f(x)=x^{\alpha}~{}(\alpha\in[1,2])$. We show that the approximation ratio of our algorithm is $2\cdot n^{(\alpha-1)(2-\alpha)}$, where the worst-case performance of $2\cdot n^{1/4}$ is attained at $\alpha=1.5$. The second example is $f(x)=\lambda x+(1-\lambda)x^{2}~{}(\lambda\in[0,1))$. For this case, the approximation ratio of our algorithm is $(2-\lambda)/(1-\lambda)$, which is a constant for a fixed $\lambda$. The third example is $f(x)=x^{2}/(\lambda x+(1-\lambda))~{}(\lambda\in[0,1])$. Note that this size function is derived by density function $\lambda\frac{w(S)}{|S|}+(1-\lambda)\frac{w(S)}{|S|^{2}}$. The approximation ratio of our algorithm is $4/(1+\lambda)$, which is at most $4$.

The case where $f$ is concave.

We next describe our results for the case where $f$ is concave. A function $f:\mathbb{Z}_{\geq 0}\rightarrow\mathbb{R}_{\geq 0}$ is said to be concave if $f(x)-2f(x+1)+f(x+2)\leq 0$ holds for any $x\in\mathbb{Z}_{\geq 0}$. Unlike the above convex case, $f$-DS in this case can be solved exactly in polynomial time.

In fact, we present an LP-based exact algorithm, which extends Charikar’s exact algorithm for the densest subgraph problem [7] and Khuller and Saha’s $2$-approximation algorithm for Dal$k$S [17]. It should be emphasized that our LP-based algorithm obtains not only an optimal solution to $f$-DS but also some attractive subsets of vertices. Let us see an example in Figure 1. The graph consists of 8 vertices and 11 unweighted edges (i.e., $w(e)=1$ for every $e\in E$). For this graph, we plotted all the points contained in $\mathcal{P}=\{(|S|,w(S))\mid S\subseteq V\}$. We refer to the extreme points of the upper convex hull of $\mathcal{P}$ as the dense frontier points. The (smallest) densest subgraph is a typical subset of vertices corresponding to a dense frontier point. Our LP-based algorithm obtains a corresponding subset of vertices for every dense frontier point. It should be noted that the algorithm SSM designed by Nagano, Kawahara, and Aihara [18] can also be used to obtain a corresponding subset of vertices for every dense frontier point. The difference between their algorithm and ours is that their algorithm is based on the computation of a minimum norm base, whereas ours solves linear programming problems.

Moreover, in this concave case, we design a combinatorial exact algorithm for unweighted graphs. Our algorithm is based on the standard technique for fractional programming. By using the technique, we can reduce $f$-DS to a sequence of submodular function minimizations. However, the direct application of a submodular function minimization algorithm leads to a computationally expensive algorithm that runs in $O(n^{5}(m+n)\cdot\log n)$ time. To reduce the computation time, we replace a submodular function minimization algorithm with a much faster flow-based algorithm that substantially extends a technique of Goldberg’s flow-based algorithm for the densest subgraph problem [14]. The total running time of our algorithm is $O(n^{3})$. Modifying this algorithm, we also present an $O\left(\frac{n^{3}}{\log n}\cdot\log\left(\frac{\log n}{\epsilon}\right)\right)$-time $(1+\epsilon)$-approximation algorithm for weighted graphs.

Although our flow-based algorithm is much faster than the reduction-based algorithm, the running time is still long for large-sized graphs. To design an algorithm with much higher scalability, we adopt the greedy peeling. As mentioned above, this algorithm runs in $O(m+n\log n)$ time for weighted graphs and $O(m+n)$ time for unweighted graphs. We prove that the algorithm yields a $3$-approximate solution for any instance.

1.2 Related Work

Tsourakakis et al. [20] introduced a general optimization problem to find dense subgraphs, which is referred to as the optimal $(g,h,\alpha)$-edge-surplus problem. In this problem, given an unweighted undirected graph $G=(V,E)$, we are asked to find $S\subseteq V$ that maximizes $\textsf{edge-surplus}_{\alpha}(S)=g(|E(S)|)-\alpha h(|S|)$, where $g$ and $h$ are strictly monotonically increasing functions, and $\alpha>0$ is a constant. The intuition behind this optimization problem is the same as that of $f$-DS. In fact, the first term $g(|E(S)|)$ prefers $S\subseteq V$ that has a large number of edges, whereas the second term $-\alpha h(|S|)$ penalizes $S\subseteq V$ with a large size. Tsourakakis et al. [20] were motivated by finding near-cliques (i.e., relatively small dense subgraphs), and they derived the function $\textsf{OQC}_{\alpha}(S)=|E(S)|-\alpha\binom{|S|}{2}$, which is called the OQC function, by setting $g(x)=x$ and $h(x)=x(x-1)/2$. For OQC function maximization, they adopted the greedy peeling and a simple local search heuristic.

Recently, Yanagisawa and Hara [21] introduced density function $|E(S)|/|S|^{\alpha}$ for $\alpha\in(1,2]$, which they called the discounted average degree. For discounted average degree maximization, they designed an integer-programming-based exact algorithm, which is applicable only to graphs with a maximum of a few thousand edges. They also designed a local search heuristic, which is applicable to web-scale graphs but has no provable approximation guarantee. As mentioned above, our algorithm for $f$-DS with convex function $f$ runs in $O(m+n\log n)$ time, and has an approximation ratio of $2\cdot n^{(\alpha-1)(2-\alpha)}$ for $f(x)=x^{\alpha}$ ($\alpha\in[1,2]$).

2 Convex Case

In this section, we investigate $f$-DS with convex function $f$. A function $f:\mathbb{Z}_{\geq 0}\rightarrow\mathbb{R}_{\geq 0}$ is said to be convex if $f(x)-2f(x+1)+f(x+2)\geq 0$ holds for any $x\in\mathbb{Z}_{\geq 0}$. We remark that $f(x)/x$ is monotonically non-decreasing for $x$ since we assume that $f(0)=0$. It should be emphasized that any optimal solution to $f$-DS with convex function $f$ has a size smaller than or equal to that of any densest subgraph. To see this, let $S^{*}\subseteq V$ be any optimal solution to $f$-DS and $S^{*}_{\text{DS}}\subseteq V$ be any densest subgraph. Then we have

[TABLE]

This implies that $|S^{*}|\leq|S^{*}_{\text{DS}}|$ holds because $f(x)/x$ is monotonically non-decreasing.

2.1 Hardness

We first prove that $f$-DS with convex function $f$ contains Dam$k$S as a special case.

Theorem 1**.**

For any integer $k\in[2,n]$, $S\subseteq V$ is optimal to Dam$k$S if and only if $S$ is optimal to $f$-DS with (convex) function $f(x)=\max\left\{x,~{}\frac{w(V)}{w(e)/2}(x-k)+k\right\}$, where $e$ is an arbitrary edge.

Proof.

Since the maximum of linear functions is convex, the function $f$ is convex. We remark that

[TABLE]

For any $S\subseteq V$ with $|S|\leq k$, we have $w(S)/f(|S|)=w(S)/|S|$. On the other hand, for any $S\subseteq V$ with $|S|>k$, we have

[TABLE]

which implies that $S$ is not optimal to $f$-DS. Thus, we have the theorem. ∎

2.2 Our Algorithm

In this subsection, we provide an algorithm for $f$-DS with convex function $f$. Our algorithm consists of the following two procedures, and outputs the better solution found by them. Let $S^{*}\subseteq V$ be an optimal solution to the problem. The first one is based on the brute-force search, which obtains an $\frac{f(2)/2}{f(|S^{*}|)/|S^{*}|^{2}}$-approximate solution in $O(m+n)$ time. The second one adopts the greedy peeling [3], which obtains a $\frac{2f(n)/n}{f(|S^{*}|)-f(|S^{*}|-1)}$-approximate solution in $O(m+n\log n)$ time. Combining these results, both of which will be proved later, we have the following theorem.

Theorem 2**.**

Let $S^{*}\subseteq V$ be an optimal solution to $f$-DS with convex function $f$. For the problem, our algorithm runs in $O(m+n\log n)$ time, and has an approximation ratio of

[TABLE]

2.2.1 Brute-Force Search

As will be shown below, to obtain an $\frac{f(2)/2}{f(|S^{*}|)/|S^{*}|^{2}}$-approximate solution, it suffices to find the heaviest edge (i.e., $\mathop{\rm argmax}\{w(e)\mid e\in E\}$), which can be done in $O(m+n)$ time. Here we present a more general algorithm, which is useful for some case. Our algorithm examines all the subsets of vertices of size at most $k$, and then returns an optimal subset among them, where $k$ is a constant that satisfies $k\geq 2$. For reference, we describe the procedure in Algorithm 1. This algorithm can be implemented to run in $O((m+n)n^{k})$ time because the number of subsets with at most $k$ vertices is $\sum_{i=0}^{k}\binom{n}{i}=O(n^{k})$ and the value of $w(S)/f(|S|)$ for each $S\subseteq V$ can be computed in $O(m+n)$ time.

We analyze the approximation ratio of the algorithm. Let $S_{i}^{*}\subseteq V$ denote a maximum weight subset of size $i\geq 2$, i.e., $S_{i}^{*}\in\mathop{\rm argmax}\{w(S)\mid S\subseteq V,~{}|S|=i\}$. We refer to $w(S_{i}^{*})/\binom{i}{2}$ as the edge density of $i$ vertices. The following lemma gives a fundamental property of the edge density.

Lemma 1**.**

The edge density is monotonically non-increasing for the number of vertices, i.e., $w(S_{i}^{*})/\binom{i}{2}\geq w(S_{j}^{*})/\binom{j}{2}$ holds for any $2\leq i\leq j\leq n$.

Proof.

It suffices to show that $w(S^{*}_{i})/\binom{i}{2}\geq w(S^{*}_{i+1})/\binom{i+1}{2}$ holds for any positive integer $i\in[2,n-1]$. For $S\subseteq V$ and $v\in S$, let $d_{S}(v)$ denote the weighted degree of $v$ in the induced subgraph $G[S]$, i.e., $d_{S}(v)=\sum_{u\in V:\,\{u,v\}\in E(S)}w(\{u,v\})$. Take a vertex $u\in\mathop{\rm argmin}\{d_{S^{*}_{i+1}}(v)\mid v\in S^{*}_{i+1}\}$. Then we obtain $d_{S^{*}_{i+1}}(u)\leq\frac{1}{i+1}\sum_{v\in S^{*}_{i+1}}d_{S^{*}_{i+1}}(v)=\frac{2}{i+1}\cdot w(S^{*}_{i+1})$. Hence, we have

[TABLE]

as desired. ∎

Using the above lemma, we can provide the approximation ratio.

Lemma 2**.**

Let $S^{*}\subseteq V$ be an optimal solution to $f$-DS with convex function $f$. If $|S^{*}|\leq k$, then Algorithm 1 obtains an optimal solution. If $|S^{*}|\geq k$, then it holds that

[TABLE]

Proof.

If $|S^{*}|\leq k$, then Algorithm 1 obtains an optimal solution because $S^{*}\in\{S^{*}_{2},\dots,S^{*}_{k}\}$. If $|S^{*}|\geq k$, then we have

[TABLE]

where the first inequality follows from Lemma 1, and the last inequality follows from $k\geq 2$. ∎

From this lemma, we see that Algorithm 1 with $k=2$ has an approximation ratio of $\frac{f(2)/2}{f(|S^{*}|)/|S^{*}|^{2}}$.

2.2.2 Greedy Peeling

Here we adopt the greedy peeling. For $S\subseteq V$ and $v\in S$, let $d_{S}(v)$ denote the weighted degree of $v$ in the induced subgraph $G[S]$, i.e., $d_{S}(v)=\sum_{u\in V:\,\{u,v\}\in E(S)}w(\{u,v\})$. Our algorithm iteratively removes the vertex with the smallest weighted degree in the currently remaining graph, and then returns $S\subseteq V$ with maximum $w(S)/f(|S|)$ over the iterations. For reference, we describe the procedure in Algorithm 2. This algorithm can be implemented to run in $O(m+n\log n)$ time for weighted graphs and $O(m+n)$ time for unweighted graphs.

The following lemma provides the approximation ratio.

Lemma 3**.**

Let $S^{*}\subseteq V$ be an optimal solution to $f$-DS with convex function $f$. Algorithm 2 returns $S\subseteq V$ that satisfies

[TABLE]

Proof.

Choose an arbitrary vertex $v\in S^{*}$. By the optimality of $S^{*}$, we have

[TABLE]

By using the fact that $w(S^{*}\setminus\{v\})=w(S^{*})-d_{S^{*}}(v)$, the above inequality can be transformed to

[TABLE]

Let $l$ be the smallest index that satisfies $S_{l}\supseteq S^{*}$, where $S_{l}$ is the subset of vertices of size $l$ appeared in Algorithm 2. Note that $v_{l}\ (\in\mathop{\rm argmin}_{v\in S_{l}}d_{S_{l}}(v))$ is contained in $S^{*}$. Then we have

[TABLE]

where the first inequality follows from the greedy choice of $v_{l}$, the second inequality follows from $S_{l}\supseteq S^{*}$, the third inequality follows from inequality (2), and the last inequality follows from the monotonicity of $f(x)/x$. Since Algorithm 2 considers $S_{l}$ as a candidate subset of the output, we have the lemma. ∎

2.3 Examples

Here we observe the behavior of the approximation ratio of our algorithm for three concrete convex size functions. We consider size functions between linear and quadratic because $f$-DS with any super-quadratic size function is a trivial problem; in fact, it only produces constant-size optimal solutions. This follows from the inequality $\frac{f(2)/2}{f(|S^{*}|)/|S^{*}|^{2}}\geq 1$ (i.e., $f(2)/2\geq f(|S^{*}|)/|S^{*}|^{2}$) by Lemma 2.

(i) $\bm{f(x)=x^{\alpha}\ (\alpha\in[1,2])}$.

The following corollary provides an approximation ratio of our algorithm.

Corollary 1**.**

For $f$-DS with $f(x)=x^{\alpha}~{}(\alpha\in[1,2])$, our algorithm has an approximation ratio of $2\cdot n^{(\alpha-1)(2-\alpha)}$.

Proof.

Let $s=|S^{*}|$. By Theorem 2, the approximation ratio is

[TABLE]

The first inequality follows from the fact that $s^{\alpha}-(s-1)^{\alpha}=s^{\alpha}-(s-1)^{\alpha-1}(s-1)\geq s^{\alpha}-s^{\alpha-1}(s-1)=s^{\alpha-1}$. The last inequality follows from the fact that the first term and the second term of the minimum function are, respectively, monotonically non-decreasing and non-increasing for $s$, and they have the same value at $s=n^{\alpha-1}$. ∎

Note that an upper bound on $2\cdot n^{(\alpha-1)(2-\alpha)}$ is $2\cdot n^{1/4}$, which is attained at $\alpha=1.5$.

(ii) $\bm{f(x)=\lambda x+(1-\lambda)x^{2}\ (\lambda\in[0,1))}$.

The following corollary provides an approximation ratio of Algorithm 1, which is a constant for a fixed $\lambda$.

Corollary 2**.**

For $f$-DS with $f(x)=\lambda x+(1-\lambda)x^{2}~{}(\lambda\in[0,1))$, Algorithm 1 with $k=2$ has an approximation ratio of $(2-\lambda)/(1-\lambda)$. Furthermore, for any $\epsilon>0$, Algorithm 1 with $k\geq\frac{2}{\epsilon}\cdot\frac{\lambda}{1-\lambda}$ has an approximation ratio of $2+\epsilon$.

Proof.

Let $s=|S^{*}|$. By Lemma 2, the approximation ratio is

[TABLE]

Thus, by choosing $k=2$, the approximation ratio is at most $(2-\lambda)/(1-\lambda)$. For any $\epsilon>0$, by choosing $k\geq\frac{2}{\epsilon}\cdot\frac{\lambda}{1-\lambda}$, the approximation ratio is at most $2+\epsilon$. ∎

(iii) $\bm{f(x)=x^{2}/(\lambda x+(1-\lambda))\ (\lambda\in[0,1])}$.

This size function is derived by density function $\lambda\frac{w(S)}{|S|}+(1-\lambda)\frac{w(S)}{|S|^{2}}$. The following corollary provides an approximation ratio of our algorithm, which is at most $4$.

Corollary 3**.**

For $f$-DS with $f(x)=x^{2}/(\lambda x+(1-\lambda))~{}(\lambda\in[0,1))$, our algorithm has an approximation ratio of $4/(1+\lambda)$.

Proof.

Let $s=|S^{*}|$. By Theorem 2, the approximation ratio is

[TABLE]

where the second inequality follows from the fact that the first term and the second term of the minimum function are, respectively, monotonically non-decreasing and non-increasing for $s$, and they have the same value at $s=\frac{(1+\lambda)n}{\lambda n+(1-\lambda)}$. ∎

3 Concave Case

In this section, we investigate $f$-DS with concave function $f$. A function $f:\mathbb{Z}_{\geq 0}\rightarrow\mathbb{R}_{\geq 0}$ is said to be concave if $f(x)-2f(x+1)+f(x+2)\leq 0$ holds for any $x\in\mathbb{Z}_{\geq 0}$. We remark that $f(x)/x$ is monotonically non-increasing for $x$ since we assume that $f(0)=0$. It should be emphasized that any optimal solution to $f$-DS with concave function $f$ has a size larger than or equal to that of any densest subgraph. This follows from inequality (1) and the monotonicity of $f(x)/x$.

3.1 Dense Frontier Points

Here we define the dense frontier points and prove some basic properties. We denote by $\mathcal{P}$ the set $\{(|S|,w(S))\mid S\subseteq V\}$. A point in $\mathcal{P}$ is called a dense frontier point if it is a unique maximizer of $y-\lambda x$ over $(x,y)\in\mathcal{P}$ for some $\lambda>0$. In other words, the extreme points of the upper convex hull of $\mathcal{P}$ are dense frontier points. The (smallest) densest subgraph is a typical subset of vertices corresponding to a dense frontier point. We prove that (i) for any dense frontier point, there exists some concave function $f$ such that any optimal solution to $f$-DS with the function $f$ corresponds to the dense frontier point, and conversely, (ii) for any strictly concave function $f$ (i.e., $f$ that satisfies $f(x)-2f(x+1)+f(x+2)<0$ for any $x\in\mathbb{Z}_{\geq 0}$), any optimal solution to $f$-DS with the function $f$ corresponds to a dense frontier point.

We first prove (i). Note that each dense frontier point can be written as $(i,w(S_{i}^{*}))$ for some $i\in\{0,1,\dots,n\}$, where $S_{i}^{*}\subseteq V$ is a maximum weight subset of size $i$. Let $(k,w(S_{k}^{*}))$ be a dense frontier point and assume that it is a unique maximizer of $y-\hat{\lambda}x$ over $(x,y)\in\mathcal{P}$ for $\hat{\lambda}>0$. Consider the concave function $f$ such that $f(x)=\hat{\lambda}(x-k)+w(S_{k}^{*})$ for $x>0$ and $f(0)=0$ (see Figure 2). The concavity of $f$ follows from $w(S_{k}^{*})-\hat{\lambda}k\geq w(S_{0}^{*})-\hat{\lambda}\cdot 0=0=f(0)$. Then, any optimal solution $S^{*}\subseteq V$ to $f$-DS with the function $f$ corresponds to the dense frontier point (i.e., $(|S^{*}|,w(S^{*}))=(k,w(S_{k}^{*}))$ holds) because $w(S)/f(|S|)$ is greater than or equal to $1$ if and only if $w(S)-\hat{\lambda}|S|\geq w(S_{k}^{*})-\hat{\lambda}k$ holds.

We next prove (ii). Let $f$ be any strictly concave function. Let $S_{k}^{*}\subseteq V$ be any optimal solution to $f$-DS with the function $f$, and take $\hat{\lambda}$ that satisfies $(f(k)-f(k-1))\cdot\frac{w(S_{k}^{*})}{f(k)}>\hat{\lambda}>(f(k+1)-f(k))\cdot\frac{w(S_{k}^{*})}{f(k)}$ (see Figure 2). Note that the strict concavity of $f$ guarantees the existence of such $\hat{\lambda}$. Since $f$ is strictly concave, we have

[TABLE]

for any $S\subseteq V$, and the inequalities hold as equalities only when $(|S|,w(S))=(k,w(S_{k}^{*}))$. Thus, $(k,w(S_{k}^{*}))$ is a unique maximizer of $y-\hat{\lambda}x$ over $(x,y)\in\mathcal{P}$, and hence is a dense frontier point.

3.2 LP-Based Algorithm

We provide an LP-based polynomial-time exact algorithm. We introduce a variable $x_{e}$ for each $e\in E$ and a variable $y_{v}$ for each $v\in V$. For $k=1,\dots,n$, we construct the following linear programming problem:

[TABLE]

For an optimal solution $(\bm{x}^{k},\bm{y}^{k})$ to $\mathrm{LP}_{k}$ and a real parameter $r$, we define a sequence of subsets $S^{k}(r)=\{v\in V\mid y^{k}_{v}\geq r\}$. For $k=1,\dots,n$, our algorithm first solves $\mathrm{LP}_{k}$ to obtain an optimal solution $(\bm{x}^{k},\bm{y}^{k})$, and then computes $r^{*}_{k}$ that maximizes $w(S^{k}(r))/f(|S^{k}(r)|)$. Note here that to find such $r^{*}_{k}$, it suffices to check all the distinct sets $S^{k}(r)$ by simply setting $r=y^{k}_{v}$ for every $v\in V$. The algorithm returns $S\in\{S^{1}(r^{*}_{1}),\dots,S^{n}(r^{*}_{n})\}$ that maximizes $w(S)/f(|S|)$. For reference, we describe the procedure in Algorithm 3. Clearly, the algorithm runs in polynomial time.

In what follows, we demonstrate that Algorithm 3 obtains an optimal solution to $f$-DS with concave function $f$. The following lemma provides a lower bound on the optimal value of $\mathrm{LP}_{k}$.

Lemma 4**.**

For any $S\subseteq V$, the optimal value of $\mathrm{LP}_{|S|}$ is at least $w(S)$.

Proof.

For $S\subseteq V$, we construct a solution $(\bm{x},\bm{y})$ of $\mathrm{LP}_{|S|}$ as follows:

[TABLE]

Then we can easily check that $(\bm{x},\bm{y})$ is feasible for $\mathrm{LP}_{|S|}$ and its objective value is $w(S)$. Thus, we have the lemma. ∎

We prove the following key lemma.

Lemma 5**.**

Let $S^{*}\subseteq V$ be an optimal solution to $f$-DS with concave function $f$, and let $k^{*}=|S^{*}|$. Furthermore, let $(\bm{x}^{*},\bm{y}^{*})$ be an optimal solution to $\mathrm{LP}_{k^{*}}$. Then, there exists a real number $r$ such that $S^{k^{*}}(r)$ is optimal to $f$-DS with concave function $f$.

Proof.

For each $e=\{u,v\}\in E$, we have $x^{*}_{e}=\min\{y^{*}_{u},y^{*}_{v}\}$ from the optimality of $(x^{*},y^{*})$. Without loss of generality, we relabel the indices of $(\bm{x}^{*},\bm{y}^{*})$ so that $y^{*}_{1}\geq\dots\geq y^{*}_{n}$. Then we have

[TABLE]

where $[y^{*}_{u}\geq r~{}\text{and}~{}y^{*}_{v}\geq r]$ is the function of $r$ that takes 1 if the condition in the square bracket is satisfied and 0 otherwise, and the last inequality follows from Lemma 4. Moreover, we have

[TABLE]

where $y_{n+1}^{*}$ is defined to be [math] for convenience, and the inequality holds by the concavity of $f$ (i.e., $f(h+2)-f(h+1)\leq f(h+1)-f(h)$), $\sum_{h=1}^{n}y^{*}_{h}=k^{*}$, and $y^{*}_{h}\leq 1$.

Let $r^{*}$ be a real number that maximizes $w(S^{k^{*}}(r))/f(|S^{k^{*}}(r)|)$ in $[0,y_{1}^{*}]$. Using inequalities (3.2) and (3.2), we have

[TABLE]

This completes the proof. ∎

Algorithm 3 considers $S^{k^{*}}(r^{*})$ as a candidate subset of the output. Therefore, we have the desired result.

Theorem 3**.**

Algorithm 3 is a polynomial-time exact algorithm for $f$-DS with concave function $f$.

By Lemma 5, for any concave function $f$, an optimal solution to $f$-DS with the function $f$ is contained in $\{S^{k}(r)\mid k=1,\dots,n,~{}r\in[0,1]\}$ whose cardinality is at most $n^{2}$. As shown above, for any dense frontier point, there exists some concave function $f$ such that any optimal solution to $f$-DS with the function $f$ corresponds to the dense frontier point. Thus, we have the following result.

Theorem 4**.**

We can find a corresponding subset of vertices for every dense frontier point in polynomial time.

3.3 Flow-Based Algorithm

We provide a combinatorial exact algorithm for unweighted graphs (i.e., $w(e)=1$ for every $e\in E$). We first show that using the standard technique for fractional programming, we can reduce $f$-DS with concave function $f$ to a sequence of submodular function minimizations. The critical fact is that $\max_{S\subseteq V}w(S)/f(|S|)$ is at least $\beta$ if and only if $\min_{S\subseteq V}(\beta\cdot f(|S|)-w(S))$ is at most [math]. Note that for $\beta\geq 0$, the function $\beta\cdot f(|S|)-w(S)$ is submodular because $\beta\cdot f(|S|)$ and $-w(S)$ are submodular [12]. Thus, we can calculate $\min_{S\subseteq V}(\beta\cdot f(|S|)-w(S))$ in $O(n^{5}(m+n))$ time using Orlin’s algorithm [19], which implies that we can determine $\max_{S\subseteq V}w(S)/f(|S|)\geq\beta$ or not in $O(n^{5}(m+n))$ time. Hence, we can obtain the value of $\max_{S\subseteq V}w(S)/f(|S|)$ by binary search. Note that the objective function of $f$-DS on unweighted graphs may have at most $O(mn)$ distinct values since $w(S)$ is a nonnegative integer at most $m$. Thus, the procedure yields an optimal solution in $O(\log(mn))=O(\log n)$ iterations. The total running time is $O(n^{5}(m+n)\cdot\log n)$.

To reduce the computation time, we replace Orlin’s algorithm with a much faster flow-based algorithm that substantially extends a technique of Goldberg’s flow-based algorithm for the densest subgraph problem [14]. The key technique is to represent the value of $\min_{S\subseteq V}(\beta\cdot f(|S|)-w(S))$ using the cost of minimum cut of a certain directed network constructed from $G$ and $\beta\geq 0$.

For a given unweighted undirected graph $G=(V,E,w)$ (i.e., $w(e)=1$ for every $e\in E$) and a real number $\beta\geq 0$, we construct a directed network $(U,A,w_{\beta})$ as follows. Note that for later convenience, we discuss the procedure on weighted graphs. The vertex set $U$ is defined by $U=V\cup P\cup\{s,t\}$, where $P=\{p_{1},\dots,p_{n}\}$. The edge set $A$ is given by $A=A_{s}\cup A_{t}\cup A_{1}\cup A_{2}$, where

[TABLE]

The edge weight $w_{\beta}:A\to\mathbb{R}_{\geq 0}$ is defined by

[TABLE]

where $d(v)$ is the (weighted) degree of vertex $v$, and

[TABLE]

Note that $a_{k}\geq 0$ holds since $f$ is a monotonically non-decreasing concave function. For reference, Figure 3 depicts the network $(U,A,w_{\beta})$.

The following lemma reveals the relationship between a minimum $s$–$t$ cut in $(U,A,w_{\beta})$ and the value of $\min_{S\subseteq V}(\beta\cdot f(|S|)-w(S))$. Note that an $s$–$t$ cut in $(U,A,w_{\beta})$ is a partition $(X,Y)$ of $U$ (i.e., $X\cup Y=U$ and $X\cap Y=\emptyset$) such that $s\in X$ and $t\in Y$, and the cost of $(X,Y)$ is defined to be $\sum_{(u,v)\in A:\,u\in X,v\in Y}w_{\beta}(u,v)$.

Lemma 6**.**

Let $(X,Y)$ be any minimum $s$–$t$ cut in the network $(U,A,w_{\beta})$, and let $S=X\cap V$. Then, the cost of $(X,Y)$ is equal to $w(V)+\beta\cdot f(|S|)-w(S)$.

Proof.

We first show that for any positive integer $s~{}(\leq n)$, it holds that

[TABLE]

By the definition of $a_{k}$, we get

[TABLE]

Thus, we have

[TABLE]

We are now ready to prove the lemma. Note that $p_{k}\in X$ if $|S|>k$ and $p_{k}\in Y$ if $|S|<k$. Therefore, the cost of the minimum cut $(X,Y)$ is

[TABLE]

where the first equality follows from equality (5). ∎

From this lemma, we see that the cost of a minimum $s$–$t$ cut is $w(V)+\min_{S\subseteq V}(\beta\cdot f(|S|)-w(S))$. Therefore, for a given value $\beta\geq 0$, we can determine whether there exists $S\subseteq V$ that satisfies $w(S)/f(|S|)\geq\beta$ by checking the cost of a minimum $s$–$t$ cut is at most $w(V)$ or not. Our algorithm applies binary search for $\beta$ within the possible objective values of $f$-DS (i.e., $\{p/f(q)\mid p=0,1,\dots,m,~{}q=2,3,\dots,n\}$). For reference, we describe the procedure in Algorithm 4. The minimum $s$–$t$ cut problem can be solved in $O(N^{3}/\log N)$ time for a network with $N$ vertices [8]. Thus, the running time of our algorithm is $O(\frac{n^{3}}{\log n}\cdot\log(mn))=O(n^{3})$ since $|U|=2n+2$. We summarize the result in the following theorem.

Theorem 5**.**

Algorithm 4 is an $O(n^{3})$-time exact algorithm for $f$-DS with concave function $f$ on unweighted graphs.

For $f$-DS with concave function $f$ on weighted graphs, the binary search used in Algorithm 4 is not applicable because there may be exponentially many possible objective values in the weighted setting. Alternatively, we present an algorithm that employs another binary search strategy (Algorithm 5). We have the following theorem.

Theorem 6**.**

Algorithm 5 is an $O\left(\frac{n^{3}}{\log n}\cdot\log\left(\frac{\log n}{\epsilon}\right)\right)$-time $(1+\epsilon)$-approximation algorithm for $f$-DS with concave function $f$.

Proof.

Let $i^{*}$ be the number of iterations executed by Algorithm 5, and let $\hat{S}=X^{(i^{*})}\cap V$. Then we have $\max_{S\subseteq V}w(S)/f(|S|)\leq\beta_{\mathrm{ub}}^{(i^{*})}$, $\beta_{\mathrm{ub}}^{(i^{*})}<(1+\epsilon)\cdot\beta_{\mathrm{lb}}^{(i^{*})}$, and $\beta_{\mathrm{lb}}^{(i^{*})}\leq w(\hat{S})/f(|\hat{S}|)$. Combining these inequalities, we have $\max_{S\subseteq V}w(S)/f(|S|)<(1+\epsilon)\cdot w(\hat{S})/f(|\hat{S}|)$, which means that $\hat{S}$ is a $(1+\epsilon)$-approximate solution.

In what follows, we analyze the time complexity of the algorithm. For each $i\in\{0,1,\dots,i^{*}-1\}$, it holds that

[TABLE]

Since $\beta_{\mathrm{ub}}^{(0)}/\beta_{\mathrm{lb}}^{(0)}=w(V)/w(S_{2}^{*})\leq m$, we have $\beta_{\mathrm{ub}}^{(i)}/\beta_{\mathrm{lb}}^{(i)}\leq m^{1/2^{i}}$ for $i=1,\dots,i^{*}$. Note that $i^{*}$ is the minimum index $i$ that satisfies $\beta_{\mathrm{ub}}^{(i)}/\beta_{\mathrm{lb}}^{(i)}\leq 1+\epsilon$. Thus, we see that $i^{*}$ is upper bounded by $O\left(\log\left(\frac{\log n}{\log(1+\epsilon)}\right)\right)$. Therefore, the total running time of the algorithm is $O\left(\frac{n^{3}}{\log n}\cdot\log\left(\frac{\log n}{\log(1+\epsilon)}\right)\right)=O\left(\frac{n^{3}}{\log n}\cdot\log\left(\frac{\log n}{\epsilon}\right)\right)$, where the equality follows from the fact that $\lim_{\epsilon\to+0}\frac{\epsilon}{\log(1+\epsilon)}=1$ holds. ∎

3.4 Greedy Peeling

Finally, we provide an approximation algorithm with much higher scalability. Specifically, we prove that the greedy peeling (Algorithm 2) has an approximation ratio of 3 for $f$-DS with concave function $f$. As mentioned above, the algorithm runs in $O(m+n\log n)$ time for weighted graphs and $O(m+n)$ time for unweighted graphs.

We prove the approximation ratio. Recall that $S_{n},\dots,S_{1}$ are the subsets of vertices produced by the greedy peeling. We use the following fact, which implies that there exists a $3$-approximate solution for Dal$k$S in $S_{n},\dots,S_{k}$.

Fact 1** (Theorem 1 in Andersen and Chellapilla [1]).**

For any integer $k~{}(\leq n)$, it holds that

[TABLE]

Theorem 7**.**

The greedy peeling (Algorithm 2) has an approximation ratio of $3$ for $f$-DS with concave function $f$.

Proof.

Let $S^{*}\subseteq V$ be an optimal solution to $f$-DS with concave function $f$, and let $s^{*}=|S^{*}|$. Let $S\subseteq V$ be the output of the greedy peeling for the problem. Then we have

[TABLE]

where the second inequality follows from Fact 1, and the third inequality follows from the monotonicity of $f(x)/x$. ∎

Acknowledgments

The authors would like to thank Yoshio Okamoto for pointing out the reference [18]. The first author is supported by a Grant-in-Aid for Young Scientists (B) (No. 16K16005). The second author is supported by a Grant-in-Aid for JSPS Fellows (No. 26-11908).