Title: A CLOSER LOOK AT THE APPROXIMATION CAPABILI-TIES OF NEURAL NETWORKS

Abstract: The universal approximation theorem, in one of its most general versions, says that if we consider only continuous activation functions σ, then a standard feedforward neural network with one hidden layer is able to approximate any continuous multivariate function f to any given approximation threshold ε, if and only if σ is non-polynomial. In this paper, we give a direct algebraic proof of the theorem. Furthermore we shall explicitly quantify the number of hidden units required for approximation. Specifically, if X ⊆ R n is compact, then a neural network with n input units, m output units, and a single hidden layer with n+d d hidden units (independent of m and ε), can uniformly approximate any polynomial function f : X → R m whose total degree is at most d for each of its m coordinate functions. In the general case that f is any continuous function, we show there exists some N ∈ O(ε -n ) (independent of m), such that N hidden units would suffice to approximate f . We also show that this uniform approximation property (UAP) still holds even under seemingly strong conditions imposed on the weights. We highlight several consequences: (i) For any δ > 0, the UAP still holds if we restrict all non-bias weights w in the last layer to satisfy |w| < δ. (ii) There exists some λ > 0 (depending only on f and σ), such that the UAP still holds if we restrict all non-bias weights w in the first layer to satisfy |w| > λ. (iii) If the non-bias weights in the first layer are fixed and randomly chosen from a suitable range, then the UAP holds with probability 1.

1. INTRODUCTION AND OVERVIEW:

A standard (feedforward) neural network with n input units, m output units, and with one or more hidden layers, refers to a computational model N that can compute a certain class of functions ρ : R n → R m , where ρ = ρ W is parametrized by W (called the weights of N ). Implicitly, the definition of ρ depends on a choice of some fixed function σ : R → R, called the activation function of N . Typically, σ is assumed to be continuous, and historically, the earliest commonly used activation functions were sigmoidal.

A key fundamental result justifying the use of sigmoidal activation functions was due to Cybenko (1989) , Hornik et al. (1989) , and Funahashi (1989) , who independently proved the first version of what is now famously called the universal approximation theorem. This first version says that if σ is sigmoidal, then a standard neural network with one hidden layer would be able to uniformly approximate any continuous function f : X → R m whose domain X ⊆ R n is compact. Hornik (1991) extended the theorem to the case when σ is any continuous bounded non-constant activation function. Subsequently, Leshno et al. (1993) proved that for the class of continuous activation functions, a standard neural network with one hidden layer is able to uniformly approximate any continuous function f : X → R m on any compact X ⊆ R n , if and only if σ is non-polynomial.

Although a single hidden layer is sufficient for the uniform approximation property (UAP) to hold, the number of hidden units required could be arbitrarily large. Given a subclass F of real-valued continuous functions on a compact set X ⊆ R n , a fixed activation function σ, and some ε > 0, let N = N (F, σ, ε) be the minimum number of hidden units required for a single-hidden-layer neural network to be able to uniformly approximate every f ∈ F within an approximation error threshold of ε. If σ is the rectified linear unit (ReLU) x → max(0, x), then N is at least Ω( 1 √ ε ) when F is the class of C 2 non-linear functions (Yarotsky, 2017) , or the class of strongly convex differentiable functions (Liang & Srikant, 2016) ; see also (Arora et al., 2018) . If σ is any smooth non-polynomial function, then N is at most O(ε -n ) for the class of C 1 functions with bounded Sobolev norm (Mhaskar, 1996) ; cf. (Pinkus, 1999, Thm. 6.8) , (Maiorov & Pinkus, 1999) . As a key highlight of this paper, we show that if σ is an arbitrary continuous non-polynomial function, then N is at most O(ε -n ) for the entire class of continuous functions. In fact, we give an explicit upper bound for N in terms of ε and the modulus of continuity of f , so better bounds could be obtained for certain subclasses F, which we discuss further in Section 4. Furthermore, even for the wider class F of all continuous functions f : X → R m , the bound is still O(ε -n ), independent of m.

To prove this bound, we shall give a direct algebraic proof of the universal approximation theorem, in its general version as stated by Leshno et al. (1993) (i.e. σ is continuous and non-polynomial). An important advantage of our algebraic approach is that we are able to glean additional information on sufficient conditions that would imply the UAP. Another key highlight we have is that if F is the subclass of polynomial functions f : X → R m with total degree at most d for each coordinate function, then n+d d hidden units would suffice. In particular, notice that our bound N ≤ n+d d does not depend on the approximation error threshold ε or the output dimension m.

We shall also show that the UAP holds even under strong conditions on the weights. Given any δ > 0, we can always choose the non-bias weights in the last layer to have small magnitudes no larger than δ. Furthermore, we show that there exists some λ > 0 (depending only on σ and the function f to be approximated), such that the non-bias weights in the first layer can always be chosen to have magnitudes greater than λ. Even with these seemingly strong restrictions on the weights, we show that the UAP still holds. Thus, our main results can be collectively interpreted as a quantitative refinement of the universal approximation theorem, with extensions to restricted weight values.

Outline: Section 2 covers the preliminaries, including relevant details on arguments involving dense sets. Section 3 gives precise statements of our results, while Section 4 discusses the consequences of our results. Section 5 introduces our algebraic approach and includes most details of the proofs of our results; details omitted from Section 5 can be found in the appendix. Finally, Section 6 concludes our paper with further remarks.

2. PRELIMINARIES:



2.1. NOTATION AND DEFINITIONS:

Let N be the set of non-negative integers, let 0 n be the zero vector in R n , and let Mat(k, ) be the vector space of all k-bymatrices with real entries. For any function f : R n → R m , let f [t] denote the t-th coordinate function of f (for each 1 ≤ t ≤ m). Given α = (α 1 , . . . , α n ) ∈ N n and any n-tuple x = (x 1 , . . . , x n ), we write x α to mean

x α1 1 • • • x αn n . If x ∈ R n , then x α is a real number, while if x is a sequence of variables, then x α is a monomial, i.e. an n-variate polynomial with a single term. Let W n,m N := {W ∈ Mat(n + 1, N ) × Mat(N + 1, m)} for each N ≥ 1, and define W n,m = N ≥1 W n,m N . If the context is clear, we supress the superscripts n, m in W n,m N and W n,m .

Given any X ⊆ R n , let C(X) be the vector space of all continuous functions f : X → R. We use the convention that every f ∈ C(X) is a function f (x 1 , . . . , x n ) in terms of the variables x 1 , . . . , x n , unless n = 1, in which case f is in terms of a single variable x (or y). We say f is non-zero if f is not identically the zero function on X. Let P(X) be the subspace of all polynomial functions in C(X).

For each d ∈ N, let P ≤d (X) (resp. P d (X)) be the subspace consisting of all polynomial functions of total degree ≤ d (resp. exactly d). More generally, let C(X, R m ) be the vector space of all continuous functions f : X → R m , and define P(X, R m ), P ≤d (X, R m ), P d (X, R m ) analogously.

Throughout, we assume that σ ∈ C(R). For every W = (W (1) , W (2) ) ∈ W, let w (k) j be the j-th column vector of W (k) , and let w (k) i,j be the (i, j)-th entry of W (k) (for k = 1, 2). The index i begins at i = 0, while the indices j, k begin at j = 1, k = 1 respectively. For convenience, let w W is given by the map

x → w (2) 0,j + N i=1 w (2) i,j σ(w (1) i • (1, x)),

