On the Use of Entire-Domain Basis Functions and Fast Factorizations for the Design of Modulated Metasurface