Как идея. По формуле Эйлера x/a(mod p) = x*a^(p-2) (mod p). Можно перебирать x от l до r и пытаться делить в попытке прийти к a. Если получилось — печатаем. Если нет (зациклились), возможно запоминаем все значения как неудачные чтобы в будущем быстро отбрасываться.