where "•" denotes dot product, and (1, x) denotes a column vector in R n+1 formed by concatenating 1 before x. The class of functions that neural networks N with one hidden layer can compute is precisely {ρ σ W : W ∈ W}, where σ is called the activation function of N (or of ρ σ W ). Functions ρ σ W satisfying W ∈ W N correspond to neural networks with N hidden units (in its single hidden layer). Every w (k) i,j is called a weight in the k-th layer, where w

[k] i,j is called a bias weight (resp. non-bias weight) if i = 0 (resp. i = 0).

Notice that we do not apply the activation function σ to the output layer. This is consistent with previous approximation results for neural networks. The reason is simple:

σ • ρ σ [j] W (restricted to domain X ⊆ R n ) cannot possibly approximate f : X → R if there exists some x 0 ∈ X such that σ(X) is bounded away from f (x 0 ). If instead f (X) is contained in the closure of σ(X), then applying σ to ρ σ [j]

W has essentially the same effect as allowing for bias weights w

(2) 0,j . Although some authors, e.g. (Leshno et al., 1993) , do not explicitly include bias weights in the output layer, the reader should check that if σ is not identically zero, say σ(y 0 ) = 0, then having a bias weight w

(2) 0,j = c is equivalent to setting w

(2) 0,j = 0 (i.e. no bias weight in the output layer) and introducing an (N + 1)-th hidden unit, with corresponding weights w (1) 0,N +1 = y 0 , w (1) i,N +1 = 0 for all 1 ≤ i ≤ n, and w (2) N +1,j = c σ(y0) ; this means our results also apply to neural networks without bias weights in the output layer (but with one additional hidden unit).

2.2. ARGUMENTS INVOLVING DENSE SUBSETS:

A key theme in this paper is the use of dense subsets of metric spaces. We shall consider several notions of "dense". First, recall that a metric on a set S is any function d : S × S → R such that for all x, y, z ∈ S, the following conditions hold:

(i) d(x, y) ≥ 0, with equality holding if and only if x = y;

(ii) d(x, y) = d(y, x); (iii) d(x, z) ≤ d(x, y) + d(y, z).

The set S, together with a metric on S, is called a metric space. For example, the usual Euclidean norm for vectors in R n gives the Euclidean metric (u, v) → u -v 2 , hence R n is a metric space.

In particular, every pair in W N can be identified with a vector in R (m+n+1)N , so W N , together with the Euclidean metric, is a metric space.

Given a metric space X (with metric d), and some subset U ⊆ X, we say that U is dense in X (w.r.t. d) if for all ε > 0 and all x ∈ X, there exists some u ∈ U such that d(x, u) < ε. Arbitrary unions of dense subsets are dense. If U ⊆ U ⊆ X and U is dense in X, then U must also be dense in X.

A basic result in algebraic geometry says that if p ∈ P(R n ) is non-zero, then {x ∈ R n : p(x) = 0} is a dense subset of R n (w.r.t. the Euclidean metric). This subset is in fact an open set in the Zariski topology, hence any finite intersection of such Zariski-dense open sets is dense; see (Eisenbud, 1995) . More generally, the following is true: Let p 1 , . . . , p k ∈ P(R n ), and suppose that X := {x ∈ R n : p i (x) = 0 for all 1 ≤ i ≤ k}. If p ∈ P(X) is non-zero, then {x ∈ X : p(x) = 0} is a dense subset of X (w.r.t. the Euclidean metric). In subsequent sections, we shall frequently use these facts.

Let X ⊆ R n be a compact set. (Recall that X is compact if it is bounded and contains all of its limit points.) For any real-valued function f whose domain contains X, the uniform norm of

f on X is f ∞,X := sup{|f (x)| : x ∈ X}. More generally, if f : X → R m , then we define f ∞,X := max{ f [j] ∞,X : 1 ≤ j ≤ m}.

The uniform norm of functions on X gives the uniform metric (f, g) → f -g ∞,X , hence C(X) is a metric space.

2.3. BACKGROUND ON APPROXIMATION THEORY:

Theorem 2.1 (Stone-Weirstrass theorem). Let X ⊆ R n be compact. For any f ∈ C(X), there exists a sequence {p k } k∈N of polynomial functions in P(X) such that lim k→∞ f -p k ∞,X = 0.

Let X ⊆ R be compact. For all d ∈ N and f ∈ C(X), define

E d (f ) := inf{ f -p ∞,X : p ∈ P ≤d (X)}.

(1)

A central result in approximation theory, due to Chebyshev, says that for fixed d, f , the infimum in (1) is attained by some unique p * ∈ P ≤d (R); see (Rivlin, 1981, Chap. 1) . (Notice here that we define p * to have domain R.) This unique polynomial p * is called the best polynomial approximant to f of degree d.

Given a metric space X with metric d, and any uniformly continuous function f : X → R, the modulus of continuity of f is a function

ω f : [0, ∞] → [0, ∞] defined by ω f (δ) := sup{|f (x) -f (y)| : x, y ∈ X, d(x, y) ≤ δ}.

