Until now, the analysis of long wave run-up on a plane beach has been focused on finding its maximum value, failing to capture the existence of resonant regimes. One-dimensional numerical simulations in the framework of the nonlinear shallow water equations are used to investigate the boundary value problem for plane and nontrivial beaches. Monochromatic waves, as well as virtual wave-gage recordings from real tsunami simulations, are used as forcing conditions to the boundary value problem. Resonant phenomena between the incident wavelength and the beach slope are found to occur, which result in enhanced run-up of nonleading waves. The evolution of energy reveals the existence of a quasiperiodic state for the case of sinusoidal waves. Dispersion is found to slightly reduce the value of maximum run-up but not to change the overall picture. Run-up amplification occurs for both leading elevation and depression waves.