теория-чисел - Квадратичные сравнения по простому модулю

Решаю сравнение $%x^2=3022(\mod 7001)$%. В данном примере $% p=3022$%, а $%p=7001$%
Получаю следующие значения переменных:
Раскладываю $%р-1$%, по следующей формуле: $%p-1=q \cdot 2^r$%;
$%7001 - 1 = 875 \cdot 2^3$%
Следовательно $%r=3$%; $%q=875$%;
Найду $% b $% по формуле $% b=a^q(\mod 7001) $%
Ищу наименьший квадратичный невычет:
$%L(2;7001)=1$%,
$%L(3;7001)=-1$%, следовательно $%f=3$%;
$%g=f^q=3477$%;
Ищу наименьшее неотрицательное $%m$% при котором остаток от $%b^{2^m}$% будет равен 1 по модулю 7001. Такое значение получается при $%m=2$%. Рассчитаю значение $%k$%: $%k=2^{r-m}=2^{3-2}=2$%;
По формуле Шенкса $%x = a^{(q-1)/2}g^{k/2} (\mod 7001)$% получаю ответ $%3664$%, что не есть верно. wolfram дает ответ $%+(-)155$%, что является правильным ответом.
Подскажите, пожалуйста, где допущена ошибка, для большинства других значений, написанная программа работает, при расчете вручную получаю то же самое неверное решение.