By the Heine-Cantor theorem, any continuous f with a compact domain is uniformly continuous. Theorem 2.2 (Jackson's theorem; see (Rivlin, 1981, Cor. 1.4.1) ). Let d ≥ 1 be an integer, and let Y ⊆ R be a closed interval of length r ≥ 0. Suppose f ∈ C(Y ), and let p * be the best polynomial approximant to f of degree d.

Then f -p * ∞,Y = E d (f ) ≤ 6ω f ( r 2d ).

3. MAIN RESULTS:

Throughout this section, let X ⊆ R n be a compact set.

Theorem 3.1. Let d ≥ 2 be an integer, and let f ∈ P ≤d (X, R m ) (i.e. each coordinate function f [t] has total degree ≤ d). If σ ∈ C(R)\P ≤d-1 (R), then for every ε > 0, there exists some W ∈ W ( n+d d ) such that f -ρ σ W ∞,X < ε. Furthermore, the following holds: (i) Given any λ > 0, we can choose this W to satisfy the condition that |w

(2) i,j | < λ for all non-bias weights w

(2) i,j (i.e. i = 0) in the second layer.

(ii) There exists some λ > 0, depending only on f and σ, such that we could choose the weights of W in the first layer to satisfy the condition that w

(1) j 2 > λ for all j. Theorem 3.2. Let f ∈ C(X, R m ), and suppose σ ∈ C(R)\P(R). Then for every ε > 0, there exists an integer N ∈ O(ε -n ) (independent of m), and some W ∈ W N , such that f -ρ σ W ∞,X < ε. In particular, if we let D := sup{ x -y 2 : x, y ∈ X} be the diameter of X, then we can set N = n+dε dε , where

d ε := min{d ∈ Z : d ≥ 2, ω f [t] ( D 2d ) < ε 6 for all 1 ≤ t ≤ m}. (Note that d ε is well-defined, since lim δ→0 + ω f [t] (δ) = 0 for each t.

) Furthermore, we could choose this W to satisfy either (i) or (ii), where (i), (ii) are conditions on W as described in Theorem 3.1. Theorem 3.3. Let f ∈ C(X, R m ), and suppose that σ ∈ C(R)\P(R). Then there exists λ > 0 (which depends only on f and σ) such that for every ε > 0, there exists an integer N (independent of m) such that the following holds:

Let W ∈ W N such that each w

(1) j ∈ R n (for 1 ≤ j ≤ N ) is chosen uniformly at random from the set {u ∈ R n : u 2 > λ}. Then, with probability 1, there exist choices for the bias weights w

(1) 0,j (for 1 ≤ j ≤ N ) in the first layer, and (both bias and non-bias) weights w

(2) i,j in the second layer, such that f -ρ σ W ∞,X < ε. Moreover, N ∈ O(ε -n ) for general f ∈ C(X, R m ), and we can let N = n+d d if f ∈ P ≤d (X, R m ).

4. DISCUSSION:

The universal approximation theorem (version of Leshno et al. (1993) ) is an immediate consequence of Theorem 3.2 and the observation that σ must be non-polynomial for the UAP to hold, which follows from the fact that the uniform closure of P ≤d (X) is P ≤d (X) itself, for every integer d ≥ 1. Alternatively, we could infer the universal approximation theorem by applying the Stone-Weirstrass theorem (Theorem 2.1) to Theorem 3.1.

Given fixed n, m, d, a compact set X ⊆ R n , and σ ∈ C(R)\P ≤d-1 (R), Theorem 3.1 says that we could use a fixed number N of hidden units (independent of ε) and still be able to approximate any function f ∈ P ≤d (X, R m ) to any desired approximation error threshold ε. Our ε-free bound, although possibly surprising to some readers, is not the first instance of an ε-free bound: Neural networks with two hidden layers of sizes 2n + 1 and 4n + 3 respectively are able to uniformly approximate any f ∈ C(X), provided that we use a (somewhat pathological) activation function (Maiorov & Pinkus, 1999) ; cf. (Pinkus, 1999) . Lin et al. (2017) showed that for fixed n, d, and a fixed smooth non-linear σ, there is a fixed N (i.e. ε-free), such that a neural network with N hidden units is able to approximate any f ∈ P ≤d (X). An explicit expression for N is not given, but we were able to infer from their constructive proof that N = 4 n+d+1 d -4 hidden units are required, over d -1 hidden layers (for d ≥ 2). In comparison, we require less hidden units and a single hidden layer.

Our proof of Theorem 3.2 is an application of Jackson's theorem (Theorem 2.2) to Theorem 3.1, which gives an explicit bound in terms of the values of the modulus of continuity ω f of the function f to be approximated. The moduli of continuity of several classes of continuous functions have explicit characterizations. For example, given constants k > 0 and 0 < α ≤ 1, recall that a continuous function f :

X → R (for compact X ⊆ R n ) is called k-Lipschitz if |f (x) -f (y)| ≤ k x -y for all x, y ∈ X, and it is called α-Hölder if there is some constant c such that |f (x)-f (y)| ≤ c x-y α for all x, y ∈ X. The modulus of continuity of a k-Lipschitz (resp. α-Hölder) continuous function f is ω f (t) = kt (resp. ω f (t) = ct α ), hence Theorem 3.2 implies the following corollary. Corollary 4.1. Suppse σ ∈ C(R)\P(R). (i) If f : [0, 1] n → R is k-Lipschitz continuous, then for every ε > 0, there exists some W ∈ W N that satisfies f -ρ σ W ∞,X < ε, where N = n+ 3k ε n . (ii) If f : [0, 1] n → R is α-Hölder continuous, then there is a constant k such that for every ε > 0, there exists some W ∈ W N that satisfies f -ρ σ W ∞,X < ε, where N = n+d d , and d = 1 2 ( k ε ) 1/α .

An interesting consequence of Theorem 3.3 is the following: The freezing of lower layers of a neural network, even in the extreme case that all frozen layers are randomly initialized and the last layer is the only "non-frozen" layer, does not necessarily reduce the representability of the resulting model. Specifically, in the single-hidden-layer case, we have shown that if the non-bias weights in the first layer are fixed and randomly chosen from some suitable fixed range, then the UAP holds with probability 1, provided that there are sufficiently many hidden units. Of course, this representability does not reveal anything about the learnability of such a model. In practice, layers are already pre-trained before being frozen. It would be interesting to understand quantitatively the difference between having pre-trained frozen layers and having randomly initialized frozen layers.

Theorem 3.3 can be viewed as a result on random features, which were formally studied in relation to kernel methods (Rahimi & Recht, 2007) . In the case of ReLU activation functions, Sun et al.

(2019) proved an analog of Theorem 3.3 for the approximation of functions in a reproducing kernel Hilbert space; cf. (Rahimi & Recht, 2008) . For a good discussion on the role of random features in the representability of neural networks, see (Yehudai & Shamir, 2019) .

The UAP is also studied in other contexts, most notably in relation to the depth and width of neural networks. Lu et al. (2017) proved the UAP for neural networks with hidden layers of bounded width, under the assumption that ReLU is used as the activation function. Soon after, Hanin (2017) strengthened the bounded-width UAP result by considering the approximation of continuous convex functions. Recently, the role of depth in the expressive power of neural networks has gathered much interest (Delalleau & Bengio, 2011; Eldan & Shamir, 2016; Mhaskar et al., 2017; Montúfar et al., 2014; Telgarsky, 2016) . We do not address depth in this paper, but we believe it is possible that our results can be applied iteratively to deeper neural networks, perhaps in particular for the approximation of compositional functions; cf. (Mhaskar et al., 2017) .

5. AN ALGEBRAIC APPROACH FOR PROVING UAP:

We begin with a "warm-up" result. Subsequent results, even if they seem complicated, are actually multivariate extensions of this "warm-up" result, using very similar ideas. Theorem 5.1. Let p(x) be a real polynomial of degree d with all-non-zero coefficients, and let a 1 , . . . , a d+1 be real numbers. For each 1 ≤ j ≤ d + 1, define f j : R → R by x → p(a j x). Then f 1 , . . . , f d+1 are linearly independent if and only if a 1 , . . . , a d+1 are distinct.

Proof. For each 0 ≤ i, k ≤ d and each 1 ≤ j ≤ d + 1, let f (i) j (resp. p (i) ) be the i-th derivative of f j (resp. p), and let α (i) k be the coefficient of x k in p (i) (x). Recall that the Wronskian of (f 1 , . . . , f d+1 ) is defined to be the determinant of the matrix M (x) := [f (i-1) j (x)] 1≤i,j≤d+1 . Since f 1 , . . . , f d+1 are polynomial functions, it follows that (f 1 , . . . , f d+1 ) is a sequence of linearly independent functions if and only if its Wronskian is not the zero function (LeVeque, 1956, Thm. 4.7(a) ). Clearly, if a i = a j for i = j, then det M (x) is identically zero. Thus, it suffices to show that if a 1 , . . . , a d+1 are distinct, then the evaluation det M (1) of this Wronskian at x = 1 gives a non-zero value. Now, the (i, j)-th entry of M (1) equals a i-1 j p (i-1) (a j ), so M (1) = M M , where M is an upper triangular matrix whose (i, j)-th entry equals α (i-1) j-i , and M = [a i-1 j ] 1≤i,j≤d+1 is the transpose of a Vandermonde matrix, whose determinant is

det(M ) = 1≤i<j≤d+1 (a j -a i ). Note that the k-th diagonal entry of M is α (k-1) 0 = (k -1)!α (0)

k-1 , which is non-zero by assumption, so det(M ) = 0. Thus, if a 1 , . . . , a d+1 are distinct, then det

M (1) = det(M ) det(M ) = 0. Definition 5.2. Given N ≥ 1, W ∈ W n,m N , x 0 ∈ R n , and any function g : R → R, let F g,x0 (W ) denote the sequence of functions (f 1 , . . . , f N ), such that each f j : R n → R is defined by the map x → g( w (1) j • (x -x 0 )). Also, define the set g W ind n,N ;x0 := {W ∈ W n,m N : F g,x0 (W ) is linearly independent}.

Note that the value of m is irrelevant for defining g W ind n,N ;x0 . Remark 5.3. Given a = (a 1 , . . . , a n ) ∈ R n , consider the ring automorphism ϕ :

P(R n ) → P(R n ) induced by x i → x i -a i for all 1 ≤ i ≤ n. For any f 1 , . . . , f k ∈ P(R n ) and scalars α 1 , . . . , α k ∈ R, note that α 1 f 1 + • • • + α k f k = 0 if and only if α 1 ϕ(f 1 ) + • • • + α k ϕ(f k ) = 0, thus linear independence is preserved under ϕ. Consequently, if the function g in Definition 5.2 is polynomial, then g W ind n,N ;x0 = g W ind n,N ;0n for all x 0 ∈ R n . Corollary 5.4. Let m be arbitrary. If p ∈ P d (R) has all-non-zero coefficients, then p W ind 1,d+1;0 is a dense subset of W 1,m

d+1 (in the Euclidean metric).

Proof. For all 1 ≤ j < j ≤ N , let A j,j := {W ∈ W 1,m d+1 : w (1) 1,j -w (1) 1,j = 0}, and note that A j,j is dense in W 1,m d+1 . So by Theorem 5.1, p W ind 1,d+1;0 = 1≤j<j ≤N A j,j is also dense in W 1,m d+1 .

As we have seen in the proof of Theorem 5.1, Vandermonde matrices play an important role. To extend this theorem (and Corollary 5.4) to the multivariate case, we need a generalization of the Vandermonde matrix as described in (D' Andrea & Tabera, 2009) . (Other generalizations of the Vandermonde matrix exist in the literature.) First, define the sets

Λ n ≤d := {(α 1 , . . . , α n ) ∈ N n : α 1 + • • • + α n ≤ d}; M n ≤d := {(x → x α ) ∈ P(R n ) : α ∈ Λ n ≤d }. It is easy to show that |Λ n ≤d | = n+d d ,

and that the set M n ≤d of monomial functions forms a basis for

P ≤d (R n ). Sort the n-tuples in Λ n ≤d in colexicographic order, i.e. (α 1 , . . . , α n ) < (α 1 , . . . , α n ) if and only if α i < α i for the largest index i such that α i = α i . Let λ 1 < • • • < λ ( n+d d ) denote all the n+d d

n-tuples in Λ n ≤d after sorting. Analogously, let q 1 , . . . q ( n+d d ) be all the monomial functions in M n ≤d in this order, i.e. each q k : R n → R is given by the map x → x λ k . Given any sequence

(v 1 , . . . , v ( n+d d ) ) of vectors in R n , the generalized Vandermonde matrix associated to it is Q = Q[v 1 , . . . , v ( n+d d ) ] := [q i (v j )] 1≤i,j≤( n+d d ) .

(2)

Definition 5.5. Given any W ∈ W n,m n+d d

, we define the non-bias Vandermonde matrix of W to be the generalized Vandermonde matrix Q

[W ] := [q i ( w (1) j )] 1≤i,j≤( n+d d ) associated to ( w (1) 1 , . . . , w (1) n+d d

).

Theorem 5.6. Let m be arbitrary, let p ∈ P d (R n ), and suppose that p has all-non-zero coefficients. Also, suppose that p 1 , . . . , p k ∈ P(W n,m n+d d

) are fixed polynomial functions on the non-bias weights of the first layer. Define the following sets:

U := {W ∈ W n,m n+d d : p i (W ) = 0 for all 1 ≤ i ≤ k}; p U ind := {W ∈ U : F p,0n (W ) is linearly independent}.

If there exists W ∈ U such that the non-bias Vandermonde matrix of W is non-singular, then p U ind is dense in U (w.r.t. the Euclidean metric).

Proof. We essentially extend the ideas in the proofs of Theorem 5.1 and Corollary 5. such that the non-bias Vandermonde matrix of W is non-singular; see Appendix A.2 for proof details.

Remark 5.8. The proof of Corollary 5.7 still holds even if we restrict every non-bias weight w

(1) i,j in the first layer to satisfy |w (1) i,j | < λ for some fixed constant λ > 0.

For the rest of this section, let {λ k } k∈N be a divergent increasing sequence of positive real numbers, and let {Y k } k∈N be a sequence of closed intervals of R, such that Y k ⊆ Y k whenever k ≤ k, and such that each interval Y k = [y k , y k ] has length λ k . Let d ≥ 1 be an integer, and suppose σ ∈ C(R).

For each k ∈ N, let σ k be the best polynomial approximant to σ| Y k of degree d. Given r > 0 and any integer N ≥ 1, define the closed ball

B N r := {x ∈ R N : x 2 ≤ r}. Lemma 5.9. If d ≥ 2, lim k→∞ E d (σ| Y k ) = ∞,

and λ k ∈ Ω(k γ ) for some γ > 0, then for every ε > 0, there is a subsequence {k t } t∈N of N, and a sequence {y kt } t∈N of real numbers, such that y kt < y kt < y kt , σ(y kt ) = σ kt (y kt ), and

min{|y kt -y kt |, |y kt -y kt |} |y kt -y kt | > 1 d + 1 -ε,

for all t ∈ N. (See Appendix B for proof details.)

The proofs of the next three lemmas can be found in Appendix C. Lemma 5.10. For any constant γ > 0,

lim k→∞ σ k -σ ∞,Y k (λ k ) 1+γ = 0.

Lemma 5.11. Let K ≥ N ≥ 1 be integers, let r 0 , . . . , r N ≥ 1 be fixed real numbers, and let S(λ) be a set {p 0 (λ), . . . , p N (λ)} of N + 1 affinely independent points in R K , parametrized by λ > 0, where each point p i (λ) has (Cartesian) coordinates (λ ri p i,1 , . . . , λ ri p i,K ) for some fixed non-zero scalars p i,1 , . . . , p i,K . Let ∆(λ) be the convex hull of S(λ), i.e. ∆(λ) is an N -simplex, and for each 0 ≤ i ≤ N , let h i (λ) be the height of ∆(λ) w.r.t. apex p i (λ). Let h(λ) := max{h i (λ) : 0 ≤ i ≤ N } and r min := min{r 1 , . . . , r N }. If r j > r min for some 0 ≤ j ≤ N , then there exists some γ > 0 such that h(λ) ∈ Ω(λ rmin+γ ). Lemma 5.12. Let M, N ≥ 1 be integers, let τ > 0, and let 0 < θ < 1. Suppose ϕ : R M → R N is a continuous open map such that ϕ(0 M ) = 0 N , and ϕ(λx) ≥ λϕ(x) for all x ∈ R M , λ > 0. Let {U k } k∈N be a sequence where each U k is a dense subspace of B M λ k \B M θλ k . Then for every δ > 0, there exists some (sufficiently large) k ∈ N, and some points u 0 , . . . , u N in U k , such that for each point p ∈ B N τ , there are scalars b

0 , . . . , b N ≥ 0 satisfying p = N i=0 b i ϕ(u i ), b 0 + • • • + b N = 1, and |b i -1 N | < δ for all 0 ≤ i ≤ N .

Outline of strategy for proving Theorem 3.1. The first crucial insight is that P ≤d (R n ), as a real vector space, has dimension n+d d . Our strategy is to consider N = n+d d hidden units. Every hidden unit represents a continuous function g j : X → R determined by its weights W and the activation function σ. If g 1 , . . . , g N can be well-approximated (on X) by linearly independent polynomial functions in P ≤d (R n ), then we can choose suitable linear combinations of these N functions to approximate all coordinate functions f [t] (independent of how large m is). To approximate each g j , we consider a suitable sequence {σ λ k } ∞ k=1 of degree d polynomial approximations to σ, so that g j is approximated by a sequence of degree d polynomial functions { g W j,k } ∞ k=1 . We shall also vary W concurrently with k, so that w (1) j 2 increases together with k. By Corollary 5.7, the weights can always be chosen so that g W 1,k , . . . , g W N,k are linearly independent. The second crucial insight is that every function in P ≤d (R n ) can be identified geometrically as a point in Euclidean n+d d -space. We shall choose the bias weights so that g W 1,k , . . . , g W N,k correspond to points on a hyperplane, and we shall consider the barycentric coordinates of the projections of both f [t] and the constant function onto this hyperplane, with respect to g W 1,k , . . . , g W N,k . As the values of k and w (1) j 2 increase, both projection points have barycentric coordinates that approach ( 1 N , . . . , 1 N ), and their difference approaches 0; cf. Lemma 5.12. This last observation, in particular, when combined with Lemma 5.9 and Lemma 5.10, is a key reason why the minimum number N of hidden units required for the UAP to hold is independent of the approximation error threshold ε.

Proof of Theorem 3.1. Fix some ε > 0, and for brevity, let N = n+d d . Theorem 3.1 is trivially true when f is constant, so assume f is non-constant. Fix a point x 0 ∈ X, and define f

0 ∈ C(X, R m ) by f [t] 0 := f [t] -f [t] (x 0 ) for all 1 ≤ t ≤ m.

Next, let r X (x 0 ) := sup{ x -x 0 2 : x ∈ X}, and note that r X (x 0 ) < ∞, since X is compact. By replacing X with a closed tubular neighborhood of X if necessary, we may assume without loss of generality that r X (x 0 ) > 0.

Define {λ k } k∈N , {Y k } k∈N and {σ k } k∈N as before, with an additional condition that λ k ∈ Ω(k τ ) for some τ > 0. Assume without loss of generality that there exists a sequence {y k } k∈N of real numbers, such that y k < y k < y k , σ(y k ) = σ k (y k ), and

EQUATION

)

for all k ∈ N. The validity of this assumption in the case

lim k→∞ E d (σ| Y k ) = ∞ is given by Lemma 5.9. If instead lim k→∞ E d (σ| Y k ) < ∞, then as k → ∞, the sequence {σ k } k∈N converges to some σ ∈ P ≤d (R).

Hence, the assumption is also valid in this case, since for any y ∈ R such that σ( y) = σ( y), we can always choose {Y k } k∈N to satisfy

y k +y k 2

= y for all k ∈ N, which then allows us to choose {y k } k∈N that satisfies lim k→∞

min{|y k -y k |,|y k -y k |} λ k = 1 2 > 1 d+2 .

By Lemma 5.10, we may further assume that

σ k -σ ∞,Y k < ε(λ k ) 1+γ C

for all k ∈ N, where C > 0 and γ > 0 are constants whose precise definitions we give later. Also, for any W ∈ W n,m N , we can choose σ ∈ C(R) that is arbitrarily close to σ in the uniform metric, such that ρ σ W -ρ σ W ∞,X is arbitrarily small. Since σ ∈ C(R)\P ≤d-1 (R) by assumption, we may hence perturb σ if necessary, and assume without loss of generality that every σ k is a polynomial of degree d with all-non-zero coefficients, such that σ k (y k ) = 0.

For every r > 0 and k ∈ N, let W r := {W ∈ W n,m N : w (1) j 2 ≤ r for all 1 ≤ j ≤ N }, and define

λ k := sup r > 0 : {y k + w (1) j • (x -x 0 ) ∈ R : x ∈ X, W ∈ W r } ⊆ Y k for all 1 ≤ j ≤ N . Each λ k is well-defined, since r X (x 0 ) < ∞. Note also that λ k r X (x 0 ) = min{|y k -y k |, |y k -y k |} by definition, hence it follows from (3) that λ k λ k < (d + 2)r X (x 0 ).

In particular, {λ k } k∈N is a divergent increasing sequence of positive real numbers.

Given any p ∈ P ≤d (R n ), let ν(p) ∈ R N denote the vector of coefficients with respect to the basis {q 1 (x -x 0 ), . . . , q N (x -x 0 )} (i.e. if ν(p) = (ν 1 , . . . , ν N ), then p(x) = 1≤i≤N ν i q i (x -x 0 )), and let ν(p) ∈ R N -1 be the truncation of ν(p) by removing the first coordinate. Note that q 1 (x) is the constant monomial, so this first coordinate ν 1 is the coefficient of the constant term. For convenience, let ν i (p) (resp. ν i (p)) be the i-th entry of ν(p) (resp. ν(p)).

For each k ∈ N, W ∈ W λ k , 1 ≤ j ≤ N , define functions g W j,k , g W j,k in C(X) by x → σ(w (1) j • (1, x)) and x → σ k (w (1) j • (1, x)

) respectively. By definition, ν i ( g W j,k ) can be treated as a function of W , and note that ν i ( g λW j,k ) = λ deg qi ν i ( g W j,k ) for any λ > 0. (Here, deg q i denotes the total degree of q i .) Since deg q i = 0 only if i = 1, it then follows that ν i ( g λW j,k ) ≥ λ ν i ( g W j,k ) for all λ > 0. For each k ∈ N, define the "shifted" function

σ k : Y k → R by y → σ k (y + y k ). Next, let W k := σ k W ind n,N ;x0 ∩ (W λ k \W 0.5λ k

), and suppose W ∈ W k . Note that in the definition of W k , we do not impose any restrictions on the bias weights. Thus, given any such W , we could choose the bias weights of W (1) to be w (1) j,0 = y k -w (1) j • x 0 for all 1 ≤ j ≤ N . This implies that each g W j,k

EQUATION

By Corollary 5.7 and Remark 5.3, W k is dense in (W λ k \W 0.5λ k ), so such a W exists (with its bias weights given as above). By the definition of σ k W ind n,N ;x0 , we infer that { g W 1,k , . . . , g W N,k } is linearly independent and hence spans P ≤d (X). Thus, for every 1 ≤ t ≤ m, there exist a

[t] 1,k , . . . , a [t] N,k ∈ R, which are uniquely determined once k is fixed, such that f [t] 0 = a [t] 1,k g W 1,k +• • •+a [t]

N,k g W N,k . Evaluating both sides of this equation at x = x 0 , we then get

a [t] 1,k + • • • + a [t]

N,k = 0.

(5)

For each ∈ R, define the hyperplane

H := {(u 1 , . . . , u N ) ∈ R N : u 1 = }. Recall that q 1 (x)

is the constant monomial, so the first coordinate of each ν( g W j,k ) equals σ(y k ), which implies that

ν( g W 1,k ), . . . , ν( g W N,k ) are N points on H σ(y k ) ∼ = R N -1 . Let c f := max{ ν(f [t] ) 2 : 1 ≤ t ≤ m}.

(This is non-zero, since f is non-constant.) Note that 0 N -1 and ν(f [t] ) (for all t) are points in

B N -1 c f

. So for any δ > 0, Lemma 5.12 implies that there exists some sufficiently large k ∈ N such that we can choose some W ∈ W k , so that there are non-negative scalars b

[t] j,k , b j,k (for 1 ≤ j ≤ N , 1 ≤ t ≤ m) contained in the interval ( 1 N -δ, 1 N + δ) that satisfy the following: b [t] 1,k + • • • + b [t] N,k = b 1,k + • • • + b N,k = 1 (for all 1 ≤ t ≤ m); 0 N -1 = N j=1 b j,k ν( g W j,k ); ν(f [t] ) = N j=1 b [t] j,k ν( g W j,k ) (for all 1 ≤ t ≤ m).

Note that ν(f

[t] 0 + σ(y k )) = b [t] 1,k ν( g W 1,k ) + • • • + b [t] N,k ν( g W N,k ) and (0 N -1 , σ(y k )) = b 1,k ν( g W 1,k ) + • • • + b N,k ν( g W N,k ), so we get f [t] 0 = b [t] 1,k -b 1,k g W 1,k + • • • + b [t] N,k -b N,k g W N,k .

Since a

[t] 1,k , . . . , a [t]

N,k are unique (for fixed k), we infer that a

[t] j,k = b [t]

j,k -b j,k for each 1 ≤ j ≤ N . Thus, for this sufficiently large k, it follows from b

EQUATION

Let S k := { ν( g W 1,k ), . . . , ν( g W N,k )}, let ∆ k be the convex hull of S k , and for each j, let h j (∆ k ) be the height of

∆ k w.r.t. apex ν( g W j,k ). Let h(∆ k ) := max{h j (∆ k ) : 1 ≤ j ≤ N }. Since ν i ( g λW j,k ) = λ deg qi ν i ( g W j,k

) for all i, and since d ≥ 2 (i.e. deg q N > 1), it follows from Lemma 5.11 that there exists some γ > 0 such that h(∆ k ) ∈ Ω((λ k ) 1+γ ). Using this particular γ > 0, we infer that there exists some constant 0 < C < ∞ such that

(λ k ) 1+γ

h(∆ k ) < C for all sufficiently large k. Note that 2δ is an upper bound of the normalized difference for each barycentric coordinate of the two points ν(f [t] ) and 0

N -1 (contained in B N -1 c f ), which satisfies 2δ ≤ c f h(∆ k ) = c f (λ k ) 1+γ • λ k λ k 1+γ • (λ k ) 1+γ h(∆ k ) < c f (λ k ) 1+γ [(d + 2)r X (x 0 )] 1+γ C . (7) Now, define C := 2N c f [(d + 2)r X (x 0 )] 1+γ C > 0.

Thus, for sufficiently large k, it follows from (5), ( 6) and (7) that

EQUATION

For this sufficiently large k, define g ∈ C(X, R m ) by g

[t] = a [t] 1,k g W 1,k + • • • + a [t]

N,k g W N,k for each t. Using ( 4) and ( 8), it follows that

f [t] 0 -g [t]

∞,X = a

[t] 1,k (g W 1,k -g W 1,k ) + • • • + a [t] N,k (g W N,k -g W N,k ) ∞,X ≤ |a [t] 1,k | • g W 1,k -g W 1,k ∞,X + • • • + |a [t] N,k | • g W N,k -g W N,k ∞,X < ε.

Finally, for all 1 ≤ t ≤ m, let w

(2) j,t = a [t] j,k for each 1 ≤ j ≤ N , and let w

(2)

0,t = f [t] (x 0 ). This gives ρ σ [t] W = g [t] + f [t] (x 0 ). Therefore, the identity f [t] = f [t] 0 + f [t] (x 0 ) implies f -ρ σ W ∞,X < ε.

Notice that for all δ > 0, we showed in (6) that there is a sufficiently large k such that a [t] j,k ≥ -2δ. A symmetric argument yields a [t] j,k ≤ 2δ. Thus, for all λ > 0, we can choose W so that all non-bias weights in W (2) are contained in the interval (-λ, λ); this proves assertion (i) of the theorem.

Note also that we do not actually require δ > 0 to be arbitrarily small. Suppose instead that we choose k ∈ N sufficiently large, so that the convex hull of S k contains 0 N -1 and ν(f [t] ) (for all t). In this case, observe that our choice of k depends only on f (via ν(f [t] )) and σ (via the definition of {λ k } k∈N ). The inequality (7) still holds for any δ satisfying b [t] j,k , b j,k ∈ ( 1 N -δ, 1 N + δ) for all j, t. Thus, our argument to show f -ρ σ W ∞,X < ε holds verbatim, which proves assertion (ii).

Proof of Theorem 3.2. Fix some ε > 0, and consider an arbitrary t ∈ {1, . . . , m}. For each integer d ≥ 1, let p [t] d be the best polynomial approximant to f [t] of degree d. By Theorem 2.2, we have

f [t] -p [t] d ∞,X ≤ 6ω f [t] ( D 2d

) for all d ≥ 1, hence it follows from the definition of d ε that

f [t] -p [t] dε ∞,X ≤ 6ω f [t] ( D 2dε ) < ε. Define ε := ε-max{6ω f [t] ( D 2dε ) : 1 ≤ t ≤ m}. Note that ε > 0, and f [t] -p [t]

dε ∞,X ≤ ε-ε (for all 1 ≤ t ≤ m). By Theorem 3.1, there exists some W ∈ W n+dε dε satisfying p [t] dε -ρ 

σ [t] W ∞,X < ε for all 1 ≤ t ≤ m, which implies f [t] -ρ σ [t] W ∞,X ≤ f [t] -p [t] dε ∞,X + p [t] dε -ρ σ [t] W ∞,X < (ε -ε ) + ε = ε, therefore f -ρ σ W ∞,X < ε. Conditions (i) and (ii) follow from Theorem 3.1. Finally, note that ω f [t] ( D 2d ) ∈ O( 1 d ) (for fixed D), i.e. d ε ∈ O( 1 ε ), hence n+dε dε = n(n-1)...(n-dε+1) n! ∈ O(ε -n ).

6. CONCLUSION AND FURTHER REMARKS:

Theorem 5.6 is rather general, and could potentially be used to prove analogs of the universal approximation theorem for other classes of neural networks, such as convolutional neural networks and recurrent neural networks. In particular, finding a single suitable set of weights (as a representative of the infinitely many possible sets of weights in the given class of neural networks), with the property that its corresponding "non-bias Vandermonde matrix" (see Definition 5.5) is non-singular, would serve as a straightforward criterion for showing that the UAP holds for the given class of neural networks (with certain weight constraints). We formulated this criterion to be as general as we could, with the hope that it would applicable to future classes of "neural-like" networks.

We believe our algebraic approach could be emulated to eventually yield a unified understanding of how depth, width, constraints on weights, and other architectural choices, would influence the approximation capabilities of arbitrary neural networks.

Finally, we end our paper with an open-ended question. The proofs of our results in Section 5 seem to suggest that non-bias weights and bias weights play very different roles. We could impose very strong restrictions on the non-bias weights and still have the UAP. What about the bias weights? K. Wolsson. Linear dependence of a function set of m variables with vanishing generalized Wronskians. Linear Algebra Appl., 117:73-80, 1989 . ISSN 0024-3795. doi: 10.1016 A GENERALIZED WRONSKIANS AND THE PROOF OF THEOREM 5.6

First, we recall the notion of generalized Wronskians as given in (LeVeque, 1956, Chap. 4.3) . Let ∆ 0 , . . . , ∆ N -1 be any N differential operators of the form

∆ k = ∂ ∂x1 α1 • • • ∂ ∂xn αn , where α 1 + • • • + α n ≤ k.

Let f 1 , . . . , f N ∈ P(R n ). The generalized Wronskian of (f 1 , . . . , f N ) associated to ∆ 0 , . . . , ∆ N -1 is defined as the determinant of the matrix M = [∆ i-1 f j (x)] 1≤i,j≤N . In general, (f 1 , . . . , f N ) has multiple generalized Wronskians, corresponding to multiple choices for ∆ 0 , . . . , ∆ N -1 .

A.1 PROOF OF THEOREM 5.6

For brevity, let N = n+d d and let x = (x 1 , . . . , x n ). Recall that λ 1 < • • • < λ N are all the n-tuples in Λ n ≤d in the colexicographic order. For each

1 ≤ i, k ≤ N , write λ k = (λ k,1 , . . . , λ k,n ), define the differential operator ∆ λ k = ∂ ∂x1 λ k,1 • • • ∂ ∂xn λ k,n , and let α (i)

λ k be the coefficient of the monomial q k (x) in ∆ λi p(x). Consider an arbitrary W ∈ U, and for each 1 ≤ j ≤ N , define f j ∈ P ≤d (R n ) by the map x → p(w (1) 1,j x 1 , . . . , w

(1) n,j x n ). Note that F p,0n (W ) = (f 1 , . . . , f N ) by definition. Next, define the matrix M W (x) := [∆ i f j (x)] 1≤i,j≤N , and note that det M W (x) is the generalized Wronskian of (f 1 , . . . , f N ) associated to ∆ 1 , . . . , ∆ N . In particular, this generalized Wronskian is well-defined, since the definition of the colexicographic order implies that λ k,1 + • • • + λ k,n ≤ k for all possible k. Similar to the univariate case, (f 1 , . . . , f N ) is linearly independent if (and only if) its generalized Wronskian is not the zero function (Wolsson, 1989) . Thus, to show that W ∈ p U ind , it suffices to show that the evaluation det M W (1 n ) of this generalized Wronskian at x = 1 n gives a non-zero value, where 1 n denotes the all-ones vector in R n .

Observe that the (i, j)-th entry of M W (1 n ) equals ( w

(1) j ) λi (∆ λi p)( w (1) j ), hence we can check that M W (1 n ) = M M ,

where M is an N -by-N matrix whose (i, j)-th entry is given by

M i,j = α (i) λj -λi , if λ j -λ i ∈ Λ n ≤d ; 0, if λ j -λ i ∈ Λ n ≤d ; and where M = Q[W ] is the non-bias Vandermonde matrix of W .

It follows from the definition of the colexicographic order that λ j -λ i necessarily contains at least one strictly negative entry whenever j < i, hence we infer that M is upper triangular. The diagonal entries of M are α 

i) 0n = (λ i,1 ! • • • λ i,n !)α (1) λi for each 1 ≤ i ≤ N , where λ i,1 ! • • • λ i,n

