  We investigate the clustering properties of 2,600 Lyman Break Galaxies (LBGs)
at z=3.5-5.2 in two large blank fields, the Subaru Deep Field and the
Subaru/XMM Deep Field (600arcmin^2 each). The angular correlation functions of
these LBGs show a clear clustering at both z~4 and 5. The correlation lengths
are r_0= 4.1^{+0.2}_{-0.2} and 5.9^{+1.3}_{-1.7} h_{100}^{-1} Mpc (r_0=
5.1^{+1.0}_{-1.1} and 5.9^{+1.3}_{-1.7} h_{100}^{-1} Mpc) for all the detected
LBGs (for L>L* LBGs) at z~4 and 5, respectively. These correlation lengths
correspond to galaxy-dark matter biases of b_g= 2.9^{+0.1}_{-0.1} and
4.6^{+0.9}_{-1.2} (b_g=3.5^{+0.6}_{-0.7} and 4.6^{+0.9}_{-1.2}), for all the
detected LBGs (for L>L^* LBGs) at z~4 and 5, respectively. These results,
combined with estimates for z~3 LBGs in the literature, show that the
correlation length of L>L^* LBGs is almost constant, ~5 h_{100}^{-1} Mpc, over
z~3-5, while the bias monotonically increases with redshift at z>3. We also
find that for LBGs at z~4 the clustering amplitude increases with UV-continuum
luminosity and with the amount of dust extinction. We estimate the mass of dark
halos hosting various kinds of high-z galaxies including LBGs with the analytic
model given by Sheth & Tormen (1999). We find that the typical mass of dark
halos hosting L>L^* LBGs is about 1x10^{12} h_{70}^{-1}Msol over z~3-5, which
is comparable to that of the Milky Way Galaxy. A single dark halo with ~10^{12}
h_{70}^{-1} Msol is found to host 0.1-0.3 LBG on average but host about four
K-band selected galaxies.
