We present an analytical integral equation theory for polyelectrolyte solutions modeled as linear freely-jointed tangent hard-sphere polyanionic chains and cationic hard-sphere monomeric counterions embedded in a continuum dielectric medium. Each hard-sphere segment on the polyelectrolyte chain and hard-sphere counterion are univalent with unit diameters. The model was formulated in the context of the multi density Ornstein-Zernike integral equation theory within the mean spherical approximation. Analytical solutions for the model were obtained using the ideal chain approximation. The contact values of the radial distribution functions, internal energy, Helmholtz energy, osmotic pressure, and activity coefficient of the system were derived as a function of chain length, density, and Bjerrum length via the energy route. Predictions from the theory were compared with computer simulation data reported in the literature, and very good agreement was found.