We present a method for computing excitation energies for molecules in solvent, based on the combination of a minimal parameter implicit solvent model and the equation-of-motion coupled-cluster singles and doubles method (EOM-CCSD). In this method, the solvent medium is represented by a smoothly varying dielectric function, constructed directly from the quantum mechanical electronic density using only two tunable parameters. The solvent-solute electrostatic interactions are computed by numerical solution of the nonhomogeneous Poisson equation and incorporated at the Hartree-Fock stage of the EOM-CCSD calculation by modification of the electrostatic potential. We demonstrate the method by computing excited state transition energies and solvent shifts for several small molecules in water. Results are presented for solvated H2O, formaldehyde, acetone, and trans-acrolein, which have low-lying n → π* transitions and associated blue shifts in aqueous solution. Comparisons are made with experimental data and other theoretical approaches, including popular implicit solvation models and QM/MM methods. We find that our approach provides surprisingly good agreement with both experiment and the other models, despite its comparative simplicity. This approach only requires modification of the Fock operator and total energy expressions at the Hartree-Fock level-solvation effects enter into the EOM-CCSD calculation only through the Hartree-Fock orbitals. Our model provides a theoretically and computationally simple route for accurate simulations of excited state spectra of molecules in solution, paving the way for studies of larger and more complex molecules.