Based on the flavor structure of four-quark effective operators, we develop an automatic computation program to calculate hadronic two-body $B$ meson decay amplitudes, and apply it to their global analysis in the perturbative QCD (PQCD) approach. Fitting the PQCD factorization formulas for $B\to PP,VP$ decays at leading order in the strong coupling $\alpha_s$ to measured branching ratios and direct CP asymmetries, we determine the Gegenbauer moments in light meson light-cone distribution amplitudes (LCDAs). It is found that most of the fitted Gegenbauer moments of the twist-2 and twist-3 LCDAs for the pseudoscalar meson $P$ ($P=\pi$, $K$) and vector meson $V$ ($V=\rho$, $K^*$) agree with those derived in QCD sum rules. The shape parameter for the $B_s$ meson distribution amplitude and the weak phase $\phi_3(\gamma)=(75.1\pm2.9)^\circ$ consistent with the value in Particle Data Group are also obtained. It is straightforward to extend our analysis to higher orders and higher powers in the PQCD approach, and to the global determination of LCDAs for other hadrons.