! denotes the product of the factorials of the entries of the n-tuple λ i . In particular, λ i,1 ! • • • λ i,n ! = 0, and α (1) λi , which is the coefficient of the monomial q i (x) in p(x), is non-zero. Thus, det(M ) = 0.

We have come to the crucial step of our proof. If we can show that det(M ) = det(Q[W ]) = 0, then det(M W (1 n )) = det(M ) det(M ) = 0, and hence we can infer that W ∈ p U ind . This means that p U ind contains the subset U ⊆ U consisting of all W such that Q[W ] is non-singular.

Note that det(Q[W ]:

) is a polynomial in terms of the non-bias weights in W (1) as its variables, so we could write this polynomial as r = r(W ). Consequently, if we can find a single W ∈ U such that Q[W ] is non-singular, then r(W ) is not identically zero on U, which then implies that U = {W ∈ U : r(W ) = 0} is dense in U (w.r.t. the Euclidean metric).

A.2 PROOF OF COROLLARY 5.7

Let N := n+d d . By Theorem 5.6, it suffices to show that there exists some W ∈ W n,m N such that the non-bias Vandermonde matrix of W is non-singular. Consider W ∈ W n,m N such that w

(1) i,j = (w (1) 1,j ) (d+1) i . Recall that the monomials in M n ≤d are arranged in colexicographic order, i.e.

