Most Tensor Problems Are NP-Hard 3: The Clique Number Is the Unique l for Which an Explicit Graph Tensor Has Spectral Norm 1Research Paper
Motivation
The spectral norm of a matrix, its largest singular value, is computable in polynomial time to any accuracy. Its analogue for 3-tensors appears in the analysis of multilinear systems, in the best rank-1 approximation of a data array, in quantum information (the geometric measure of entanglement), and in the injective norm of tensor products. Hillar and Lim, Most Tensor Problems Are NP-Hard (J. ACM 60(6), 2013, Article 45; final version arXiv:0911.1393v5), show that, in contrast with matrices, deciding whether a fixed nonzero rational number is a singular value or the spectral norm of a rational 3-tensor is NP-hard (their Theorems 6.5 and 1.10).
The proof is a reduction from the clique number of a graph. It combines three classical ingredients: the theorem of Motzkin and Straus (1965), which expresses as the value of a quadratic program over the simplex; a theorem of Banach (1938) on the spectral norm of symmetric tensors; and an observation from He, Li and Zhang (2010) on maxima of sums of squared bilinear forms. Deciding is one of Karp's (1972) NP-complete problems.
Setting
A real 3-tensor is an array with trilinear form . With the Euclidean norm , its spectral norm is
A real is a unit -singular value of if there are unit vectors with
for all : the stationarity conditions of on a product of unit spheres.
Let be a simple graph on , , with edges , , and clique number . Let be the standard simplex, the adjacency matrix and the all-ones matrix. For a positive integer set and . Let , and let be the maximum over unit , of the multilinear form
The graph tensor has as entry the coefficient of in (26).
Formalization targets
Goal: Theorems 6.5 and 1.10, encoding-free core
For every graph on vertices and every enumeration of its edges:
so is the largest for which is a unit singular value of . The paper's -query Turing reductions rest on exactly this.
Milestones (in attack order)
- Theorem 6.7 (Motzkin–Straus): .
- (22): ; so , for , for .
- Lemma 6.8: , where .
- Theorem 6.9 (Banach), (24): for a symmetric 4-tensor, .
- Proposition 6.10 (He–Li–Zhang): for symmetric , .
- Lemma 6.11: , hence iff , iff , iff .
- Proof of Theorem 6.5: , is a unit singular value of , and every unit singular value has absolute value at most .
- §7 (content of Theorem 1.13): , attained only with .
Significance
The result. The goal shows that a single comparison with of the spectral norm, or of the singular-value spectrum, of an explicit rational tensor of size encodes whether . It is the step that makes the spectral norm, the largest singular value, and (through §7) the best rank-1 approximation of a 3-tensor NP-hard, whereas all three are polynomial-time for matrices. It explains why the higher-order power method and its relatives carry no global guarantees, and it is the starting point for inapproximability results such as the paper's Theorem 1.11.
Formalizing it. The results are proved on paper; none of them is formalized. Motzkin–Straus, Banach's polarization theorem for symmetric 4-tensors, and the He–Li–Zhang identity are classical results that are not in Mathlib and not on the platform; each is reusable well beyond this mission (Turán-type bounds, polynomial optimization over spheres, injective norms). The mission also checks the paper's construction at the level of indices and constants.
Difficulty
The construction of is explicit, and the comparison with is arithmetic once is known. Everything rests on the three cited theorems. Motzkin–Straus is an optimization statement over the simplex whose maximizers may lie on the boundary and need not be unique, so a first-order computation at an interior critical point does not settle it. Banach's theorem asserts equality, not equality up to a constant: the standard polarization identity bounds the 4-linear form by its diagonal values only up to a factor. Proposition 6.10 is where the paper's own argument is loose: the 4-linear form it introduces is not symmetric, so (24) does not apply as printed. Finally, the singular-value clauses require that the maximum of a trilinear form on a product of spheres is a stationary value. This must be proved, not assumed.
Formalization scope
Tensors are arrays Fin l → Fin m → Fin n → ℝ with 0-based indices: the paper's is A (i-1) (j-1) (k-1). The Euclidean norm is Real.sqrt (∑ i, x i ^ 2). Suprema and maxima are sSup of sets that are nonempty and bounded under the stated hypotheses; where the paper's "max" carries content (Motzkin–Straus, (22), §7), attainment is stated with IsGreatest/IsLeast. The graph is a SimpleGraph (Fin v) with , clique number Mathlib's SimpleGraph.cliqueNum, and its edges are enumerated by an arbitrary Fin e ≃ G.edgeSet. Every statement holds for each enumeration, so the tensor depends on the choice and the conclusions do not. is a rational tensor, read over by casting.
Explicit readings of the paper:
- Definition 6.1 has no norm constraint. Since preserves (20), unnormalized singular values make every nonzero real a singular value of , so the clause "not for " would be false. Singular vectors are therefore unit vectors, as in the paper's derivation on p. 0:21.
- The statement of Lemma 6.11 prints . There is no in §6, and the proof uses , which is formalized.
- is written entrywise as , which equals for a loopless edge.
- " is just the maximum -singular value of " is made explicit as milestone 7.
- The goal does not claim that fails to be a singular value of for ; the paper neither claims nor uses this.
This is the mathematical core of a polynomial-time Turing reduction, not an NP-hardness theorem. The polynomial size of , the -query procedure, the scaling to a general fixed , and the NP-completeness of CLIQUE (Karp 1972, cited) are not formalized. Theorem 1.11 and Corollary 1.12 are out of scope.
Two formalizations would make the goal trivial, and both are excluded. One is dropping the unit-norm constraint, which makes every nonzero number a singular value. The other is a spectral norm taken as a supremum over all vectors without normalization, which makes it unbounded and junk-valued; here the quotients are normalized.
Contributions are welcome on each milestone. Motzkin–Straus, Banach (24) and He–Li–Zhang are self-contained and reusable. Lemma 6.11 and milestone 7 need Cauchy–Schwarz and a Lagrange-type stationarity argument on spheres, which Mathlib supports through IsLocalExtrOn and the Lagrange multiplier files.
Selected references
- C. J. Hillar, L.-H. Lim, Most Tensor Problems Are NP-Hard, J. ACM 60(6), Article 45, 2013. https://doi.org/10.1145/2512329 (final version arXiv:0911.1393v5, https://arxiv.org/abs/0911.1393v5)
- T. S. Motzkin, E. G. Straus, Maxima for graphs and a new proof of a theorem of Turán, Canadian J. Math. 17, 533–540, 1965. https://doi.org/10.4153/CJM-1965-053-6
- S. Banach, Über homogene Polynome in (L²), Studia Math. 7, 36–44, 1938. https://doi.org/10.4064/sm-7-1-36-44
- S. He, Z. Li, S. Zhang, Approximation algorithms for homogeneous polynomial optimization with quadratic constraints, Math. Program. Ser. B 125(2), 353–383, 2010. https://doi.org/10.1007/s10107-010-0409-z
- R. M. Karp, Reducibility among combinatorial problems, in Complexity of Computer Computations, Plenum, 85–103, 1972. https://doi.org/10.1007/978-1-4684-2001-2_9