Recursive solutions of rational interpolation problems via fast matrix factorization