应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler 工程 第三十七题
71/1返回列表
查看: 942  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-07-08 16:57:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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 ]
赞一下(1人) 回复此楼
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 ]
赞一下(3人) 回复此楼
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"和"[]" ...

看来昨天大脑缺氧了，哈哈。
赞一下(1人) 回复此楼
4楼2011-07-09 09:45:16
已阅   回复此楼   关注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测试的。
赞一下(1人) 回复此楼
5楼2011-07-09 10:02:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
xzhdty(金币+1): 欢迎常来程序语言 2011-07-10 08:17:05
引用回帖:
Originally posted by holmescn at 2011-07-09 10:02:16:
按说你这个pmlist相当于一个hash，应该比set要快啊。set显然是用的树，应该有二分法那个查找效率。

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

不过，为什么我的程序要比你的快一点点呢。我的大概要 1.3s， ...

我用timer测试的更快
CODE:
>>>
10000000 loops, best of 3: 0.0224 usec per loop

估计是我自己写的tic/toc函数开销太大?
赞一下(2人) 回复此楼
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
6楼2011-07-09 14:21:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lilac_c

至尊木虫 (知名作家)
xzhdty: 欢迎常来程序语言 2011-07-10 08:17:26
路过，闲看．．．．．．．．．
赞一下(1人) 回复此楼
我生活在一个经常爆发地震的年代
7楼2011-07-10 07:53:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 holmescn 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 材料与化工（0856）304求B区调剂 +7 邱gl 2026-03-10 11/550 2026-03-14 12:18 by 邱gl
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 云南财经大学信息学院计算机学硕专硕学位点 +3 zjptai 2026-03-10 5/250 2026-03-14 01:23 by 飞行琦
[考研] 求调剂，一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 工科，求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +5 18266118446 2026-03-09 5/250 2026-03-11 16:51 by 2020015
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
关闭
请填处理意见
关闭