要查看内存错误,请以limit = 100运行代码,然后打印出链.
也许您想序列化递归代码:
lengths = {1: 1}
def collatz(i):
i0 = i
acc = 0
while True:
if i in lengths:
lengths[i0] = acc + lengths[i]
return acc + lengths[i]
acc += 1
i = (i * 3 + 1) if i % 2 else i // 2
longest = 1
for i in range(1, 1000000):
c = collatz(i)
if c > longest:
longest = c
print(i, c)
当然,仍然可以通过许多方式对其进行优化,但是可以在4秒内产生预期的结果.
编辑:
您的方法将创建一个列表,该列表的长度达到有史以来的最高期限.对于limit = 100,这是9232.这不是很多.但是对于极限= 1000000,它是56991483520(以704511开头的链),这真是个菜鸟.如果仅是int32之上的数组,则该内存已经是212 GB,实际上还不止于此.
这是麻烦的链条:704511、2113534、1056767、3170302、1585151、4755454、2377727,7133182、3596591、10699774、5349887、16049662、8024831、24074494、12037247、36111742、18055871、54167614、27083807、81525122、40625711、121877134、60938567 ,182815702、91407851、274223554、137111777、411335332、205667666、102833833、308501500、154250750、77125375、231376126、115688063、347064190、173532095、520596286、260298143、780894430、390447215、11713416365、17567075、178670780 ,1976639030,988319515,2964958546,1482479273,4447437820,2223718910,1111859455,3335578366,1667789183,5003367550,2501683775,7505051326,3752525663,11257576990,5628788495,16886365486,8443182743,25329548230,12680774679 867 ,3561967720、1780983860、890491930、445245965、1335737896、667868948、333934474、166967237、500901712、250450856、125225428、62612714、31306357、93919072、46959536、234 79768,11739884,5869942,2934971,8804914,4402457,13207372,6603686,3301843,9905530,4952765,14858296,7429148,3714574,1857287,5571862,2785931,8357794,4178897,12536692,6268346,3134173,9402520,4701260,2350260 1175315,3525946,1762973,5288920,2644460,1322230,661115,1983346,991673,2975020,1487510,743755,2231266,1115633,3346900,1673450,836725,2510176,1255088,627544,313772,156886,78443,235330,117330 352996、176498、88249、264748、132374、66187、198562、99281、297844、148922、74661、223384、111692、55846、27923、83770、41885、125656、62828、31414、15707、47122、23561、70684、35342, 17671,53014,26507,79522,39761,119284,59642,29821,89464,44732,22366,11183,33550,16775,50326,25163,75490,37745,113236,56618,28309,84928,42464,21232,10616, 5308,2654,1327,3982,1991,5974,2987,8962,4481,13444,6722,3361,10084,5042,2521,7564,3782,1891,5674,2837,8512,4256,2128,1064,532, 266、133、400、200、100、50、25、76、38、19、5 8、29、88、44、22、11、34、17、52、26、13、40、20、10、5、16、8、4、2、1
使用您的确切递归思想,但使用dict(稀疏)而不是列表(运行时没有问题):
lengths = {1: 1}
def collatz(i):
if i in lengths: return lengths [i]
j = (i * 3 + 1) if i % 2 else i // 2
c = collatz (j) + 1
lengths [i] = c
return c
longest = 1
for i in range(1, 1000000):
c = collatz(i)
if c > longest:
longest = c
print(i, c)