Time Dependent Density Functional Theory offers a good accuracy/computational cost ratio among different methods used to predict the electronic structure for molecules of practical interest. The Coupled Electronic Oscillator (CEO) formalism was recently shown to accurately predict Nonlinear Optical (NLO) properties of organic chromophores when combined with Time Dependent Density Functional Theory. Unfortunately, CEO does not lend itself easily to interpretation of the structure activity relationships of chromophores. On the other hand, the Sum Over States formalism in combination with semiempirical wavefunction methods has been used in the past for the design of simplified essential states models. These models can be applied to optimization of NLO properties of interest for applications. Unfortunately, TD-DFT can not be combined directly with SOS because state-to-state transition dipoles are not defined in the linear response TD approach. In this work, a second order CEO approach to TD-DFT is simplified so that properties of double excited states and state-to-state transition dipoles may be expressed through the combination of linear response properties. This approach is termed the a posteriori Tamm-Dancoff approximation (ATDA), and validated against high-level wavefunction theory methods. Sum over States (SOS) and related Two-Photon Transition Matrix formalism are then used to predict Two-Photon Absorption (2PA) profiles and anisotropy, as well as Second Harmonic Generation (SHG) properties. Numerical results for several conjugated molecules are in excellent agreement with CEO and finite field calculations, and reproduce experimental measurements well.