Модульные уравнения в Haskell

Я хочу решить линейные и квадратичные модульные уравнения в Haskell в одной переменной. То, как я делаю это прямо сейчас, - это положить x = [1..] в уравнении один за другим и найти остаток (expr `rem` p == 0, если уравнение по модулю p (не обязательно главным), где expr имеет переменную x). Я считаю, что это очень неэффективный процесс. Есть ли лучший способ сделать это?