Mapping n Grid Points Onto a Square Forces an Arbitrarily Large Lipschitz Constant

We prove that the regular n×n\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n\times n}$$\end{document} square grid of points in the integer lattice Z2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbb{Z}^{2}}$$\end{document} cannot be recovered from an arbitrary n2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n^{2}}$$\end{document}-element subset of Z2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbb{Z}^{2}}$$\end{document} via a mapping with prescribed Lipschitz constant (independent of n). This answers negatively a question of Feige from 2002. Our resolution of Feige’s question takes place largely in a continuous setting and is based on some new results for Lipschitz mappings falling into two broad areas of interest, which we study independently. Firstly the present work contains a detailed investigation of Lipschitz regular mappings on Euclidean spaces, with emphasis on their bilipschitz decomposability in a sense comparable to that of the well known result of Jones. Secondly, we build on work of Burago and Kleiner and McMullen on non-realisable densities. We verify the existence, and further prevalence, of strongly non-realisable densities inside spaces of continuous functions.


Introduction
The main objective of this paper is to answer a question of Feige originating in the 1990s, which asks whether every subset of the integer lattice Z 2 with cardinality n 2 for some n ∈ N can be mapped via a bijection f onto the regular n × n grid {1, . . . , n} 2 in such a way that the Lipschitz constant of f may be bounded above independently of n:  [Mat02] and appeared in the technical report [MN11, Question 2.12], which lists many open problems arising from this meeting and, in several cases, their subsequent solutions. However, to the best of the authors' knowledge, there has been almost no V.K. was partially supported by the project GAČR 16-01602Y. He was also supported by the project FWF P23628-N18 during his stays at Universität Innsbruck. E.K. was partially supported by the project FWF P23628-N18.  [Kal14]; some ideas contained there were helpful in the development of the present work. Feige's motivation to ask the question stemmed from his work on the so-called 'graph bandwidth' problem. In this problem the goal is to find a bijection l : V → {1, . . . , n} for a given n-vertex graph (V, E) with the 'bandwidth' as small as possible, that is, minimising the quantity max uv∈E |l(u) − l(v)|. It is known that finding the optimal solution to this problem is NP-hard [Pap76]. In [Fei00] Feige designated a randomised approximation algorithm that produces a solution with bandwidth larger than the optimum by a factor polylogarithmic in n. His algorithm can also be adapted to a generalisation of the bandwidth problem to two dimensions; however, in this case it does not produce a bijection between the set of vertices V and the grid {1, . . . , √ n} 2 , but rather an injection into a larger grid leaving some of the grid points unused [Fei17]. Feige then asked whether one can map such a set bijectively onto the grid {1, . . . , √ n} 2 without increasing the bandwidth 'too much'.
We will prove that the answer to Feige's question is negative in all dimensions d ≥ 2: Moreover, we show how almost any positive continuous function on the unit cube [0, 1] d can be used to construct a sequence of sets S n ∈ F n verifying Theorem 1.2. It is then natural to ask how fast C n grows. Whilst it is straightforward to verify that C n ≤ √ dn, the results of the present article do not permit finer estimates of C n , either from above or below. The question of the rate of growth of the sequence (C n ) will be an interesting topic for future research.
In order to answer Feige's question, we adapt a technique developed independently by Burago and Kleiner [BK98] and McMullen [McM98], which translates questions about Lipschitz mappings on discrete sets into a continuous setting. The papers [BK98] and [McM98] present a negative answer to the question of whether every two separated nets in the plane are bilipschitz equivalent, or put differently, whether every separated net M ⊆ Z 2 admits a bilipschitz bijection f : M → Z 2 . This 'discrete' question is shown to be equivalent to the following 'continuous' one: Does there exist for every measurable function ρ : [0, 1] 2 → (0, ∞) with 0 < inf ρ < sup ρ < ∞ a bilipschitz mapping f : [0, 1] 2 → R 2 with ρ = Jac(f ) almost everywhere.
(1.1) 'pieces' whose images overlap in T . Therefore any local 'bad behaviour' of ρ on one piece of f −1 (T ) may be compensated for by its values on other pieces, allowing f to satisfy equation (1.2). So excluding Lipschitz regular solutions of (1.2) crucially requires global control of the density ρ.
We overcome these problems, in part, by using a bilipschitz decomposition result for Lipschitz regular mappings, which derives from a result of Bonk  This result turns out to be very useful. Given a Lipschitz regular mapping f : [0, 1] d → R d , it allows us to carefully choose the set T ⊆ f ([0, 1] d ), referred to in the discussion above so that its pre-image f −1 (T ) decomposes precisely as a finite union of open sets on which f is bilipschitz; see Proposition 2.15. We then extend the existing techniques for constructing bilipschitz non-realisable densities to fit this more complex situation, where there are multiple bilipschitz mappings in play instead of just one.
1.1 Structure of the paper. In Section 2 we investigate Lipschitz regular mappings on Euclidean spaces, in particular proving Theorem 2.10. Besides the material necessary for the proof of Theorem 2.10, and subsequently for the proof of Proposition 2.15, which we use to answer Question 1.1, Section 2 also contains a discussion of optimality and limits of Theorem 2.10.
The groundwork for our verification of non-realisable densities is laid in Section 3, where we develop a construction of Burago and Kleiner in order to derive certain powerful properties of bilipschitz Jacobians; see Lemma 3.1. This section is the most technical part of the present paper; readers willing to accept Lemma 3.1 may safely skip this Section during the first pass-over, and perhaps, return to it only later. Both Sections 2 and 3 can be read independently of the rest of the paper.
Section 4 is devoted to the proof of existence of non-realisable densities and Theorem 4.1. The only material from the previous sections needed in Section 4 is Proposition 2.15 and Lemma 3.1. Having gathered together all the necessary ingredients, we resolve Feige's question in Section 5.
In an effort to make the paper easier to read and to bring the more important ideas of the arguments to the forefront, we postpone the proofs of technical lemmas, or formal verifications of intuitive statements until an appendix. Many of the results in the appendix can be treated as exercises and a reader interested in the core argument can safely skip them. For the sake of completeness we include all proofs. The appendix is divided into four subsections, each corresponding to a section of the paper. Where specific notation is introduced within a given section, we adopt the same notation in its appendix. 1.2 Notation. We conclude this introduction with a summary of the notation and key definitions common to all sections of the paper: Sets and measures. For i ∈ N we write [i] for the set {1, 2, . . . , i}. Given a finite set F , we let |F | denote the cardinality of F . For an infinite set F we put |F | equal to ∞. Throughout the paper I will denote the unit interval I := [0, 1]. The closure, interior and boundary of a set A are written as A, int A and ∂A respectively. We adopt the convention of Mattila [Mat99] and do not distinguish between outer measures and measures. The symbol L will refer to the Lebesgue measure and we write a.e. instead of 'almost everywhere' or 'almost every' with respect to L. We let denote the average value of a measurable, real-valued function ρ on a measurable set S of positive Lebesgue measure. We will use the term density to refer to a nonnegative, measurable real valued function. Thus, each density ρ : I d → [0, ∞) can be associated to a measure ρL on I d as defined earlier in this introduction.
Norms and balls. We write − 2 for the Euclidean norm and − ∞ for both the supremum norm and (briefly) the L ∞ norm. In addition, the symbol T op will represent the operator norm of a linear mapping T . An open ball with centre a and radius r will be denoted by B(a, r). Most of the time balls will be in Euclidean spaces, but we sometimes consider balls in general metric spaces including spaces of functions. It will be clear from the context which norm or metric is relevant. Occasionally we will extend this notation to denote neighbourhoods of sets: Given a set A ⊆ R d and r > 0 we let B(A, r) := a∈A B(a, r). To denote the closure of a ball or set neighbourhood we write B instead of B.
Mappings. For L ≥ 1, we call a mapping f : then we say that f is (b, L)-bilipschitz. Moreover, we use the term L-bilipschitz to refer to the special case where b = 1/L. Occasionally we will consider Lipschitz and bilipschitz mappings between more general metric spaces and generalise the above notions in the standard way. All bilipschitz mappings between metric spaces are Lipschitz regular. A clasic example of a non-bilipschitz (in fact non-injective) Lipschitz regular mapping is given by a folding mapping of the plane R 2 , i.e. take the plane and fold it along the y-axis. This defines a mapping R 2 → R 2 which is Lipschitz regular with regularity constant 2.
David and Semmes studied Lipschitz regular mappings in the context of general metric spaces and Euclidean spaces. In the Euclidean space setting, David proves that Lipschitz regular mappings behave somewhat like bilipschitz mappings. More precisely, that inside any ball B in the domain of a Lipschitz regular map- We provide answers to these questions in the case that the dimension of the domain is equal to the dimension of the co-domain.
The main result of the present section can be derived quickly from a result of Bonk and Kleiner [BK02,Theorem 3.4]. Indeed an answer to the first part of the above question is implicitly present in [BK02]. Readers primarily interested in the resolution of Feige's question may therefore choose to take a short cut, beginning GAFA MAPPING N GRID POINTS ONTO A SQUARE 595 after the introduction of necessary notation and background on Lipschitz regular mappings; see the paragraph under the heading 'Short cut' on page 8. Our argument for the proof of Theorem 2.10 appears to be new. Moreover, since [BK02, Theorem 3.4] is a statement for more general mappings than we consider here, many of the difficulties which occur in the proof in [BK02] are not present in our setting. Our argument is shorter and simpler, and thus, we present it in a full detail. After the preparation of the necessary background, our argument for the proof of Theorem 2.10 occupies roughly one page. We discuss the result [BK02, Theorem 3.4] of Bonk and Kleiner and its proof further later on.
The material of the present section remaining after Proposition 2.15, which appears under the heading 'Optimality of Theorem 2.10' on page 15, can be seen as a complement of Theorem 2.10 and Proposition 2.15 and is independent of Feige's question and of [BK02].
Notation and convention. We use the term (C, L)-regular mapping to denote a Lipschitz regular mapping f with Lip(f ) ≤ L and Reg(f ) ≤ C. Let R > 0, we call a set S in a metric space (M, dist) R-separated if for every two distinct points x, y ∈ S we have dist(x, y) > R. We write H d for the d-dimensional Hausdorff measure. A ball is always assumed to be open if not said otherwise.
Before we start the exposition of the results we list general properties of Lipschitz regular mappings that will be needed later.
The upper bound in the previous lemma comes from the Lipschitz property, while the lower bound can be derived easily using Lipschitz regularity and the definition of Hausdorff measure. As a corollary, Lipschitz regular mappings possess Luzin's properties (N ) and (N −1 ) given in the next Corollary: A converse of Lemma 2.2 holds in the setting of Ahlfors regular spaces. We do not introduce the definition of Ahlfors regularity here, 2 since we will work only in the setting of Euclidean spaces, which are also Ahlfors regular.
for every x ∈ M and every r > 0, then f is Lipschitz regular.
Since the above lemma plays a key part in our resolution of Feige's question, we give a proof for the required case M = M = R d .
Proof. Without loss of generality we assume Lip(f ) = 1. Fix x ∈ R d , r > 0 and let Γ be a maximal 2r-separated subset of f −1 (B(x, r)). Then the balls (B(y, r)) y∈Γ are pairwise disjoint and Hence, up to multiplication by a fixed constant depending only on d, we have We deduce that f −1 (B(x, r)) can be covered by C · 2 d balls of radius 2r. Thus, f is Lipschitz regular with Reg(f ) ≤ max 2, C · 2 d .
We also add an easy observation, which, however, will prove useful later.
Observation 2.5. Let f : M → M be Lipschitz regular and y ∈ M . Then we have Proof. To the contrary, we assume there are pairwise distinct points x 1 , . . . , x Reg(f )+1 ∈ f −1 ({y}). Let us denote by r the minimum distance between x i and x j for 1 ≤ i < j ≤ Reg(f ) + 1. Then no ball in M of radius r 2 can contain more than one of the points x 1 , . . . , x Reg(f )+1 . Therefore, f −1 (B M (y, r 2 Reg(f ) )) cannot be covered with at most Reg(f ) balls of radius r 2 in M , a contradiction. Occasionally, we will be given an open set U and a Lipschitz regular mapping f defined on U , the closure of U . Then by N (f ) we mean the set N (f | U ) ⊆ U . Note that it is then still true that f (N (f )) has full measure in f (U )-we will use this fact several times.
Topological degree. An important tool that we use in our work, besides differentiability, is the notion of topological degree. We briefly introduce it here; for a detailed treatment of this topic, we refer to [ The degree function is defined explicitly in [Dei85, Chapter 2]. We just point out that in the special case where g ∈ C 1 (U, R d ) and for every point x ∈ g −1 ({y}) the derivative Dg(x) is invertible, then the degree function is given by the expression sign(Jac(g)(x)), (2.2) G A F A (see [Dei85, Definition 2.1]). In particular, we have that deg(g, U, y) = 0 whenever y ∈ R d \g(U ). We will require some further properties of the degree which follow easily from the properties (d1), (d2) and (d3). All of the statements of the next Proposition are contained in [Dei85, Theorem 3.1].
Proposition 2.6. Let U ⊆ R d be an open, bounded set, f ∈ C(U, R d ) and y ∈ R d \f (∂U ).
(i) If y and y belong to the same connected component of In the next Proposition, we extend the formula (2.2) to Lipschitz mappings. sign(Jac(f )(x)).
Since the proof of Proposition 2.7 is a rather technical exercise, we include it only in "Appendix A".
Before we present our version of the proof of Theorem 2.10, let us state one additional auxiliary lemma. It says that whenever a continuous mapping in R d has derivative of full rank at a point, it preserves neighbourhoods of this point. We believe that such a statement may be a folklore; however, we did not find any reference.
Bilipschitz decomposition of Lipschitz regular mappings. Our main goal in this section is to show that Lipschitz regular mappings in Euclidean spaces decompose into bilipschitz mappings in a nice way: Before we prove Theorem 2.10, let us put it briefly into context. For a general Lipschitz mapping h : R d → R d it is known that one can obtain a different bilipschitz decomposition using Sard's theorem; see e.g. [Fed96, Lemma 3.2.2]. One can start with sets defined for every k ∈ N and then cut these sets into pieces of diameter less than 1/k forming a decomposition (A n ) ∞ n=1 . Then Sard's theorem implies that L h R d \ n∈N A n = 0. When compared to the decomposition of Theorem 2.10, the difference is that the sets A n are not necessarily open, the lower bilipschitz constant of each h| An may depend on n and ∞ n=1 A n need not be a large subset of the domain in any sense.
If the decomposition that was just described is applied to a Lipschitz regular mapping, the resulting sets A n occupy almost all of the domain, since the set N (f ) has a full measure in the domain. But the sets A n still need not be open. The fact that for Lipschitz regular mappings it is possible to ensure the openness of bilipschitz pieces A n will be of crucial importance to us.
The first quantitative version of the decomposition using Sard's theorem was provided by David [Dav88, Proposition 1] for general Lipschitz mappings f : ) is large in measure, then B contains a set E large in measure such that f | E is bilipschitz. When applied to a Lipschitz regular mapping f , using the measure-preserving property expressed in Lemma 2.2, the condition that L(f (B)) is large in measure is satisfied automatically; for this version of David's result, see [DS00, Theorem 4.1].
A well-known result of Jones [Jon88] provides another quantitative version of the decomposition for Lipschitz mappings I d → R m . In the decomposition of Jones as well as that of David the bilipschitz pieces may have empty interior.

Question. Can we hope for any control of the measure of the bilipschitz pieces in a bilipschitz decomposition of Lipschitz regular mappings if one requires the pieces being open?
For example, can we hope for any control of the measure of the set ∞ n=1 A n given by the conclusion of Theorem 2.10?
The answer to the previous question is no: The decomposition from Theorem 2.10 cannot be strengthened in this way for a general Lipschitz regular mapping. A detailed discussion of these questions is contained in the subsection 'Optimality of Theorem 2.10' at the end of this section.
Our proof of Theorem 2.10 can be divided into three parts. The first one is to find, for any given open set in the domain, an open subset on which the given mapping is almost injective (this notion is formalised below). The second part is to show that a Lipschitz regular, almost injective map on an open set is injective and the third part is to prove that a Lipschitz regular, injective map on an open set with a convex image is bilipschitz. In each of these steps we rely on the Lipschitz regularity of the mapping in question.
Let us remark that the first two steps described above, which comprise of Lemmas 2.12 and 2.13 in the following, may be replaced by an application of [BK02, Theorem 3.4]. Bonk and Kleiner work in [BK02] with much more general mappings; instead of assuming that f : U ⊆ R d → R d is Lipschitz regular, they only require that f is continuous and that there is some constant C > 0 such that f −1 ({y}) ≤ C for all y ∈ R d . The latter condition is referred to as 'bounded multiplicity'. Moreover, the domain U may be replaced by any compact metric space X with the property that every non-empty open subset of X has topological dimension d.
The argument we present below is different to that of Bonk and Kleiner in [BK02]. However, a key aspect of both proofs appears to be finding points x in the domain such that f (x) is an interior point of the image f (O) for every neighbourhood O of x. The most difficult part of Bonk and Kleiner's argument is to show that such points exist. However, for Lipschitz regular mappings we can easily find many such points using almost everywhere differentiability of Lipschitz mappings, the regularity condition and Lemma 2.8. Indeed, note that all points in the set N (f ) have this property. Therefore, our argument below may be a more accessible approach to [BK02, Theorem 3.4] for the special case where the mappings considered are Lipschitz regular.
We start with the following definition.
Definition 2.11. We say that a mapping h : As advertised above, we begin by showing that a Lipschitz regular mapping is almost injective on some open set:

non-empty and open and f : U → R d be Lipschitz regular. Then there is a non-empty open set V ⊆ U such that f | V is almost injective and f (V ) is an open ball.
Proof. The proof relies heavily on the special properties of the set N (f ) (see (2.1) on page 597).
As the next step, we use the degree to show that whenever a Lipschitz regular mapping is almost injective on an open set U , it is injective on U .
Proof. Suppose not. This means we can find two points . Any such point has at least two preimages in B 1 by Proposition 2.7; again, this is a contradiction.
The third step towards the proof of Theorem 2.10 is to show that a Lipschitz regular, injective mapping with a convex image is bilipschitz.
Lemma 2.14. Let U ⊆ R d be an open set and f : U → R d be an injective, Lipschitz regular mapping such that f (U ) is convex. Then f is bilipschitz with lower bilipschitz constant at least We note that the same statement also appears in [BK02, Lemma 4.2]. For reader's convenience, we include its short proof here as well.
Proof. By Brouwer's invariance of domain [Hat02, Thm. 2B.3] the mapping f is a homeomorphism onto its image.
For every two distinct points x, y ∈ U we consider the line segment f (x)f (y) ⊂ f (U ) connecting their images. Its preimage under f , we denote it by γ(x, y) := f −1 f (x)f (y) , is a curve with endpoints x and y. By Lipschitz regularity, the curve γ(x, y) can be covered by at most Reg(f ) balls of radius Reg Consequently, the distance between x and y cannot be larger than 2 Proof of Theorem 2.10. We start with a countable basis (U n ) n∈N for the subspace topology on U . By a consecutive application of Lemmas 2.12, 2.13 and 2.14 we get a collection of open sets (V n ) n∈N such that for every n ∈ N we have V n ⊆ U n and that f | Vn is bilipschitz with lower bilipschitz constant b = 1/ 2 Reg(f ) 2 . Now we set A 1 := V 1 and inductively define A n := V n \ n−1 j=1 A j . By construction, the set ∞ n=1 A n is dense in U , and hence, also in U .
Using Theorem 2.10 we can deduce that a Lipschitz regular mapping on an open set can be expressed, on some open subset of the image, as a sum of bilipschitz homeomorphisms. Such form of a decomposition is needed for the non-realisability results contained in the next section and, ultimately, for the resolution of Feige's question 1.1.
n=1 be the open sets from the conclusion of Theorem 2.10 applied to the mapping f . Let y ∈ f (U ) be such that the number Assuming that this claim is valid we may define the desired sets . Thus the proof can be completed by verifying the earlier claim, that is, by proving . However, this contradicts the fact that f is bilipschitz on A βi , and therefore also bilipschitz on A βi . We conclude that 2.1 Optimality of Theorem 2.10. The remainder of the current section is devoted to discussion of limits and optimality of Theorem 2.10. The content here is independent of the rest of the article, so the reader interested mainly in the resolution of Feige's question 1.1 can safely skip the rest of this section.
In the above, we have raised a question of optimality of Theorem 2.10 in terms of the measure of the bilipschitz pieces A n . Theorem 2.10 does not offer any control of their measure; below we will show that this is unavoidable. However, in a special case that a Lipschitz regular mapping f has Reg(f ) ≤ 2, we can provide a stronger bilipschitz decomposition; namely, the bilipschitz pieces A n , in addition to the conclusions of Theorem 2.10, can cover almost all of the domain.
Proof. From Observation 2.5 we know that every point y ∈ f (U ) has either one or two preimages. Since L(∂U ) = 0, the set f (N (f ))\f (∂U ) has full measure in f (U ) by Luzin's property (N ). Let y ∈ f (N (f ))\f (∂U ). Using Lemma 2.8, we may choose r > 0 sufficiently small so that B(y, r) ⊆ f (U )\f (∂U ).
If deg(f, U, y) ≡ 1 (mod 2), then Proposition 2.7 implies that y has exactly one preimage. Using Proposition 2.6, part (i), we deduce that the same is true of all points y ∈ f (N (f )) ∩ B(y, r). Thus the mapping f : f −1 (B(y, r)) → B(y, r) is almost injective. We may now apply Lemma 2.13 and then Lemma 2.14 to conclude that f | f −1 (B(y,r)) is bilipschitz with lower bilipschitz constant On the other hand, if deg(f, U, y) ≡ 0 (mod 2) then y must have two distinct pre-images x 1 , x 2 ∈ N (f ). Let B 1 , B 2 be disjoint balls with x 1 ∈ B 1 and x 2 ∈ B 2 . From Lemma 2.8 we deduce that f (B 1 ) ∩ f (B 2 ) contains a non-empty open ball G containing the point y. Then every point in G has exactly one pre-image in each of the balls B 1 and B 2 . Hence f | f −1 (G)∩Bi is injective for i = 1, 2 and, applying Lemma 2.14, we conclude that these mappings are also bilipschitz with lower bilipschitz constant 1 2 Reg(f ) 2 .
In the above we established that for every point y ∈ f (N (f ))\f (∂U ) there is an open ball B containing y such that f −1 (B) decomposes precisely as the union of at most two sets on which f is bilipschitz with lower bilipschitz constant For every x ∈ X and every δ > 0 the mapping f is not injective on the ball B(x, δ). Moreover, there are disjoint, non-empty, open balls U 1 , We give a proof for the case d = 1. The example for d ≥ 1 can easily be constructed from this: If f : I → R is the example for the case d = 1 with an appropriate choice of ε, then the function h : verifies Example 2.17 for general d ≥ 1. The formal proof of this is left to the reader. Given a point a ∈ (0, 1) and c > 0 we will denote by F a,c the interval [a, a + 3c]. Next, we define a 1-Lipschitz function g(a, c) : I → I that makes two folds on F a,c in a sense; see Figure 1.
More precisely, we let We will now summarise various properties of the function g(a, c) which will be needed in the following construction. It is clear that g(a, c) is 1-Lipschitz and Denoting by J 1 , J 2 the two components of the set I\F a,c we further point out that the sets g(a, c)(J 1 ), g(a, c)(J 2 ) and g(a, c)(F a,c ) are pairwise disjoint subsets of I, and that g restricted to each J i is a translation. Therefore, for any interval U ⊆ g(a, c)(I), the preimage g(a, c) −1 (U ) is an isometric copy of U whenever U does not intersect g(a, c)(F a,c ), and g(a, c) −1 (U ) may be covered by 3 intervals of length L(U ) whenever U intersects g(a, c)(F a,c ).
Let X ⊆ I be a fat Cantor set 4 with L(X) ≥ 1 − ε and (A n ) ∞ n=1 be an enumeration of the components of I\X. In what follows we will use the fact that every neighbourhood of a given point x ∈ X contains some of the intervals (A n ) ∞ n=1 . The idea of the construction is to 'pleat' inside each of the intervals A n using mappings of the form g(a, c) defined above; see Figure 2. Now we describe the construction more formally. We start with f 0 := id. Let a n be a midpoint of the interval A n . For n ∈ N we write g n := g(f n−1 (a n ), c n ) and f n := g n • f n−1 , where c n > 0 are chosen small enough with respect to several constraints, which will be described during the course of the construction. Then we define f as the limit of f n .
The first requirement on c n is that F an,cn ⊂ A n . Second, in order for f to be well-defined, we want to choose c n so that the sequence (f n ) ∞ n=1 is Cauchy. We have already observed that g n − id ∞ ≤ 2c n . Thus, choosing c n ≤ ε 2 n+1 , we get that f n − f n−1 ∞ ≤ 2c n ≤ ε 2 n and that f is well-defined. Moreover, f clearly satisfies condition (i).
To see that f is (3, 1)-regular, note that f is 1-Lipschitz because f is the uniform limit of a sequence of 1-Lipschitz functions. However, obtaining the regularity estimate is a bit more tricky.
Given any open interval U ⊆ f (I) we have Imagine, for the time being, that the latter set can be covered by 3 intervals of length D(L(U ) + 2 f n − f ∞ ) for some D < 3. Then letting n → ∞ we deduce that f −1 (U ) can be covered by 3 closed intervals of length DL(U ), which in turn can be covered by 3 open intervals of length 3L(U ). Therefore, we fix a strictly increasing sequence of numbers D n ∈ [1, 3) such that D := sup D n < 3 and show that c n can be chosen so that the following holds true: For every n ∈ N and every open interval U ⊂ I the set f −1 n (U ) can be covered by 3 intervals of length D n L(U ).
For n = 0 the condition is clearly satisfied by any D 0 ≥ 1. For a general n ∈ N we will distinguish three cases. If U is disjoint from f n (F an,cn ), then f −1 is the preimage under f n−1 of a translation of U , which can be covered by 3 intervals of length D n−1 L(U ) by induction, which is less than D n L(U ).
If U intersects the interval f n (F an,cn ), but is disjoint from the set n−1 i=1 f n (F ai,ci ), then f −1 n (U ) is a translation of g −1 n (U ), which can be covered by 3 intervals of length L(U ), as was already noted in the discussion of the properties of g(a, c) above.
We are left with the option that U intersects f n (F an,cn ) and also the set . However, since the intervals in (A n ) ∞ n=1 are pairwise disjoint, this means that L(U ) must be quite large; namely, L(U ) ≥ L(An) 2 − 3c n , since a n is the midpoint of A n . On the other hand, the inequality f n − f n−1 ∞ ≤ 2c n implies that By induction, the latter set can be covered by 3 intervals of length D n−1 (L(U )+4c n ). The last quantity can be made smaller than D n L(U ) using the lower bound on L(U ) and choosing c n small enough. This finishes the proof that the function f is (3, 1)regular. By construction, the function f is not injective on any of the intervals F an,cn ⊂ A n , but it maps each of the three subintervals [a n + (j − 1)c n , a n + jc n ] for j ∈ [3] isometrically onto the same interval. Since every neighbourhood of any point of X contains some of the intervals (A n ) ∞ n=1 , this verifies condition (iii).
A natural question that can come to the reader's mind is where we can put Lipschitz regular mappings on the imaginary scale between bilipschitz and Lipschitz mappings? Are they closer to general Lipschitz mappings or rather to bilipschitz ones?
We can show that a typical 1-Lipschitz mapping, in the sense of the Baire Category Theorem, is not injective on any open subset of the domain, and hence, in the light of Theorem 2.10, a typical 1-Lipschitz mapping is not Lipschitz regular: Let D = B(u, r) ∈ B, g ∈ I(D) and η > 0. To verify that I(D) is nowhere dense we will find g ∈ V and s > 0 such that g − g ∞ < η and B(g , s) ∩ I(D) = ∅.
Choose ε < min {L(B(u, r/2), r/2, η} and let f : I d → I d and X ⊆ I d be given by Example 2.17.   Another question that a curious reader may ask is whether Lipschitz regular mappings can be characterised as Lipschitz mappings admitting a bilipschitz decomposition as in Theorem 2.10.
However, this turns out not to be the case. It is easy to construct an example with infinitely many overlapping images of bilipschitz pieces. But even more is true: It is possible to construct an injective 1-Lipschitz function on the unit interval that has a decomposition as in Theorem 2.10, but, at the same time, is not Lipschitz regular. An example f is given by the formula where g : I → I is any positive, bounded, measurable function which is constant and equal to one on a dense collection of open subintervals of I and not a.e. bounded away from zero.

Geometric Properties of Bilipschitz Mappings
Bilipschitz mappings of a Euclidean space R d transform volume according to the formula L(f (A)) = A |Jac(f )| dL. In this section we establish that bilipschitz mappings cannot transform volume too wildly. In some sense we show that sufficiently G A F A fine grids of cubes must witness 'continuity' of the volume transform. This in turn places rather restrictive conditions on the Jacobian of a bilipschitz mapping, which we will exploit in Section 4 in order to find non-realisable densities. Our work in this section is an interpretation of the construction of Burago and Kleiner [BK98], which we modify in various ways, leading to some extensions of the results in [BK98]. The main result of this section will be the following lemma: Statement 1 expresses that each collection of cubes S i+1 is much finer than the previous collection S i . The inequality of statement 2 can be interpreted geometrically as stating that the volume of the image of the cube S under h j is very close to the volume of the image of its neighbour S . Put differently, we may rewrite the inequality of 2 in the following form:

L(h j (S)) − L(h j (S )) ≤ ζL(S).
It is possible to assemble Lemma 3.1 using predominantly arguments contained in the article [BK98] of Burago  any version of Lemma 3.1 explicitly and to prove Lemma 3.1 it is not sufficient to just take some continuous part of their argument. One needs to inspect their whole proof in detail and work considerably to put together all of the pieces correctly. Therefore, we present a complete proof of Lemma 3.1 in which we introduce some new elements. The proof of Lemma 3.1 requires some preparation and will be given at the end of this section. Variants of the Burago-Kleiner construction with additional details have been employed in a pure discrete setting in the works [Gar09,Mag11] and [CN16].
Lying behind all of the results of the present section is a simple property of Lipschitz mappings of an interval: If [0, c] ⊆ R is an interval and a Lipschitz mapping h : [0, c] → R n stretches the endpoints 0, c almost as much as its Lipschitz constant allows, then it is intuitively clear that the mapping h is close to affine. The next dichotomy can be thought of as a 'discretised' version of this statement: Statement 2 is a discrete formulation of the condition that the Lipschitz constant of h is not almost realised by the endpoints 0, c. Statement 1 expresses in a discrete way that h is close to affine; after partitioning the interval [0, c] into N subintervals of equal length this statement asserts that h looks like an affine mapping on nearly all pairs of adjacent subintervals.
We now formulate a multi-dimensional version of Lemma 3.2; see Figure 3. We consider thin cuboids in We save the proof of both Lemma 3.2 and Lemma 3.3 until "Appendix B", since a slightly weaker version of Lemma 3.3, though not explicitly, is present in [BK98, Lemma 3.2]. Our proof of the one-dimensional statement Lemma 3.2 follows [BK98] closely, but we develop a new induction argument to deduce Lemma 3.3 from Lemma 3.2. In doing so we hope to expose clearly that the property of bilipschitz mappings established in Lemma 3.3 is of a one-dimensional nature. For now, let us demonstrate how we intend to apply Lemma 3.3. First, we show that whenever statement 1 of Lemma 3.3 holds for a bilipschitz mapping h into R d , there are adjacent cubes S i and S i+1 whose images under h have almost the same measure. Eventually this will lead to conclusion 2 of Lemma 3.1.
Lemma 3.4. Let L ≥ 1, ε ∈ (0, 1/2L), d ∈ N and N 0 = N 0 (d, L, ε) be given by the conclusion of Lemma 3.3. Let N ≥ N 0 Let the mappings f 1 : 1). These conditions imply a bound on the difference in volume of the images f 1 (S i ) and f 2 (S i ), namely For a verification see Lemma B.2 in the appendix. Since φ is a translation, this establishes the required inequality.
Given a bilipschitz mapping g : [0, c] × [0, c/N] d−1 → R n , we now seek to repetitively apply Lemma 3.3 on smaller and smaller scales in order to, in some sense, eliminate statement 2 of the dichotomy of Lemma 3.3. Consequently, we find cubes (scaled and translated copies of the sets [(i − 1)c/N, ic/N ] × [0, c/N] d−1 ) on which g satisfies inequality (3.1) of statement 1 of Lemma 3.3. This will allow us to apply Lemma 3.4.

Sketch of the elimination of statement 2 from Lemma 3.3.
Let all parameters d, L, ε, M , ϕ, N 0 , c, n and N be given by the statement of Lemma 3.3. We consider an L-bilipschitz mapping g : [0, c] × [0, c/N] d−1 → R n . If statement 2 holds for g, there is a pair of points a 1 := z, b 1 := z+ c NM e 1 which the mapping g stretches by a factor (1 + ϕ) more than it stretches the pair a 0 := 0 and b 0 := ce 1 . We may now consider the restriction of g to a rescaled copy of the original cuboid [0, c] × [0, c/N] d−1 with vertices a 1 and b 1 corresponding to 0 and ce 1 respectively. If, again, it is the case that statement 2 is valid for this mapping, then we find points a 2 , b 2 inside the new Figure 4: An illustration of the strategy to eliminate statement 2 from Lemma 3.3. The stretch factor of g on the points a 1 and b 1 is at least (1 + ϕ)-times larger then the stretch factor of g on the points 0 and ce 1 . If statement 2 applies in the next iteration, we find two points stretched by g with factor at least (1 + ϕ) 2 -times the stretch factor on 0 and ce 1 cuboid which g stretches by a factor (1 + ϕ) more than it stretches the pair a 1 and b 1 , and so a factor (1 + ϕ) 2 times more than it stretches a 0 and b 0 . The process is illustrated in Figure 4. We iterate this procedure as long as possible to obtain sequences (a i ) and (b i ) satisfying where the final bound is given by the lower bilipschitz inequality for g. It is clear now that the procedure described above cannot continue forever: Otherwise, for i sufficiently large, the inequality above contradicts the L-Lipschitz condition on g. Thus, Lemma 3.3 tells us that after at most r-iterations of the procedure, where r ∈ N is a number determined by d, L and ε, we must have that statement 1 is valid for the appropriate restriction of the mapping g.
Let us now present the conclusions of the above sketch formally. We postpone the formal proof of this statement until the appendix.
such that statement 1 of Lemma 3.3 is valid for the mapping g p : [0, We are now ready to give a proof of Lemma 3.1: Let us now verify that the above defined families S 1 , . . . , S r satisfy condition 1 in the statement of Lemma 3.1. It is immediate from Definition 3.6 that Thus, given 1 ≤ i < r and S ∈ S i , we have that Note that any cube in the collections S i+1 has the form for some w ∈ c i+1 Z d and l ∈ [N ]. Since S ∈ Q ci/N and c i /N = Mc i+1 , such a cube can only intersect S in a set of positive Lebesgue measure when w ∈ S. Therefore, the number of cubes in S i+1 that can intersect S ∈ Q ci/N in a set of positive Lebesgue measure is bounded above by It follows that where, in the above, we use c i+1 = c i /N M and L(S) = (c i /N ) d and prescribe that N is sufficiently large so that the inequality holds. Thus, statement 1 is satisfied. Turning now to statement 2, we consider a k-tuple (h 1 , . . . , h k ) of L-bilipschitz mappings h i : U → R d and define a mapping g : U → R kd co-ordinate-wise by . Therefore for each i ∈ Ω and each h t,p : [0, We fix i ∈ Ω and impose the condition ε < 1 2L on ε. Then the conditions of Lemma 3.4 are satisfied for It is clear upon reference to Definition 3.6 that S and S are e 1 -adjacent cubes belonging to the family S p . Moreover, we have h t (S) = h t,p (S i ) and h t (S ) = h t,p (S i+1 ) for all t ∈ [k]. Therefore S and S verify statement 2 of Lemma 3.1 for the k- tuple (h 1 , . . . , h k ). This completes the proof of Lemma 3.1.

Realisability in Spaces of Functions
The objective of the present section is to prove that in some sense almost all continuous functions ρ ∈ C(I d , R) do not admit a Lipschitz regular mapping f : where we view C(I d , R) as a Banach space with the supremum norm − ∞ . More precisely, we prove the following result: Then E is a σ-porous subset of C(I d , R).
Remark 4.2. To be able to work with functions ρ ∈ C(I d , R) attaining negative values as well, we extend the definition of the pushforward measure to such functions: where by ρ + , ρ − we mean the positive and the negative part of ρ. Technically speaking, the pushforward measure is no longer a measure, but a difference of two measures. 5 However, we will use it only in the form of (4.1), that is, when the result is again a measure. This is only a technical tool that helps us treat functions attaining negative values properly, but it does not bring in any additional difficulty to the present work. An alternative option would be to say that, by definition, no function with negative values satisfies (4.1), but the statement of Theorem 4.1 would be then seemingly weaker.
Burago and Kleiner [BK98] and McMullen [McM98] prove the existence of a positive function ρ : I 2 → R for which equation (4.1) has no bilipschitz solutions f : I 2 → R 2 . We point out that Theorem 4.1 extends this result in various ways. Firstly, Lipschitz regular mappings of I d into R d form a larger class than the class of bilipschitz mappings from I d to R d . Thus, Theorem 4.1 establishes the existence of a density ρ which admits no solutions f to equation (4.1) inside a larger class of mappings. Secondly, Theorem 4.1 asserts the existence of not only one such density ρ, but states that almost all continuous functions ρ ∈ C(I d , R) are not realisable in the sense of (4.1) for Lipschitz regular mappings f : I d → R d . That bilipschitz nonrealisable functions contain a dense G δ subset of both the set of positive continuous functions and the set of positive, L ∞ -bounded, measurable functions on the unit square [0, 1] 2 , was recently proved by Viera in [Vie16], but [Vie16] is completely independent from the present work. , 1] × I d−1 and write f : I d → D 1 for the mapping which 'folds D 2 onto D 1 '. More precisely, the mapping f is defined as the identity mapping on D 1 and as the reflection in the hyperplane 1 2 × R d−1 on D 2 . Let ψ ∈ C(D 1 , R) be a positive, bilipschitz non-realisable density with values in (0, 1). We impose the additional mild condition that ψ is constant with value 1 2 inside the hyperplane 1 2 × R d−1 . The existence of such a density ψ follows easily from the d-dimensional analog of [BK98, Theorem 1.2].
Set ρ = ψ on D 1 . The bilipschitz non-realisability of ρ is now already assured, no matter how we define ρ on D 2 . To make ρ Lipschitz regular realisable, we define ρ on D 2 by The function ρ : I d → R is continuous and positive, whilst the mapping f : I d → R d is Lipschitz regular and satisfies f (I d ) = D 1 . Moreover, for any measurable set where, for the penultimate equation, we use the change of variables formula and the fact that f restricted to the set D 2 is an affine isometry. This verifies the Lipschitz regular realisability of ρ.
Porous and σ-porous sets. We recall the definitions of porosity and σ-porosity according to [Zaj05, Definition 2.1], where they are referred to as 'lower porosity' and 'lower σ-porosity' respectively.  Figure 6: The diagram of a bilipschtz decomposition for a density ρ ∈ E C,L,n Definition 4.4. Let (X, − ) be a Banach space.
(i) A set P ⊆ X is called porous at a point x ∈ X if there exist ε 0 > 0 and α ∈ (0, 1) such that for every ε ∈ (0, ε 0 ) there exists y ∈ X such that

be expressed as a countable union of porous subsets of X.
The class of σ-porous subsets of a Banach space X is strictly contained in the class of subsets of X of the first category in the sense of the Baire Category Theorem. Further, the notions of porosity and σ-porosity extend to metric spaces in the natural way. For a survey on porous and σ-porous sets we refer the reader to [Zaj05]. Due to its relevance later in the paper, we point out that porosity of a set P ⊆ X is a weaker condition than requiring P to be porous at all points x ∈ X (not just at points x ∈ P ). For example, the set 1 n : n ∈ Z\ {0} is porous in R but is not porous at the point 0.
Let us begin working towards the proof of Theorem 4.1.
Porous decompositon of E. In the present paragraph we describe how to partition the set E into a countable family of porous sets (E C,L,n ). We will need the lower bilipschitz constant b(·) given by the conclusion of Proposition 2.15. Let (O n ) ∞ n=1 be a countable basis for the topology of I d . For C, L, n ∈ N we let E C,L,n denote the set of all functions ρ ∈ C( for a.e. y ∈ O n . (4.2) 618 M. DYMOND, V. KALUŽA, E. KOPECKÁ G A F A Note that the basis set O n 'generates' the diagram of bilipschitz homeomorphisms and V = f 1 (O n ); see Figure 6. However, the critical role of O n in the definition above is to prescribe the portion of the domain I d on which all functions ρ ∈ E C,L,n have the special form given by (4.2).
To explain the origins of equation (4.2), we refer the reader back to Proposition 2.15. Consider a Lipschitz regular mapping f : I d → R d and the non-empty open set T ⊆ f (I d ) given by the conclusion of Proposition 2.15. Because the preimage f −1 (T ) decomposes precisely as a union of N sets on which f defines a bilipschitz homeomorphism to T , a pushforward (signed) measure of the form f ρL with ρ ∈ C(I d , R) can be expressed on T as a sum of integrals involving ρ and Jacobians of N bilipschitz homeomorphisms f 1 , . . . , f N . Thus, whenever f ρL = L| f (I d ) we obtain some equation relating ρ to finitely many bilipschitz homeomorphisms and their Jacobians. We will see that this equation has precisely the form of (4.2).
For C, L ∈ N, let E C,L denote the subset of C(I d , R) consisting of all functions ρ for which there exists a (C, L)-regular mapping f : Clearly we have E = C,L∈N E C,L . In the next lemma we prove that E C,L is covered by the sets (E C,L,n ).
To see that these choices witness that ρ ∈ E C,L,n , it only remains to verify equation (4.2). Note that f −1 (V ) = N i=1 Y i . Therefore, for every measurable set S ⊆ V we have that We conclude that Recall that the sets Y i and V are all bilipschitz homeomorphic via the mappings f i : Y i → V . Therefore, we may make the substitution x = f 1 (y) in the above equation, after which a simple rearrangement and an application of a 'chain rule identity' for Jacobians yields (4.2).
If, for the time being, we treat the terms ρ(f −1 i • f 1 (y)) in (4.2) as constants, then, on the open set O n , functions ρ ∈ E C,L,n are linear combinations of at most C Jacobians of L/b(C)-bilipschitz mappings. The purpose of the next lemma is to provide, for given constants k and L, a function ψ ∈ C(I d , R) which is small in supremum norm, but far away from being a linear combination of k L-bilipschitz Jacobians.
Lemma 4.6. Let ε, ζ ∈ (0, 1), L ≥ 1, k ∈ N and U ⊆ I d be an open set. Then there exists a function ψ ∈ C(I d , R) such that ψ ∞ ≤ ε, supp(ψ) ⊆ U and for every k- tuple (h 1 , h 2 , . . . , h k ) To prove Lemma 4.6 it suffices to consider the families of tiled cubes S 1 , . . . , S r given by the conclusion of Lemma 3.1 applied to d, k, U , L, ζ and some very small η ∈ (0, 1), and to define ψ as a 'chessboard function' whose average value on e 1adjacent cubes makes jumps of size at least ε. From the conclusion 1 of Lemma 3.1 we may essentially regard the cubes from two different families S i , S j as pairwise disjoint; choosing η sufficiently small ensures that the values of ψ on S j have negligible impact on the average values of ψ on cubes in S i for i < j. We postpone the formal description of this construction until "Appendix C", since it is a standard argument without any deep ideas. For now, let us proceed to the key proof of the present section, namely the verification of porosity of the sets (E C,L,n ). We actually prove that the sets (E C,L,n ) possess a stronger property: Lemma 4.7. For every C, L, n ∈ N, E C,L,n is a porous subset of C(I d , R). In fact, the set E C,L,n is porous at every point φ ∈ C(I d , R).
Before we begin the proof, we will outline the strategy. For given C, L, n ∈ N, φ ∈ C(I d , R) and ε ∈ (0, 1), our task is to find a function φ ∈ C(I d , R) so that φ − φ ∞ ≤ ε and B( φ, αε) ∩ E C,L,n = ∅ for some α = α(C, L, n, φ). We will exploit the uniform continuity of φ: By prescribing at the start a sufficiently small open set U ⊆ O n we may treat φ as constant (relative to ε) on U and indeed on any L/b(C)-bilipschitz image of U . Thus, when using the condition (4.2) for functions ρ ∈ E C,L,n we will always be able to treat the terms ρ(f −1 i • f 1 (y)) as constant. In other words, on U we will have that all functions in E C,L,n are linear combinations of at most C L/b(C)-bilipschitz Jacobians. We set φ = φ + ψ where ψ is given by the conclusion of Lemma 4.6 for ζ = αε, L = L/b(C) and an appropriate choice of Proof of Lemma 4.7. Let C, L, n ∈ N, φ ∈ C(I d , R) and ε ∈ (0, 1). We will construct to be determined later in the proof.
Using that φ is uniformly continuous, we may choose δ > 0 sufficiently small so that and (4.4) holds for ψ. Using (4.2), we may now write for a.e. y ∈ O n . To complete the proof we will show that the average value of the final expression over the cube S is too close to its average value over S , that is, closer than the condition (4.4) on ψ allows.
Let i ∈ {2, 3, . . . , N}. Then we have that whenever y, z ∈ S ∪ S ⊆ U . Therefore, in light of (4.5) and the fact that for all y ∈ S ∪ S . Thus, we have (4.8) For the final inequality above we used ρ diam(S ∪ S ) ≤ diam(U ) < δ we deduce that the average value of the right hand side of (4.6) over the cube S differs from its average value over S by at most However, with the setting ζ = ε 2(4 + C( φ ∞ + 1 + 6(L/b(C)) d )) this number is strictly less than ε, contrary to (4.4). Thus, we conclude that B φ, ε 2(4 + C( φ ∞ + 1 + 6(L/b(C)) d )) ∩ E C,L,n = ∅, which demonstrates the porosity of E C,L,n at φ.
It is now a simple task to combine the previous Lemmas for a proof of Theorem 4.1.
Proof of Theorem 4.1. From Lemma 4.5 we have whilst Lemma 4.7 asserts that each of the sets in the union on the right hand side is porous.
Readers interested in the resolution of Feige's question may proceed immediately to Section 5. In the remainder of the present section we discuss an independent topic of interest.

Realisability in L ∞ spaces.
Until now we have only studied realisability in spaces of continuous functions. However, functions ρ admitting a bilipschitz or Lipschitz regular solution f : (4.9) need not be continuous. Therefore, it is natural to study the set of realisable functions in the less restrictive setting of L ∞ (I d ), the space of all Lebesgue measurable, realvalued functions ρ defined on I d , which are bounded with respect to the L ∞ -norm We will prove that the set of all bilipschitz realisable functions in L ∞ (I d ) is a σporous set. For bilipschitz mappings f , (4.9) is equivalent to the equation (4.10) The question of whether Lipschitz regular realisable densities are also σ-porous, or in some sense negligible, in L ∞ spaces remains open.
Theorem 4.8. Let Then G is a σ-porous subset of L ∞ (I d ). In fact G may be decomposed as a countable union of sets (G L ) ∞ L=1 so that each G L is porous at every point ρ ∈ L ∞ (I d ).
Remark. For 1 ≤ p < ∞, the question of whether the set of bilipschitz realisable densities is small in L p (I d ) is not interesting because the set of all a.e. bounded functions is already σ-porous in this space.
The proof of Theorem 4.8 will require the following lemma, for which we recall the notation of Section 3. We postpone the proof of the lemma until "Appendix C", since it is based on a slightly more delicate version of the construction employed in the proof of Lemma 4.6. Fix L ≥ 1, ρ ∈ L ∞ (I d ) and ε > 0. We will find ρ ∈ L ∞ (I d ) with ρ − ρ ∞ ≤ ε and B( ρ, ε/16) ∩ G L = ∅. This will verify the porosity of the set G L and complete the proof of the theorem.

GAFA MAPPING N GRID POINTS ONTO A SQUARE 623
Let U ⊆ I d be an arbitrary, non-empty, open set, ζ = ε/2 and let η ∈ (0, 1) be a parameter to be determined later in the proof. Let r ∈ N and the tiled families S 1 , . . . , S r of cubes contained in U be given by the conclusion of Lemma 3.1 applied to d, k = 1, L, η and ζ. We define a sequence of functions ( where ψ(S i , ρ, ε) is given by the conclusion of Lemma 4.9. Now let ρ ∈ L ∞ (I d ) be defined by when we set η = 1 16 . Together with Lemma 3.1, part 2 and the setting ζ = ε/2, this implies that equation (4.10) with ρ = φ has no L-bilipschitz solutions f : We provide a negative solution to Question 5.1 in a strong sense in all dimensions d ∈ N, d ≥ 2. For d = 1 the answer is trivially positive. First, we transform Feige's question into a question about densities of measures supported on I d . In order to do so, we adapt a construction of Burago and Kleiner [BK98] that encodes densities into separated nets in R d . Then we prove that if we plug in any of the positive non-realisable functions whose existence is ensured by Theorem 4.1, the sequence of discrete sets that arise from the chosen function provide a negative solution to Question 5.1.

Feige's Question
We provide an equivalent version of Question 5.1 that fits better the tools we have.
for every x, y ∈ X i . Thus, we need to examine the behaviour of the sequence From the definition of n i,k we have the following bounds on Using the identity − I d ρ = 1 and the fact that Jac(ϕ i ) = l d i we immediately obtain that Since mi li → 0, the sequence ni li i∈N is bounded and converges to 1, and hence, li ni → 1 as well. Consequently, for any L > L we can find i 0 ∈ N such that for every i ≥ i 0 the mappings g i are L -Lipschitz. By trimming off the initial segment of the sequence (g i ) i∈N up to i 0 , we can assume that all mappings g i are L -Lipschitz for any chosen L > L. We extend each g i , by Kirszbraun's extension theorem [Kir34], [Fed96,2.10.43], to a mappingḡ i : I d → R d such that Lip(ḡ i ) = Lip(g i ). By the Arzelà-Ascoli theorem we know that the sequence (ḡ i ) ∞ i=1 subconverges to a limit f , which is also L -Lipschitz. 8 By passing to a convergent subsequence, we may assume thatḡ i ⇒ f .
Before we prove that f transforms measure according to the equation we have to establish additional notation and present two lemmas about weak convergence of measures. 9 For i ≥ 1 we define a measure μ i on I d by In order to show that f (ρL) = L| f (I d ) we first prove that μ i converges weakly to ρL on I d . Moreover, this will be shown to imply that (ḡ i ) (μ i ) converges weakly to f (ρL). Finally, we prove that (ḡ i ) (μ i ) also converges weakly to L| f (I d ) , and hence, f (ρL) and L| f (I d ) must be the same by the uniqueness of weak limits.

GAFA MAPPING N GRID POINTS ONTO A SQUARE 629
To this end we use the following two lemmas, which are probably a part of a common knowledge in measure theory, but the authors were unable to find a proper reference. Since their proofs are straightforward, we include them only in "Appendix D".
Lemma 5.5. Let ν and (ν n ) ∞ n=1 be finite Borel measures on a compact metric space K. Moreover, assume that there is, for each n ∈ N, a finite collection Q n of Borel subsets of K that cover ν-almost all of K and, at the same time, Then ν n converges weakly to ν.
Lemma 5.6. Let K be a compact space and (ν n ) n∈N be a sequence of finite, Borel measures on K converging weakly to a finite Borel measure ν. Let X be a metric space and h n : K → X be a sequence of continuous mappings converging uniformly to h. Then (h n ) (ν n ) converges weakly to h (ν).
Equipped with the two lemmas above, we resume proving Theorem 5.3.
Claim 5.3.1. The sequence of measures (μ i ) i∈N converges weakly to ρL.
non-empty open balls and f : U 1 ∪ U 2 → R d be a continuous mapping such that f | Ui is a homeomorphism for i = 1, 2 and f (U 1 ) = f (U 2 ). Then there exists s > 0 such that for any continuous mapping h :

B Appendix to Section 3: Geometric Properties of Bilipschitz Mappings
Proof. Let M ∈ N and ϕ ∈ (0, 1) be parameters to be determined later in the proof. Let c > 0, n ∈ N, N ∈ N and h : [0, c] → R n be an L-bilipschitz mapping.
, t = t(L, ε) ∈ (ϕ, 1) be some parameter to be determined later in the proof and For x ∈ P we have This inequality follows from the definition of P , the inequality (B.2) and t > ϕ. For the remaining co-ordinate functions we have Combining the two inequalities above we deduce Rearranging this inequality, we obtain where, for the last inequality, we apply N ≥ 2. It follows that the set c NM Z ∩ [0, c − c N ]\P can intersect at most 6ϕ ϕ+t (N − 1) intervals S i . Letting we deduce that |Ω| ≥ (1 − 6ϕ ϕ+t )(N − 1). Moreover for any i ∈ Ω and x ∈ S i , we can find x ∈ P with |x − x| ≤ c/N M . This allows us to apply (B.3) to get We are now ready to specify the parameters t, M and ϕ, so that the inequalities obtained above verify statement 1. First, we prescribe that t ∈ (0, 1) is sufficiently small and M ∈ N is sufficiently large so that √ 5tL+ 2L M < ε. Finally we demand that ϕ ∈ (0, t) is small enough so that 6ϕ ϕ+t < ε. (3.1) Proof. In this proof we will sometimes add the superscript d or d − 1 to objects such as the Lebesgue measure L or vectors e i , 0 in order to emphasise the dimension of the Euclidean space to which they correspond. For d ≥ 2, we will express points in R d in the form x = (x 1 , x 2 , . . . , x d ). Given x = (x 1 , . . . , x d ) ∈ R d and s ∈ R we let  (x 1 , x 2 , . . . , x d−1 , s).
Thus, we get that for each s ∈ [0, c/N] at least one of the following statements holds: Therefore, by Fubini's theorem, For each i ∈ [N − 1] we let S i : and observe that Combining the two inequalities derived above for L d (A), we deduce Moreover, for any i ∈ Ω and any cube Q ⊆ S i with side length (2 √ θL d (S i )) 1 d we have A ∩ Q = ∅. Therefore, for any i ∈ Ω and any x ∈ S i we can find x ∈ A ∩ S i with Using this approximation, we obtain Thus, statement 1 is verified when we prescribe that θ > 0 is sufficiently small and N 0 is sufficiently large so that For the inequality we use 2Lε ∈ (0, 1) and apply Bernoulli's inequality. We conclude that L(f 1 (S)) − L(f 2 (S)) ≤ 2dL d+1 εL(S).
Since the above argument is completely symmetric with respect to f 1 and f 2 , we also have L(f 2 (S)) − L(f 1 (S)) ≤ 2dL d+1 εL(S).