Paper:
Minimizing a Submodular Function from Samples
In this paper we consider the problem of minimizing a submodular function from training data. Submodular functions can be efficiently minimized and are consequently heavily applied in machine learning. There are many cases, however, in which we do not know the function we aim to optimize, but rather have access to training data that is used to learn it. In this paper we consider the question of whether submodular functions can be minimized when given access to its training data. We show that even learnable submodular functions cannot be minimized within any non-trivial approximation when given access to polynomially-many samples. Specifically, we show that there is a class of submodular functions with range in [0, 1] such that, despite being PAC-learnable and minimizable in polynomial-time, no algorithm can obtain an approximation strictly better than 1/2 o(1) using polynomially-many samples drawn from any distribution. Furthermore, we show that this bound is tight via a trivial algorithm that obtains an approximation of 1/2.

Review:
Informally stated, the paper states the following.
There exists a family F of submodular functions valued in [0, 1], so that if after polynomially many samples any algorithm estimates some set S as the minimizer, there exists some f in F that is consistent with the samples and on which S is suboptimal optimum by at least 1/2.
The argument is based on the probabilistic method, i.e., by showing that such functions exist with non-vanishing probability.
Moreover, the class that is considered is shown to be PAC-learnable, effectively showing that for submodular functions PAC-learnability does not imply "minimizability" from samples.
However, one has to be careful about this argument.
Namely, PAC learnability guarantees low error only on samples coming from the training distribution, while to minimize the function it is necessary to approximate it well near the optimum (where we might lack samples) up to an additive factor.
Thus, it is not very surprising to me that the former does not imply the latter and that one can construct function classes to fool any distribution.
This is a purely theoretical result that both answers an interesting question, and also provides with a construction of intuitive class of submodular functions which are (as proven) hard to minimize from samples.
I only read the (sketched) proofs in the main text, which seem to be very rigorously executed.
The paper is very dense, which I do not think can be easily improved on, but it is nicely written.
Based on this grounds, I would recommend this paper for acceptance.
I think it would be more interesting if the authors present some natural settings where the learning setting they consider is used in practice.
In other words, when does learn a submodular functions from training data of the form (S_i, F(S_i))?
This setting is different from how submodular functions are learned in the vision domain (where one might argue is the primary application field of submodular minimization).
There, one is given a set of instances, which we want to have a low energy, either using a min-margin approach, or via maximum likelihood.
Questions / Remarks:
1. The argument l129-132 is unclear to me. Do you want to say that as n -> oo both P_{f\in F}(f(S) >= 1/2) and |F'|/|F| go to one, there has to exist an f in |F'| satisfying f(S) >= 1/2 - o(1)?
2. In prop. 2 you do not define F_A(B), which I guess it means F(A \cup B) - F(B), which is confusing as it is typically defined the other way round. You could also use F(A | B), which is I guess the more commonly used notation.
3. I think the discussion after prop 2 and before section 3 should be moved somewhere else, as it does not following from the proposition itself.
4. Maybe add a note that U(A) means drawn uniformly at random from the set A.
5. l212: It is not motivated why you need one predictor for each possible masking element.
6. When you minimize a symmetric submodular function it is to minimize it over all sets different from the empty set and the ground set. This is exactly what is done using Queyranne's algorithm in cubic time.
7. Shouldn't the 1 be 1/2 in the last equation of the statement of lemma 4?
8. I find the notation in the two inequalities in Lemma 5 (after "we get") confusing. Does it mean that it holds for any h(n) = o(1) function with the given probability? Also, maybe add an "and" after the comma there.
9. l98: a word is missing after "desirable optimization"