Humble Number 丑数和丑代码

这道题目的经历让我感触很深，我做的题目是浙大的，链接在这里；

Humble Number 的原题可能是北大的，北大的题目更广义一些，是自定义的基数，浙大的只有（2,3,5,7）四个基。

 

题目看了一遍以后第一感觉就是2000000，n多个零的那个数，让我感觉遍历肯定是会傻的。所以一定是构造。

于是我写了下面第一个代码，思路是这样的，假设我们从某一段开始，已有了[1,2,3,4,5,6,7]这几个结果，

那么接下去的一个数，一定是前面的这些数乘以2,3,5,7得到的。

 

而对于1, 1×2=2, 1×3=3, 1×5=5, 1×7=7, 可见1的机会已经用完了。

 

那么对于2, 2×2=4, 2×3=6, 意味着还有5和7可以作为因子相乘。

 

于是我定义了一个结构：这个结构定义了with*，表示*已经被乘过了。

class node: value=None; with2=None; with3=None; with5=None; with7=None; def __init__(self,value,with2,with3,with5,with7): self.with2=with2; self.with3=with3; self.with5=with5; self.with7=with7; self.value=value; 

并且可以找出如下的结果，假设到12的时候，12以后这个数只会在2×7, 3×5, 4×3, 5×3, 6×2

中间产生。而我们发现这些数都是大于等于12， 所以要是等于12，就吧with* 标记成1。如果大于12，就找大于12最小数。

 

而找这些可能的结果也很简单，发现下图第二个矩阵中1排列是个左上三角，沿着边找就可以了。

 

 

这样看来这是一个极好的思路，我以为应该就能解决掉这个问题了。

于是有了下面的代码：

import time class node: value=None; with2=None; with3=None; with5=None; with7=None; def __init__(self,value,with2,with3,with5,with7): self.with2=with2; self.with3=with3; self.with5=with5; self.with7=with7; self.value=value; #------------------------------------------------------------------------------- def initial(): #n1=node(1,1,1,1,1); n2=node(2,1,1,0,0); n3=node(3,1,0,0,0); n4=node(4,0,0,0,0); n5=node(5,0,0,0,0); n6=node(6,0,0,0,0); n7=node(7,0,0,0,0); #base.append(n1); base.append(n2); base.append(n3); base.append(n4); base.append(n5); base.append(n6); base.append(n7); #------------------------------------------------------------------------------- def printbase(): for i in range(len(base)): print base[i].value; #------------------------------------------------------------------------------- def printbasedetail(base): for i in range(len(base)): out=str(base[i].value)+" :"+str(base[i].with2); out=out+","+str(base[i].with3); out=out+","+str(base[i].with5); out=out+","+str(base[i].with7); print out; #------------------------------------------------------------------------------- def findResonableBase(): resonableBase=[]; for i in range(len(base)): if base[i].with2!=0 and base[i].with7==0: resonableBase.append(base[i]); else: resonableBase.append(base[i]); break; return resonableBase; #------------------------------------------------------------------------------- def findTheMultipler(node): if node.with2==0: return 2; elif node.with3==0: return 3; elif node.with5==0: return 5; elif node.with7==0: return 7; else: return None; #------------------------------------------------------------------------------- def MarkWithValue(i): if base[i].with2==0: base[i].with2=1; elif base[i].with3==0: base[i].with3=1; elif base[i].with5==0: base[i].with5=1; elif base[i].with7==0: base[i].with7=1; #------------------------------------------------------------------------------- def findRightBase(resonableBase): minList=[]; minIndexList=[]; min=resonableBase[0].value*findTheMultipler(resonableBase[0]); index=0; minList.append(min); minIndexList.append(index); for i in range(1,len(resonableBase)): temp=resonableBase[i].value*findTheMultipler(resonableBase[i]); if temp<minList[0]: minList=[]; minIndexList=[]; minList.append(temp); minIndexList.append(i); elif temp==minList[0]: minList.append(temp); minIndexList.append(i); result.append(minList[0]); n=node(minList[0],0,0,0,0); base.append(n); for i in range(len(minIndexList)): MarkWithValue(minIndexList[i]); if base[index].with2==1 and base[index].with3==1 and base[index].with5==1 and base[index].with7==1: base.remove(base[index]); return resonableBase[index]; #------------------------------------------------------------------------------- result=[1,2,3,4,5,6,7]; base=[]; initial(); t1=time.time(); while len(result)<5000: resonableBase=findResonableBase(); ## printbasedetail(resonableBase); ## print '============' rightBase=findRightBase(resonableBase); ## printbasedetail(base); ## print '............' ## print result; ## print '************' t2=time.time(); print result[99]; print t2-t1  

而事实上是结果出来了，但是非常的慢，后来我用time测试了一下，计算要求的5000个，话了6秒，晕，肯定过不了了。

 

 

后来的思路来自于这片文章，看了他的程序我豁然开朗。

并且找到了慢的原因，基数太多了！

 

人家用的大致思路是一样的，探索找下一个，然后填满5000个。

但是人家用了4个质数做了基，那么比较的次数就大大降低了，

而我用的是对角线数做基数，并且还有标记的工作，加上使用了类有开销，所以慢了100倍！！！我去……

 

我用python重写了一遍：

import time t1=time.time(); pindex=[0,0,0,0]; prim=[2,3,5,7]; result=[1]; for i in range(5842): for j in range(4): while result[pindex[j]]*prim[j]<=result[i]: pindex[j]=pindex[j]+1; min=2000000001; for j in range(4): if result[pindex[j]]*prim[j]<min: min=result[pindex[j]]*prim[j]; result.append(min); print result[99]; t2=time.time(); print t2-t1 

 

两种结果的比较，一个6秒一个0.06秒

我于是在想，习惯了用方便自己的思路做事，苦了我的电脑，哈哈！

算法的力量真大啊！