Computing the multiplicity structure in solving polynomial systems