The effective linear and nonlinear optical properties of metal/dielectric composite media, in which ellipsoidal metal inclusions are distributed in shape, are investigated. The shape distribution function P(L x, L y) is assumed to be 2Δ-2θ(L x - 1/3 + Δ/3)θ(L y - 1/3 + Δ/3)θ(2/3 + Δ/3 - L x - L y), where θ( . . . ) is the Heaviside function, Δ is the shape variance and Li are the depolarization factors of the ellipsoidal inclusions along i-symmetric axes (i = x, y). Within the spectral representation, we adopt Maxwell-Garnett type approximation to study the effect of shape variance Δ on the effective nonlinear optical properties. Numerical results show that both the effective linear optical absorption α ∼ ωIm( ) and the modulus of the effective third-order optical nonlinearity enhancement |χ(3) e|/χ(3) 1 exhibit the nonmonotonic behavior with Δ. Moreover, with increasing Δ, the optical absorption and the nonlinearity enhancement bands become broad, accompanied with the decrease of their peaks. The adjustment of Δ from 0 to 1 allows us to examine the crossover behavior from no separation to large separation between optical absorption and nonlinearity enhancement peaks. As Δ → 0, i.e., the ellipsoidal shape deviates slightly from the spherical one, the dependence of |χ(3) e|/χ(3) 1 on Δ becomes strong first and then weak with increasing the imaginary part of inclusions' dielectric constant. In the dilute limit, the exact formula for the effective optical nonlinearity is derived, and the present approximation characterizes the exact results better than old mean field one does.