While reading a paper by Owens (arXiv:1906.05913) about embeddings of rational homology balls in the complex projective plane, I found out the following somewhat unexpected number theory corollary (here $F(n)$ denotes the $n$-th Fibonacci number, with $F(1)=1$ and $F(2)=1$).
If $n>1$, the square of the odd Fibonacci number $F(2n+1)$ can be written as the sum of exactly $F(2n+1)+1$ nonzero squares.
I will put an answer with the proof. However, I am curious about whether a proof of this result can be obtained by purely numer-theoretic methods.
It follows from a result in Conway's little book The Sensual Quadratic Form, and a bit of computation, that every number $n \geq 34$ is the sum of five nonzero squares. That missing number $33$ is a one-off. There is a proof in Niven and Zuckerman where $34$ is replaced by the easier $170,$ pages 318-319 in the Fifth edition (with Montgomery, 1991).
With $k \geq 6:$ any number $n \geq k + 14$ can be expressed as the sum of $k$ nonzero squares.
The numbers that are not the sum of five nonzero squares are: $ 1,2, 3, 4, 6, 7, 9, 10, 12, 15, 18, 33 $ This list is at OEIS with references
not six: $1, 2, 3, 4, 5, 7, 8, 10, 11, 13, 16, 19.$
The induction step: if $n \geq n_k$ means that $n$ is the sum of $k$ nonzero squares, then $ n \geq 1 + n_k$ means that $n$ is the sum of $1+k$ nonzero squares