古代數學巨著《九章算數》中有這么一道題叫“五家共井,甲二綆(汲水用的井繩)不足,如(接上)乙一綆;乙三綆不足,如丙一綆;
丙四綆不足,如丁一綆;丁五綆不足,如戊一綆;戊六綆不足,如甲一綆,皆及。
意思就是說五家人共用一口井,甲家的繩子用兩條不夠,還要再用乙家的繩子一條才能打到井水;乙家的繩子用三條不夠,還要再用丙家的繩子
一條才能打到井水;丙家的繩子用四條不夠,還要再用丁家的繩子一條才能打到井水;丁家的繩子用五條不夠,還要再用戊家的繩子一條才能打
到井水;戊家的繩子用六條不夠,還要再用甲家的繩子一條才能打到井水。
最后問:井有多深?每家的繩子各有多長?
分析:同樣這套題也是屬於不定方程,拿這個題目的目地就是讓大家能夠在不定方程組這種范疇問題上做到“舉一反三”,根據題意
我們設井深為h,各家分別為a,b,c,d,e,則可以列出如下方程組:
2a+b=h ①
3b+c=h ②
4c+d=h ③
5d+e=h ④
6e+a=h ⑤
首先我們看下普通青年的想法,他們的想法是找a,b,c,d,e之間的對應關系。
依次將②代入①,③代入②,④代入③,⑤代入④可得如下方程組:
a=b+c/2
b=c+d/3
c=d+e/4
d=e+a/5
從計算機的角度來說,我不希望有小數的出現,所以我可推斷: c一定是2的倍數,d一定是3的倍數,e一定是4的倍數,a一定是5的倍數,根據這種關系我們
就可以有如下代碼:
1 namespace Test 2 { 3 class Program 4 { 5 static void Main(string[] args) 6 { 7 int a, b, c, d, e, h; 8 9 a = b = c = d = e = h = 0; 10 11 bool flag = true; 12 13 while (flag) 14 { 15 //4的倍數 16 e += 4; 17 18 a = 0; 19 20 while (flag) 21 { 22 //5的倍數 23 a += 5; 24 25 d = e + a / 5; 26 27 c = d + e / 4; 28 29 if (c % 2 != 0) 30 continue; 31 32 if (d % 3 != 0) 33 continue; 34 35 b = c + d / 3; 36 37 if (b + c / 2 < a) 38 break; 39 40 if (b + c / 2 == a) 41 flag = false; 42 } 43 } 44 45 h = 2 * a + b; 46 47 Console.WriteLine("a={0},b={1},c={2},d={3},e={4} ------h={5}\n", a, b, c, d, e, h); 48 49 Console.Read(); 50 } 51 } 52 }
同樣我們的時間復雜度是O(N2),急需優化。
我們再來看看文藝青年的想法,他們的想法是找a,b,c,d,e中的某個數與h的對應關系。
比如我就找c與h的對應關系,上面的①②③④⑤可寫成如下方程組:
b=h-2a ⑥
c=h-3b ⑦
d=h-4c ⑧
e=h-5d ⑨
a=h-6e ⑩
將⑥,⑧,⑨,⑩分別代入⑦,一陣痙攣后可知:
c=(148/721)h
上面的公式也就表明了c和h的比例關系,我們令 h=721k,則 c=148k,將其代入⑥,⑦,⑧,⑨,⑩可得如下方程組
a=265k
b=191k
c=148k
d=129k
e=76k
x=721k
又因為k>0,所以題目有無數個解。這里我就取0<k<5,否則繩子已經到達極限了,需要用蛟龍號去深潛了。
1 namespace Test 2 { 3 class Program 4 { 5 static void Main(string[] args) 6 { 7 for (int k = 1; k < 5; k++) 8 { 9 int h = 721 * k; 10 11 int a = 265 * k; 12 13 int b = 191 * k; 14 15 int c = 148 * k; 16 17 int d = 129 * k; 18 19 int e = 76 * k; 20 21 Console.WriteLine("a={0},b={1},c={2},d={3},e={4} ------h={5}\n", a, b, c, d, e, h); 22 } 23 24 Console.Read(); 25 } 26 } 27 }
相信大家以后遇到類似的問題,應該會胸有成竹了。