We present numerical simulations on the propagation of ultra-high-frequency protons with energies of 10 19.5 -10 22 eV in extragalactic magnetic fields over 1 Gpc. We use the Optical Redshift Survey (ORS) galaxy sample, which allows us to accurately quantify the contribution of nearby sources to the energy spectrum and the arrival distribution, as a source model. The sample is corrected, taking the selection effect and absence of galaxies in the zone of avoidance (jbj < 20 ) into account. We calculate three observable quantities-cosmicray spectrum, harmonic amplitude, and two-point correlation function-from our data of numerical simulations. With these quantities, we compare the results of our numerical calculations with the observation. We find that the arrival distribution of ultra-high-energy cosmic rays (UHECRs) becomes most isotropic as sources are restricted to to luminous galaxies ðM lim ¼ À20:5Þ. However, it is not isotropic enough to be consistent with the Akeno Giant Air Shower Array (AGASA) observation, even for M lim ¼ À20:5. In order to obtain a sufficiently isotropic arrival distribution, we randomly select sources more luminous than À20.5 mag from the ORS sample, which contribute to the observed cosmic-ray flux, and investigate the dependence of the results on their number. We show that three observable quantities, including the GreisenZatsepin-Kuz'min (GZK) cutoff of the energy spectrum, can be reproduced in the case that the number fraction $10 À1.7 of the ORS galaxies more luminous than À20.5 mag is selected as UHECR sources. In terms of the source number density, this constraint corresponds to $10 À6 Mpc À3 . However, since the mean number of sources within the GZK sphere is only $0.5 in this case, the eight AGASA events above 10 20.0 eV, which do not together constitute clustered events, cannot be reproduced. On the other hand, if the cosmic-ray flux measured by the High Resolution Fly's Eye EHE Cosmic-Ray Detector (HiRes), which is consistent with the GZK cutoff, is correct and observational features about the arrival distribution of UHECRs are same as the AGASA, our source model can explain both the arrival distribution and the flux at the same time. Thus, we conclude that a large fraction of the eight AGASA events above 10 20 eV might originate in the top-down scenarios or that the cosmic-ray flux measured by the HiRes experiment might be better. We also discuss the origin of UHECRs below 10 20.0 eV through comparisons between the number density of astrophysical source candidates and our result ($10 À6 Mpc À3 ).