A mixed finite element method with exactly divergence-free velocities for incompressible magnetohydrodynamics