歷史考證
埃及分數埃及同中國一樣，也是世界上著名的文明古國。人們在考察古埃及歷史時注意到象阿基米德這樣的數學巨匠，居然也研究過埃及分數。本世紀一些最偉大的數學家也研究埃及分數，例如，沃爾夫數學獎得主，保羅-歐德斯，他提出了著名的猜想 4/n=1/x+1/y+1/z. 難倒了世界上第一流的數學家。當9個麵包要平均分給 10個人的時候，古埃及人不知道每個人可以取得 9/10，而是說每人1/3+1/4+1/5+1/12+1/30。真叫人難以想像，你連9/10都搞不清楚，怎么知道9/10=1/3+1/4+1/5+1/12+1/30。所以幾千年來，數學史家一直堅持認為，
古埃及人不會使用分數。 1858年，蘇格蘭考古學家萊登買到了一份古埃及草紙檔案，經過鑑定這是繁生於尼羅河泛濫形成的池塘和沼澤地里的草製成的紙，成文年代約在公元前1700年。
那么，古埃及的人們，是怎么算的呢？首先，把 2 個物品分成 4 個 1/2，先給每個人 1 個 1/2，剩下的 1 個1/2 再分成 3 等分，均分結果，每人分到 1/2 加 1/2 的 1/3，也就是 1/2 + 1/6 = 2/3。這份至今保存在大英博物館的“萊登”草紙，用很大的篇幅記載著將真分數分解成單分子分數，這種運算方式，遭到現代數學家們紛紛責難，認為埃及人之所以未能把算術和代數發展到較高水平，其分數運算之繁雜也是原因之一。
相關傳說
埃及金字塔是舉世聞名的，表明古埃及人具有高超的建築技巧和超凡的智力，難道最簡單的現代分數也不懂？金字塔所蘊含的難道是一篇粗劣的作品？
現代數學已經發展到十分抽象和複雜的程度，而埃及分數卻是這樣粗糙，在人們的記憶里早該煙消雲散了，然而，它產生的問題直到今天仍然引起人們的重視。
四川大學已故老校長柯召寫道：“埃及分數所產生的問題有的已成為至今尚未解決的難題和猜想，他們難住了許多當代數學家”。柯召本人至死都沒有能夠證明這個猜想。
一個古老的傳說是：
老人彌留之際，將家中11匹馬分給3個兒子，老大1/2，老二1/4，老三1/6。二分之一是5匹半馬，總不能把馬殺了吧，正在無奈之際，鄰居把自己家的馬牽來，老大二分之一，牽走了6匹；老二四分之一，牽走了3匹；老三六分之一，牽走了2匹。一共11匹，分完後，鄰居把自己的馬牽了回去。即11/12=1/2+1/4+1/6。
奇妙的埃及分數終於調動自己的潛在難度擊敗了敢於輕視他們的人們。並且給與嘲笑他的人以難堪的回答。
現代探索
相關發現
兩千多年後的數學家終於發現：2/n=1/[(n+1)/2]+1/[(n+1)n/2]； 1/n=1/(n+1)+1/[n(n+1)]；1=1/2+1/3+1/6。此時才大夢初醒。埃及分數以旺盛的生命力屹立在世界數壇，使三千年後的數學家也自嘆弗如。例如，分馬問題，能否設計出（n-1)/n=1/x+1/y+1/z .。經過2000多年的努力，終於揭開其中的奧秘：有6種可能，共7種分法。7/8=1/2+1/4+1/8；11/12=1/2+1/4+1/6=1/2+1/3+1/12；17/18=1/2+1/3+1/9；19/20=1/2+1/4+1/5；23/24=1/2+1/3+1/8；41/42=1/2+1/3+1/7。原先人們以為，這樣的情況大概有無窮多個，可是，繼續追擊卻一無所獲，真是難以預料。黑龍江的關春河發現共有43種情況。這是正確的。
求解過程
當限定分母為奇數時，把“1”分解為埃及分數，項數限定為9項，共有5組解：
1=1/3+1/5+1/7+1/9+1/11+1/15+1/35+1/45+1/231。
1=1/3+1/5+1/7+1/9+1/11+1/15+1/21+1/135+1/10395。
1=1/3+1/5+1/7+1/9+1/11+1/15+1/21+1/165+1/693。
1=1/3+1/5+1/7+1/9+1/11+1/15+1/21+1/231+1/315。
1=1/3+1/5+1/7+1/9+1/11+1/15+1/33+1/45+1/385。
以上5組解是在1976年才找到。限定為11項時，發現了1組解 最小分母是105。若大於105則有很多的解。
1/n型分數還可以表示成為級數分解式：
1/n=1/(n+1)+1/(n+1)^2+1/(n+1)^3+1/(n+1)^4+....+1/(n+1)^k+1/n(n+1)^k.
埃及分數成為不定方程中一顆耀眼的明珠。
埃及分數最著名的猜想是Erods猜想：1950年Erods猜想，對於n〉1的正整數，
總有：
4/n=1/x+1/y+1/z. （1）
其中，x，y，z。都是正整數。
Stralss進一步猜想，當n≥2時，方程的解x，y，z滿足x≠y，y≠z，z≠x。x〈y〈z。
1963年柯召，孫奇，張先覺證明了Erods猜想stralss猜想等價。幾年後yamanot又把結果發展到10的7次方。以後一些數學家又把結果推向前去，始終未獲根本解決。對於4/n=1/x+1/y+1/z，只需要考慮n=p為素數的情況，因為若（1）式成立，則對於任何整數m，m<1，
4/pm=1/xm+1/ym+1/zm，（2）
也成立。
一切奇素數都可以表示為4R+1與4R+3型。對於p=4R+3型，（參見《單位分數》人民教育出版社1962年）：（1）式是顯然的。
2002年王曉明提出：
如果設X=AB,Y=AC,Z=ABCP,
即：
4/P=1/AB+1/AC+1/ABCP.（3）
對於p=4R+3型，（3）式是顯然的。
因為這時A=(p+1)/4 ，B=1。C=P+1.。
即：
4 /P = { 1/ [(P+1)/4] } + { 1 / [(P+1)(p+1)/4] } + { 1/ [p(p+1)(p+1)/4] }。 (4)
例如：4/7=1/2+1/16+1/112
對於p=4R+1 型的素數，把（3)式整理成 ：
4ABC=PC+PB+1 (5)
A = (PC+PB+1)/4BC (6)
在（6）式中，若要 B|(PC+PB+1)，需使得B|(PC+1)，設PC+1=TB；若要C|(PC+PB+1)，需使得C|(PB+1)，設PB+1=SC；對於P=4R+1形，若要4|p(C+B)+1],需C+B=4K-1，對於P=4R+3形，若要4|[P（C+B）+1]，需C+B=4K+1。於是，形成一個二元一次不定方程組：
-PC+TB=1 (7)
SC+(-P)B=1 (8)
例如p=17時，A=3，B=2，C=5，T=43，S=7，k=2 。
4 /17=[1/(2×3)]+[1/(3×5)]+[1/(3×2×5×17 )]
即4/17=1/6+1 /15+1/510.
等價於下面的式子:
（-17）×5+43×2=1
7×5+（-17）×2=1
注意：P=(4ABC-1)/(B+C). (9)
由於4ABC-1是4R+3型，所以，當P=4R+1型時，B+C=4K-1型；P=4R+3,B+C=4K+1型。.
因為對於二元一次不定方程組，我們有得是辦法。根據《代數學辭典》上海教育出版社1985年（376頁）：“
方程組:ax+by=c
a'x+b'y=c'
公共解（整數解）x,y的充分必要條件是(ab'-a'b)不等於0，並且 (ab'-a'b) | (bc'-b'c) 和 (ab'-a'b) | (ca'-c'a)。”
我們把（7）（8）式的C與B當成上面的x,y. 在（7）式中，只要（P,T）=1；就有無窮多組B和C整數解;在（7）中，只要（P,S）=1，就有B和C的整數解。根據已知的定理（柯召，孫奇《談談不定方程》）13 至17頁，聯立二元一次不定方程，就知道（7）（8）式必然有公共整數解（用到矩陣，單位模變換等知識）。即ST-P×P≠0，(ST-P×P) | (P+T)； (ST-P×P) | (P+S)。為什麼說是必然有解，只要有一個素數有解，其它素數必然有解。在中國象棋中，“馬”從起點可以跳到所有的點，那么，馬在任何一個點就可以跳到任何點。因為馬可以從任何一個點退回的起點。
下面是一些p值的解：
--p---|---A---|---B---|----C-----|------T-----|------S-------|-------K-----|
------------------------------------------------------------------------------|
--5---|--2----|---1----|---2------|-----11-----|----3---------|------1------|
-29--|---2----|---4----|---39----|----283----|----3---------|------11-----|
-37--|---2----|---5----|--62-----|---459-----|----3---------|-------17----|
-53--|---2----|---7----|--124----|---939-----|----3--------|-------33----|
-61--|---2----|---8----|--163----|---1243----|----3--------|-------43----|
-173-|--2----|----22--|--1269---|--9979----|----3--------|------323----|
-----------------------------------------------------------------------------------------
以上是P=4R+1，R為奇數時的解，此時，A=2；S=3。
---------------------------------------------------------------------------------
-17--|--3-----|---2----|-----5------|----43-----|-----7--------|-----2-------|
-41--|--12----|---1----|----6-------|---247----|----7---------|-----2-------|
-41--|--6------|---3----|----4-------|---55-----|-----31-------|-----2-------|
-73--|---10----|---2----|---21------|----767--|-----7---------|-----6-------|
- 97--|---17---|---2----|----5-------|---243---|----39--------|-----2-------|
-113-|--5------|---6----|---97------|--1827---|----7---------|----26-------|
-409-|--59-----|---2---|----13------|--2659---|----63-------|----4--------|
-409-|--22-----|---5---|-----66-----|--5399---|----31-------|-----18-----|
-409-|--11-----|---11--|----60-----|---2231--|----75-------|-----18-----|
---------------------------------------------------------------------------------------
以上是p=4R+1，R是偶數時的解。
41有兩組解；409有三組解。就是說4/41=1/（12×1）+1/（12×6）+1/（12×1×6×41）=1/12+1/72+1/2952
4/41=1/（6×3）+1/（6×4）+1/（6×3×4×41）=1/18+1/24+1/2952。
-41×6+247×1=1
7×6+（-41）×1=1
和第二組解；
-41×4+55×3=1
31×4+（-41×3）=1
（2）式是對於所有的p值都有解，但不是全部解。(例如，4/41有7組解，而（2）式只求證4/p=1/AB+1/AC+1/ABCP
的形式解。請注意普遍解與全部解的區別。
在七十年代，人們又提出了5/P的情況，所有的素數P都可以表示成5R+1；5R+2；5R+3；5R+4形。
對於P= 5R+4形，5/（5R+4）=1/（R+1）+1/[（5R+4）（R+1）]
其中任何一個：1/N=1/（N+1）+1/[N（N+1）]。
例如，5/9=1/2+1/18，而1/2=1/3+1/6；或者1/18=1/19+1/（18×19）。
對於P=5R+3形，5/（5R+3）=1/（R+1）+2/[（5R+3）（R+1）]
其中任何一個：2/N=1/[（N+1）/2]+1/[N（N+1）/2]
例如，5/13=1/3+2/39，而2/39=1/[（39+1）/2]+1/[39×（39+1）/2]。
對於P=5R+2形，5/（5R+2）=1/（R+1）+3/[（5R+2）（R+1）]
R必然是奇數，（R+1）必然是偶數。
而：3/[（5R+2）（R+1）]=1/[（5R+2）（R+1）]+1/[（5R+2）（R+1）/2]
例如，5/37=1/8+3/（37×8）；而3/（37×8）=1/（37×8）+1/（37×4）。
對於P=5R+1形，
設5/P=1/AB+1/AC+1/ABCP （8）。
5ABC=PC+PB+1 （9）
A=（PC+PB+1）/5BC （10）。
同樣可以整理成（6）（7）式，同樣有解。B+C=5K-1形。
下面是一些p=5R+1形的素數的解。
5/11=1/3+1/9+1/99，A=3，B=1，C=3，T=34，S=4；
5/31=1/7+1/56+1/1736，A=7，B=1，C=8，T=248，S=4；
5/41=1/9+1/93+1/11439，A=3，B=3，C=31，T=424，S=4；
5/61=1/14+1/95+1/81130，A=1，B=14，C=95，T=414，S=9；
5/71=1/15+1/267+1/94785，A=3，B=5，C=89，T=1264，S=4；
5/101=1/21+1/531+1/375417，A=3，B=7，C=177，T=2554，S=4；
5/131=1/27+1/885+1/1043415，A=3，B=9，C=295，T=4294，S=4；
方法同4/P一樣。請讀者自己完成。
為什麼（6）（7）式可以必然有解？
兩聯二元一次不定方程：
a1x+b1y=1
a2x+b2y=1.
有解的充分條件是(a1b2-a2b1)|(a1-a2);(a1b2-a2b1)|(b2-b1).
我們考察一聯二元一次不定方程：
ax+by=1.(14)
根據已知定理，只要（a,b)=1,(14)式就有整數x,y的解。並且是有無窮多組解。
例如，5x-2y=1.
x; y
-----------------
1, 2;
3, 7;
5, 12;
7, 17;
9, 22;
11，27；
13，32；
15，37；
17, 42;
19, 47;
...........
換句話說，（14）式中，x與y也互素。這就是聯立方程組有公共解的基礎。我們把a,b與x,y互換，
以上例為例子，5x-2y=1換成5a-2b=1,x=5,y=2.
3x-7y=1
17x-42y=1
形成二聯二元一次不定方程。
5x-12y=1
19x-47y=1
7x-17y=1
形成三聯二元一次不定方程。
（4）式可以表示成一個素數的式子：
p=(4ABC-1)/(C+B)。例如p=41時，41=（4x6x3x4-1）/（4+3）；41=（5x3x3x31-1）/（31+3）；
41=(6x1x8x47-1)/(8+47)；41=(7x1x7x36-1)/(7+36)；41=(8x6x1x6-1)/(1+6)；41=(9x1x6x19-1)/(6+19)；
41=(10x1x6x13-1)/(6+13)；41=(11x1x4x55-1)/(4+55);；41=(12x4x1x6-1)/(1+6);；41=(13x1x4x15-1)/(4+15)；
41=(14x1x3x124-1)/(3+124).。到n=15就沒有了：41= (nABC-1)/(B+C)都有效。
人們於是問：是否一切n
p=(nABC-1)/(B+C).
有三個未知變數的素數公式，可以求得一切素數：
P=(4ABC-1)/(B+C).(15)。
（15）式對於一切p=4r+1形式的素數都可以。
例如，17.：17=（4x3x2x5-1)/(2+5)。
（15）式對於一切p=4r+3形式的素數，A=(P+1)/4，,B=1,，C=P+1。例如11=（4x3x1x12-1)/(1+12).。
對於合數n=4r+3形式。n=（4xBXC－1）/（B+C）．
例如51=（4x13x664－1）/（13+664）。B=（P+1）/4，C=n（n+1）/4+1．
算法解決
題目
埃及分數
在古埃及，人們使用單位分數的和（形如 1/a 的，a 是自然數）表示一切有理數。 如：
2/3=1/2+1/6，但不允許 2/3=1/3+1/3，因為加數中有相同的。 對於一個分數 a/b，表示方
法有很多種，但是哪種最好呢？
首先，加數少的比加數多的好，其次，加數個數相同的，最小的分數越大越好。 如：
最好的是最後一種，因為 1/18 比 1/180、1/45、1/30、1/180 都大。
【輸入檔案】
給出兩個正整數 a、b（0 < a < b < 1000），編程計算對於分數 a/b 最好的表達方式。
【輸出檔案】
若干個數，自小到大排列，依次是單位分數的分母。
【樣例輸入】
19 45
【樣例輸出】
5 6 18
pascal 代碼（疊代加深，有更好的算法）
var
temp,ans:array[1..20]of longint;
flag:boolean;
aim:extended;
a,b,te,maxd:longint;
function gcd(a,b:longint):int64;
begin
if b=0 then exit(a);
exit(gcd(b,a mod b));
end;
function lcm(a,b:longint):int64;
var
t:longint;
begin
if a
begin
t:=a;a:=b;b:=t;
end;
exit(a div gcd(a,b)*b);
end;
procedure sum(var s1,s2:int64;m:longint);
var
t:int64;
begin
t:=lcm(s2,m);
if t>100000 then
begin
s1:=10000;
s2:=1;
exit;
end;
s1:=t div s2*s1+t div m;
s2:=t;
t:=gcd(s1,s2);
s1:=s1 div t;
s2:=s2 div t;
end;
procedure fc(s1,s2:longint);
var
i:longint;
begin
if (s1=a)and(s2=b) then
begin
flag:=true;
if ans[maxd]>temp[maxd] then
ans:=temp;
end;
end;
procedure dfs(s:extended;s1,s2,m,i:int64);
var
j:longint;
up,down,t1,t2:int64;
begin
if s1/s2>aim then exit;
if i>maxd then begin fc(s1,s2); exit;end;
up:=trunc(1/(aim-s));
if up
down:=trunc((maxd-i+1)/(aim-s))+100;
for j:=up to down do
begin
temp[i]:=j;
t1:=s1;t2:=s2;
sum(t1,t2,j);
dfs(s+1/j,t1,t2,j+1,i+1);
end;
end;
procedure print;
var
i:longint;
begin
for i:=1 to maxd-1 do
write(ans[i],' ');
writeln(ans[maxd]);
end;
begin
fillchar(ans,sizeof(ans),$3f);
read(a,b);
te:=gcd(a,b);
a:=a div te;
b:=b div te;
aim:=a/b;
maxd:=1;
flag:=false;
while not flag do
begin
inc(maxd);
dfs(0,0,1,2,1);
end;
print;
end.
還有更基礎的解法，初學者可用：
var a,b,c:integer;
begin
write('a,b=");readln(a,b);
write(a,"/',b,'=");
repeat
c:=b div a +1;
a:=a*c-b;
b:=b*c;
write("1/',c);
write('+');
until (a=1) or (b mod a=0);
if (b mod a=0)and(a<>1)
then write('1/',b div a);
if a=1 then write('1/',b);
readln;
end.
未來發展
實際上這個問題還遠遠沒有解決。但是已經給出了前進的方向。
.埃及分數，一個曾被人瞧不起的，古老的課題，它隱含了何等豐富的內容，許多新奇的謎等待人們去揭開。
