解二次同余

Question

我想在R中解决这个等式x ^ 2 = 4 mod 3，我尝试了很多命令，但这些命令与这个等式无关。请帮我找到命令。 例如，我试过

modlin(2,4,3)
# [1] 2

modpower(2,4,3)
# [1] 1

以及其他五个命令。

Answer 1

这是一个自制的函数，它解决了x ^ 2 = n mod m：

f <- function(n,m){which(((0:(m-1))^2 %% m)==(n %% m))[1] - 1}

> f(4,3)
[1] 1
> sapply(-10:10,f,5)
 [1]  0  1 NA NA  2  0  1 NA NA  2  0  1 NA NA  2  0  1 NA NA  2  0
>