Conjugated oligofluorenes are chemically reduced up to five charges in tetrahydrofuran solvent and confirmed with clear spectroscopic evidences. Stimulated by these experimental results, we have conducted a comprehensive computational study of the electronic structure and the solvation structure of representative oligofluorene anions with a focus on the pairing between sodium ions abd these multi-anions. Using density functional theory (DFT) methods and a solvation model of both explicit solvent molecules and implicit polarizable continuum, we first elucidate the structure of tightly solvated free sodium ions, and then explore the pairing of sodium ions either in contact with reduced oligofluorenes or as solvent-separated ion pairs. Computed time-dependent-DFT absorption spectra are compared with experiments to assign the dominant ion pairing structure for each multi-anion. Computed ion pair binding energies further support our assignment. The availability of different length and reducing level of oligofluorenes enables us to investigate the effects of total charge and charge density on the binding with sodium ions, and our results suggest both factors play important roles in ion pairing for small molecules. However, as the oligofluorene size grows, its charge density determines the binding strength with the sodium ion.