1, x 1 , x 2 1 , . . . , x d 1 , x 2 , x 1 x 2 , x 2 1 x 2 , . . . , x 2 2 , x 1 x 2 2 , . . . , x d n . Thus, there are fixed integers 0

= β 1 < β 2 < • • • < β N , such that the (i, j)-th entry of Q[W ] is (w (1)

1,j ) βi . Such matrices are well-studied in algebraic combinatorics, and the determinant of Q[W ] is a Schur polynomial; see (Stanley, 1999) . In particular, if we choose positive pairwise distinct values for w

(1) 1,j (for 1 ≤ j ≤ N ), then Q[W ] is non-singular, since a Schur polynomial can be expressed as a (non-negative) sum of certain monomials; see (Stanley, 1999, Sec. 7 .10) for details.

B AN ANALOG OF KADEC'S THEOREM AND THE PROOF OF LEMMA 5.9

Throughout this section, suppose σ ∈ C(R) and let d ≥ 1 be an integer. We shall use the same definitions for {λ k } k∈N , {Y k } k∈N and {σ k } k∈N as given immediately after Remark 5.8. Our goal for this section is to prove Theorem B.1 below, so that we can infer Lemma 5.9 as a consequence of Theorem B.1. Note that Theorem B.1 is an analog of the well-known Kadec's theorem (Kadec, 1960) from approximation theory. To prove Theorem B.1, we shall essentially follow the proof of Kadec's theorem as given in (Kadec, 1963) .

