HDP算法