Submodular set function

In mathematics, a submodular set function (also known as a submodular function) is a set function whose value, informally, has the property that the difference in the incremental value of the function that a single element makes when added to an input set decreases as the size of the input set increases. Submodular functions have a natural diminishing returns property which makes them suitable for many applications, including approximation algorithms, game theory (as functions modeling user preferences) and electrical networks. Recently, submodular functions have also found immense utility in several real world problems in machine learning and artificial intelligence, including automatic summarization, multi-document summarization, feature selection, active learning, sensor placement, image collection summarization and many other domains.^[1][2][3][4]

Definition[]

If ${\displaystyle \Omega }$ is a finite set, a submodular function is a set function ${\displaystyle f:2^{\Omega }\rightarrow \mathbb {R} }$, where ${\displaystyle 2^{\Omega }}$ denotes the power set of ${\displaystyle \Omega }$, which satisfies one of the following equivalent conditions.^[5]

1. For every ${\displaystyle X,Y\subseteq \Omega }$ with ${\displaystyle X\subseteq Y}$ and every ${\displaystyle x\in \Omega \setminus Y}$ we have that ${\displaystyle f(X\cup \{x\})-f(X)\geq f(Y\cup \{x\})-f(Y)}$.
2. For every ${\displaystyle S,T\subseteq \Omega }$ we have that ${\displaystyle f(S)+f(T)\geq f(S\cup T)+f(S\cap T)}$.
3. For every ${\displaystyle X\subseteq \Omega }$ and ${\displaystyle x_{1},x_{2}\in \Omega \backslash X}$ such that ${\displaystyle x_{1}\neq x_{2}}$ we have that ${\displaystyle f(X\cup \{x_{1}\})+f(X\cup \{x_{2}\})\geq f(X\cup \{x_{1},x_{2}\})+f(X)}$.

A nonnegative submodular function is also a subadditive function, but a subadditive function need not be submodular. If ${\displaystyle \Omega }$ is not assumed finite, then the above conditions are not equivalent. In particular a function ${\displaystyle f}$ defined by ${\displaystyle f(S)=1}$ if ${\displaystyle S}$ is finite and ${\displaystyle f(S)=0}$ if ${\displaystyle S}$ is infinite satisfies the first condition above, but the second condition fails when ${\displaystyle S}$ and ${\displaystyle T}$ are infinite sets with finite intersection.

Types of submodular functions[]

Monotone[]

A submodular function ${\displaystyle f}$ is monotone if for every ${\displaystyle T\subseteq S}$ we have that ${\displaystyle f(T)\leq f(S)}$. Examples of monotone submodular functions include:

Linear (Modular) functions

Any function of the form ${\displaystyle f(S)=\sum _{i\in S}w_{i}}$ is called a linear function. Additionally if ${\displaystyle \forall i,w_{i}\geq 0}$ then f is monotone.

Budget-additive functions

Any function of the form ${\displaystyle f(S)=\min \left\{B,~\sum _{i\in S}w_{i}\right\}}$ for each ${\displaystyle w_{i}\geq 0}$ and ${\displaystyle B\geq 0}$ is called budget additive.^{[citation needed]}

Coverage functions