We begin with a crucial observation. For every best polynomial approximant σ k to σ| Y k of degree d, it is known that there are (at least) d + 2 values

y k ≤ a (k) 0 < a (k) 1 < • • • < a (k) d+1 ≤ y k ,

and some sign δ k ∈ {±1}, such that σ(a (Rivlin, 1981, Thm. 1.7) . Define Note that e k ∈ C(R) and φ k+1 ∈ P ≤d (R). Since y k+1 ≤ a (k) i ≤ y k+1 by assumption, it follows from the definition of σ k+1 that -E d (σ| Y k+1 ) ≤ e k+1 (a (k) i ) ≤ E d (σ| Y k+1 ). By the definition of a (k) i , we have e k (a (k) i ) = (-1 .

(k) i ) -σ k (a (k) i ) = (-1) i δ k E d (σ| Y k ) for all 0 ≤ i ≤ d + 1; see

∆ k := max a (k) i -y k y k -y k - i d + 1 : 0 ≤ i ≤ d + 1 .

) i δ k E d (σ| Y k ). Consequently, E d (σ| Y k ) -E d (σ| Y k+1 ) ≤ (-1) i δ k (e k -e k+1 )(a (k) i ) ≤ E d (σ| Y k ) + E d (σ| Y k+1 ), or equivalently, -E d (σ| Y k ) -E d (σ| Y k+1 ) ≤ (-1) i δ k φ k+1 (a (k) i ) ≤ E d (σ| Y k+1 ) -E d (σ| Y k ). Since Y k ⊆ Y k+1 implies E d (σ| Y k ) ≤ E d (σ| Y k+1 ), it follows that a 2i-1 ≤ a (k)

