We have made semi-analytical studies to investigate the configurations of caustics and the probability distribution of the flux factor K for the binary microlensing including external shears. A parametric equation of critical curve is derived in a 4th order complex polynomial. We present the topological dependencies of the caustics for selected gamma parameters (0, 0.3, 0.6, 1.3, 2.0, and 2.5) and convergence terms (0., 0.4, 0.8, 1.2, 1.6, and 2.0). For the purpose of analyzing the efficiency of High Amplification Event (HAE) on each caustics, we examine the probability distribution of the flux factor by a Monte Carlo method. Changing the separation of the binary system from 0.8 to 1.3 (in normalied unit), we examine the probability distribution of the K-values in various gamma parameters. The relationship between gamma parameters, seperations and their probabilties of the flux factor K have been studied. Our results show that the relatively higher K values (K>1.5) are increased as increasing the separation of the binary system. We therfore conclude that, in the N-body microlensing, the probabilities of higher HAEs are inversely proportional to the star density as well. We also point out that the present research might be used as a preliminary step toward investigating heavy N-body microlensing simulations.