I have not seen AR described in terms of shocks, only MA. Whether the past errors are an "actual" error term or not is a fascinating question. It is not an
observable error term but population standard deviations are not observable (knowable) either usually - that does not mean they don't exist
I think they solve MA indirectly through working through AR calculations perhaps although it has been a very long time since I read that. MA and AR are related a MA of one form can represent a AR of another form which is why in ARIMA you may use a MA of 1 rather than an AR of 2 in some situations- effectively they are the same and the rule of simplicity comes into play.
Thankfully I don't worry about things such as how algorithms work - not smart enough to understand it anyway. I found the Hyndman book on transfer functions essentially incomprehensible.
If you are learning ARIMA I strongly recommend the Duke website on it. It is truly awesome although designed for practice not theory.