Discrete symmetries in the three-Higgs-doublet model

N-Higgs-doublet models (NHDM) are among the most popular examples of electroweak symmetry breaking mechanisms beyond the Standard Model. Discrete symmetries imposed on the NHDM scalar potential play a pivotal role in shaping the phenomenology of the model, and various symmetry groups have been studied so far. However, in spite of all efforts, the classification of finite Higgs-family symmetry groups realizable in NHDM for any N>2 is still missing. Here, we solve this problem for the three-Higgs-doublet model by making use of Burnside's theorem and other results from pure finite group theory which are rarely exploited in physics. Our method and results can be also used beyond high-energy physics, for example, in study of possible symmetries in three-band superconductors.