Superfast Iterative Solvers for Linear Matrix Equations