It then follows from (Kadec, 1963, Lem. 2) that

EQUATION

)

where θ is an arbitrary real number satisfying 0 < θ < 1 2 . Therefore, the assertion follows by letting γ = 1+τ D .

Proof of Lemma 5.9. Fix ε > 0. By Theorem B.1, we have lim inf k→∞ ∆ k λ k k γ = 0 for any γ > 0. Thus, by the definition of lim inf, there exists a subsequence {k t } t∈N of N such that

∆ k t λ k t (k t ) γ < ε

for all t ∈ N (given any γ > 0). Since λ k is at least Ω(k γ ) for some γ > 0, we can use this particular γ to get that lim inf t→∞ λ k t (k t ) γ > 0. Consequently, there is a subsequence {k t } t∈N of {k t } t∈N such that |∆ kt | < ε for all t ∈ N. Since d ≥ 2 by assumption, it then follows that ). Recall that any modulus of continuity ω f is subadditive (i.e. ω f (x + y) ≤ ω f (x) + ω f (y) for all x, y); see (Rivlin, 1981, Chap. 1) . Thus for fixed d, we have

1 d + 1 -ε < a (kt

ω σ| Y k ( λ k 2d ) ∈ O(λ k ), which implies k → σ k -σ ∞,Y k ∈ o(λ 1+γ k ).

C.2 PROOF OF LEMMA 5.11

Our proof of Lemma 5.11 is a straightforward application of both the Cayley-Menger determinant formula and the Leibniz determinant formula. For each 0 ≤ i ≤ N , let S i (λ) := S(λ)\{p i (λ)}, and let ∆ i (λ) be the convex hull of S i (λ). Let V(∆(λ)) (resp. V( ∆ i (λ))) denote the N -dimensional (resp. (N -1)-dimensional) volume of ∆(δ) (resp. ∆ i (λ)). Define the (N + 2)-by-(N + 2) matrix M (λ) = [M i,j (λ)] 0≤i,j≤N +1 as follows: M i,j (λ) = p i (λ) -p j (λ) 2 2 for all 0 ≤ i, j ≤ N ; M N +1,i (λ) = M i,N +1 (λ) = 1 for all 0 ≤ i ≤ N ; and M N +1,N +1 (λ) = 0.

The Cayley-Menger determinant formula gives [V(∆(λ))] 2 = (-1) N +1

(N !) 2 2 N det(M (λ)). Analogously, if we let M (λ) be the square submatrix of M (λ) obtained by deleting the first row and column from

:



Figures and Tables:

(k) j denote the truncation of w (k) j obtained by removing the first entry w (k) 0,j . Define the function ρ σ W : R n → R m so that for each 1 ≤ j ≤ m, the j-th coordinate function ρ σ [j]

Proof of Theorem 3.3. Most of the work has already been done earlier in the proofs of Theorem 3.1 and Theorem 3.2. The key observation is that det(Q[W ]) is a non-zero polynomial in terms of the weights W , hence {det(Q[W ]) = 0 : W ∈ W n+d d } is dense in W n+d d , or equivalently, its complement has Lebesgue measure zero.



If lim k→∞ E d (σ| Y k ) = ∞, then for any γ > 0, we have lim inf k→∞ ∆ k λ k k γ = 0.Proof. For every k ∈ N, define the functions e k := σ -σ k and φ k+1 := σ k -σ k+1 = e k+1 -e k .

a 2i (for each 0 ≤ i ≤ d + 1), where a 2i-1 and a 2i are the roots of the equation|φ k+1 (y)| = E d (σ| Y k+1 ) -E d (σ| Y k ). If E d (σ| Y k+1 ) = E d (σ| Y k ), then σ k+1 = σ k by definition, so we could set a (k+1) i = a (k) i for all i, i.e. there is nothing to prove in this case. Henceforth, assume E d (σ| Y k+1 ) = E d (σ| Y k ), and consider the polynomial function φ(y) := φ k+1 (y -y k ) E d (σ| Y k+1 ) -E d (σ| Y k )

λ k ) D = ∞ for any D > 1. If we compare the divergent series for any τ > 0), we thus getlim inf k→∞ ∆ k λ k k (1+τ )/D = 0.

