24小时热门版块排行榜    

查看: 950  |  回复: 6
本帖产生 2 个 程序强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

holmescn

金虫 (正式写手)

[交流] Euler 工程 第三十七题 已有2人参与

3797这个数很有意思,它本身是一个质数。同时,无论你是从左边删除数字,还是从右边删除数字,它还是个质数。
比如从左边删除数字:3797,797,97,7
从右边删除数字:3797,379,37,3

这样的质数一共有11个。请给出他们的和。

注意:2,3,5,7不算在内。

[ Last edited by holmescn on 2011-7-9 at 09:47 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)


dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-12 16:15:16
引用回帖:
Originally posted by libralibra at 2011-07-09 00:18:33:
11个啊兄弟,你来个7个,我以为我搞错了
Find the sum of the only eleven primes that are both truncatable from left to right and right to left.

你结束值太小了,不过不确定"in"和"[]" ...

按说你这个pmlist相当于一个hash,应该比set要快啊。set显然是用的树,应该有二分法那个查找效率。

主要是生成1000万以下的质数用时太长了。

不过,为什么我的程序要比你的快一点点呢。我的大概要 1.3s,用time测试的。
5楼2011-07-09 10:02:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

holmescn

金虫 (正式写手)

★ ★ ★
余泽成(金币+3, 程序强帖+1): 鼓励交流! 2011-07-11 20:42:48
Python穷举
CODE:
# Project Euler Problem 37
#

# Gen a prime list
from math import sqrt
UBound = 100000
primes = range(2, UBound)
for x in primes:
    if x > 0:
        for n in xrange(2*x, UBound, x):
            primes[n-2] = 0

# change to a set to make "not in" faster
primes = set([x for x in primes if x > 0])

# Check each prime
for x in primes:
    if x > 10:
        s = str(x)
        flag = True
        for i in xrange(1, len(s)):
            # if any truncated number is not
            # a prime, break and set fail
            if int(s[i:]) not in primes or int(s[:len(s)-i]) not in primes:
                flag = False
                break
        if flag:
            print x

23, 37, 53, 73, 313, 317, 373, 797, 3137, 3797, 739397

和是748317

[ Last edited by holmescn on 2011-7-9 at 09:52 ]
2楼2011-07-08 23:22:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
xzhdty(金币+1): 欢迎常来程序语言看看 2011-07-09 07:01:26
余泽成(金币+2, 程序强帖+1): 欢迎常来程序语言版! 2011-07-11 20:43:27
11个啊兄弟,你来个7个,我以为我搞错了
Find the sum of the only eleven primes that are both truncatable from left to right and right to left.

你结束值太小了,不过不确定"in"和"[]"哪个快?
生成素数列表我没有删除非素数,直接置0,下面用[]>0判断的
发现一个可怕的事实:vim运行python竟然比IDLE快,这个程序快1s
IDLE需要2.9s左右
CODE:
11 748317
[23, 37, 53, 73, 313, 317, 373, 797, 3137, 3797, 739397]
Elapsed time: 1.97395228 seconds

CODE:
#! /usr/bin/env python

from mytictoc import tic, toc

tic()

# produce prime list
stop = 1000000
pmlist = range(stop+1)
pmlist[:2] = [0,0]
i = 2
while i*i<=stop:
    if pmlist[i]>0:
        for j in xrange(2*i,stop+1,i):
            pmlist[j] = 0
    i += 1

nums = []
for i in xrange(23,stop+1):
    if pmlist[i]>0 and '0' not in str(i): # isprime and don't contain 0
        # get all sub number
        allnum = []
        for j in xrange(1,len(str(i))):
            allnum.extend([int(str(i)[:j]),int(str(i)[-j:])])
        # check all sub number
        if all([pmlist[x]>0 for x in allnum]):
            nums.append(i)

# print result
print len(nums),sum(nums)
print nums
        
toc()

[ Last edited by libralibra on 2011-7-9 at 01:05 ]
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
3楼2011-07-09 00:18:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)


dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-12 16:15:10
引用回帖:
Originally posted by libralibra at 2011-07-09 00:18:33:
11个啊兄弟,你来个7个,我以为我搞错了
Find the sum of the only eleven primes that are both truncatable from left to right and right to left.

你结束值太小了,不过不确定"in"和"[]" ...

看来昨天大脑缺氧了,哈哈。
4楼2011-07-09 09:45:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +8 Ncdx123456 2026-03-13 9/450 2026-03-18 14:40 by haxia
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 0703化学调剂 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
信息提示
请填处理意见