Efficient design methods for large-scale metalenses are crucial for various applications. The conventional phase-mapping method shows a weak performance under large phase gradients, thus limiting the efficiency and quality of large-scale, high-numerical-aperture metalenses. While inverse design methods can partially address this issue, existing solutions either accommodate only small-scale metalenses due to high computational demands or compromise on focusing performance. We propose an efficient large-scale design method based on an optimization approach combined with the adjoint-based method and the level-set method, which first forms a one-dimensional metalens and then extends it to two dimensions. Taking fabrication constraints into account, our optimization method for large-area metalenses with a near-unity numerical aperture (NA = 0.99) has improved the focusing efficiency from 42% to 60% in simulations compared to the conventional design method. Additionally, it has reduced the deformation of the focusing spot caused by the ultrahigh numerical aperture. This approach retains the benefits of the adjoint-based method while significantly reducing the computational burden, thereby advancing the development of large-scale metalenses design. It can also be extended to other large-scale metasurface designs.