  Superconductivity in the Hubbard model on a square lattice near half filling
is studied using an optimization (or correlated) variational Monte Carlo
method. Second-order processes of the strong-coupling expansion are considered
in the wave functions beyond the Gutzwiller factor. Superconductivity of
d_x^2-y^2-wave is widely stable, and exhibits a crossover around U=U_co\sim 12t
from a BCS type to a new type. For U\gsim U_co (U\lsim U_co), the energy gain
in the superconducting state is derived from the kinetic (potential) energy.
Condensation energy is large and \propto exp(-t/J) [tiny] on the strong [weak]
coupling side of U_co. Cuprates belong to the strong-coupling regime.
