Vibronic structure of nonlinear optical (NLO) coefficients is developed within the Condon approximation, displaced harmonic oscillators, and crude adiabatic states. The displacements of backbone modes of conjugated polymers are taken from vibrational data on the ground and 1B excited state. NLO resonances are modeled by three excitations and transition moments taken from Pariser-Parr-Pople (PPP) theory and optimized to polydiacetylene (PDA) spectra in crystals and films, with blue-shifted 1B exciton. The joint analysis of third-harmonic-generation, two-photon absorption, and nondegenerate four-wave-mixing spectra of PDA crystals and films shows weak two-photon absorption to 2A below 1B, leading to overlapping resonances in the THG spectrum, strong two-photon absorption to an nA state some 35% above 1B, and weak Raman resonances in nondegenerate FWM spectra. The full pi-pi* spectrum contributes to Stark shifts and field-induced transitions, as shown by PPP results for PDA oligomers. The Stark shift dominates high-resolution electroabsorption (EA) spectra of PDA crystals below 10 K. The close correspondence between EA and the first-derivative I’(omega) of the linear absorption above the 1B exciton in PDA crystals provides an experimental separation of vibrational and electronic contributions that limits any even-parity state in this 0.5 eV interval. An oscillator-strength sum rule is applied to the convergence of PDA oligomers with increasing length, N, and the crystal oscillator strengths are obtained without adjustable parameters. The sum rule for the 1B exciton implies large transition moments to higher-energy A(g) states, whose locations in recent models are contrasted to PPP results. Joint analysis of NLO and EA spectra clarifies when a few electronic excitations are sufficient, distinguishes between vibrational and electronic contributions, and supports similar pi-electron interactions in conjugated molecules and polymers.