We study the luminosity function and the correlation function of about 1200 z~4 Lyman break galaxies (LBGs) with i'<26 that are photometrically selected from deep BRi' imaging data of a 618 arcmin^2 area in the Subaru/XMM-Newton Deep Field taken with Subaru Prime Focus Camera. The contamination and completeness of our LBG sample are evaluated, on the basis of the Hubble Deep Field-North (HDF-N) objects, to be 17% and 45%, respectively. We derive the UV (rest 1700A) luminosity functions (LFs) and find a large population of UV-luminous galaxies at z~4. The LFs of the red and blue subsamples imply that the bright LBGs are redder in the UV continuum than the average color of the LBGs. Then we calculate the correlation function over theta = 2''-1000'' and find that it is fitted fairly well by a power law, omega(theta)=A_omega theta^(-0.8), with A_omega=0.71 +/- 0.26. We estimate the correlation length r_0 (in comoving units) of the two-point spatial correlation function xi(r) = (r/r_0)^(-1.8) to be r_0=2.7 +0.5/-0.6 h^(-1) Mpc (Omega_m=0.3 and Omega_Lambda=0.7). The correlation function shows an excess of omega (theta) on small scales (theta < 5''), departing from the power-law fit at > 3 sigma significance level. Interpreting this as being due to galaxy mergers, we evaluate the fraction of galaxies undergoing mergers to be 3.0 +/- 0.9%, which is significantly smaller than those of galaxies at intermediate redshifts.