4, using the notion of generalized Wronskians; see Appendix A.1 for proof details. Corollary 5.7. Let m be arbitrary. If p ∈ P(R) is a fixed polynomial function of degree d with all-non-zero coefficients, then p W ind

Table content: <table><tr><td>n, n+d d</td><td>d ;0n is a dense subset of W n,m n+d</td><td>.</td></tr><tr><td colspan="3">Proof. By Theorem 5.6, it suffices to show that there is some W ∈ W n,m n+d</td></tr><tr><td/><td/><td>d</td></tr></table>

/0024-3795(89) 90548-X. URL https://doi-org.library.sutd.edu.sg:2443/10.1016/0024-3795(89)90548-X.

Table content: <table><tr><td>Dmitry Yarotsky. Error bounds for approximations with deep relu networks. Neural Networks,</td></tr><tr><td>94:103 -114, 2017. ISSN 0893-6080. doi: https://doi.org/10.1016/j.neunet.2017.07.002. URL</td></tr><tr><td>http://www.sciencedirect.com/science/article/pii/S0893608017301545.</td></tr><tr><td>Gilad Yehudai and Ohad Shamir. On the power and limitations of random features for understanding</td></tr><tr><td>neural networks. Preprint arXiv:1904.00687v2 [stat.ML], June 2019.</td></tr></table>

(σ| Y k ) = ∞ by assumption, the infinite product (σ| Y k+1 ) -E d (σ| Y k ) E d (σ| Y k+1 ) + E d (σ| Y k )also diverges. It then follows from (9) that

Table content: <table><tr><td>Since lim the series</td><td>∞</td><td>E d</td><td>∞ k=0</td><td>E d (σ| Y k+1 ) E d (σ| Y k )</td><td>diverges, and thus</td></tr><tr><td colspan="2">k=0</td><td/><td/><td/><td/></tr></table>

Now σ -σ kt is continuous, so by the definition of a kt ) = σ kt (y kt ). From (10), we thus infer thatmin{|y k t -y k t |,|y k t -y k t |} λ k t Theorem 2.2 gives σ k -σ ∞,Y k = E d (σ| Y k ) ≤ 6ω σ| Y k ( λ k 2d

Table content: <table><tr><td>1</td><td>) λ kt -y kt</td><td>&lt;</td><td colspan="2">a (kt) 2 λ kt -y kt</td><td>&lt;</td><td>d d + 1</td><td colspan="2">+ ε.</td><td>(10)</td></tr><tr><td/><td/><td/><td>(kt) i</td><td colspan="4">, there is some a (kt) 1</td><td>&lt; y kt &lt; a (kt) 2</td><td>such that</td></tr><tr><td colspan="9">σ(y &gt; 1 d+1 -ε as desired.</td></tr><tr><td colspan="2">C PROOFS OF REMAINING LEMMAS</td><td/><td/><td/><td/><td/><td/></tr><tr><td>C.1 PROOF OF LEMMA 5.10</td><td/><td/><td/><td/><td/><td/><td/></tr></table>

Peer review (iclr2020):
This this exciting submission presents a new proof of Leshno*s version of the universal approximation property (UAP) for neural networks -- one of the foundational pillars of our understanding of neural networks. The new proof provides new insights into the universal approximation property. I consider these the main contribution of the paper. Specifically, the authors - provide an upper bound on the required width for the neural network - show that the approximation property still holds even if strong further requirements are imposed on the weights of the first or last layer. I rate this submission a weak accept. It’s a very good paper. The work makes useful contributions that should and will be of interest to many in the field. The paper is generally well-written. Some remarks: - Being somewhat long, the “Proof of Theorem 3.1” would be a much better read if the authors prefixed it with an outline of the strategy that the proof takes. - The authors point out that the lack of dependence of Theorem 3.1 on epsilon is surprising, and cite Lin’s work from 2017 who previously found such an independence. Lin’s derivation of the epsilon-independent UAP is much more intuitive than that of this submission, in which the epsilon independence really pops out somewhat magically and for me only made sense when I read the paper again. I would encourage the authors to add to Lin’s paper’s citation sentence that this paper motivates the epsilon independence well. Alternatively, the authors could add a few sentences to their paper to provide intuition on how the epsilon-independence comes about in their line of argument.