Let ${\displaystyle \Omega =\{E_{1},E_{2},\ldots ,E_{n}\}}$ be a collection of subsets of some ground set ${\displaystyle \Omega '}$. The function ${\displaystyle f(S)=\left|\bigcup _{E_{i}\in S}E_{i}\right|}$ for ${\displaystyle S\subseteq \Omega }$ is called a coverage function. This can be generalized by adding non-negative weights to the elements.

Entropy

Let ${\displaystyle \Omega =\{X_{1},X_{2},\ldots ,X_{n}\}}$ be a set of random variables. Then for any ${\displaystyle S\subseteq \Omega }$ we have that ${\displaystyle H(S)}$ is a submodular function, where ${\displaystyle H(S)}$ is the entropy of the set of random variables ${\displaystyle S}$, a fact known as Shannon's inequality.^[6] Further inequalities for the entropy function are known to hold, see entropic vector.

Matroid rank functions

Let ${\displaystyle \Omega =\{e_{1},e_{2},\dots ,e_{n}\}}$ be the ground set on which a matroid is defined. Then the rank function of the matroid is a submodular function.^[7]

Non-monotone[]

A submodular function which is not monotone is called non-monotone.

Symmetric[]

A non-monotone submodular function ${\displaystyle f}$ is called symmetric if for every ${\displaystyle S\subseteq \Omega }$ we have that ${\displaystyle f(S)=f(\Omega -S)}$. Examples of symmetric non-monotone submodular functions include:

Graph cuts

Let ${\displaystyle \Omega =\{v_{1},v_{2},\dots ,v_{n}\}}$ be the vertices of a graph. For any set of vertices ${\displaystyle S\subseteq \Omega }$ let ${\displaystyle f(S)}$ denote the number of edges ${\displaystyle e=(u,v)}$ such that ${\displaystyle u\in S}$ and ${\displaystyle v\in \Omega -S}$. This can be generalized by adding non-negative weights to the edges.

Mutual information

Let ${\displaystyle \Omega =\{X_{1},X_{2},\ldots ,X_{n}\}}$ be a set of random variables. Then for any ${\displaystyle S\subseteq \Omega }$ we have that ${\displaystyle f(S)=I(S;\Omega -S)}$ is a submodular function, where ${\displaystyle I(S;\Omega -S)}$ is the mutual information.

Asymmetric[]

A non-monotone submodular function which is not symmetric is called asymmetric.

Directed cuts

Let ${\displaystyle \Omega =\{v_{1},v_{2},\dots ,v_{n}\}}$ be the vertices of a directed graph. For any set of vertices ${\displaystyle S\subseteq \Omega }$ let ${\displaystyle f(S)}$ denote the number of edges ${\displaystyle e=(u,v)}$ such that ${\displaystyle u\in S}$ and ${\displaystyle v\in \Omega -S}$. This can be generalized by adding non-negative weights to the directed edges.

Continuous extensions[]

Lovász extension[]

This extension is named after mathematician László Lovász. Consider any vector ${\displaystyle \mathbf {x} =\{x_{1},x_{2},\dots ,x_{n}\}}$ such that each ${\displaystyle 0\leq x_{i}\leq 1}$. Then the Lovász extension is defined as ${\displaystyle f^{L}(\mathbf {x} )=\mathbb {E} (f(\{i|x_{i}\geq \lambda \}))}$ where the expectation is over ${\displaystyle \lambda }$ chosen from the uniform distribution on the interval ${\displaystyle [0,1]}$. The Lovász extension is a convex function if and only if ${\displaystyle f}$ is a submodular function.

Multilinear extension[]

Consider any vector ${\displaystyle \mathbf {x} =\{x_{1},x_{2},\ldots ,x_{n}\}}$ such that each ${\displaystyle 0\leq x_{i}\leq 1}$. Then the multilinear extension is defined as ${\displaystyle F(\mathbf {x} )=\sum _{S\subseteq \Omega }f(S)\prod _{i\in S}x_{i}\prod _{i\notin S}(1-x_{i})}$.

Convex closure[]

Consider any vector ${\displaystyle \mathbf {x} =\{x_{1},x_{2},\dots ,x_{n}\}}$ such that each ${\displaystyle 0\leq x_{i}\leq 1}$. Then the convex closure is defined as ${\displaystyle f^{-}(\mathbf {x} )=\min \left(\sum _{S}\alpha _{S}f(S):\sum _{S}\alpha _{S}1_{S}=\mathbf {x} ,\sum _{S}\alpha _{S}=1,\alpha _{S}\geq 0\right)}$. The convex closure of any set function is convex over ${\displaystyle [0,1]^{n}}$. It can be shown that ${\displaystyle f^{L}(\mathbf {x} )=f^{-}(\mathbf {x} )}$ for submodular functions.

Concave closure[]

Consider any vector ${\displaystyle \mathbf {x} =\{x_{1},x_{2},\dots ,x_{n}\}}$ such that each ${\displaystyle 0\leq x_{i}\leq 1}$. Then the concave closure is defined as ${\displaystyle f^{+}(\mathbf {x} )=\max \left(\sum _{S}\alpha _{S}f(S):\sum _{S}\alpha _{S}1_{S}=\mathbf {x} ,\sum _{S}\alpha _{S}=1,\alpha _{S}\geq 0\right)}$.

Properties[]

1. The class of submodular functions is closed under non-negative linear combinations. Consider any submodular function ${\displaystyle f_{1},f_{2},\ldots ,f_{k}}$ and non-negative numbers ${\displaystyle \alpha _{1},\alpha _{2},\ldots ,\alpha _{k}}$. Then the function ${\displaystyle g}$ defined by ${\displaystyle g(S)=\sum _{i=1}^{k}\alpha _{i}f_{i}(S)}$ is submodular.
2. For any submodular function ${\displaystyle f}$, the function defined by ${\displaystyle g(S)=f(\Omega \setminus S)}$ is submodular.
3. The function ${\displaystyle g(S)=\min(f(S),c)}$, where ${\displaystyle c}$ is a real number, is submodular whenever ${\displaystyle f}$ is monotone submodular. More generally, ${\displaystyle g(S)=h(f(S))}$ is submodular, for any non decreasing concave function ${\displaystyle h}$.
4. Consider a random process where a set ${\displaystyle T}$ is chosen with each element in ${\displaystyle \Omega }$ being included in ${\displaystyle T}$ independently with probability ${\displaystyle p}$. Then the following inequality is true ${\displaystyle \mathbb {E} [f(T)]\geq pf(\Omega )+(1-p)f(\varnothing )}$ where ${\displaystyle \varnothing }$ is the empty set. More generally consider the following random process where a set ${\displaystyle S}$ is constructed as follows. For each of ${\displaystyle 1\leq i\leq l,A_{i}\subseteq \Omega }$ construct ${\displaystyle S_{i}}$ by including each element in ${\displaystyle A_{i}}$ independently into ${\displaystyle S_{i}}$ with probability ${\displaystyle p_{i}}$. Furthermore let ${\displaystyle S=\cup _{i=1}^{l}S_{i}}$. Then the following inequality is true ${\displaystyle \mathbb {E} [f(S)]\geq \sum _{R\subseteq [l]}\Pi _{i\in R}p_{i}\Pi _{i\notin R}(1-p_{i})f(\cup _{i\in R}A_{i})}$.^{[citation needed]}

Optimization problems[]

Submodular functions have properties which are very similar to convex and concave functions. For this reason, an optimization problem which concerns optimizing a convex or concave function can also be described as the problem of maximizing or minimizing a submodular function subject to some constraints.

Submodular set function minimization[]

The simplest minimization problem is to find a set ${\displaystyle S\subseteq \Omega }$ which minimizes a submodular function; this is the unconstrained problem. This problem is computable in (strongly)^[8][9] polynomial time.^[10][11] Computing the minimum cut in a graph is a special case of this general minimization problem. However, adding even a simple constraint such as a cardinality lower bound makes the minimization problem NP hard, with polynomial factor lower bounds on the approximation factor.^[12][13]

Submodular set function maximization[]

Unlike the case of minimization, maximizing a submodular functions is NP-hard even in the unconstrained setting. Theory and enumeration algorithms for finding local and global maxima (minima) of submodular (supermodular) functions can be found in B. Goldengorin. European Journal of Operational Research 198(1):102-112, DOI: 10.1016/j.ejor.2008.08.022. For instance max cut is a special case even when the function is required only to be non-negative. The unconstrained problem can be shown to be inapproximable if it is allowed to be negative. There has been extensive work on constrained submodular function maximization when the functions are non-negative. Typically, the approximation algorithms for these problems are based on either greedy algorithms or local search algorithms. The problem of maximizing a non-negative symmetric submodular function admits a 1/2 approximation algorithm.^[14] Computing the maximum cut of a graph is a special case of this problem. The more general problem of maximizing a non-negative submodular function also admits a 1/2 approximation algorithm.^[15] The problem of maximizing a monotone submodular function subject to a cardinality constraint admits a ${\displaystyle 1-1/e}$ approximation algorithm.^{[16][page needed][17]} The maximum coverage problem is a special case of this problem. The more general problem of maximizing a monotone submodular function subject to a matroid constraint also admits a ${\displaystyle 1-1/e}$ approximation algorithm.^[18][19][20] Many of these algorithms can be unified within a semi-differential based framework of algorithms.^[13]

Related optimization problems[]

Apart from submodular minimization and maximization, another natural problem is Difference of Submodular Optimization.^[21][22] Unfortunately, this problem is not only NP hard, but also inapproximable.^[22] A related optimization problem is minimize or maximize a submodular function, subject to a submodular level set constraint (also called submodular optimization subject to submodular cover or submodular knapsack constraint). This problem admits bounded approximation guarantees.^[23] Another optimization problem involves partitioning data based on a submodular function, so as to maximize the average welfare. This problem is called the submodular welfare problem.^[24]

Applications[]

Submodular functions naturally occur in several real world applications, in economics, game theory, machine learning and computer vision. Owing to the diminishing returns property, submodular functions naturally model costs of items, since there is often a larger discount, with an increase in the items one buys. Submodular functions model notions of complexity, similarity and cooperation when they appear in minimization problems. In maximization problems, on the other hand, they model notions of diversity, information and coverage. For more information on applications of submodularity, particularly in machine learning, see ^[4][25][26]

See also[]

• Supermodular function
• Matroid, Polymatroid
• Utility functions on indivisible goods

Citations[]

References[]

• Schrijver, Alexander (2003), Combinatorial Optimization, Springer, ISBN 3-540-44389-4
• Lee, Jon (2004), A First Course in Combinatorial Optimization, Cambridge University Press, ISBN 0-521-01012-8
• Fujishige, Satoru (2005), Submodular Functions and Optimization, Elsevier, ISBN 0-444-52086-4
• Narayanan, H. (1997), Submodular Functions and Electrical Networks, ISBN 0-444-82523-1
• Oxley, James G. (1992), Matroid theory, Oxford Science Publications, Oxford: Oxford University Press, ISBN 0-19-853563-5, Zbl 0784.05002

External links[]

• http://www.cs.berkeley.edu/~stefje/references.html has a longer bibliography