雞兔同籠和丟番圖

Review:

Last time, we went over the basic theorems of Linear Diophantine Equation. This time, we will focus on how to solve the problem.

丟番圖方程的用處是可以解 "雞兔同籠".

假設(shè)一些小兔子和小雞 在一個籠里一共有100只腿液样,那可以有多少兔子眶明,和多少小雞呢阴颖?
image

2*c + 4*r =100
c: # of chickens
r: # of rabbits</pre>

Solution:

Steps for Slove Linear Diophantine Equation: mx+ ny= w

1. Check the existence of Solution:
gcd (2,4) =2 (2x+4y=100)
2| 100 => the equation has at least one solution.

2. Find the solution of x, y such that 2*x+ 4*y= gcd(m,n)=2
x=-1 , y=1
2*(-1) + 4*1=2

3. Multiply w/gcd(m,n) to the equation
100/2=50
50*2*(-1)+ 50*4*1 = -100 +200= 100

4. Select solutions within the constraints
c=-50, r=200 would be a theoretical solution.
but we are talking about rabbits and chicken, so we need positive numbers.
We can do it by trial or follow the formula, we have a set of solution:
Paired Solution:
( x_{initial}+ (m/gcd(n,m))*z , y_{initial} - (n/gcd(n,m)*z))

Answer:

(2 chicken, 24 rabbits), (4 chicken,23 rabbits), (6 chicken , 22 rabbits )....
Now we have figured out the set of chicken and rabbits will make 100 legs in total: (2 ??, 24 ??), (4 ??, 23 ??) , (6 ??, 22 ??) ....

Python Solution:

Now we can try a more technical solution:

>#linear diophantine
def isolve(a,b,c):
    q,r=divmod(a,b)
    if r==0:
        return([0,c/b])
    else:
        sol=isolve(b,r,c)
        u=sol[0]
        v=sol[1]
        return([v,u-q*v])

Out[14]: [50.0, 0.0]
#Note: this is one solution, we need to find the rest.

# Alternative Solution
from sympy.solvers.diophantine import diop_linear
diop_linear(2*a + 4*b-100)

Out[13]: {(2*t_0 + 50, -t_0)}

Let t_0= -24, then we have
(2, 24) ....

As we can see, python 只給我們提供一個初始解顷蟀,和公式,我們要自己帶數(shù)字.

Summary + Story:

Although Python is very fast with the initial solutions of Linear Diophantine equation, it is always a good idea we understand the Key Concept: GCD behind the solutions.

I still remember how much fun just to play around the numbers.
我小時候最喜歡的物理老師胡寧說他最討嫌"雞兔同籠" 的題目, 因為哪里會有人蠢到真的把它們放到一個籠子里? 十多年過去了, Renee 告訴我 Highland Park, NJ 的動物園里面真的把兔子和小雞放在一個籠子里绍赛,我想我應(yīng)該去照個照片蔓纠,希望有機會再見到胡老師 ! :)

Puzzle :

Diophantine Equation came from Diophantus, a Hellenistic mathematician from Egypt, He had a lot of great work done in arithmetics. It was said his age was engraved as a puzzle, could you figure out his age?

'Here lies Diophantus,' the wonder behold.
Through art algebraic, the stone tells how old:
'God gave him his boyhood one-sixth of his life,
One twelfth more as a youth while whiskers grew rife;
And then yet one-seventh ere marriage began;
In five years there came a bouncing new son.
Alas, the dear child of master and sage After attaining half the measure of his father's life chill fate took him.
After consoling his fate by the science of numbers for four years,
he ended his life.'

Happy Studying! ??

Reference:

https://brilliant.org/wiki/linear-diophantine-equations-one-equation/

https://www.math.utah.edu/~carlson/hsp2004/PythonShortCourse.pdf

https://docs.sympy.org/latest/modules/solvers/diophantine.html

https://en.wikipedia.org/wiki/Diophantus

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末吗蚌,一起剝皮案震驚了整個濱河市腿倚,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌褪测,老刑警劉巖猴誊,帶你破解...
    沈念sama閱讀 206,482評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件潦刃,死亡現(xiàn)場離奇詭異,居然都是意外死亡懈叹,警方通過查閱死者的電腦和手機乖杠,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,377評論 2 382
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來澄成,“玉大人胧洒,你說我怎么就攤上這事∧矗” “怎么了卫漫?”我有些...
    開封第一講書人閱讀 152,762評論 0 342
  • 文/不壞的土叔 我叫張陵,是天一觀的道長肾砂。 經(jīng)常有香客問我列赎,道長,這世上最難降的妖魔是什么镐确? 我笑而不...
    開封第一講書人閱讀 55,273評論 1 279
  • 正文 為了忘掉前任包吝,我火速辦了婚禮,結(jié)果婚禮上源葫,老公的妹妹穿的比我還像新娘诗越。我一直安慰自己,他們只是感情好息堂,可當我...
    茶點故事閱讀 64,289評論 5 373
  • 文/花漫 我一把揭開白布嚷狞。 她就那樣靜靜地躺著,像睡著了一般荣堰。 火紅的嫁衣襯著肌膚如雪床未。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 49,046評論 1 285
  • 那天持隧,我揣著相機與錄音即硼,去河邊找鬼。 笑死屡拨,一個胖子當著我的面吹牛,可吹牛的內(nèi)容都是我干的褥实。 我是一名探鬼主播呀狼,決...
    沈念sama閱讀 38,351評論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼损离!你這毒婦竟也來了哥艇?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 36,988評論 0 259
  • 序言:老撾萬榮一對情侶失蹤僻澎,失蹤者是張志新(化名)和其女友劉穎貌踏,沒想到半個月后十饥,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,476評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡祖乳,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 35,948評論 2 324
  • 正文 我和宋清朗相戀三年逗堵,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片眷昆。...
    茶點故事閱讀 38,064評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡蜒秤,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出亚斋,到底是詐尸還是另有隱情作媚,我是刑警寧澤,帶...
    沈念sama閱讀 33,712評論 4 323
  • 正文 年R本政府宣布帅刊,位于F島的核電站纸泡,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏赖瞒。R本人自食惡果不足惜弟灼,卻給世界環(huán)境...
    茶點故事閱讀 39,261評論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望冒黑。 院中可真熱鬧田绑,春花似錦、人聲如沸抡爹。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,264評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽冬竟。三九已至欧穴,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間泵殴,已是汗流浹背涮帘。 一陣腳步聲響...
    開封第一講書人閱讀 31,486評論 1 262
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留笑诅,地道東北人调缨。 一個月前我還...
    沈念sama閱讀 45,511評論 2 354
  • 正文 我出身青樓,卻偏偏與公主長得像吆你,于是被迫代替她去往敵國和親弦叶。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 42,802評論 2 345

推薦閱讀更多精彩內(nèi)容