Most tensor problems are NP-hard
Christopher Hillar, Lek-Heng Lim
TL;DR
The paper demonstrates that a broad class of multilinear tensor problems—including tensor eigenvalues, singular values, spectral norms, rank decompositions, and hyperdeterminants—are NP-hard, even when restricted to symmetric tensors. By constructing reductions from classical NP-hard problems like graph 3-colorability, the authors show that these problems remain hard in real, complex, and even certain finite-field settings, and they discuss related #P- and VNP-hardness results for eigenvector counting and hyperdeterminants. They also explore approximation limits, showing there is no PTAS for key tensor tasks like eigenvector approximation and spectral-norm evaluation unless P=NP, and they discuss the limits of computability, including undecidability results for certain bilinear and bivariate matrix-function problems. The work places tensor problems at a boundary between tractable linear/convex computations and intractable nonlinear/nonconvex problems, influencing both theoretical understanding and practical use of tensor methods. Overall, the article underscores the substantial computational obstacles in extending linear-algebraic tools to higher-order multilinear settings.
Abstract
We prove that multilinear (tensor) analogues of many efficiently computable problems in numerical linear algebra are NP-hard. Our list here includes: determining the feasibility of a system of bilinear equations, deciding whether a 3-tensor possesses a given eigenvalue, singular value, or spectral norm; approximating an eigenvalue, eigenvector, singular vector, or the spectral norm; and determining the rank or best rank-1 approximation of a 3-tensor. Furthermore, we show that restricting these problems to symmetric tensors does not alleviate their NP-hardness. We also explain how deciding nonnegative definiteness of a symmetric 4-tensor is NP-hard and how computing the combinatorial hyperdeterminant of a 4-tensor is NP-, #P-, and VNP-hard. We shall argue that our results provide another view of the boundary separating the computational tractability of linear/convex problems from the intractability of nonlinear/nonconvex ones.
