求解二次同餘

Question

我想在R中求解這個方程x^2 = 4 mod 3，我嘗試了許多命令，但是這個命令與這個方程無關。請幫我找到命令。 例如，我試過

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

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

和其他五個命令。

Answer 1

這裏是一個自制功能，解決了X^2 =正模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 
>