In the scenario that a solid-fuel launch vehicle maneuvers in outer space at high angles of attack and sideslip for energy management, Approximate Analytical Solutions (AAS) for the three-dimensional (3D) ascent flight states are derived, which are the only solutions capable of considering time-varying Mass Flow Rate (MFR) at present. The uneven MFR makes the thrust vary nonlinearly and thus increases the difficulty of the problem greatly. The AAS are derived based on a 3D Generalized Ascent Dynamics Model (GADM) with a normalized mass as the independent variable. To simplify some highly nonlinear terms in the GADM, several approximate functions are introduced carefully, while the errors of the approximations relative to the original terms are regarded as minor perturbations. Notably, a finite series with positive and negative exponents, called Exponent-Symmetry Series (ESS), is proposed for function approximation to decrease the highest exponent in the AAS so as to reduce computer round-off errors. To calculate the ESS coefficients, a method of seeking the Optimal Interpolation Points (OIP) is proposed using the least-squares-approximation theory. Due to the artful design of the approximations, the GADM can be decomposed into two analytically solvable subsystems by a perturbation method, and thus the AAS are obtained successfully. Finally, to help implement the AAS, two indirect methods for measuring the remaining mass and predicting the burnout time in flight are put forward using information from accelerometers. Simulation results verify the superiority of the AAS under the condition of time-varying MFR.