Inverse design of large-area metasurfaces