Introductory textbooks have the responsibility of conveying a lot of knowledge to the beginner in a reasonably short amount of time. In this endeavor they have to make use of all kinds of metaphors in order to make the facts memorable. Only a small fraction of students would stay with a book that builds everything on axioms and lists all the preassumptions in tedious detail like a math book (actually, the math guys often wish that to be true, because they have been socialized differently, and they are confused by all the fuzzy assumptions they never heard about).
The problem begins if one of the mnemonics is taken too seriously. Then the learner is going to get stuck because he thinks that there has to be some sense in what this or that authority has claimed to be true.
Two things in your question strike the eye with respect to what I have said above:
- the claim that imaginary quantities are "unphysical"
- the expression $p=\sqrt{\dots}$ which is claimed to be the "classical formula for momentum"
As to 1), the imaginary unit is nothing special when it comes to physics. The mathematicians can do all sorts of interesting and meaningful stuff with complex numbers. But for a physicist what counts is whether the mathematical quantities are able to represent measurements (at least indirectly, e.g. via the probability amplitude in QM) in a way that allow predictions or at least correlations. You can simply represent the imaginary unit as a real matrix
$$i=\left( \begin{array}{cc} 0 & 1 \\ -1 & 0 \end{array}\right)$$
As you can easily check, $i^2=-1$ where $1$ denotes the unit matrix in 2D. Would you say that matrices indicate something "physically illegal"? Certainly no. And, as it turns out, the above matrix is the infinitesimal generator (the tangent, so to say) of rotations, and so you can most naturally describe two-dimensional rotations by unit complex numbers (you can even describe 3D rotations by a generalization of complex numbers, called quaternions!). The harmonic oscillator is such a kind of rotation: in a sense, energy "rotates" between kinetic and potential energy, or more precisely, the system rotates between maximum speed and maximum displacement. The fact that you have two degrees of freedom for such a rotation reflects the fact that the oscillator is second order, linear and has real coefficients (making the complex conjugate oscillator equivalent to the original one) and therefore, its solutions can be decomposed into sine and cosine solutions without losing information.
That is just an example of how the imaginary unit can represent a physical quantity. On the other hand, if you start out with the assumption of real coordinates and the solution of your differential equation or whatever yields complex or imaginary coordinates, you are in the trouble of explaining what that means in terms of measurements. Usually, one is not inclined to assume that complex position means the existence of some sort of extra-dimension. Ockham's razor demands that we first check the simplest explanation. And in the case of the harmonic oscillator it is simply that we are always able to choose real solutions and superpose them, so why arbitrarily choose complex solutions? In that sense, complex solutions of the harmonic oscillator are unphysical: we don't need them to explain reality, but we may use them if we explain what using them means (namely superposition of two orthogonal solutions in a very practical complex exponential expression). With quantum mechanics it is completely different: we can't explain "quantum reality" by a single-valued real wave function. We need to use $i$ (or an equivalent matrix) to make everything fit together. So, the imaginary unit is not unphysical in quantum mechanics, at least not in the same sense as in classical mechanics. It is, however, partially redundant because it is related to the gauge symmetry of the electromagnetic field. If we gauge the electromagnetic potentials differently, we can partly change the phase of the wave function (which determines its complex nature), but that is only up to $2\pi$ rotations.
As to 2), there are lots of relations in quantum mechanics which look suspiciously "classical", and there are also a lot of relations in quantum mechanics that are in striking contradiction to classical mechanics. Just because an arbitrary formula looks like a classical one if you are using the right letters of the alphabet doesn't mean that the quantities they represent are the same that you can measure in an experiment. After all, they are just letters. Replace $p$ by $R$ and it doesn't look so intriguing anymore. The reason why the authors have chosen to present the stuff that way was to wake up the reader and tell him: "hey man, there is pot of gold lying in that quantum mechanics, and if you dig this up, you are going to be rich".
Actually, there is no such thing as "classical momentum" in quantum mechanics, and particularly an imaginary momentum is meaningless in quantum mechanics. Maybe you can use that imaginary classical momentum to describe other things nicely (e.g. the tunnel effect if "classical momentum" becomes imaginary), but these descriptions are only useful in a limited context (otherwise we wouldn't need QM in the first place). Therefore, the formula you have given does not represent momentum. In QM, momentum of a given state of the system can only be defined as a set of probability amplitudes $\langle p|\psi\rangle$ (the Fourier components of the wave function, represented by the eigenstates $|p\rangle$ of the momentum operator) for the possible momentum measurements. If you measure the system, the probability of getting a certain momentum value is the probability amplitude squared. Ehrenfest's theorem connects that to classical mechanics by dealing with the expectation values, but that is only valid for a huge number of repeated measurements on the same system (or measuring an ensemble of equivalent systems). Momentum of a single measurement is undetermined (unless the system is in a momentum eigenstate, which you have ruled out by referring to the time-independent Schrödinger equation with a potential).
So, summarizing the above and short answer to your question: the imaginary unit has no predefined meaning in classical mechanics. It can be anything you want it to be and that it is able to mathematically represent (like rotations in phase space, rotations in physical space, etc.). Hence, it is not at all unphysical per-se, unless you have constructed the "unphysicalness" into it by using a bare mathematical formalism.