Optimal lower bounds for universal relation, and for samplers and finding duplicates in streams
Michael Kapralov, Jelani Nelson, Jakub Pachocki, Zhengyu Wang, David, P. Woodruff, Mobin Yahyazadeh

TL;DR
This paper establishes tight lower bounds for the universal relation problem in communication complexity, leading to optimal bounds for sampling and duplicate detection in data streams, using novel proof techniques involving encoding and reductions.
Contribution
It provides the exact randomized one-way communication complexity of the universal relation problem and introduces two innovative proofs, including a new reduction from Augmented Indexing.
Findings
Lower bounds match upper bounds for the problem.
Optimal bounds for $ ext{ell}_p$-sampling in turnstile streams.
Efficient duplicate detection in streaming models.
Abstract
In the communication problem (universal relation) [KRW95], Alice and Bob respectively receive with the promise that . The last player to receive a message must output an index such that . We prove that the randomized one-way communication complexity of this problem in the public coin model is exactly for failure probability . Our lower bound holds even if promised . As a corollary, we obtain optimal lower bounds for -sampling in strict turnstile streams for , as well as for the problem of finding duplicates in a stream. Our lower bounds do not need to use large weights, and hold even if promised at all points in the stream. We give two different proofs of our main…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Optimal lower bounds for universal relation, and for samplers and finding duplicates in streams111This paper is a merger of [NPW17],
and of work of Kapralov, Woodruff, and Yahyazadeh.
Michael Kapralov EPFL. [email protected].
Jelani Nelson Harvard University. [email protected]. Supported by NSF grant IIS-1447471 and CAREER award CCF-1350670, ONR Young Investigator award N00014-15-1-2388, and a Google Faculty Research Award.
Jakub Pachocki OpenAI. [email protected]. Work done while affiliated with Harvard University, under the support of ONR grant N00014-15-1-2388.
Zhengyu Wang Harvard University. [email protected]. Supported by NSF grant CCF-1350670.
David P. Woodruff IBM Research Almaden. [email protected].
Mobin Yahyazadeh Sharif University of Technology. [email protected]. Work done while an intern at EPFL.
Abstract
In the communication problem (universal relation) [KRW95], Alice and Bob respectively receive with the promise that . The last player to receive a message must output an index such that . We prove that the randomized one-way communication complexity of this problem in the public coin model is exactly for failure probability . Our lower bound holds even if promised . As a corollary, we obtain optimal lower bounds for -sampling in strict turnstile streams for , as well as for the problem of finding duplicates in a stream. Our lower bounds do not need to use large weights, and hold even if promised at all points in the stream.
We give two different proofs of our main result. The first proof demonstrates that any algorithm solving sampling problems in turnstile streams in low memory can be used to encode subsets of of certain sizes into a number of bits below the information theoretic minimum. Our encoder makes adaptive queries to throughout its execution, but done carefully so as to not violate correctness. This is accomplished by injecting random noise into the encoder’s interactions with , which is loosely motivated by techniques in differential privacy. Our correctness analysis involves understanding the ability of to correctly answer adaptive queries which have positive but bounded mutual information with ’s internal randomness, and may be of independent interest in the newly emerging area of adaptive data analysis with a theoretical computer science lens. Our second proof is via a novel randomized reduction from Augmented Indexing [MNSW98] which needs to interact with adaptively. To handle the adaptivity we identify certain likely interaction patterns and union bound over them to guarantee correct interaction on all of them. To guarantee correctness, it is important that the interaction hides some of its randomness from in the reduction.
1 Introduction
In turnstile -sampling, a vector starts as the zero vector and receives coordinate-wise updates of the form “” for . During a query, one must return a uniformly random element from . The problem was first defined in [FIS08], where a data structure (or “sketch”) for solving it was used to estimate the Euclidean minimum spanning tree, and to provide -approximations of a point set in a geometric space (that is, one wants to maintain a subset such that for any set in a family of bounded VC-dimension, such as the set of all axis-parallel rectangles, ). Sketches for -sampling were also used to solve various dynamic graph streaming problems in [AGM12a] and since then have been crucially used in almost all known dynamic graph streaming algorithms222The spectral sparsification algorithm of [KLM*+*14] is a notable exception., such as for: connectivity, -connectivity, bipartiteness, and minimum spanning tree [AGM12a], subgraph counting, minimum cut, and cut-sparsifier and spanner computation [AGM12b], spectral sparsifiers [AGM13], maximal matching [CCHM15], maximum matching [AGM12a, BS15, Kon15, AKLY16, CCE*+*16, AKL17], vertex cover [CCHM15, CCE*+*16], hitting set, -matching, disjoint paths, -colorable subgraph, and several other maximum subgraph problems [CCE*+*16], densest subgraph [BHNT15, MTVV15, EHW16], vertex and hyperedge connectivity [GMT15], and graph degeneracy [FT16]. For an introduction to the power of -sketches in designing dynamic graph stream algorithms, see the recent survey of McGregor [McG14, Section 3]. Such sketches have also been used outside streaming, such as in distributed algorithms [HPP*+*15, PRS16] and data structures for dynamic connectivity [KKM13, Wan15, GKKT15].
Given the rising importance of -sampling in algorithm design, a clear task is to understand the exact complexity of this problem. The work [JST11] gave an -bit space lower bound for data structures solving even the case which fail with constant probability, and otherwise whose query responses are -close to uniform in statistical distance. They also gave an upper bound for with failure probability , which in fact gave uniform samples from the support of , using space (here denotes ). Thus we say their data structure actually solves the harder problem of -samplingk for with failure probability , where in -samplingk the goal is to recover uniformly random elements, without replacement, from . The upper and lower bounds in [JST11] thus match up to a constant factor for and a constant. We note though in many settings, even if the final application desires constant failure probability, -samplingk with either failure probability or (or both) is needed as a subroutine (see Figure 1).
Universal relation.
The work of [JST11] obtains its lower bound for -sampling (and some other problems) via reductions from universal relation (). The problem was first defined in [KRW95] and arose in connection with work of Karchmer and Wigderson on circuit depth lower bounds [KW90]. For , is the minimum depth of a fan-in circuit over the basis computing . Meanwhile, the (deterministic) communication complexity is defined as the minimum number of bits that need to be communicated in a correct protocol for Alice and Bob to solve the following communication problem: Alice receives and Bob receives (and hence in particular ), and they must both agree on an index such that . It is shown in [KW90] that , where they then used this correspondence to show a tight depth lower bound for monotone circuits solving undirected - connectivity. The work of [KRW95] then proposed a strategy to separate the complexity classes and : start with a function on bits requiring depth , then “compose” it with itself times (see [KW90] for a precise definition of composition). If one could prove a strong enough direct sum theorem for communication complexity after composition, even for a random , such a -fold composition would yield a function that is provably in (and in fact, even in ), but not in . Proving such a direct sum theorem is still wide open, and the statement that it is true is known as the “KRW conjecture”; see for example the recent works [GMWW14, DM16] toward resolving this conjecture. As a toy problem en route to resolving it, [KRW95] suggested proving a direct sum theorem for -fold composition of a particular function that they defined. That task was positively resolved in [EIRS91] (see also [HW90]).
The problem abstracts away the function , and Alice and Bob are simply given with the promise that . The players must then agree on any index with . The deterministic communication complexity of is nearly completely understood, with upper and lower bounds that match up to an additive bits, even if one imposes an upper bound on the number of rounds of communication [TZ97]. Henceforth we also consider a generalized problem , where the output must be distinct indices on which differ. We also use to denote the variants when promised , and also Bob knows . Clearly can only be harder than , respectively.
More than twenty years after its initial introduction in connection with circuit depth lower bounds, Jowhari et al. in [JST11] demonstrated the relevance of in the randomized one-way communication model for obtaining space lower bounds for certain streaming problems, such as various sampling problems and finding duplicates in streams. In the one-way version, Bob simply needs to find such an index after a single message from Alice, and we only charge Alice’s single message’s length as the communication cost. If denotes the randomized one-way communication complexity of in the public coin model with failure probability , [JST11] showed that the space complexity of FindDuplicate with failure probability is at least . In FindDuplicate, one is given a length- stream of integers in , and the algorithm must output some element which appeared at least twice in the stream (note that at least one such element must exist, by the pigeonhole principle). The work [JST11] then showed a reduction demonstrating that any solution to -sampling with failure probability in turnstile streams immediately implies a solution to FindDuplicate with failure probability at most in the same space (and thus the space must be at least ). The same result is shown for -sampling for any , in which the output index should equal with probability , and a similar result is shown even if the distribution on only has to be close to this -distribution in variational distance (namely, the distance should be bounded away from ). It is then shown in [JST11] that for any bounded away from . The approach used though unfortunately does not provide an improved lower bound for .
Seemingly unnoticed in [JST11], we first point out here that the lower bound proof for in that work actually proves the same lower bound for the promise problem . This observation has several advantages. First, it makes the reductions to the streaming problems trivial (they were already quite simple when reducing from , but now they are even simpler). Second, a simple reduction from to sampling problems provides space lower bounds even in the strict turnstile model, and even for the simpler support-finding streaming problem for which when queried is allowed to return any element of , without any requirement on the distribution of the index output. Both of these differences are important for the meaningfulness of the lower bound. This is because in dynamic graph streaming applications, typically is indexed by for some graph on vertices, and is the number of copies of edge in some underlying multigraph. Edges then are not deleted unless they had previously been inserted, thus only requiring correctness for strict turnstile streams. Also, for every single application mentioned in the first paragraph of Section 1 (except for the two applications in [FIS08]), the known algorithmic solutions which we cited as using -sampling as a subroutine actually only need a subroutine for the easier support-finding problem. Finally, third and most relevant to our current work’s main focus, the straightforward reductions from to the streaming problems we are considering here do not suffer any increase in failure probability, allowing us to transfer lower bounds on for small to lower bounds on various streaming problems for small . The work [JST11] could not provide lower bounds for the streaming problems considered there in terms of for small .
We now show simple reductions from to FindDuplicate and from to support-findingk. In support-findingk we must report elements in . In the claims below, is the failure probability for the considered streaming problem.
Claim 1**.**
Any one-pass streaming algorithm for FindDuplicate must use space.
Proof.
We reduce from . Suppose there were a space- algorithm for FindDuplicate. Alice creates a stream consisting of all elements of and runs on those elements, then sends the memory contents of to Bob. Bob then continues running on arbitrarily chosen elements of . Then there must be a duplicate in the resulting concatenated stream, satisfies iff is a duplicate. ∎
Claim 2**.**
Any one-pass streaming algorithm for support-findingk in the strict turnstile model must use bits of space, even if promised that at all points in the stream.
Proof.
This is again via reduction from . Let be a space- algorithm for support-findingk in the strict turnstile model. For each , Alice sends the update to . Alice then sends the memory contents of to Bob. Bob then for each sends the update to . Now note that is exactly the indicator vector of the set . ∎
Claim 3**.**
Any one-pass streaming algorithm for -sampling for any in the strict turnstile model must use bits of space, even if promised at all points in the stream.
Proof.
This is via straightforward reduction from support-findingk, since reporting elements of satisfying some distributional requirements is only a harder problem than finding any elements of . ∎
The reductions above thus raise the question: what is the asymptotic behavior of ?
Our main contribution:
We prove for any bounded away from and , where . Given known upper bounds in [JST11], our lower bounds are optimal for FindDuplicate, support-finding, and -sampling for any for nearly the full range of (namely, for ). Also given an upper bound of [JST11], our lower bound is optimal for -samplingk for nearly the full range of parameters (namely, for ). Previously no lower bounds were known in terms of (or ). Our main theorem:
Theorem 1**.**
For any bounded away from and , .
We give two different proofs of Theorem 1 (in Sections 3 and 4). Our upper bound is also new, though follows by minor modifications of the upper bound in [JST11] and thus we describe it in the appendix. The previous upper bound was . We also mention here that it is known that the upper bound for both and -samplingk in two rounds (respectively, two passes) is only [JST11]. Thus, one cannot hope to extend our new lower bound to two or more passes, since it simply is not true.
1.1 Related work
The question of whether -sampling is possible in low memory in turnstile streams was first asked in [CMR05, FIS08]. The work [FIS08] applied -sampling as a subroutine in approximating the cost of the Euclidean minimum spanning tree of a subset of a discrete geometric space subject to insertions and deletions. The algorithm given there used space bits to achieve failure probability (though it is likely that the space could be improved to with a worse failure probability, by replacing a subroutine used there with a more recent -estimation algorithm of [KNW10]). As mentioned, the currently best known upper bound solves -samplingk using bits [JST11], which Theorem 1 shows is tight.
For -sampling, conditioned on not failing, the data structure should output with probability . The first work to realize its importance came even earlier than for -sampling: [CK04] showed that an -sampler using small memory would lead to a nearly space-optimal streaming algorithm for multiplicatively estimating in the turnstile model, but did not know how to implement such a data structure. The first implementation was given in [MW10], achieving space with . . For the space was improved to bits for constant [AKO11]. In [JST11] this bound was improved to bits for failure probability when and . For the space bound achieved by [JST11] was a factor worse: bits.
For finding a duplicate item in a stream, the question of whether a space-efficient randomized algorithm exists was asked in [Mut05, Tar07]. The question was positively resolved in [GR09], which gave an -space algorithm with constant failure probability. An improved algorithm was given in [JST11], using bits of space for failure probability .
2 Overview of techniques
We now describe our two proofs of Theorem 1. For the upper bound, [JST11] achieved , but in the appendix we show that slight modifications to their approach yield . Our main contribution is in proving an improved lower bound. Assume for some sufficiently small constant (since otherwise we already obtain an lower bound). In both our lower bound proofs in this regime, the proof is split into two parts: we show and separately. We give an overview the former here, which is the more technically challenging half. Our two proofs of the latter are in Sections 3.2 and 4.2.
2.1 Lower bound proof via encoding subsets and an adaptivity lemma
Our first proof of the lower bound on is via an encoding argument. Fix . A randomized encoder is given a set with and must output an encoding , and a decoder sharing public randomness with the encoder must be able to recover given only . We consider such schemes in which the decoder must succeed with probability , and the encoding length is a random variable. Any such encoding must use bits in expectation for some .
There is a natural, but sub-optimal approach to using a public-coin one-way protocol for to devise such an encoding/decoding scheme. The encoder pretends to be Alice with input being the indicator set of , then lets be the message Alice would have sent to Bob. The decoder attempts to recover by iteratively pretending to be Bob times, initially pretending to have input , then iteratively adding elements found in to ’s support. Henceforth let denote the indicator vector of a set .
One might hope to say that if the original failure probability were , then by a union bound, with constant probability every iteration succeeds in finding a new element of (or one could even first apply some error-correction to so that the decoder could recover even if only a constant fraction of iterations succeeded). The problem with such thinking though is that this decoder chooses ’s adaptively! To be specific, being a correct protocol means
[TABLE]
where is the public random string that both Alice and Bob have access to. The issue is that even in the second iteration (when ), Bob’s “input” depends on , since depends on the outcome of the first iteration! Thus the guarantee of (1) does not apply.
One way around the above issue is to realize that as long as every iteration succeeds, is always a subset of . Thus it suffices for the following event to occur: . Then by a union bound, which is at most for . We have thus just shown that .
Our improvement is as follows. Our new decoder again iteratively tries to recover elements of as before. We will give up though on having iterations and hoping for all (or even most) of them to succeed. Instead, we will only have iterations, and our aim is for the decoder to succeed in finding a new element in for at least a constant fraction of these iterations. Simplifying things for a moment, let us pretend for now that all iterations do succeed in finding a new element. will then be Alice’s message , together with the set of size not recovered during the rounds, explicitly written using bits. If the decoder can then recover these remaining elements, this then implies the decoder has recovered , and thus we must have . The decoder proceeds as follows. Just as before, initially the decoder starts with and lets be the output of Bob on and adds it to . Then in iteration , before proceeding to the next iteration, the decoder randomly picks some elements from and adds them into , so that the number of elements left to be uncovered is some fixed number . These extra elements being added to should be viewed as “random noise” to mask information about the random string used by , an idea very loosely inspired by ideas in differential privacy. For intuition, as an example suppose the iteration succeeds in finding some . If the decoder were then to add to , as well as random elements from to , then the resulting reveals only bit of information about (and hence about ). This is as opposed to the bits could have revealed if the masking were not performed. Thus the next query in round , although correlated with , has very weak correlation after masking and we thus might hope for it to succeed. This intuition is captured in the following lemma, which we prove in Section 3.1:
Lemma 1**.**
Consider : and uniformly random. If where , then for any random variable supported on ,
[TABLE]
where is the mutual information between and , and is the binary entropy function.
Fix some . One should imagine here that is iff fails when Alice has input and Bob has input in a instance, and the public random string is . Then the lemma states that if is not arbitrary, but rather random (and correlated with ), then the failure probability of the protocol is still bounded as long as the mutual information between and is bounded. It is also not hard to see that this lemma is sharp up to small additive terms. Consider the case , and iff . Then if is uniform, for all we have . Now consider the case where is random and equal to with probability and is uniform in with probability . Then in expectation reveals bits of , so that . It is also not hard to see that .
In light of the strategy stated so far and Lemma 1, the path forward is clear: at each iteration , we should add enough random masking elements to to keep the mutual information between and all previously added elements below, say, . Then we expect a constant fraction of iterations to succeed. The encoder knows which iterations do not succeed since it shares public randomness with the decoder (and can thus simulate it), so it can simply tell the decoder which rounds are the failed ones, then explicitly include in correct new elements of for the decoder to use in the place of Bob’s wrong output in those rounds. A calculation shows that if one adds a fraction of the remaining items in to after drawing one more support element from Bob, the mutual information between the next query to Bob and the randomness used by will be (see Lemma 5). Thus we do this for a sufficiently small constant times . We will then have . Note that we cannot continue in this way once (since the number of “random noise” elements we inject should at least be one). Thus we are forced to stop after iterations. We then set , so that as desired.
The argument for lower bounding is a bit simpler, and in particular does not need rely on Lemma 1. Both the idea and rigorous argument can be found in Section 3.2, but again the idea is to use a protocol for this problem to encode appropriately sized subsets of .
As mentioned above, our lower bounds use protocols for and to establish protocols for encoding subsets of some fixed size of . These encoders always consist of some message Alice would have sent in a or protocol, together with a random subset (using bits, to represent both and the set itself). Here is a random variable. These encoders are thus Las Vegas: the length of the encoding is a random variable, but the encoder/decoder always succeed in compressing and recovering the subset. The final lower bounds then come from the following simple lemma, which follows from the source coding theorem.
Lemma 2**.**
Let s denote the number of bits used by the or protocol, and let denote the expected number of bits to represent . Then . In particular, .
Section 3.1 provides our first proof that . We extend our results in Section 3.2 to for , proving a lower bound of communication even for constant failure probability.
2.2 Lower bound proof via reduction from
Our second lower bound proof for is via a randomized reduction from [MNSW98]. In this problem, Charlie receives and Diane receives together with for , and Diane must output . It is shown in [MNSW98] that for any bounded away from . In our reduction, .
For , we can also think of the problem as Alice being given and Bob being given , and Bob must output some element of . In , Charlie views his input as blocks of bits of nearly equal size, where the th block represents a subset of in some collection of sets, for some carefully chosen universe sizes per block. Here is a collection of subsets of of size of maximal size such any two sets in the collection have intersection size strictly less than . Furthermore, Diane’s index is in some particular block of bits corresponding to some set , and Diane also knows for .
Now we explain the reduction. We assume some protocol for , and we give a protocol for . First, we define the universe , which has size . Charlie then defines . Charlie and Diane use public randomness to define a uniformly random permutation on . Charlie can compute . Also, since Diane knows for , she can define and compute . Then and are the inputs to Alice and Bob in the protocol for . Charlie sends Diane the message Alice would have sent Bob in if her input had been , and Diane simulates Bob to recover an element in . Importantly, Alice and Bob do not know anything about at this point other than that and . Thus, the protocol for , if it succeeds, outputs an arbitrary element , which is a deterministic function of the labels of elements in and and the randomness that Alice and Bob share, which is independent from the randomness in . Since is still a uniformly random map conditioned on and for each , and , it follows that is a uniformly random element of . After receiving , if , then Charlie and Diane reveal the pairs for each to Alice and Bob and Bob updates his set to include for each . One can show that at each step in this process, if Alice and Bob succeed in outputting an arbitrary item from , then this is a uniformly random item from . The fact that this item is uniformly random is crucial for arguing the number of computation paths of the protocol of Alice and Bob is with good probability, over , so that one can argue (see below) that with good probability on every such computation path Alice and Bob succeed on that path, over their randomness . Although the idea of using a random permutation appeared in [JST11] to show that any public coin protocol can be made into one in which a uniformly random element of is output, here we must use this idea adaptively, slowly revealing information about and arguing that this property is maintained for each of Bob’s successive queries.
Due to geometrically increasing repetitions of items for increasing , a uniformly random element in is roughly times more likely to correspond to an item in than in for . Thus if Diane simulates Bob to recover a random element in , it is most likely to recover an element of . She can then tell Bob to include and its redundant copies to and iterate.
There are several obstacles to overcome to make this work. First, iterating means using adaptively, which was the same issue that arose in Section 2.1. Second, a constant fraction of the time (), we expect to obtain an element not in , but rather from some for . If this happened too often, then Diane would need to execute many queries to recover a sufficiently large number of elements from in order to solve . This would then require a union bound over too many possible computation paths, which would not be possible as Alice likely would fail on one of them (over the choice of ). However, since the random permutation argument above ensures that at each step we receive a uniformly random item from the current set , if we continue for iterations, we can argue that with large probability, our sequence of inputs over the iterations with which Diane invokes Bob’s output are all likely to come from a family of size at most . Here we need to carefully construct this family to contain a smaller number of sets from levels for which is larger so that the overall number of sets is small. Given this, we can union bound over all such , for total failure probability . Furthermore, we can also argue that after iterations, it is likely that we have recovered at least of the elements from , which is enough to uniquely identify by the limited intersection property of .
3 Lower bounds via the adaptivity lemma
3.1 Communication Lower Bound for
Consider a protocol for with failure probability , operating in the one-way public coin model. When Alice’s input is and Bob’s is , Alice sends to Bob, and Bob outputs , which with probability at least is in . As mentioned in Section 2, we use as a subroutine in a scheme for encoding/decoding elements of for . We assume , since for larger we have an lower bound.
3.1.1 Encoding/decoding scheme
We now describe our encoding/decoding scheme for elements in , which uses in a black-box way. The parameters shared by ENC and DEC are given in Algorithm 2.
As discussed in Section 2, on input , ENC computes as part of its output. Moreover, ENC also outputs a subset computed as follows. Initially and . ENC proceeds in rounds. In round , ENC computes . Let denote a binary string of length , where records whether succeeds in round . ENC also outputs . If , i.e. succeeds, ENC sets and removes from (since the decoder can recover from the -protocol, ENC does not need to include it in ); otherwise ENC sets . At the end of round , ENC picks a uniformly random set in . In particular, ENC uses its shared randomness with DEC to generate in such a way that agree on the sets (DEC will actually iteratively construct ). We present ENC in Algorithm 3.
The decoding process is symmetric. Let and . DEC proceeds in rounds. On round , DEC obtains by invoking . By construction of (to be described later), it is guaranteed that . Therefore, DEC recovers exactly the same as ENC. DEC initially assigns . If , DEC adds to both and . At the end of round , DEC inserts many random items from into so that . DEC can achieve this because of the shared random permutation when constructing . In the end, DEC outputs . We present DEC in Algorithm 4.
3.1.2 Analysis
We have two random objects in our encoding/decoding scheme: (1) the random source used by , denoted by , and (2) the random permutation . These are independent.
First, we can prove that . That is, for any fixing of the randomness in and , DEC will always decode successfully. It is because ENC and DEC share and , so that DEC essentially simulates ENC. We formally prove this by induction in Lemma 3.
Now our goal is to prove that by using the -protocol, the number of bits that ENC saves in expectation over the naive -bit encoding is bits. Intuitively, it is equivalent to prove the number of elements that ENC saves is . We formalize this in Lemma 4. Note that ENC also needs to output (i.e., whether the succeeds on rounds), which takes bits. By our setting of parameters, we can afford the loss of bits. Thus it is sufficient to prove .
We have . In Lemma 1, we prove the probability that fails on round is upper bounded by , where is the mutual information between and . Furthermore, we will show in Lemma 5 that is upper bounded by . By our setting of parameters, we have and thus .
Lemma 3**.**
.
Proof.
We claim that for , is a partition of ( is defined in Algorithm 3, and in Algorithm 4). We prove the claim by induction on . Our base case is , for which the claim holds since , .
Assume the claim holds for (), and we consider round . On round , by induction , the index obtained by both ENC and DEC are the same. Initially and , and so is a partition of . If is a valid sample (i.e. ), then , and ENC removes from and in the meanwhile DEC inserts into , so that remains a partition of . Next, ENC repeats removing the from with the smallest value until . Symmetrically, DEC repeats inserting the into with the smallest value among , until . In the end we have , so ENC and DEC execute repetition the same number of times. Moreover, we can prove that during the same iteration of this repeated insertion, the element removed from is exactly the same element inserted to . This is because in the beginning of a repetition is a partition of . We have . Let denote that minimizes . Then (since will be removed from , it has no chance to be included in in ENC, so that contains ), and is also the smallest among . Thus both ENC and DEC will take (for ENC, to remove from , and for DEC, to insert into ). Therefore, remains a partition of .
Given the fact that is a partition of , the are the same in ENC and DEC. Furthermore, are the same in ENC and DEC. We know . Since ENC outputs , and DEC outputs , we have . ∎
Lemma 4**.**
Let be a random variable with and . Then .
Proof.
[TABLE]
Taking expectation on both sides, we have . ∎
Lemma 1 (restated). Consider : and uniformly random. If where , then for any r.v. supported on ,
[TABLE]
where is the mutual information between and , and is the binary entropy function.
Proof.
It is equivalent to prove
[TABLE]
By definition of mutual entropy , where and we must show
[TABLE]
The upper bound for is obtained by considering the following one-way communication problem: Alice knows both and while Bob only knows , and Alice must send a single message to Bob so that Bob can recover . The expected message length in an optimal protocol is exactly . Thus, any protocol gives an upper bound for , and we simply take the following protocol: Alice prepends a bit to her message iff (taking bits in expectation). Then if , Alice sends directly (taking bits). Otherwise, when , Alice sends the index of in (taking bits). ∎
Corollary 1**.**
Let denote the random source used by the -protocol with failure probability at most . If is a fixed set and , .
Lemma 5**.**
, for .
Proof.
Note that . Since and , . Here is the main idea to lower bound : By definition of conditional entropy, . We fix an arbitrary . If we can prove that for any where , , then by definition of entropy we have .
First we can prove for any fixed ,
[TABLE]
We have . On round (), ENC removes elements (at least of which are chosen all at random) from to obtain . Conditioned on the event that , the probability that is at most , where the equation achieves when , and ENC takes a uniformly random subset of of size , so that the subset does not intersect with .
Next we can prove
[TABLE]
For notational simplicity, let denote . We have
[TABLE]
By telescoping,
[TABLE]
Moreover,
[TABLE]
By our setting of parameters
[TABLE]
Therefore, for ,
[TABLE]
By Taylor series , and thus , for . For , we have .
By Lemma 6, we have . Therefore, the right hand side of (7) is upper bounded by . Together with (6), we prove (4) holds.
Finally, let , we have and thus . Therefore, and so . ∎
Lemma 6**.**
Let and . We have .
Proof.
First, we bound the product of first terms. Note that for . Therefore,
[TABLE]
Then, we bound the product of the rest terms
[TABLE]
Multiplying two parts proves the lemma. ∎
Theorem 2**.**
, given that .
Proof.
By Lemma 3, the success probability of protocol is . By Lemma 2, we have , where . The size of is . By Corollary 1, conditioned on , . By Lemma 5, (Note that when , ). Therefore, . By the setting of parameters (see Algorithm 2) we have . Therefore, . By Lemma 4, . Furthermore, . Thus we obtain . ∎
3.2 Communication Lower Bound for
In this section, we prove the lower bound . In fact, our lower bound holds for any failure probability bounded away from . Let denote a -protocol where Alice sends to Bob, and Bob outputs . We consider the following encoding/decoding scheme for . computes as part of its message. In addition, includes constructed as follows, spending bits. Initially , and proceeds in rounds. Let where is generated by sub-sampling each element in with probability . In round (), tries to obtain elements from by invoking , denoted by , and removes (whose expected size is ) from . Note that is able to recover the elements in . For each round the failure probability of is at most . Thus we have . Furthermore, each element contains bits of information, thus yielding a lower bound of bits.
In this section we assume , since for larger we have an lower bound.
3.2.1 Encoding/decoding scheme
3.2.2 Analysis
Theorem 3**.**
, given that and .
Proof.
Let . Let SUCC denote the event that . Note that . By the Chernoff bound, . In the following, we argue conditioned on SUCC. Namely, in each round , there are at least items in .
Similar to Lemma 3, we can prove the protocol always succeeds. By Lemma 2, we have , where . The size of is . The randomness used by is independent from for every . Therefore, , and is independent from . We have , and thus . By Lemma 4, . Moreover, and . Thus we have . ∎
4 Lower bounds proofs via augmented indexing
Here we show another route to proving via reduction from augmented indexing. We again separately prove lower bounds for and . Both proofs make use of the following standard lemma. The proof can be found in the appendix (see Section A.2).
Lemma 7**.**
For any integers and , there exists a collection with such that for all , .
Both our lower bounds in Sections 4.1 and 4.2 reduce from augmented indexing (henceforth ) to either with low failure probability, or with constant failure probability, in the public coin one-way model of communication. We remind the reader of the setup for the problem. There are two players, Charlie and Diane. Charlie receives and Diane receives together with . Charlie must send a single message to Diane such that Diane can then output . The following theorem is known.
Theorem 4**.**
[MNSW98]** .
We show that if there is an -bit communication protocol for on -bit vectors with failure probability (or for with constant failure probability), that implies the existence of an -bit protocol for for some (or for ). The lower bound on then follows from Theorem 4.
4.1 Communication Lower Bound for
Set . In this section we assume and show . This implies a lower bound of for all bounded away from .
As mentioned, we assume we have an -bit protocol for with failure probability , with players Alice and Bob.We use to give an -bit protocol for , which has players Charlie and Diane, for .
The protocol operates as follows. Without loss of generality we may assume that, using the notation of Lemma 7, is a power of for as in the lemma statement. This is accomplished by simply rounding down to the nearest power of by removing elements arbitrarily. Also, define for some sufficiently small constant to be determined later. Now, Charlie partitions the bits of his input into consecutive sequences of bits such that the th chunk of bits for each can be viewed as specifying an element for and . Lemma 7 gives , which is for . Thus . Given these sets , we now discuss how Charlie generates a vector . Charlie then simulates Alice on to generate the message Alice would have send to Bob in protocol , then sends that same message to Diane.
To generate , assume Charlie and Diane have sampled a bijection from
[TABLE]
to uniformly at random. We denote this bijection by . This is possible since . Then Charlie defines to be the indicator vector , where
[TABLE]
then sends a message to Diane, equal to Alice’s message with input . This completes the description of Charlie’s behavior in the protocol .
We describe how Diane uses to solve . Diane’s input lies in some chunk . We now show how Diane can use to recover with probability (and thus in particular recover ). Since Diane knows for , she knows for . She can then execute the following algorithm.
In Algorithm 8 Diane is building up a subset of . Once , Diane can uniquely recover by the limited intersection property of guaranteed by Lemma 7. Until then, she uses to recover elements of , which, as we now show, are chosen uniformly at random from .
Claim 4**.**
For every protocol for Alice and Bob that uses shared randomness with Bob’s behaviour given by , for every choice of shared random string of Alice and Bob, for every , the following conditions hold. If is a uniformly random permutation, the success or failure of is determined by and the image of under . Conditioned on a choice of , and such that succeeds, one has that is a uniformly random element of .
Proof.
The first claim follows by noting that the message that Alice sends to Bob is solely a function of and . The behaviour of Bob is determined by and (and the latter is determined by ).
Now condition on the values of , and such that succeeds, and let denote the output. Note that by our conditioning is a fixed quantity. The only randomness left is the exact mapping of to . This mapping is independent of and and uniformly random, so is a uniformly random element of , as required. ∎
Fix any protocol (not necessarily the one that Charlie and Diane use; see analysis of the idealized process below). Now fix together with values of , and such that succeeds.
Elements in are unlikely to be recovered.
Given Claim 4, since the elements of appear with frequency in , they are less likely to be returned by when is small. More precisely, as long as , for any
[TABLE]
Here again the probability is over the choice of (recall that we condition on the image under , but not on the actual mapping).
We now define the set of typical intermediate sets, which plays a crucial role in our analysis. Let for denote . Let be the collection of all such that (1) for all , and (2) for each , . The following claim will be useful:
Claim 5**.**
For the set defined above one has .
Proof.
[TABLE]
∎
We will show that for most choices of and shared random string Algorithm 8 (a) never leaves the set and (b) successfully terminates. Note that Algorithm 8 is a random process whose sample space is the product of the set of all possible permutations and shared random strings . As before, we denote this process by . It is useful for analysis purposes to define another process , which is an idealized version of . In this process instead of running Alice runs , which is guaranteed to output an element of for every choice of , shared random string , , and . The proof proceeds in three steps.
Step 1: proving that succeeds in recovering and never leaves with high probability. Choose uniformly at random. By (11), as long as , the expected number of items recovered by from for in the first iterations is at most . Thus the probability of recovering more than items from is at most by Markov’s inequality. Note that the probability is over the choice of only, as is assumed to succeed with probability by definition of . Thus
[TABLE]
In particular this means that with probability at least at most items from are recovered in the first (or fewer, if the algorithm terminates earlier) iterations. This also implies that with probability at least if the algorithm proceeds for the entire iterations, it recovers at least elements of and hence terminates. We thus get that succeeds at least with probability .
Step 2: coupling to on most of the probability space. For every and every let be the probabilistic event (over the choice of ’s random string ) that succeeds in returning an element in . Note that is a subset of the probability space of shared random strings , and depends on . We let
[TABLE]
to simplify notation. Using Claim 5 and the union bound we have for every
[TABLE]
as long as for for a sufficiently small constant.
Now recall that is an idealized protocol, which is guaranteed to output an element of for every choice of , shared random string , , and . We have just shown that for every the event occurs with probability at least over the choice of . Now define as equal to for all (the typical set of intermediate sets) and such that , and extend to return an arbitrary element of for remaining tuples . Note that defined in this way is a deterministic function once , , and are fixed.
Note that with probability at least over the choice of and one has for all , as required.
Step 3: arguing that succeeds with high probability. Choose uniformly at random. By Step 2 we have that with probability at least over this choice for all . At the same time we have by Step 1 that with probability at least over the choice of the idealized process succeeds in recovering and never leaves . Putting the two bounds together, we get that succeeds with probability at least , showing the following theorem.
Theorem 5**.**
For any and integer with , for .
Corollary 2**.**
For any and integer , .
4.2 Communication Lower Bound for
The idea for lower bounding is as in Section 4.1, but slightly simpler. That is because for the protocol for , Diane will not make adaptive queries to Bob in the protocol for . Rather, she will only make one query using Bob and will be able to decide with good probability from that single query. We make use of the following lemma from [JST11], whose proof is similar to our analysis in Section 4.1.
Lemma 8**.**
[JST11]** Any public coin protocol for can be turned into one that outputs every index with with the same probability. The number of bits sent, failure probability, and number of rounds do not change. Similarly, any protocol can be turned into one in which all subsets of of size on which differ are equally likely to be output.
Henceforth we assume outputs random differing indices, which is without loss of generality by Lemma 8.
Again Charlie receives and Diane receives and and they want to solve . Charlie views his input as consisting of blocks for for a sufficiently small constant , and the th block for specifies a set for and . As before, for sufficiently small we have . The bijection and set are defined exactly as in Section 4.1, and Charlie simulates Alice to send the message to Diane that Alice would have sent to Bob on input . Again, Diane knows for , where lies in the th block of bits. Diane’s algorithm to produce her output is then described in Algorithm 9.
Recall Bob, when he succeeds, returns uniformly random elements from . Meanwhile, only has elements for some small constant . As in Section 4.1, almost all of the support of comes from items in block , and hence we expect almost all our samples to come from (and be uniform in) items corresponding to elements of .
We now provide a formal analysis. Henceforth we condition on Bob succeeding, which happens with probability . The number of elements in corresponding to an element of is , whereas the number of elements corresponding to an element of for is
[TABLE]
Thus, we expect at most elements in to correspond to elements in for , and the probability that we have at least such elements in is less than by Markov’s inequality. We henceforth condition on having less than such elements in . Now we know contains at least elements corresponding to , chosen uniformly from . For any given element , the probability that none of the elements in from correspond to is for sufficiently small (where ). Thus the expected number of not covered by is less than . Thus the probability that fewer than elements are covered by is a most by Markov’s inequality (and otherwise, Diane succeeds). Thus, the probability that Diane succeeds is at least . We have thus shown the following theorem.
Theorem 6**.**
For any integers , for .
Corollary 3**.**
For any integers , .
Remark 1**.**
One may wish to understand for near (or at least, larger than ). Such a lower bound is given in Theorem 3. The proof given above as written would yield no lower bound in this regime for since is in fact easy when the failure probability is allowed to be least (Charlie can send no message at all, and Diane can simply guess via a coin flip). One can however get a handle on by instead directly reducing from the following variant of augmented indexing: Charlie receives and Diane receives and and must output , where the are as above. One can show that unless Charlie sends almost his entire input, Diane cannot have success probability significantly better than random guessing (which has success probability ). The proof is nearly identical to the analysis of augmented indexing over large domains [EJS10, JW13]. Indeed, the problem is even almost identical, except that here we consider Charlie receiving a vector whose entries come from different alphabet sizes (since the are different), whereas in [EJS10, JW13] all the entries come from the same alphabet. **
Acknowledgments
Initially the authors were focused on proving optimal lower bounds for samplers, but we thank Vasileios Nakos for pointing out that our lower bound immediately implies a tight lower bound for finding a duplicate in data streams as well. Also, initially our proof of Lemma 1 incurred an additive in the numerator of the right hand side of (2). This is clearly suboptimal for small (for example, consider , in which case the right hand side should be and not )). We thank T.S. Jayram for pointing out that a slight modification of our proof could actually replace the additive with the binary entropy function (and also for showing us a different proof of this lemma, which resembles the standard proof of Fano’s inequality).
Appendix A Appendix
A.1 A tight upper bound for
In [JST11, Proposition 1] it is shown that for . Here we show that a minor modification of their protocol in fact shows the correct complexity , which given our new lower bound, is optimal up to a constant factor for the full range of as long as is bounded away from .
Recall Alice and Bob receive , respectively, and share a public random string. Alice must send a single message to Bob, from which Bob must recover indices for which . Bob is allowed to fail with probability . The fact that is obvious: Alice can simply send the message , and Bob can then succeed with failure probability [math]. We thus now show for some constant , which completes the proof of the upper bound. We assume (otherwise, Alice sends explicitly).
As mentioned, the protocol we describe is nearly identical to one in [JST11] (see also [CF14]). We will describe the new protocol here, then point out the two minor modifications that improve the bound to in Remark 2. We first need the following lemma.
Lemma 9**.**
Let be a finite field and an integer. Then for any , there exists for s.t. for any with , .
Proof.
The proof is via the probabilistic method. iff . Note has . Thus it suffices to show that such a exists with no -sparse vector in its kernel. The number of vectors with is at most . For any fixed , . Thus
[TABLE]
by a union bound. The above is strictly less than for , yielding the claim. ∎
Corollary 4**.**
Let be a finite field and an integer. Then for any , there exists for together with an algorithm such that for any with , .
Proof.
Given Lemma 9, a simple such is as follows. Given some with , loops over all in with and outputs the first one it finds for which . ∎
The protocol for is now as follows. Alice and Bob use public randomness to pick commonly known random functions for , such that for any and for any , . They also agree on a matrix and as described in Corollary 4 for a sufficiently large constant to be determined later, with . Thus has rows. Alice then computes for where , and her message to Bob is . For and an -dimensional vector, denotes the -dimensional vector with for , and for . Note Alice’s message is bits, as desired. Bob then executes the following algorithm and outputs the returned values.
The correctness analysis is then as follows, which is nearly the same as the -sampler of [JST11]. If Alice’s input is and Bob’s is , let , so that can be viewed as an element of . Also let . Then , and since , there either (1) exists a unique such that , or (2) (in which case we define ). Let be the event that simultaneously for all . Let be the event that either we are in case (2), or we are in case (1) and holds. Note that conditioned on both occurring, Bob succeeds by Corollary 4.
We now just need to show . We use the union bound. First, consider . If , then . If , then , which is by the Chernoff bound since . Next we bound . For , we know . Thus, letting denote ,
[TABLE]
for some constant by the Chernoff bound and the fact that . Recall that the Chernoff bound states that for a sum of independent Bernoullis,
[TABLE]
Then by a union bound over and applying (12),
[TABLE]
Remark 2**.**
As already mentioned, the protocol given above and the one described in [JST11] using bits differ in minor points. First: the protocol there used different hash functions , but as seen above, only are needed. This already improves one factor to . The other improvement comes from replacing the -sparse recovery structure with rows used in [JST11] with our Corollary 4. Note the matrix in our corollary has even more rows, but the key point is that the bit complexity is improved. Whereas using a -sparse recovery scheme as described in [JST11] would use linear measurements of a -sparse vector with bits per measurement (for a total of bits), we use measurements with only bits per measurement. The key insight is that we can work over instead of when the entries of are in , which leads to our slight improvement. **
A.2 Proof of the existence of the desired
Lemma 7 (restated). For any integers and , there exists a collection with such that for all , .
Proof.
The proof is via the probabilistic method. We pick independently, each one uniformly at random from . Fix . Imagine being fixed and picking the elements of one by one. Let denote the indicator random variable for the event that the th element picked is also in . Then , and we set , which is by linearity of expectation. We have for . The are not independent, but they are negatively dependent. Thus the Chernoff bound yields
[TABLE]
Setting so that , by a union bound with positive probability for all , simultaneously, as desired. Note for this choice of , we have . ∎
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[AGM 12a] Kook Jin Ahn, Sudipto Guha, and Andrew Mc Gregor. Analyzing graph structure via linear measurements. In Proceedings of the 23 rd ACM-SIAM Symposium on Discrete Algorithms (SODA) , pages 459–467, 2012.
- 2[AGM 12b] Kook Jin Ahn, Sudipto Guha, and Andrew Mc Gregor. Graph sketches: sparsification, spanners, and subgraphs. In Proceedings of the 31 st ACM SIGMOD-SIGACT-SIGART Symposium on Principles of Database Systems (PODS) , pages 5–14, 2012.
- 3[AGM 13] Kook Jin Ahn, Sudipto Guha, and Andrew Mc Gregor. Spectral sparsification in dynamic graph streams. In Proceedings of the 16 th International Workshop on Approximation Algorithms for Combinatorial Optimization Problems (APPROX) , pages 1–10, 2013.
- 4[AKL 17] Sepehr Assadi, Sanjeev Khanna, and Yang Li. On estimating maximum matching size in graph streams. In Proceedings of the 28 th Annual ACM-SIAM Symposium on Discrete Algorithms (SODA) , pages 1723–1742, 2017.
- 5[AKLY 16] Sepehr Assadi, Sanjeev Khanna, Yang Li, and Grigory Yaroslavtsev. Maximum matchings in dynamic graph streams and the simultaneous communication model. In Proceedings of the 27 th Annual ACM-SIAM Symposium on Discrete Algorithms (SODA) , pages 1345–1364, 2016.
- 6[AKO 11] Alexandr Andoni, Robert Krauthgamer, and Krzysztof Onak. Streaming algorithms via precision sampling. In Proceedings of the 52 nd Annual IEEE Symposium on Foundations of Computer Science (FOCS) , pages 363–372, 2011.
- 7[BHNT 15] Sayan Bhattacharya, Monika Henzinger, Danupon Nanongkai, and Charalampos E. Tsourakakis. Space- and time-efficient algorithm for maintaining dense subgraphs on one-pass dynamic streams. In Proceedings of the 47 th Annual ACM on Symposium on Theory of Computing (STOC) , pages 173–182, 2015.
- 8[BS 15] Marc Bury and Chris Schwiegelshohn. Sublinear estimation of weighted matchings in dynamic data streams. In Proceedings of the 23 rd Annual European Symposium on Algorithms (ESA) , pages 263–274, 2015.
