鱼C论坛

 找回密码
 立即注册
查看: 3650|回复: 15

[技术交流] 小练习结果揭晓:第一个拥有超过500个约数的三角形数是多少?

[复制链接]
发表于 2016-5-12 21:41:43 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
首先给大家道个歉,这两天事情较多,没能抽出整块时间看大家的答案,一直拖到今天晚上。
根据题目要求
需在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宝贝归来 运算时间过长



我写的程序如下,供参考:
from math import sqrt
from itertools import combinations

def tri(n): #三角形数
    return n * (n + 1) //  2

def length(mylist):
    len1 = 0
    for i in range(len(mylist)):
        a =combinations(mylist, i + 1)
        b = list(set(a))
        len1 += len(b)
    return len1 + 2

#生成质数列表
primeList = [2]
item = 101
n = 3
while len(primeList) < item:
    for i in primeList:
        if n % i == 0: break
        if i > sqrt(n):
            primeList.append(n)
            break
    n += 2

n = 15
while True:
    list1 = []
    num = tri(n)
    num1 = num
    for i in primeList:
        while True:
            if num1 % i == 0:
                num1 //= i
                list1.append(i)
            else:
                break
        if num1 == 1:
            break
    if length(list1) > 500:
        print(num)
        break
    n += 1


本帖被以下淘专辑推荐:

想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复

使用道具 举报

发表于 2016-5-12 21:49:24 | 显示全部楼层
目测版主的电脑使用pp模块1秒左右能出结果。。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2016-5-12 22:05:39 | 显示全部楼层
zooo 发表于 2016-5-12 21:49
目测版主的电脑使用pp模块1秒左右能出结果。。

我下载pp-1.6.5.exe,可是安装出错(python3.5.1),等有时间把python再重新装一下试试。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2016-5-12 22:23:34 | 显示全部楼层
又下载pp-1.6.4.4.zip,安装成功。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2016-5-12 22:28:38 | 显示全部楼层
程序老是败在算法。。。。一直是不查资料直接做题,看样子还得多补充下基本的数学算法
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2016-5-12 23:16:00 | 显示全部楼层
冬雪雪冬 发表于 2016-5-12 22:23
又下载pp-1.6.4.4.zip,安装成功。

使用pp模块多少时间算完的
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2016-5-12 23:22:08 | 显示全部楼层
本帖最后由 冬雪雪冬 于 2016-5-13 13:36 编辑
zooo 发表于 2016-5-12 23:16
使用pp模块多少时间算完的


加速不是很明显:
76576500.0 5.24573826789856
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2016-5-13 12:06:50 | 显示全部楼层
版主 你的电脑好快啊...
QQ截图20160513120526.jpg
我这里需要2分多时间
你那里10秒钟就算完了....
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2016-5-13 13:14:55 | 显示全部楼层
晕,后面发现定义了一个函数但主程序没有使用,所以修改了一下,把那个函数删除了,结果没注意到却是在限定时间之后了
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

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

嗯,截至时间后会取消回帖仅作者可见,为公平起见只好忍痛割爱了。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2016-5-13 15:11:28 | 显示全部楼层
冬雪雪冬 发表于 2016-5-13 15:01
嗯,截至时间后会取消回帖仅作者可见,为公平起见只好忍痛割爱了。

嗯,支持,重在参与,动动手学到了知识就是最好的奖励了
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2016-5-13 15:20:40 | 显示全部楼层
老忘 发表于 2016-5-13 15:11
嗯,支持,重在参与,动动手学到了知识就是最好的奖励了

是的,在这个过程中自己思考把题目做好,公布结果后又能看看别人的思路,这样对提高编程能力大有帮助。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2016-5-13 22:05:23 | 显示全部楼层
学了太久c敲python时老是不自觉声明变量,不自觉加分号,明显写代码变得严谨,感觉对不起python啊。。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2016-5-13 23:27:45 | 显示全部楼层
哎,我没有做出来
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2016-5-14 12:57:54 | 显示全部楼层
我大python威武雄壮啊
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2016-5-16 07:40:16 | 显示全部楼层
xuexi
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2024-11-24 10:57

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表