求解二次同余

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
>