I'll show the direct computation, because I think it is the best way to see this. Work in this map $\phi$, which is either defined on the upper hemisphere $\{x_{n+1} >0\}$, or lower hemisphere $\{x_{n+1}<0\}$,
$$ \phi(x_1,\dots,x_{n+1}) = (x_1,\dots,x_n,\sum_{i=1}^{n+1}x_i^2 -R^2) = (y_1,\dots,y_{n+1}) $$
Write $(\vec{e}_1,\dots,\vec{e}_{n+1})$ the canonical basis of $\mathbb{R}^{n+1}$. Then compute the matrices of $d\phi$ and $d\phi^{-1}$ to see that for all $i\leq n$,
$$ \partial y_i = \frac{\partial}{\partial y_i} = \vec{e}_i - \frac{x_i}{x_{n+1}}\vec{e}_{n+1} $$
Differentiate the $\partial y_i$ in euclidean $\mathbb{R}^{n+1}$,
$$ \tilde{\nabla}_{\partial y_i}\partial y_j = \left(-\frac{x_ix_j}{x_{n+1}^3}-\frac{\delta_{ij}}{x_{n+1}}\right) \vec{e}_{n+1} $$
Then we need to project $\vec{e}_{n+1}$ orthogonally on the $n$-sphere. For that, check that the metric tensor is $g_{ij} = \delta_{ij} + \frac{x_ix_j}{x_{n+1}^2}$ and its inverse is $g^{ij} = \delta_{ij} - \frac{x_ix_j}{R^2}$. The orthogonal projection of $\vec{e}_{n+1}$ is therefore $-\sum_{i=1}^n \frac{x_ix_{n+1}}{R^2}\partial y_i$. This yields the Christoffel symbols of $\nabla$ in map $\phi$ :
$$ \nabla_{\partial y_i}\partial y_j = \left(\frac{x_ix_j}{x_{n+1}^2}+\delta_{ij}\right)\sum_{k=1}^n\frac{x_k}{R^2}\partial y_k, \quad\quad \text{i.e.} \quad\quad \Gamma_{ij}^k = \left(\frac{x_ix_j}{x_{n+1}^2}+\delta_{ij}\right)\frac{x_k}{R^2} $$
As you can see, the Christoffels are fairly simple. Continue to compute the right part of the Riemann curvature tensor
$$ \sum_{\lambda=1}^n \Gamma^\rho_{\mu\lambda}\Gamma^\lambda_{\nu\sigma} - \Gamma^\rho_{\nu\lambda}\Gamma^\lambda_{\mu\sigma} = \frac{x_\rho(x_\mu\delta_{\nu\sigma}-x_\nu\delta_{\mu\sigma})}{R^2 x^2_{n+1}} $$
and the left part
$$ \partial y_\mu\Gamma^\rho_{\nu\sigma} - \partial y_\nu\Gamma^\rho_{\mu\sigma} = \left(\frac{x_\nu x_\sigma}{x^2_{n+1}}+\delta_{\nu\sigma} \right)\frac{\delta_{\rho\mu}}{R^2} - \left(\frac{x_\mu x_\sigma}{x^2_{n+1}}+\delta_{\mu\sigma} \right)\frac{\delta_{\rho\nu}}{R^2} +\frac{x_\rho x_\nu\delta_{\sigma\mu} - x_\rho x_\mu\delta_{\sigma\nu}}{R^2 x^2_{n+1}} $$
Finally, compute the traces to see that the Ricci tensor $R_{\sigma\nu}=\sum_{\rho=1}^n R^\rho_{\sigma\rho\nu}$ is equal to $\frac{n-1}{R^2}$ times the metric tensor. That means the scalar curvature is $\frac{n(n-1)}{R^2}$, as requested.
You see there is no magic : just do it. In the process you got a lot more insight into the geometry of the $n$-sphere than just its scalar curvature : its Ricci tensor, Riemann curvature tensor and the Christoffels.