一道初等数论题

一道初等数论题

问题

解一次同余方程
(1)$3x\equiv 5(\mathrm{m}\mathrm{o}\mathrm{d}7)$　　(2)$103x\equiv 57(\mathrm{m}\mathrm{o}\mathrm{d}211)$

解答

问题一是十分简单的因为显然$3x\equiv 12\equiv 5(\mathrm{m}\mathrm{o}\mathrm{d}7)$,而$(3,7)=1$，所以可以除去3,即$x\equiv 4(\mathrm{m}\mathrm{o}\mathrm{d}7)$.

而问题二就比问题一稍稍棘手.因为首先211和103都是比较大的质数,所以利用欧拉定理将其化为首1的方程是不现实的.那么我们当然考虑和问题一一样,通过将余数加上模数进行等价变形（即更换代表元）来将其化为首1的方程.可惜的是,这个数不是特别好找的.我们可以看以下的寻找过程
$$\begin{array}{ccccc}& ⋮& & & \\ (103,57+3\times 211)& =& (103,690)& =& 1\\ (103,57+2\times 211)& =& (103,479)& =& 1\\ (103,57+1\times 211)& =& (103,268)& =& 1\\ (103,57+0\times 211)& =& (103,57)& =& 1\\ (103,57-1\times 211)& =& (103,-154)& =& 1\\ (103,57-2\times 211)& =& (103,-365)& =& 1\\ (103,57-3\times 211)& =& (103,-576)& =& 1\\ & ⋮& & & \end{array}$$
可以发现直接寻找是比较困难的,所以我们不妨使用带余除法来寻找代表元.即考虑方程$103\mid 57+211a^{\prime }$,或者将其写成
$$57+211a^{\prime }=103b^{\prime },a^{\prime },b^{\prime }\in \mathbb{Z}\text{\hspace{1em}(1)}$$
这其实是一个一次不定方程,我们将其写成标准形式
$$211a+103b=57,a,b\in \mathbb{Z}$$
由于$(211,103)=1\mid 57$,利用带余除法可以解出$a,b$
$$\begin{array}{ccccccc}211& =& 103& \times & 2& +& 5\\ 103& =& 5& \times & 20& +& 3\\ 5& =& 3& \times & 1& +& 2\\ 3& =& 2& \times & 1& +& 1\end{array}$$
从下往上回代可得到$a=-41\times 57=-2366,b=84\times 57=4788$.
4788就是我们要找的一个代表元,所以,
$$x\equiv 4788\equiv 146(\mathrm{m}\mathrm{o}\mathrm{d}211)\text{\hspace{1em}(2)}$$
从另一个角度看,如果将$(1)$式写成
$$103x=211q+57,x,q\in \mathbb{Z}$$
那么我们可以将此法更直接地理解为用定义改写$103x\equiv 57(\mathrm{m}\mathrm{o}\mathrm{d}211)$,然后利用带余除法求解.利用一次不定方程的解法可得到其所有解
$$\{\begin{array}{l}x=4788+211t\\ q=-2366+103t\end{array},t\in \mathbb{Z}\text{\hspace{1em}(3)}$$
事实上$(3)$中$x$的表达式与$(2)$是等价的.

暴力点的话…

能暴力解决的事，为什么要动脑子呢？　　　　——我说的

其实这题完全可以编程暴力解决【滑稽】.显然在$[1,211]$上必有一整数满足方程,所以我们的循环甚至都不需要写成定数的,直接让循环指标从1开始自增就行.

Talk is cheap, show me the code.

#include <stdio.h>

int main(){
    int i = 0;
    int state = 1;
    while(state){
        i += 1;
        if((103*i)%211==57) state = 0;
    }
    printf("%d", i);
}

输出为

146

整挺好,跟我算出来的一样.