Abstract: A general method for classifying the possible quark models of a multi-Higgs-doublet model, in the presence of Abelian symmetries, is presented. All the possible sets of textures that can be present in a given sector are shown, thus turning the determination of the flavor models into a combinatorial problem. Several symmetry implementations are studied for two and three Higgs doublet models. Some models' implementations are explored in great detail, with a particular emphasis on models known as Branco-Grimus-Lavoura and nearest-neighbor-interaction. Several considerations on the flavor changing neutral currents of multi-Higgs models are also given.