※ 引述《sweetycool (tina)》之銘言：
: 已知 a,b,c,d皆為正整數
: a+b=c+d
: 1/13 = ac/bd
: 令s=a+b+c+d
: 求min s ?
: ans: 170
: 想問一下這題有辦法跑出來嗎?感恩
嗯，雖然數學版已經有人回說答案不對了，不過還是來跑一下
順便介紹幾個函式
首先是本題 (合併成一行之後的) 程式
Catch[
Do[
If[Length[#] != 0, Throw[Append[#, n]]] &[
Position[Outer[Times, #, #], 1/13] &[
Table[a/(n - a), {a, 1, n - 1}]
]
], {n, 1, 85}
]; {{}, -1}
]
以下由內而外解說
最內層是一個 Table，裡面列出分子分母的和是固定值 n 的所有分數
主要是看到條件裡的 ac/bd 可以看成 (a/b)*(c/d)
而這兩個分數由於有 a+b = c+d 的條件，分子分母的和是相同的
所以就先列出所有這些分數
外面一層是 Position 包 Outer，還有一個純函式的代值
先講純函式代值，這招我很常用在 one-liner 裡取代重覆的東西
基本模式是一個純函式 ...& 後面馬上用 [] 代值進去
用處是在當我要運算的式子裡有重覆的東西時可以只寫一次
這一層的函式是 Position[Outer[Times, #, #], 1/13]&，代的值是內層的 Table
其實這跟先把內層收進變數裡再重覆運用變數是一樣的
也就是這裡相當於
list = Table[...]
Position[Outer[Times, list, list], 1/13]
再來講 Outer，這個是廣義外積，簡單說明即是:
Outer[f, {a, b}, {x, y, z}] 會得到
{{f[a, x], f[a, y], f[a, z]}, {f[b, x], f[b, y], f[b, z]}}
於是在這裡 Outer[Times, list, list] 就會列出所有 list 裡任意兩個元素的乘積
這正好是我要尋找 1/13 的對象，因此就用 Position 去找出 1/13 在哪裡了
再外一層也是純函式代值，函式是
If[Length[#] != 0, Throw[Append[#, n]]]&
代的值是內層 Position 的結果
由於 Position 在找不到東西時回傳的是 {}
而找到東西時回傳的則是 {{位置一座標}, {位置二座標}, ...}
因此判斷它的長度是不是 0，如果不是的話就表示找到東西了
就在這個結果後面塞進 n 值之後 Throw 出來
這裡要講解的是 Throw / Catch 這一對函式
(jurian0101 在 #1JTHjnoX 有提過, 這裡再多講解一點)
這一對函式是屬於外包內結構，Catch 在外，Throw 在內
執行 Catch 時會先執行它的參數
如果這當中碰到 Throw 被執行了，則 Catch 的參數執行就到此中斷
Catch 的結果就是剛剛執行的 Throw 的參數
否則 (當中都沒有執行到 Throw) 則是參數執行結果
這例子中的 Catch 是在最外層
也就是說只要找到東西之後就會馬上停止計算，把結果丟出來給最外層
反之如果沒找到則繼續計算
再往外一層是一個 Do，這就只是很簡單的讓 n 從 1 開始試而已
上界填一個足夠大的數就好，這裡填 85 的原因是因為原 PO 貼的答案是 170 的關係
(不能填 Infinity，它會說這不能跑 XD)
當然萬一跑完範圍都找不到我們需要一個方式知道才行
因此分號後面就是 Do 跑完之後要執行的東西
如果真的跑到那裡的話會得到 {{}, -1} 這東西
然後因為這中間都沒碰到 Throw，所以外面的 Catch 得到的結果就會是這個
(其實不用擺也行, 因為 Do 的執行結果永遠是 Null，
只要跑完了卻沒有結果就是找不到)
最外層就是收答案用的 Catch 了
執行結果得到
{{1, 7}, {7, 1}, 14}
這表示在 n = 14 時，那個 Outer 的大陣列裡的 {1, 7} 及 {7, 1} 兩個位置是 1/13
但那個 Outer 陣列的元素是由內層的 Table 乘起來的
它的 {1, 7} 位置就是內層 Table 第一個元素跟第七個元素相乘的結果
而內層 Table 的第 a 個元素是分子為 a, 分子分母和為 n 的分數
因此從這個結果我們就能讀出 1/13 跟 7/7 的乘積是 1/13
也就是所求最小 s 的 (a, b, c, d) 是 (1, 13, 7, 7)