Advancing Hierarchical Equations of Motion for Efficient Evaluation of Coherent Two-dimensional Spectroscopy