You're correct that a spherical charge density can be considered as a series of spherical shells. The field generated by a spherical shell of charge falls off at 1/r^2 outside of the sphere, meaning you can simplify it to a point charge. But inside the sphere, the generated field is 0. So to calculate the field at any point, you can indeed model the charge distribution as a series of shells and integrate, but you only integrate up to the radius of the point of interest, not all the way out. That means the integral will be different at each point, and so the resulting field will not fall off at 1/r^2 inside of the charge distribution.
Usually, the reason that you simplify a distribution to a point charge is because it allows you to calculate a single equivalent charge, and then easily determine the field at any point away from the distribution, because the field will just fall off at 1/r^2. If you're still inside of the distribution, and the field isn't falling off at 1/r^2, then this simplification doesn't do anything for you--you still have to do a separate integral for every point, so the math will be exactly the same as it was before.
I think the easiest way to solve this is just to ignore these kinds of simplifications for the moment, and go back to the relevant equations for electrostatics. Once they've been solved for this problem, the simplifications we've discussed here will emerge as a consequence of them (in fact I suspect that might have been the point of this problem in the first place.)