Fast spline smoothing via spectral factorization concepts