A sharp estimate for the decreasing rearrangement of the length of the gradient of solutions to a class of nonlinear Dirichlet and Neumann elliptic boundary value problems is established under weak regularity assumptions on the domain. As a consequence, the problem of gradient bounds in norms depending on global integrability properties is reduced to one-dimensional Hardy-type inequalities. Applications to gradient estimates in Lebesgue, Lorentz, Zygmund, and Orlicz spaces are presented.