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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler 工程 第三十七题
71/1返回列表
查看: 943  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +9 大米饭！ 2026-03-15 12/600 2026-03-16 21:59 by 大米饭！
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 9/450 2026-03-16 17:35 by 沐霖12138
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 材料工程专硕，一志愿中国矿业大学，总分314，求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 一志愿华中农业大学071010，总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
关闭
请填处理意见
关闭