  We address a particular set of extended Su-Schrieffer-Heeger models with $2n$
sites in the unit cell [SSH($2n$)], that we designate by Sine-Cosine models
[SC$(n)$], with hopping terms defined as a sequence of $n$ sine-cosine pairs of
the form $\{\sin(\theta_j),\cos(\theta_j)\}$, $j=1, \cdots,n$. These models,
when squared, generate a block-diagonal matrix representation with one of the
blocks corresponding to a chain with uniform local potentials. We further focus
our study on the subset of SC$(2^{n-1})$ chains that, when squared an arbitrary
number of times (up to $n$), always generate a block which is again a
Sine-Cosine model, if an energy shift is applied and if the energy unit is
renormalized. We show that these $n$-times squarable models [SSC$(n)$] and
their band structure are uniquely determined by the sequence of energy unit
renormalizations and by the energy shifts associated to each step of the
squaring process. Chiral symmetry is present in all Sine-Cosine chains and edge
states levels at the respective central gaps are protected by it. Zero-energy
edge states in a SSC$(j)$ chain (with $j<n$) of the Matryoshka sequence
obtained squaring the SSC$(n)$ chain with open boundary conditions (OBC),
become finite energy edge states in non-central band gaps of the SSC$(n)$
chain. The extension to higher dimensions is discussed.
