cs.DSOct 23, 2026

Learning Positive Linear Combinations of ReLUs in Fully Polynomial Time

Authors: Ilias Diakonikolas, Jelena Diakonikolas, Puqian Wang, Nikos Zarifis

Organizations: UW Madison · CSAIL, MIT

Abstract

We study the fundamental problem of PAC learning positive linear combinations of kk arbitrary ReLU activations on Rd\mathbb{R}^d with respect to the Gaussian distribution, in L2L_2 loss. We give the first fully-polynomial time algorithm for this learning task. Specifically, our algorithm draws poly(d,k,1/ε)\mathrm{poly}(d, k, 1/ε) samples, runs in poly(d,k,1/ε)\mathrm{poly}(d, k, 1/ε) time, and outputs a hypothesis that is εε-close to the target concept with high probability. All prior algorithms for this problem require super-polynomial time in kk or 1/ε1/ε, or leverage strong well-conditioned assumptions on the weight matrix. On a technical level, our approach develops a novel boosting procedure of a weak'' learner that recovers the hidden weight vectors to constant error, in a list-decodable manner. To develop this weak'' learning procedure, we leverage implicit tensor computation techniques to construct a low-degree filtering procedure based on Hermite tensors of the target function. This procedure amplifies separation among (clusters of) weight vectors to learn each of them to constant error. We then argue that the same procedure can be used for improving the estimation of hidden weight vectors, using gradient descent-like updates on the Euclidean sphere. A crucial ingredient of our algorithm is a pruning approach that ensures only relevant vectors that approximate the hidden weight vectors at target accuracy are maintained by the algorithm.

Figures & tables

Explore similar work

Oct 18, 2026

Polynomial Learning of Positive ReLU Networks under Gaussian Inputs

We study positive, bias-free, one-hidden-layer ReLU networks under standard Gaussian inputs. We obtain algorithms with polynomial sample and real-arithmetic complexity, without separation or conditioning assumptions on the hidden directions. When the conditional label mean is a width-kk network with total output weight at most one, the learner achieves squared prediction error ε\varepsilon with output width at most k+min⁡{k,d}k+\min\{k,d\} and complexity linear in one plus the residual second moment. For arbitrary labels with bounded second moment, we obtain a positive ReLU network of width O(klog⁡(2/ε))O(k\log(2/\varepsilon)). Its risk is bounded by the irreducible noise variance plus Ck3log⁡9/2(2k/ε)Ck^3\log^{9/2}(2k/\varepsilon) times the squared distance of the conditional mean to the target class, plus ε\varepsilon. The proofs combine implicit Hermite moments, orthogonal tensor-train compression, and a direction decoder whose probability of returning a direction far from every atom is bounded linearly in the compression and estimation errors. Exact Gaussian localization refines the resulting direction covers. In the agnostic setting, a cone-constrained convex program with a scale-dependent penalty produces a positive network; in the conditional-mean setting, exact sparsification yields the smaller width. We also establish projective direction recovery and confidence amplification with sample complexity proportional to the logarithm of the reciprocal failure probability.
Jun 9, 2026

Robust Regression of General ReLUs with Queries

We study the task of agnostically learning general (as opposed to homogeneous) ReLUs under the Gaussian distribution with respect to the squared loss. In the passive learning setting, recent work gave a computationally efficient algorithm that uses poly(d,1/ε)poly(d,1/ε) labeled examples and outputs a hypothesis with error O(opt)+εO(opt)+ε, where optopt is the squared loss of the best fit ReLU. Here we focus on the interactive setting, where the learner has some form of query access to the labels of unlabeled examples. Our main result is the first computationally efficient learner that uses dpolylog(1/ε)+O~(min⁡{1/p,1/ε})d polylog(1/ε)+\tilde{O}(\min\{1/p, 1/ε\}) black-box label queries, where pp is the bias of the target function, and achieves error O(opt)+εO(opt)+ε. We complement our algorithmic result by showing that its query complexity bound is qualitatively near-optimal, even ignoring computational constraints. Finally, we establish that query access is essentially necessary to improve on the label complexity of passive learning. Specifically, for pool-based active learning, any active learner requires Ω~(d/ε)\tildeΩ(d/ε) labels, unless it draws a super-polynomial number of unlabeled examples.
Oct 16, 2026

"Deep lifts": Neural networks can generalize even with the dumbest optimizer

Given a labeled dataset {(xi,yi)}i=1n\{(x_i, y_i)\}_{i = 1}^n, consider drawing a parameter vector θ^\hatθ uniformly at random from those that interpolate the data within a width-WW depth-DD ReLU architecture with continuous weights bounded within [−M,M][-M, M]. We study the generalization properties of θ^\hatθ in the case where the dataset can, in fact, be explained by a smaller ww-by-dd ReLU network. Despite the obliviousness of the selection rule for θ^\hatθ, it adaptively generalizes at the rate O~(wdWDn)\tilde O\left(\frac{wdWD}{n}\right), rather than being governed by the VC dimension of the full WW-by-DD hypothesis class. The reason is due to a particular inductive bias in ReLU networks: the ww-by-dd parameter witnesses an exponential volume of WW-by-DD interpolating parameters, which dwarfs the volume of poorly generalizing parameters (conditioned on interpolation). We explicitly construct this set of near-equivalent WW-by-DD parameters by a proof technique called deep lifting, and then tie the achievable volume to stability properties of the ww-by-dd witness. Examples show that these techniques allow us to convert approximation-theoretic results into generalization guarantees for nonparametric classes. We further show that generalization of uniform-random interpolators can survive label corruptions.