故事虚构,是从一个真的游戏再综合新闻组的内容而来.
缘起
缘起
这是一个晴朗的星期六下午,你悠闲地在网上浏览.忽然间你看到一个留言板上的小游戏.它很简单,
问题是:
把五个数字56789,放到[][][] * [][],令结果最大.
你最先对自己说: "这有什么难,把最大的放到最大位数那里就行了."你再心算了一下,也许不对.每个结果要看其他位置上放了什么数才行.你开始觉得有些兴趣了,反正你正在学一种好玩的编程语言,何不练习一下呢?
於是你开出你心爱的Python,开始思考: "其实我要的是一个程式,我给它各种数字的组合,然后它自动帮我找出最大的一个.如果我传入1,1,1,1,1和1,1,1,1,2,它会知道要算111 * 11和111 * 12,求出较大的是111 * 12并输出这个组合以及其乘积.这个程式并不难嘛."
切换行号显示
1# calc.py
2defcalc(seq):
3maximum =0
4max_item = []
5foriinseq:
6product = (i[0]*100+ i[1]*10+ i[2]) * (i[3]*10+ i[4])
7ifproduct > maximum:
8maximum = product
9max_item = i
10elifproduct == maximum:
11max_item +=','+i
12returnmax_item, maximum
13
14seq = [ [5,6,7,8,9], [5,6,7,9,8] ]
15max_item, maximum =calc(seq)
16print"Maximum
at", max_item,",product", maximum
你试了一下,
$python calc.py
Maximum at [5, 6, 7, 9, 8] ,product 90160
没问题.现在你只要给出所有的排列即可.你打了几行,觉得[5,6,8,7,9]这样打太辛苦了,而且用i[0]*100 + i[1]*10 ...的方法好像太难看了,因此你有必要做一次修改.好,用字串吧. "56789",这样输入时较方便,而且int("567") *
int("89")要好看得多,也应该快些.另外你也把程式改得更短,看起来像是个有经验的人所写.
切换行号显示
1# calc.py
2defcalc(seq, where):
3maximum, max_item =0, []
4foriinseq:
5product =int(i[:where]) *int(i[where:])
6ifproduct > maximum:
7maximum, max_item = product, i
8elifproduct == maximum:
9max_item +=','+ i
10print"Maximum at", max_item,",product", maximum
11
12if__name__ =="__main__":
13seq = ["56789","56798"]
14where =3
15calc(seq,where)
嗯,好些了.那句ifname== "main"是为了将来你把calc.py做为模组来用时而设的.在别的程式中用import calc的话那几行就不会运行了.
现在你可以随自己的意来解更普遍的问题,比如123放在[]*[][],或是1234567放在[][][][]*[][][]这样的问法.现在你开始输入排列了. "56789",
"56798", "56879", "56897", ..........没多久你又觉得自己太愚蠢了,为什么不叫电脑程式自动产生这些无聊的排列呢? 56789一共有5!也就是120种排列方法呢!如果你想算123456789的话,用手输入可能要用去一生的时间!!
於是你开始想如何产生排列的算法了.用循环就可以了,不过循环会产生重覆的组合,譬如555*55.但我们可以加些条件式进去把重覆项拿走.於是你有了第一个程式解.
切换行号显示
1#permute1.py
2defpermute(seq):
3result = []
4forainseq:
5forbinseq:
6forcinseq:
7fordinseq:
8foreinseq:
9ifa!=banda!=canda!=danda!=eand\
10b!=candb!=dandb!=eand\
11c!=dandc!=eandd!=e:
12result.append(''.join([a,b,c,d,e]))
13returnresult
14
15seq =list("56789")
16where =3
17thelist =permute(seq)
18importcalc
19calc.calc(thelist,where)
你小心地记著用''.join()的方法产生字串要比用a+b+c+d快,同时也认为用import calc的方式会让你更容易为不同的地方做些速度的微调.你开始运行程式了:
%python permute1.py
Maxmum at 87596 ,product 84000
你成功了.啊哈,你认为可以贴到留言板上去领赏了.经过一些考虑后,你觉得还是要做到更普遍的功能,就是让用户输入排列多少个数字,怎样分割.研究了一下程式,你觉得用循环好像无法达到要求,因为你事前并不知道要排多少个数字,因此你不知道要写多少个循环才够.面对这种情况,你好像只能用递归的方法了.
你知道如何求得,例如, 5个数字的排列:先挑一个数,有五种选择;当选定了一个数之后挑第二个数时只剩下四个选择,依此类推.因此五个数共有5*4*3*2*1共120个排列.当你面对"56789"这五个数的排列问题时,你先挑出一个数,例如6,那剩下的便是一个四个数字的排列问题了.就是说, 56789的排列可以简化(或是简单复杂化:p)成字头为5的所有排列加上字头为6的所有排列加字头为7的所有排列加字头为8的所有排列再加字头为9的所有排列.想通了这点,你决定用递归函数来写程式,先依次在56789中选出5,然后把剩下的6789当做是一个新的求四位数排列的问题再次调用函数,以得到所有以5为字头的排列;再选出6,剩下的5789调用函数.而每次求6789或是5789的排列时再把它简化成另一个求3个数字的排列问题,直到要求的排列只剩下一个数.
以下就是你的函数,不过你还不知道它到底是不是正确的,因为写递归函数很易出错,因此你要先试一下.
切换行号显示
1#permute2.py
2defpermute(seq):
3l =len(seq)
4ifl ==1:
5return[seq]
6else:
7res=[]
8foriinrange(len(seq)):
9rest = seq[:i] + seq[i+1:]
10forxinpermute(rest):
11res.append(seq[i:i+1] + x)
12returnres
13
14seq =list("1234")
15thelist =permute(seq)
16thelist = [''.join(x)forxinthelist ]
17printthelist
你运行后得到以下的结果:
$ python permute2.py
['1234', '1243', '1324', '1342', '1423', '1432','2134', '2143', '2314',
'2341', '2413', '2431', '3124', '3142', '3214','3241', '3412', '3421',
'4123', '4132', '4213', '4231', '4312', '4321']
看来是正确的.但有没有办法再快一些呢?你想了半天,终於发现其实你不必等到l = 1的时候才有所动作,你可以在l = 2的时候就干些事了.因为你知道l = 2的话,排列一定是[ [0,1], [1,0] ]的,这样你起码可以用些力气帮电脑一把.当然如果你把l = 3的排列也写出来更好,但写到l = 4或以上大可不必了.这种帮它一下的做法在程式优化中的学名是unroll,你隐约记得是学过的.好,现在你有另一个程式了.
切换行号显示
1#permute3.py
2defpermute(seq):
3l =len(seq)
4ifl <=2:
5ifl ==2:
6return[ seq, [seq[1], seq[0]] ]
7else:
8return[seq]
9else:
10res=[]
11foriinrange(len(seq)):
12rest = seq[:i] + seq[i+1:]
13forxinpermute(rest):
14res.append(seq[i:i+1] + x)
15returnres
16
17seq =list("12345")
18thelist =permute(seq)
19thelist = [''.join(x)forxinthelist ]
20printthelist
现在你可以正式测试了.你把permute3.py改了一下,以便可以从命令行取得数字以及分割方法.程式变成下面的样子,同时你也对permute2.py做了相同的修改.
切换行号显示
1#permute3.py
2defpermute(seq):
3l =len(seq)
4ifl <=2:
5ifl ==2:
6return[ seq, [seq[1], seq[0]] ]
7else:
8return[seq]
9else:
10res=[]
11foriinrange(len(seq)):
12rest = seq[:i] + seq[i+1:]
13forxinpermute(rest):
14res.append(seq[i:i+1] + x)
15returnres
16
17importsys, calc
18seq =list(sys.argv[1])
19where =int(sys.argv[2])
20thelist = [''.join(x)forxinpermute(seq) ]
21Print'Got',len(thelist),'items.'
22calc.calc(thelist,where)
你开始试行了.用time方式来看程式到底运行了多长时间吧.
$ time python permute2.py 56789 3
Got 120 items.
Maximum at 87596 ,product 84000
real0m0.057s
user0m0.050s
sys0m0.000s
$ time python permute3.py 56789 3
Got 120 items.
Maximum at 87596 ,product 84000
real0m0.040s
user0m0.030s
sys0m0.010s
呵,不错.修改了的就是快些.到了这个地步,你开始觉得好奇了.像求排列这样一个常见的问题,不知道别人都是怎样做的呢.也许应该到网上去找找看,或者有一两个已经写好的程式片断可以抄的.你可不想弄错一个原来己经有标准答案的问题.把permute2.py贴上留言板或者会令自己看起来像一个三流的程式设计员,这可是你最不想见到的.於是你在网上到处搜寻.不过似乎都是以递归算法为主的,直至用了一些时间,你终於在ASPN:的网上程式码收集站上看到了这一个片断:
Getting th permutation
of a sequence « Python recipes « ActiveState Code
切换行号显示
1# permute4.py
2defpermute(seq, index):
3seqc = seq[:]
4seqn = [seqc.pop()]
5divider =2
6whileseqc:
7index, new_index =divmod(index,divider)
8seqn.insert(new_index, seqc.pop())
9divider +=1
10return''.join(seqn)
作者声称这个算法的量级是O(n)的.你觉得难以置信,但不妨一试.由於你理解到这个函数每次只传回排列中的某一项,因此你写了个小程式去验算它.
切换行号显示
1# test.py
2frompermute4.pyimportpermute
3
4seq =list("1234")
5foriinrange(30):
6printpermute(seq, i),
试验一下:
$ python test.py
1234 1243 1324 1423 1342 1432 2134 2143 3124 41233142 4132 2314
2413 3214 4213 3412 4312 2341 2431 3241 4231 34214321 1234 1243
1324 1423 1342 1432
测试显示这个函数没问题.但它怎样做到的呢?你研究了一下,每个不同的index值都传回唯一的排列,而且大过n!的index会从头来算起, divider每次都要增加一,列表的排法又是按商余数来重整.唉,不得要领.嗨!管它呢.反正能用就行了.於是你修改permute4.py,加入一个新的函数求factorial,这样就可以调用permute得到所有的排列.并进行计时.你用了更多的数字,以便速度的差别更明显些.
切换行号显示
1# permute4.py
2defpermute(seq, index):
3seqc = seq[:]
4seqn = [seqc.pop()]
5divider =2
6whileseqc:
7index, new_index =divmod(index,divider)
8seqn.insert(new_index, seqc.pop())
9divider +=1
10return''.join(seqn)
11
12deffact(x):
13f =1
14foriinrange(1,x+1):
15f *= i
16returnf
17
18importsys, calc
19seq =list(sys.argv[1])
20where =int(sys.argv[2])
21n = fact(len(seq))
22thelist = [
permute(seq, i)foriinrange(n) ]
23print'Got',len(thelist),'items.'
24calc.calc(thelist,where)
$ time cpython permute3.py 1234567 4
Got 5040 items.
Maximum at 6531742 ,product 4846002
real0m0.461s
user0m0.440s
sys0m0.020s
$ time cpython permute4.py 1234567 4
Got 5040 items.
Maximum at 6531742 ,product 4846002
real0m0.389s
user0m0.370s
sys0m0.010s
哇!的确不是盖的.很好,而且现在你知道了别人不知的新答案.就把它贴上去罢.就在你决定的时候按钮之际,你到底犹豫了: "我对这个算法不是很了解,如果别人问起的话怎样回答呢?这会让我像个东抄西抄的小偷呢!不,要不我要明白它的原理,不然就自己做一个比它更好的."你觉得壮志无限.
但是现在已经很晚,你要去睡了.无奈你在床上反覆地思考著更好的方法,你整个晚上都没睡好.
待续......
第二天
第二天
你醒来第一件事,洗脸刷牙.编程爱好者并不一定和终日蓬头垢面同义.然后呢,看看电视报纸,做些公益活动,今天是礼拜天嘛.废话少说,终於你在电脑前坐下,登入了你喜爱的Slackware /RedHat/ Redflag / Mandrake / Debian /WindowsXP / Chinese2000 / DOS / Solaris/ AIX / Unicos / OSX [作者按:请依实际情况增删,但千万拜托不要把SCO也加进来],这些都是Python能够运行的平台.
你记起你以前学到递归时听过的话:任何递归/回溯函数都可以还原成非递归形式的.於是你决定用你自己的方式一试.你默念著求排列的方法, 5个数取一个,剩下4个,再取一个,剩下3个....於是你写出一个新的程式,和最初的一个很相像:
切换行号显示
1# permute5.py
2defpermute(seq):
3result = []
4foriinseq:
5seq1 = seq[:]
6seq1.remove(i)
7forjinseq1:
8seq2 = seq1[:]
9seq2.remove(j)
10forlinseq2:
11seq3 = seq2[:]
12seq3.remove(l)
13forminseq3:
14seq4 = seq3[:]
15seq4.remove(m)
16result.append(''.join([i,j,l,m,seq4[0]]))
17returnresult
18
19printpermute(list("12345"))
这个程式依次创建5, 4, 3, 2, 1个数的列表,每个都不包括之前被选的数字,然后把5个数合起来完成一种排列.当然,你还有空间做unroll.但现在问题在於,你对程式的要求是事先并不知道要求多少个数字的排列,就是你不知道要写几个for才够.但现在你认为有一个好办法:既然Python是动态的,它可以执行自己写出来的编码,为什么不叫它自己帮自己来写这个循环程式以完成工作呢?你觉得这种让程式来为自己写程式的想法很科幻也很妙,而且让你记起了以前听到很多资深程式员爱用的m4宏语言.於是你赶紧试了一个.你认为可以用counter0, counter1,
counter2...来代替i, j, l, m...的循环子命名法.
切换行号显示
1# permute5.py
2defgenfunc(n):
3head ="""
4def permute(seq0):
5result = [] """
6boiler ="""
7for counter%i in
seq%i:
8seq%i = seq%i[:]
9seq%i.remove(counter%i)"""
10foriinrange(1,n):
11space =''*i
12head = head + boiler.replace('\n','\n'+space)%(i,i-1,i,i-1,i,i)
13temp =','.join(['counter%i'%(x)forxinrange(1,n) ] + ["seq%i[0]"%(n-1)])
14head = head +'\n'+ space +"result.append(''.join([%s]))"%(temp)
15returnhead +'\nreturn result\n'
16
17importsys
18functext = genfunc(len(sys.argv[1]))
19printfunctext
20exec(functext)
21printdir()
22thelist = permute(list(sys.argv[1]))
23print'Got',len(thelist),'items.'
运行一下,
sh-2.05b$ python permute5.py 12345 3
def permute(seq0):
result = []
for counter1in seq0:
seq1 =seq0[:]
seq1.remove(counter1)
forcounter2 in seq1:
seq2 =seq1[:]
seq2.remove(counter2)
forcounter3 in seq2:
seq3 = seq2[:]
seq3.remove(counter3)
forcounter4 in seq3:
seq4= seq3[:]
seq4.remove(counter4)
result.append(''.join([counter1,counter2,counter3,counter4,seq4[0]]))
return result
['__builtins__', '__doc__', '__name__', 'calc','functext', 'genfunc',
'permute', 'sys']
Got 120 items.
看来格式是弄对了.现在算算运行时间,会不会好些呢?也许会比以前更快,也许因为要再执行自己产生的程式而更慢,一切都要看实际的数据才能呢!你修改了permute5.py以便它能标准化地计算时间.你开始觉得import calc实在是挺聪明的设计.
切换行号显示
1# permute5.py
2defgenfunc(n):
3head ="""
4def permute(seq0):
5result = [] """
6boiler ="""
7for counter%i in
seq%i:
8seq%i = seq%i[:]
9seq%i.remove(counter%i)"""
10foriinrange(1,n):
11space =''*i
12head = head + boiler.replace('\n','\n'+space)%(i,i-1,i,i-1,i,i)
13temp =','.join(['counter%i'%(x)forxinrange(1,n) ] + ["seq%i[0]"%(n-1)])
14head = head +'\n'+ space +"result.append(''.join([%s]))"%(temp)
15returnhead +'\nreturn result\n'
16
17importsys, calc
18functext = genfunc(len(sys.argv[1]))
19#print functext
20exec(functext)
21thelist = permute(list(sys.argv[1]))
22print'Got',len(thelist),'items.'
23calc.calc(thelist,int(sys.argv[2]))
开始计时:
$ time cpython permute5.py 1234567 4
Got 5040 items.
Maximum at 6531742 ,product 4846002
real0m0.213s
user0m0.200s
sys0m0.010s
啊哈!那个什么量级O(n)的也被你打败!!你觉得它的量级其实不是O(n),那只是对找一整个排列其中一个的时候才有用,要把整个排列都算出来的话还是要回到n!上的.
你非常自豪.但这也许是适当的时候提醒自己谦虚的妤处.既然都到了这个地步了,何不再走多一步,翻一下书看看,也许你找到的方法已经早有别人发现了.果真是这样的话,你,一个无知的白痴,到处大吹大擂自己的发明是会被笑话的.
於是你找出了封尘的电脑和数学教科书.找到了排列组合一章,开始细读.终於你找到了这样的一幅图画:
[4321]
[3421]
[321]< [3241]
[21] <[231]... [3214]
[213]...
[1] <
[321]...
[21] <[231]...
[213]...
书中写到,要产生一个排列其实可以用这样一个方法:
"先选一个数1,然后第二个数2可以放在1的前面或是后面.而每一个放法都会产生一个2位数,对於每一个这样的两位数,第三个数3,都可以放在它的前面,中间,或是最后;如此产生一个3位数;而每一个3位数,第4位数都可以插入到这3个数的任何一个空位中,如此类推.书中还列出了一个程式范例呢!并声这个方法要和已知的最快的算排列的方法速度相若.
你急不可待地开始把书中的描述实现.用Python,你很快又得到了一个全新的程式:
切换行号显示
1# permute6.py
2defpermute(seq):
3seqn = [seq.pop()]
4whileseq:
5newseq = []
6new = seq.pop()
7#print "seq:",seq,'seqn', seqn ,'new', new
8foriinrange(len(seqn)):
9item = seqn[i]
10forjinrange(len(item)+1):
11newseq.append(''.join([item[:j],new,item[j:]]))
12seqn = newseq
13#print 'newseq',newseq
14returnseqn
15
16importsys, calc
17seq =list(sys.argv[1])
18where =int(sys.argv[2])
19thelist =permute(seq)
20print'Got',len(thelist),'items.'
21calc.calc(thelist,where)
测试结果如下:
$ time cpython permute6.py 1234567 4
Got 5040 items.
Maximum at 6531742 ,product 4846002
real0m0.167s
user0m0.150s
sys0m0.020s
哇塞!书中自有黄金屋咧!想不到这个才是最快的算法.你开始感到要击败这次的对手不是不件容易的事,而且现在已经很晚了,你身心也都被疲倦所包围著.你绝望地看著这个新的程式码和它那美妙的结构,作出最后的尝试:
待续...
守夜人
守夜人
Got 24 items.
[
'1234', '2134', '2314', '2341',
'1324', '3124', '3214', '3241',
'1342', '3142', '3412', '3421',
'1243', '2143', '2413', '2431',
'1423', '4123', '4213', '4231',
'1432', '4132', '4312', '4321'
]
上面就是permute6.py产生的四位数字排列结果,你细心地反覆观看,终於看出了一些端倪:其实所产生的排列是有一种对称性的,第一个和最后一个是完全次序相反的,而第二个又和倒数第二个完全相反.利用这些对称性,也许你可以把计算时间打个对折哟.而你研究了一下程式的实现方法后你发现只要改一行!就可以实现这样的功能:就是第一行seqn = [ seq.pop() ]改成seqn = [
seq.pop()+seq.pop() ].这样你就实现了只产生其中一半的排列,尔后你只要把这个列表中的元素都掉个就完成了整个排列.程式如下
切换行号显示
1# permute7.py
2defpermute(seq):
3seqn = [ seq.pop()+seq.pop() ]
4whileseq:
5newseq = []
6new = seq.pop()
7#print
"seq:",seq,'seqn', seqn ,'new', new
8foriinrange(len(seqn)):
9item = seqn[i]
10forjinrange(len(item)+1):
11newseq.append(''.join([item[:j],new,item[j:]]))
12seqn = newseq
13#print 'newseq',newseq
14returnseqn
15
16importsys, calc
17seq =list(sys.argv[1])
18where =int(sys.argv[2])
19thelist =permute(seq)
20print'Got',len(thelist),'items.'
21printthelist
22#这个等一下再探讨
23#calc.calc2(thelist,
where)
测试数据表示果然这个改良了的程式要比原来快了刚好一倍.这次应该足够了.但是要得到整个排列的话要把这半个列表重抄一次而且每个元素都要反过来: "1234" -> "4321".但是在Python之中的字串并没有反序的函数,因此你必须先把字串变成列表,再反过来,然而list.reverse()这个函数偏偏很讨厌的不会传回任何值(因为它是in-place的缘故),所以你要用i = list(item);
i.reverse; i =''.join(i);这个复杂的方法.你想了想,这个做法大概会把原来只做一半排列所省下来的时间都浪费掉了.你思考半天,终於决定重写calc.py部份以便直接利用已知的半个列表.
切换行号显示
1# calc.py
2defcalc(seq, where):
3maximum, max_item =0, []
4foriinseq:
5product =int(i[:where]) *int(i[where:])
6ifproduct > maximum:
7maximum, max_item = product, i
8elifproduct == maximum:
9max_item +=','+i
10print"Maximum at", max_item,",product", maximum
11
12defcalc2(seq, where):
13maximum, max_item =0, []
14foriinseq:
15product =int(i[:where]) *int(i[where:])
16ifproduct > maximum:
17maximum, max_item = product, i
18elifproduct == maximum:
19max_item +=','+i
20rev =list(i)
21rev.reverse()
22i =''.join(rev)
23product =int(i[:where]) *int(i[where:])
24ifproduct > maximum:
25maximum, max_item = product, i
26elifproduct == maximum:
27max_item +=','+i
28print"Maximum at", max_item,",product", maximum
当然你保留了以前的函数calc而只是新加一个专门给permute7.py调用的calc2函数.你试了一下速度,成功了比permute6.py快了一丁点.虽然只是这一点儿点儿,你也觉得快活无比.因为又一次,你为一个大家都觉得好的方法做出了改良.
虽然你知道在这个阶段如果你把calc.py整合到排列产生器中也许会得更好的提速效果,但你觉得现在这样已经可以很多人都认同你的能力了.而且能把一个高效的排列产生器独开来也许是聪明的做法,因为将来你一定会再用它的.你看了一下所有的程式,从permute1.py到permute7.py,再做了一次速度的检定.反正是最后一次了,干脆做个大的吧.
$ time python permute2.py 123456789 5
Got 362880 items.
Maximum at 875319642 ,product 843973902
real0m46.478s
user0m46.020s
sys0m0.430s
$ time python permute3.py 123456789 5
Got 362880 items.
Maximum at 875319642 ,product 843973902
real0m38.997s
user0m38.600s
sys0m0.400s
$ time python permute4.py 123456789 5
Got 362880 items.
Maximum at 875319642 ,product 843973902
real0m33.845s
user0m33.460s
sys0m0.380s
$ time python permute5.py 123456789 5
Got 362880 items.
Maximum at 875319642 ,product 843973902
real0m10.681s
user0m10.530s
sys0m0.150s
$ time python permute6.py 123456789 5
Got 362880 items.
Maximum at 875319642 ,product 843973902
real0m8.279s
user0m8.110s
sys0m0.170s
$ time cpython permute7.py 123456789 5
Got 181440 items.
Maximum at 875319642 ,product 843973902
real0m7.827s
user0m7.650s
sys0m0.180s
嗯,很不错.最快的要比原先快上近七倍呢!於是你打算明天一有空便把这个最好的公式贴到网上去,让更多人分享.你很放心地去睡觉了.
但是不知为了什么,你总觉得有些事烦扰著你,还有什么不妥的地方呢?你实在想不到了,迷迷糊糊地抱著迷惑不安的心情入梦.
终於你忍不住爬了起床,再一次回到电脑屏幕前.你想到了一个致命的问题,对於很大很大的排列, permute7.py还是会尝试一下子把所有的排列都做出来,不用多久电脑资源就会被耗光的.你也许要另想一个方法,每次只做一个排列.但你是否可以把所有的排列用1, 2, 3, 4的方法数出来呢?
你看著教科书上的那幅图画,这样的树状结构应该有办法的,你对自己说.
慢慢读著书上的文字.设想有n个数字,先取第一个数字.再取第二个数字,第二个数可以放在第一个数的左或右面,就是有0, 1两个选择.再取第三个数,放到前面选好的两个数字中,可以放在最左,中间,最右,就是有0, 1, 2三个选择.嗯,很自然吗.忽然你想到了二进位,八进位那些数系转换关系. "我可以设计这样一个数, ...xyz,其中个位数z是二进位的,也就是放第二个数的两个位置;十位数y是三进位的,化表放第三个数字的三个位子,然后百位数是四进位,千位数是五进位的,依以类推."没错,这样设计的话,如果0表示放於最左面的话,则"2021"这个数就代表了排列五个元素(abcde),取一个a,然后第二个b放在a的右面成ab,取c放到最右面成为abc,取d放到最左面成dabc;最后e放到中间去成为daebc.至於"2021"这个特别的设计的数可以用..+ x*4 + y*3 + z*2这样的计算来映对到自然数的数列上去.
没错了,如求4个数的4! = 24个排列,第18个排列可以这样求得, 18除2,余数是0,所以第二个数放在第一个数的左面;然后商9再除3,余数0,所以第三个数於在头两个数的最左;最后3除以4,余数是3,因此第四个数要放在前三个数的第4个空位,也就是最右面.利用这个方法,你就不必先求出整个排列才能开始计算.尽管这好像牺牲了时间,但省下了大量的空间.你完全可以算出1000个数的最大排列方法,纵使那可能要用去几个月的运算.你更高兴的是用这个方法,你可以把整个计算拆开成为一个一个的部份:譬如说求10! = 3628800个排列,你大可以把1到1000000让一部电脑做, 1000001到2000001让另一部来做...大家的工作并不重覆,这等於实现并行运算了!啊哈!妙极了!
忽然你灵光一闪,对了,这个不正是permute4.py的算法吗!你昨天还久思不得其解呢,现在你已经完全明白了.呜,虽然这么好的算法原来不是你原创的,但是你仍然异常兴奋.因为你的思路竟和那些大牛们互相吻合.你渐渐记起了当你还在玩DOS游戏机的年代,曾有个古怪的人向你吹嘘过某个电脑扑克游戏中用了一个威力很大的洗牌法,多么的快而且公平,不必怕黑客们用已知的随机数表来出千.现在你猜到了,那个算法很可能就是眼下这一个了.有52张牌,如果要预先算出52!个排列才能洗牌可真要命呢.
你觉得舒服多了,你整理了一下程式,打算把结果告诉大家.然而刚才的不安感又重新来袭.你再看了一次最后也应该是最棒的程式,心中安慰自己道: "千万不要跌进低等程式员的的陷阱,他们疯也似的不断追求,郤从来不知道自己的目标,也不知道什么才是好.完美的设计不在于你无法添加新的功能,完美的设计是再也不能精简现有的功能."你觉得permute7.py已迫近了这一个极限.但你内心深处并没有因此而舒坦开来,一种悬空的感觉自足下升起.也许是你太累了,于是者你决定闭目养神数分钟再开始上网,可惜你很快地迷迷糊糊地进入了梦境.
待续...
终篇
终篇
你做了一个梦,梦中你看到阿凡提骑著他那出名的毛驴来到你面前并向你提出挑战: "除非你解答了我的难题,不然我的驴子会不停在你耳边嘶叫令你无法睡好!问题是:
把数字56789放到
[][][]*[][]里
得出最大的的乘积....
你发出会心的微笑,正想祭出你的permute7.py之时忽然想起阿凡提是不可能懂得电脑编程的!你心中登时凉了一大截:阿凡提的方法一定不必用电脑算出所有的排列方法,并很快的得知答案的.随著一声震天的驴嘶,你惊醒了,发现原来你伏在电脑桌上睡去了,不小心压著了键盘上的方向键而令你的电脑发出了痛苦的BEEP声.
回想梦境,你打算暂时离开电脑,回到问题本身上来:怎样才能"看"出最大的乘积呢?
你拿出纸笔,开始计算:
假设五个数为[a][b][c]*[d][e],展开的话成为
a * 100 * d *10
+ a * 100 * e * 1
+ b * 10* d *10
+ b * 10* e *1
+ c * 1* d *10
+ c * 1* e *1
这个可以写成一个矩阵:
de
a1000100
b10010
c101
你这样想到:在整个答案中, a带来的贡献是一个百位数加上一个十位数,而d的贡献是一个百位数加十位数加个位数,因此d要比a更重要.要取得最大的积则一定要把56789中最大的9放在d的位置,然后是a,如此类推.
为了方便计算,你干脆用对数来记数100 = 10e2,用2来代表好了,因此:
d e
a3 2
b2 1
c1 0
计算每一行或列的和,把它称为该数的基值,我们得到
a = 5, b = 3, c = 1, d = 6, e = 3.
咦? b和e的基值是一样的,怎么办!
你思考著: "啊哈!因为我们用了对数,而log(1) = 0因此把b和e之间的微小分别给忽略了!"好吧,试把每个数都加大一个,得到:
d e
a4 3
b3 2
c2 1
这样基数变成了: a = 7, b = 5, c = 3, d
= 9, e = 6.这些基数代表了该位置的重要性,可以按大小分予不同的数字.
好,按基数的大小来分配数字你得到了865 * 97.一对答案,哟!不一样!正确解是875 * 96.哪里不对了呢?仔细分析下来,你发现b和e互换了.换的原因是这样的:
b的贡献: b * d * 100 + b * e *
10 e的贡献: e * a * 100 + e * b * 10 + e * c
粗看的话e的贡献要大一些,但因为我们把9配给了d而8配给了a,因此最终的结果是b的实际贡献要比e大.由於e和b的基数只相差在e * c这个个位数乘积上, d和a之间的分配结果把这个小的差异覆盖掉了.
你考虑著: "要把这样的覆盖也算上去的话,也许可以做一个二阶基数.如b * d的基数是100,但是由於d的基数为9,那b的二阶基数可以算成9,代表和b相关的是一个较为重要的数;同样e * a的基数是也是100但由於a的基数只是7,因此e的二阶基数只是7而已.这样就可以得出b要比e更重要了."
於是你有了一个想法:先写出相关矩阵,然后计算每个数的基数和二阶基数,再进行排序,当两个基数很接近时就用二阶基数来判别哪个较重要.嗯,你觉得自己聪明极了,於是开始验算,但很快又发现其实b和e的二阶基数原来也是一样的!!大家都是15.也许我们要用一个三阶基数才能分辨他们.
你又想了一些新的二阶基数的定义,有些的确给出正确的答案,但你渐渐觉得这一切并不很妥当.因为就算能计出56789,但是在更多的排列,如7位数甚至9位数的排列你怎样保证答案也一定准确呢,而两个基数到底怎样才算是比较接近呢?仔细审视你的方法,用对数来表示乃至直接相加来求所谓的基数统统都是即兴的,毫不严谨.或者要真正解决他们必需要把每一种情况都算进来,而那样做的话必然要计算n!那么多次!说到底还是要算排列的.
你有些失望,但暗中觉得松了一口气,因为到底是permute7.py得到最后的胜利.你伸了一下懒腰,原来天都快亮了.这时你感到有些饿,便去拿了半个凉馒头,冲了一些可可.你对著空空的萤光屏,静静地坐了大概十分钟,然后答案进入你的脑海,谜团被解开了.
你的方法是求出所有位置的"重要性"(用你的语言就是求基数),然后依次把最大的数字分配给最重要的位置.但是位置的重要性是和其他位置纠缠在一起的,因此一次过算出所有位置的重要性必须考虑大量不同的组合排列,并不实际.
但是,我们其实可以每次只求第一个最大的基数的位置(abc*de的情况下最大的基数是d),这个最大的基数是没有争议的.当求得这个位置时,干脆把最大的数字放到该位子上,然后再求一次基数,找出接下来最大的位子,这个位子也是无可争议的.如此一来,原来求5个数字的全排列成就简化为5次简单的回圈.一个求Factorial(n)的问题变成了n次循环!
啊哈!
你精神大好,从另一个角度切入:
假如5个数字一样, 11111,最大的乘积只能是111 * 11,现在容许改大一个数,改哪个会使结果最大?
211 * 11, 121 * 11, 112 * 11, 111 *
21, 111 * 12 ?答案是111 * 21,也就是d的位置.好,把d替换成9.
问题变成5个数字, 111 * 91,改一个数(除了9),改哪一个?
211 * 91, 121 * 91, 112 * 91, 111 *
19 ?答案是211 * 91,也就是a的位置.好,替换成8.
依此类推,答案正正是875 * 96.
你重开电脑,很快地把新方法输入程式,并改名为wise.py.
切换行号显示
1defsolve(seq,where):
2n =len(seq)
3seq.sort()
4seq.reverse()
5table = [ []foriinrange(n) ]
6left, right = where, n - where
7leftr =long('1'*left)
8rightr =long('1'*right)
9flag=[]
10foritemin[int(x)forxinseq]:
11foriinrange(left):
12table[left-i-1] = (leftr +10**i) * rightr
13foriinrange(right):
14table[right-i+where-1] = leftr * (rightr +10**i)
15foriinflag:
16table[i] =0
17tablesorted = table[:]
18tablesorted.sort()
19maxindex = table.index(tablesorted[-1])
20ifmaxindex >= where:
21rightr = rightr + (item-1) *10**(right-maxindex+where-1)
22else:
23leftr = leftr + (item-1) *10**(left-maxindex-1)
24flag.append(maxindex)
25#print maxindex, leftr, rightr
26returnleftr, rightr
27
28importsys
29leftr, rightr =
solve(list(sys.argv[1]),int(sys.argv[2]))
30print"Maximum
at", leftr,rightr,',product', leftr*rightr
你验算了一下结果,完全正确!这时你好奇地再次time了一下程式的速度
$time python permute7.py 123456789 5
Got 181440 items.
Maximum at 875319642 ,product 843973902
real0m7.827s
user0m7.650s
sys0m0.180s
$ time python wise2.py 123456789 5
Maximum at 87531 9642 ,product 843973902
real0m0.042s
user0m0.010s
sys0m0.030s
哇!快了近两百倍!当然了.如果算更多位的排列会快更多,因为wise.py跳离了n!的限制.
你现在觉得舒服多了.你真的解了这个问题.你不再怕有人会写出更快10倍的程式了.你既有了"聪明"的答案(软解)来对付阿凡提和他的驴子,而在硬解方面,你也自信有世界第一流的排列产生器.你完全满足了,你不再感到疲累,心中疑犹一扫而空.这时你身体感到一阵震栗但心中却喜乐无穷,你第一次感受到了编程之道的洗礼.并且,你学会了所有程式大师都有的态度:我没法用中文来形容,这种态度叫做"to hack".你知道只要你熟练并保持这种态度来面对生命中的难题,你很快就可以满师出山了.
你最后一次浏览了一下你的程式码,发现在wise.py中,其实每一个循环完成后,最重要的位置和最次要的位子都是不容争议的,因此大可放心地替换两个数字而不是一个,那程式可以再快一倍.不过你觉得现在己经很够了,你颇有禅机地自言自语道: "我已经找到明月,再纠缠只下去只是妄执於指月的手而已."你熟练地登出系统并关上电脑,你知道这次你可以真正安心地睡一觉了.
哎哟!天已亮了,今天是礼拜一,你要上班的.喔!又要被老板说上班一条虫,下班一条龙......惨.......
全篇完.
检讨总结
检讨总结
一)在上面的故事中,我们看到了解决编程问题的五个方法.
把问题规范成一个普遍的形式,这样更容易和别人交流以及找相关资料.
自己尝试找答案.
在网上搜寻更好的答案.
想一个方法来打败这个更好的答案.
翻查教科书或是文献,从基本开始思考有没有最好的解.这些书能被选成教本一定有它的原因.
研究问题的特殊情况,也许会有别出心裁的巧妙方法.
二)故事中每个程式都只有二三十行大小,说明Python语言表达力强且语意很浓缩,做为快速开发或是测算自己的想法都是非常好的.
三) Python程式浓缩之余,它的语言也异常的清晰.回看上面的程式,你会发现它们全都不难明白.这说明Python程式更加容易维护.
四)在故事中,我们有很大的篇幅是在讨论方法而只有小部份是在描述Python的语言特性.这证明Python更适合用来教授编程的概念.事实上, Python的作者Guido和很多人都认为Python是电脑教育的首选语言.教师可以让学生静静地思考,学通运算的法则;而不是上来就疯狂地敲打键盘,并要记住一大堆电脑语言的古怪特徵.
五)整个故事围绕於算法的改进而较少碰到Python程式的优化问题.也许在续集中(如果有的话),我们要尝试一下在固定的算法及尽量少改动程式码的条件下,提高Python程式的效率.我暂时想到的方法包括:
利用较新和较快的语法.如yield, generator.
用Python的自带优化选项以及内建模组.
用第三方的扩展模组,如Numpy, Scipy.
用编译方式代替解释,如freeze, py2exe.
用JIT类的方法,如Psyco.
用Thread,在多CPU的机器上平行运算.
最后一样要大改程式了,用C来做扩展.
更有'to hack'感觉的,修改Python主干程式,加入像string.reverse()这样的辅助函数.
六)文中所用的测试硬件:
CPU: Pentium III 866 RAM: 128 MB
文中所用的测试软件:
Slackware Linux: 9.0 Linux Kernel:2.4.2GCC: 3.2.2
Python:修改过的2.1.3
七)啃凉馒头对脑筋有帮助.
八)如果你能想到更好的方法,欢迎联络本人: