  The mechanism of the flux-line-lattice (FLL) melting in anisotropic high-T_c
superconductors in ${\bf B}\parallel {\bf \hat{c}}$ is clarified by Monte Carlo
simulations of the 3D frustrated XY model. The percentage of entangled flux
lines abruptly changes at the melting temperature T_m, while no sharp change
can be found in the number and size distribution of vortex loops around T_m.
Therefore, the origin of this melting transition is the entanglement of flux
lines. Scaling behaviors of physical quantities are consistent with the above
mechanism of the FLL melting. The Lindemann number is also evaluated without
any phenomenological arguments.
