In this paper, the non-leptonic three-body decays [Formula: see text], [Formula: see text], [Formula: see text] are studied by introducing two-meson distribution amplitude for the [Formula: see text], [Formula: see text] and [Formula: see text] pairs in naive and QCD factorization (QCDF) approaches, such that the analysis is simplified into quasi-two body decays. By considering that the vector meson is being ejected in factorization, the resonant and non-resonant contributions are analyzed by using intermediate mesons in Breit–Wigner resonance formalism and the heavy meson chiral perturbation theory (HMChPT), respectively. The calculated values of the resonant and non-resonant branching ratio of [Formula: see text], [Formula: see text] and [Formula: see text] decay modes are compared with the experimental data. For [Formula: see text] and [Formula: see text], the non-resonant contributions are about 70–80% of experimental data, for which the total results by considering resonant contributions are in good agreement with the experiment.