Loading [MathJax]/jax/output/HTML-CSS/jax.js

2014年1月3日 星期五

[數論] 不定方程

參考資料:
不定方程,凡異出版社,1996年1月初版,ISBN:9576942225

二元一次不定方程
二元一次不定方程式 ax+by=c, 其中 a,b,c 都是已知整數, 且 a,b 不全為 0.
 ax+by=c 有整數解 (a,b)c
實際求解的方法,有「誤試法(Trial and Error)」、「歐基里德算法(輾轉相除法)」,有時候解題方法不只一種。

誤試法(Trial and Error):當 |a|,|b| 不大時可以試試看, e.g. 7x+15y=1989.
歐基里德算法:e.g. 5767x+4453y=1679.
(a,b)=1, x=x0, y=y0 是方程式 ax+by=c 的一組解(特解), 則全部的整數解(通解)為 {x=x0+bty=y0at, t 是任何整數.
<pf> 設 (x,y) 是方程式的一組解, 則 ax+by=c. 再由 ax0+by0=c, 兩式相減得到 a(x0x)+b(y0y)=0, 由此知 ab(y0y). 但 (a,b)=1, 所以 a(y0y). 可寫成 y0y=at, 其中 t 為整數, 即 y=y0at. 把 y=y0at 代回 a(x0x)+b(y0y)=0, 得到 x=x0+bt.


例10 設 a,b,n 都是正整數, (a,b)=1. 當 n>abab 時, 方程 ax+by=n 有非負整數解; 當 n=abab 時, 則無非負整數解. 換句話說, 凡大於 abab 的整數必可表示為 ax+by (x0,y0) 的形式, 但 abab 不能表示為此種形式.

<pf1> 因為 (a,b)=1, 所以必有整數 x,y 使得 ax+by=n 成立. 我們總可以選擇 0x<b(思考一下,why?).
n>abab 時, by=nax>(abab)axababa(b1)=b. 因此 y>1, 證明了 y 也是非負整數.

n=abab 時, 則 ax+by=n=abab, x=ababbya=b1ba(y+1). 若要使 x0, 因為(a,b)=1, 所以 (y+1) 須為 a 的倍數, 且 ba(y+1)<b, 得到 (y+1)<a.
且假設 y 也為非負整數, 所以 (y+1)>0. 整理得到 0<(y+1)<a, 亦即 a0<(y+1)<a1, 顯然 (y+1) 不是 a 的倍數, 與假設矛盾. 故 x,y 無非負整數解.

<pf2>
首先證明方程 ax+by=nn>abab 時有非負整數解.
方程式的一般解為 {x=x0+bty=y0at, tZ. 取 適當的 t 使得 0y=y0ata1. 對於此 t, 有 ax=a(x0+bt)=nb(y0at)>ababb(a1)=a. 而 a>0, 故 x=x0+bt>1, 亦即 x0.

其次證明 n=abab 時, ax+by=n 無非負整數解.
假設方程有解, x0,y0. 由 ax+by=n=ababa(x+1)+b(y+1)=ab. 因為 (a,b)=1, 所以 a(y+1),b(x+1). 於是有 y+1a,x+1b. 所以 ab=a(x+1)+b(y+1)ab+ab=2ab, 但 a>0,b>0, 所以上式矛盾.

例4 若 a>0, b>0, 且 (a,b)=1, 試證:
1. 方程 ax+by=n 的非負整數解之個數為 [nab][nab]+1.
2. 在 0nabab 中, 恰有 (a+1)(b+1)2 個整數 n 能(或不能)表示成 ax+by 的形式, 這裡的 x,y 都是非負整數.


三元一次不定方程
例13 求不定方程 6x+20y15z=23 的全部整數解.(hint: 從係數小的著手)


一次不定方程組
消元法
例2 《張丘建算經》的百雞問題. 公雞5元一隻, 母雞3元一隻, 小雞1元一隻. 現在用100元買100隻雞, 問公雞, 母雞, 小雞應各買多少隻?
解 設公雞, 母雞, 小雞的隻數分別為 x,y,z. 由題意, 得 {x+y+z=1005z+3y+z3=100. 相應的 (x,y,z)(0,25,75), (4,18,78), (8,11,81)(12,4,84).


沒有留言:

張貼留言