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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler 工程 第三十七题
北京石油化工学院2026年研究生招生接收调剂公告
71/1返回列表
查看: 971  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 总分293求调剂 +7 加一一九 2026-03-25 10/500 2026-03-29 16:53 by 杨杨杨紫
[考研] 085600 286分 材料求调剂 +8 麻辣鱿鱼 2026-03-27 9/450 2026-03-29 15:36 by Serene1974
[考研] 0856求调剂 +7 楒桉 2026-03-28 7/350 2026-03-29 08:28 by fmesaito
[考研] 0703化学 +11 妮妮ninicgb 2026-03-27 11/550 2026-03-29 06:45 by 544594351
[考研] 数一英一271专硕（085401）求调剂，可跨 +7 前行必有光 2026-03-28 8/400 2026-03-28 23:22 by 小木虫tim
[考研] 322求调剂 +7 宋明欣 2026-03-27 7/350 2026-03-28 21:27 by sanrepian
[考研] 求调剂 +3 QiMing7 2026-03-25 4/200 2026-03-28 14:30 by QiMing7
[考研] 308求调剂 +7 墨墨漠 2026-03-27 7/350 2026-03-28 07:43 by 热情沙漠
[考研] 085404求调剂，总分309，本科经历较为丰富 +4 来财aa 2026-03-25 4/200 2026-03-28 07:41 by 棒棒球手
[考研] 315分求调剂 +7 26考研上岸版26 2026-03-26 7/350 2026-03-28 04:05 by fmesaito
[考研] 考研调剂 +4 Sanmu-124 2026-03-26 4/200 2026-03-27 17:49 by kiokin
[考研] 一志愿南师大0703化学 275求调剂 +4 Ripcord上岸 2026-03-27 4/200 2026-03-27 17:00 by zhyzzh
[考研] 材料292调剂 +12 橘颂思美人 2026-03-23 12/600 2026-03-27 15:44 by caszguilin
[考研] 348求调剂 +4 小懒虫不懒了 2026-03-27 5/250 2026-03-27 12:47 by 果果妈咪
[考研] 求调剂 +3 刘柯@ 2026-03-24 4/200 2026-03-27 11:28 by shangxh
[考研] 一志愿吉大071010，316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
[考研] 351求调剂 +4 麦克阿磊 2026-03-24 4/200 2026-03-27 00:32 by wxiongid
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 材料调剂 +3 iwinso 2026-03-23 3/150 2026-03-25 11:29 by greychen00
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
信息提示
关闭
请填处理意见
关闭