文章目录
- 二元一次不定方程|整数解
- 定理1整数解存在充要条件
- 定理2 通解
- 特解
- 知识回顾
- 利用辗转相除法求
- 例
- 使用表达式凑出通解
二元一次不定方程|整数解
二元一次不定方程的一般形式为
a
x
+
b
y
=
c
ax + by = c
ax+by=c(1)
其中 a a a、 b b b、 c c c 是整数,且 a a a、 b b b 都不为零。
解不定方程 (1) 就是要求出 (1) 式的所有整数解。
定理1整数解存在充要条件
方程 (1) 有整数解的充分必要条件是 ( a , b ) ∣ c (a,b)|c (a,b)∣c。
证明:
条件: ( a , b ) ∣ c (a,b)|c (a,b)∣c,
结论:方程(1): a x + b y = c ax + by = c ax+by=c
必要性是明显的。
由 ( a , b ) ∣ a (a,b)|a (a,b)∣a, ( a , b ) ∣ b (a,b)|b (a,b)∣b,且 x , y ∈ Z x,y\in\mathbb{Z} x,y∈Z,则由整除的性质有 ( a , b ) ∣ ( a x + b y ) (a,b)|(ax+by) (a,b)∣(ax+by),即 ( a ∣ b ) ∣ c (a|b)|c (a∣b)∣c
充分性:
由裴蜀恒等式知存在整数 u 0 u_0 u0、 v 0 v_0 v0,使
a
u
0
+
b
v
0
=
(
a
,
b
)
.
au_0 + bv_0 = (a,b).
au0+bv0=(a,b).(1-1)
,将方程(1-1)两边同时乘以
c
(
a
,
b
)
\frac{c}{(a,b)}
(a,b)c,得
a
u
0
c
(
a
,
b
)
au_{0}\frac{c}{(a,b)}
au0(a,b)c+
b
v
0
c
(
a
,
b
)
bv_{0}\frac{c}{(a,b)}
bv0(a,b)c=
c
c
c
由于 ( a , b ) ∣ c (a,b)|c (a,b)∣c,于是 c ( a , b ) ∈ Z \frac{c}{(a,b)}\in\mathbb{Z} (a,b)c∈Z, x 0 = u 0 c ( a , b ) x_0 = u_0 \frac{c}{(a,b)} x0=u0(a,b)c, y 0 = v 0 c ( a , b ) y_0 = v_0 \frac{c}{(a,b)} y0=v0(a,b)c 就是 (1) 式的一组整数解。
定理2 通解
若知道 (1) 式的一个特解 x 0 ( x 0 , y 0 ) \mathbf{x}_0(x_0, y_0) x0(x0,y0),则立刻就可以写出 (1) 式的全部解。
已知
(
x
0
,
y
0
)
(x_0, y_0)
(x0,y0)是(1)式的一个解,即
a
x
0
+
b
y
0
=
c
ax_{0}+by_{0}=c
ax0+by0=c,则(1)式的全部解为以下式(2)
{
x
=
x
0
+
b
(
a
,
b
)
t
,
y
=
y
0
−
a
(
a
,
b
)
t
,
t
=
0
,
±
1
,
±
2
,
⋯
\begin{cases} x = x_0 + \frac{b}{(a,b)}t, \\ y = y_0 - \frac{a}{(a,b)}t, \end{cases} \quad t=0,\pm 1, \pm 2, \cdots
{x=x0+(a,b)bt,y=y0−(a,b)at,t=0,±1,±2,⋯
证明:
易知(2)式给出的所有 ( x , y ) (x,y) (x,y)都满足方程(1) a x + b y = c ax + by = c ax+by=c。
反之,设 ( x , y ) (x,y) (x,y)是(1)式的任意一个解。证明其能够被表示为式(2)的形式
由
a
x
+
b
y
=
c
ax+by=c
ax+by=c及
a
x
0
+
b
y
0
=
c
,
ax_0+by_0=c,
ax0+by0=c,(2-0)
可得
a
(
x
−
x
0
)
+
b
(
y
−
y
0
)
=
0
,
a(x-x_0)+b(y-y_0)=0,
a(x−x0)+b(y−y0)=0,(2-1)
,对其两边乘以
1
(
a
,
b
)
\frac{1}{(a,b)}
(a,b)1并移项,
进而有
a
(
a
,
b
)
(
x
−
x
0
)
=
−
b
(
a
,
b
)
(
y
−
y
0
)
.
\frac{a}{(a,b)}(x-x_0)=-\frac{b}{(a,b)}(y-y_0).
(a,b)a(x−x0)=−(a,b)b(y−y0).(3)
可见 b ( a , b ) ∣ ( a ( a , b ) ( x − x 0 ) ) \frac{b}{(a,b)}|(\frac{a}{(a,b)}(x-x_{0})) (a,b)b∣((a,b)a(x−x0))
由 ( a ( a , b ) , b ( a , b ) ) = 1 \left(\frac{a}{(a,b)}, \frac{b}{(a,b)}\right)=1 ((a,b)a,(a,b)b)=1,知 b ( a , b ) ∣ ( x 0 − x ) \frac{b}{(a,b)}|(x_0-x) (a,b)b∣(x0−x)。
(注意,整除的性质:若 m ∣ a 1 a 2 m|a_{1}a_{2} m∣a1a2, ( m , a 1 ) = 1 (m,a_{1})=1 (m,a1)=1,则 m ∣ a 2 m|a_{2} m∣a2;在这里 m = b ( a , b ) m=\frac{b}{(a,b)} m=(a,b)b, a 1 a_{1} a1= a ( a , b ) \frac{a}{(a,b)} (a,b)a, a 2 = ( x − x 0 ) a_{2}=(x-x_{0}) a2=(x−x0))
因此存在整数
t
t
t,使得
x
0
−
x
x_{0}-x
x0−x=
b
(
a
,
b
)
t
\frac{b}{(a,b)}{t}
(a,b)bt,即
x
=
x
0
+
b
(
a
,
b
)
t
.
x=x_0+\frac{b}{(a,b)}t.
x=x0+(a,b)bt.(3-1)
代入(3)式(或带入(2-1)更好算),即得
y
=
y
0
−
a
(
a
,
b
)
t
y=y_0-\frac{a}{(a,b)}t
y=y0−(a,b)at(3-2)
特解
知识回顾
对整数
a
>
0
a>0
a>0,
b
>
0
b>0
b>0反复运用带余除法,可得下列等式
a
=
b
q
1
+
r
1
,
0
<
r
1
<
b
,
b
=
r
1
q
2
+
r
2
,
0
<
r
2
<
r
1
,
r
1
=
r
2
q
3
+
r
3
,
0
<
r
3
<
r
2
,
r
2
=
r
3
q
4
+
r
4
,
0
<
r
4
<
r
3
,
.
.
.
r
n
−
2
=
r
n
−
1
q
n
+
r
n
,
0
<
r
n
<
r
n
−
1
,
r
n
−
1
=
r
n
q
n
+
1
+
r
n
+
1
(
r
n
+
1
=
0
)
.
(1)
\begin{aligned} \\&a=bq_{1}+r_{1},&&0<r_{1}<b, \\&b=r_{1}q_{2}+r_{2},&&0<r_{2}<r_{1}, \\&r_{1}=r_{2}q_{3}+r_{3},&&0<r_{3}<r_{2}, \\&r_{2}=r_{3}q_{4}+r_{4},&&0<r_{4}<r_{3}, \\&... \\&r_{n-2}=r_{n-1}q_{n}+r_{n},&&0<r_{n}<r_{n-1}, \\&r_{n-1}=r_{n}q_{n+1}+r_{n+1}&&(r_{n+1}=0).\end{aligned}\tag{1}
a=bq1+r1,b=r1q2+r2,r1=r2q3+r3,r2=r3q4+r4,...rn−2=rn−1qn+rn,rn−1=rnqn+1+rn+10<r1<b,0<r2<r1,0<r3<r2,0<r4<r3,0<rn<rn−1,(rn+1=0).(1)
设 a a a, b b b是任意两个正整数,并进行了式(1)中的辗转相除法,
则有
Q
k
a
−
P
k
b
=
(
−
1
)
k
−
1
r
k
,
1
⩽
k
⩽
n
.
Q_{k}a-P_{k}b=(-1)^{k-1}r_{k},\quad 1\leqslant k\leqslant n.
Qka−Pkb=(−1)k−1rk,1⩽k⩽n. (2)
其中
{
P
0
=
1
,
P
1
=
q
1
,
P
k
=
q
k
P
k
−
1
+
P
k
−
2
,
Q
0
=
0
,
Q
1
=
1
,
Q
k
=
q
k
Q
k
−
1
+
Q
k
−
2
,
2
⩽
k
⩽
n
.
(3)
\left\{\begin{matrix} P_{0}=1,P_{1}=q_{1},P_{k}=q_kP_{k-1}+P_{k-2}, \\Q_{0}=0,Q_{1}=1,Q_{k}=q_kQ_{k-1}+Q_{k-2}, \end{matrix}\right. 2\leqslant k\leqslant n. \tag{3}
{P0=1,P1=q1,Pk=qkPk−1+Pk−2,Q0=0,Q1=1,Qk=qkQk−1+Qk−2,2⩽k⩽n.(3)
利用辗转相除法求
不妨设(1)式中的 a , b a, b a,b互素,不然去解方程 a ( a , b ) x + b ( a , b ) y = c ( a , b ) \frac{a}{(a,b)}x+\frac{b}{(a,b)}y=\frac{c}{(a,b)} (a,b)ax+(a,b)by=(a,b)c即可。
首先画出如下的框图,然后将由辗转相除法求出的 q 1 , q 2 , ⋯ , q n q_1, q_2, \cdots, q_n q1,q2,⋯,qn依次填在框图的第二行里。
接着在第三行填上 P 0 = 1 , P 1 = q 1 P_0 = 1, P_1 = q_1 P0=1,P1=q1,在第四行填上 Q 0 = 0 , Q 1 = 1 Q_0 = 0, Q_1 = 1 Q0=0,Q1=1。
由递推公式依次求出 P 2 , Q 2 , P 3 , Q 3 , ⋯ , P n , Q n P_2, Q_2, P_3, Q_3, \cdots, P_n, Q_n P2,Q2,P3,Q3,⋯,Pn,Qn。
由 辗转相除法基本定理 (取 k = n k=n k=n) 知
Q n a − P n b = ( − 1 ) n − 1 , Q_na - P_nb = (-1)^{n-1}, Qna−Pnb=(−1)n−1,
即
a [ ( − 1 ) n − 1 Q n c ] + b [ ( − 1 ) n P n c ] = c . a[(-1)^{n-1}Q_nc] + b[(-1)^nP_nc] = c. a[(−1)n−1Qnc]+b[(−1)nPnc]=c.
因此将以上求得的 P n , Q n P_n, Q_n Pn,Qn 代入下式
{
x
0
=
(
−
1
)
n
−
1
Q
n
c
,
y
0
=
(
−
1
)
n
P
n
c
,
\begin{cases} x_0 = (-1)^{n-1}Q_nc, \\ y_0 = (-1)^nP_nc, \end{cases}
{x0=(−1)n−1Qnc,y0=(−1)nPnc, (4)
即得(1)式的一个特解。
例
例:求不定方程 7 x + 4 y = 189 7x + 4y = 189 7x+4y=189 的一个特解。
( 7 , 4 ) = 1 (7,4)=1 (7,4)=1,可以用辗转相除法求一个特解
由(4)式可得方程的一个特解为
{ x 0 = − 100 , y 0 = 200. \begin{cases} x_0 = -100, \\ y_0 = 200. \end{cases} {x0=−100,y0=200.
由定理2 即得方程的全部解为
{ x = − 100 + 4 t , t = 0 , ± 1 , ± 2 , … y = 200 − 7 t . \begin{cases} x = -100 + 4t, & t = 0, \pm 1, \pm 2, \ldots \\ y = 200 - 7t. \end{cases} {x=−100+4t,y=200−7t.t=0,±1,±2,…
使用表达式凑出通解
也可以用如下方法求这个不定方程的一个特解:
由 7 x + 4 y = 100 7x + 4y = 100 7x+4y=100,推出 x = 100 − 4 y 7 = 14 + 2 − 4 y 7 x = \frac{100-4y}{7} = 14 + \frac{2-4y}{7} x=7100−4y=14+72−4y,
令 2 − 4 y 7 ∈ Z \frac{2-4y}{7}\in\mathbb{Z} 72−4y∈Z, y = 2 − 7 k 4 , k ∈ Z y=\frac{2-7k}{4},k\in\mathbb{Z} y=42−7k,k∈Z,显然 k k k是偶数才能使 y y y为整数,可取 k = 2 k=2 k=2,此时 y y y可取特值 y 0 = 4 y_0 = 4 y0=4,则 x 0 = 12 x_0 = 12 x0=12。
带入通解公式
{
x
=
x
0
+
b
(
a
,
b
)
t
,
y
=
y
0
−
a
(
a
,
b
)
t
,
t
=
0
,
±
1
,
±
2
,
⋯
\begin{cases} x = x_0 + \frac{b}{(a,b)}t, \\ y = y_0 - \frac{a}{(a,b)}t, \end{cases} \quad t=0,\pm 1, \pm 2, \cdots
{x=x0+(a,b)bt,y=y0−(a,b)at,t=0,±1,±2,⋯
不定方程的全部解为
{
x
=
12
+
4
t
,
t
=
0
,
±
1
,
±
2
,
…
y
=
4
−
7
t
,
\begin{cases} x = 12 + 4t, & t = 0, \pm 1, \pm 2, \ldots \\ y = 4 - 7t, \end{cases}
{x=12+4t,y=4−7t,t=0,±1,±2,…
虽然用以上两种方法得到的解的表达式形式上不同,但本质上是相同的。