A density matrix divide-and-conquer method is proposed for electronic structure calculation of large molecules. It is based on partition of density matrix and thus applicable to both density-functional and Hartree-Fock method. Compared to the original formulation with electron density, the present method is more efficient and as accurate.