小练习结果揭晓：第一个拥有超过500个约数的三角形数是多少？

冬雪雪冬

首先给大家道个歉，这两天事情较多，没能抽出整块时间看大家的答案，一直拖到今天晚上。
根据题目要求需在5.10 10：00之前完成。
原帖子：http://bbs.fishc.com/thread-71586-1-1.html

题目：

三角形数序列是由对自然数的连加构造而成的。所以第七个三角形数是 1 + 2 + 3 + 4 + 5 + 6 + 7 = 28. 那么三角形数序列中的前十个是：

1, 3, 6, 10, 15, 21, 28, 36, 45, 55, ...

下面我们列出前七个三角形数的约数：

1: 1
3: 1,3
6: 1,2,3,6
10: 1,2,5,10
15: 1,3,5,15
21: 1,3,7,21
28: 1,2,4,7,14,28
可以看出 28 是第一个拥有超过 5 个约数的三角形数。

那么第一个拥有超过 500 个约数的三角形数是多少？


本题的正确答案是：76576500
答题正确的鱼油是：
小火木  6.93s
EggyBruch  16.8s
DingRan  30.9s
bacon6581  10.4s    10行的逗号是中文标点！
lingtingliushen  5.8s
holdme  8.56s

优秀的鱼油是 lingtingliushen

以下鱼油为限定完成时间结束后完成或修改的：
老忘  1.14s
zooo  6.19s  pp我没安装成功，所以没测试
挥舞乾坤  3.89s
小剑剑  13.27s


此外 vpo1 和 宝贝归来 运算时间过长



我写的程序如下，供参考：

1. from math import sqrt
   
2. from itertools import combinations
   
3. 
   
4. def tri(n): #三角形数
   
5.     return n * (n + 1) //  2
   
6. 
   
7. def length(mylist):
   
8.     len1 = 0
   
9.     for i in range(len(mylist)):
   
10.         a =combinations(mylist, i + 1)
    
11.         b = list(set(a))
    
12.         len1 += len(b)
    
13.     return len1 + 2
    
14. 
    
15. #生成质数列表
    
16. primeList = [2]
    
17. item = 101
    
18. n = 3
    
19. while len(primeList) < item:
    
20.     for i in primeList:
    
21.         if n % i == 0: break
    
22.         if i > sqrt(n):
    
23.             primeList.append(n)
    
24.             break
    
25.     n += 2
    
26. 
    
27. n = 15
    
28. while True:
    
29.     list1 = []
    
30.     num = tri(n)
    
31.     num1 = num
    
32.     for i in primeList:
    
33.         while True:
    
34.             if num1 % i == 0:
    
35.                 num1 //= i
    
36.                 list1.append(i)
    
37.             else:
    
38.                 break
    
39.         if num1 == 1:
    
40.             break
    
41.     if length(list1) > 500:
    
42.         print(num)
    
43.         break
    
44.     n += 1
    

复制代码

zooo

目测版主的电脑使用pp模块1秒左右能出结果。。

冬雪雪冬

zooo 发表于 2016-5-12 21:49
目测版主的电脑使用pp模块1秒左右能出结果。。

我下载pp-1.6.5.exe，可是安装出错（python3.5.1），等有时间把python再重新装一下试试。

冬雪雪冬

又下载pp-1.6.4.4.zip，安装成功。

小火木

程序老是败在算法。。。。一直是不查资料直接做题，看样子还得多补充下基本的数学算法

zooo

冬雪雪冬 发表于 2016-5-12 22:23
又下载pp-1.6.4.4.zip，安装成功。

使用pp模块多少时间算完的

冬雪雪冬

zooo 发表于 2016-5-12 23:16
使用pp模块多少时间算完的

加速不是很明显：
76576500.0 5.24573826789856

bacon6581

版主 你的电脑好快啊...

我这里需要2分多时间
你那里10秒钟就算完了....

老忘

晕，后面发现定义了一个函数但主程序没有使用，所以修改了一下，把那个函数删除了，结果没注意到却是在限定时间之后了

冬雪雪冬

老忘 发表于 2016-5-13 13:14
晕，后面发现定义了一个函数但主程序没有使用，所以修改了一下，把那个函数删除了，结果没注意到却是在限定 ...

嗯，截至时间后会取消回帖仅作者可见，为公平起见只好忍痛割爱了。

老忘

冬雪雪冬 发表于 2016-5-13 15:01
嗯，截至时间后会取消回帖仅作者可见，为公平起见只好忍痛割爱了。

嗯，支持，重在参与，动动手学到了知识就是最好的奖励了

冬雪雪冬

老忘 发表于 2016-5-13 15:11
嗯，支持，重在参与，动动手学到了知识就是最好的奖励了

是的，在这个过程中自己思考把题目做好，公布结果后又能看看别人的思路，这样对提高编程能力大有帮助。

小剑剑

学了太久c敲python时老是不自觉声明变量，不自觉加分号，明显写代码变得严谨，感觉对不起python啊。。

mdmouse

哎，我没有做出来

上山砍柴去

我大python威武雄壮啊

lorraine

xuexi