As a continuation of a previous paper [J. Phys. Soc. Jpn. 56 (1987) 1490], the variational Monte-Carlo method is extended to include the antiferromagnetic long-range order. The theory is based on the Gutzwiller-type correlation factor and its effect is exactly taken into account by the Monte-Carlo procedure. An application is made to the half-filled-band case of one-dimensional lattice (50 sites), two-dimensional square lattice (up to 20x20 sites) and three-dimensional simple cubic lattice (6 x 6 x 6 sites). The result is qualitatively different from previous studies relying on the random-phase-type “Gutzwiller approximation.” The variational energy for two and three dimensions is favorably compared with Hirsch's quantum Monte-Carlo data.