We analyze charmless two-body non-leptonic B decays under the framework of factorization assisted topological amplitude approach. Compared with the conventional flavor diagram approach, we consider flavor SU(3) breaking effect assisted by factorization hypothesis for topological diagram amplitudes of different decay modes, factorizing out the corresponding decay constants and form factors. The non-perturbative parameters of topology diagram magnitudes χ and strong phase ϕ are universal that can be extracted by χ2 fit from current abundant experimental data of charmless B decays. The number of free parameters and the χ2 per degree of freedom are both reduced comparing with previous analysis. With these best fitted parameters, we predict branching fractions and CPasymmetry parameters of nearly 100 Bu,d and Bs decay modes. The long-standing ππ and πK-CP puzzles are resolved simultaneoulsy.