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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler 工程 第三十七题
71/1返回列表
查看: 987  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +15 熊二想上岸 2026-04-04 15/750 2026-04-08 22:52 by 猪会飞
[考研] 化学调剂求助 +13 LULONG1 2026-04-03 18/900 2026-04-08 21:54 by LULONG1
[考研] 269求调剂 +6 啊啊我我 2026-04-07 6/300 2026-04-08 20:04 by 我减肥1
[考研] 考研调剂 +11 冰冰，，， 2026-04-07 11/550 2026-04-08 15:08 by Delta2012
[考博] 申博 +8 IQwQl 2026-04-04 8/400 2026-04-08 09:43 by 0608104024
[考研] 362求调剂 +8 我要考大 2026-04-06 12/600 2026-04-08 00:32 by JourneyLucky
[考研] 生物工程求调剂 +13 喜欢还是不甘心 2026-04-05 13/650 2026-04-07 16:55 by Ecowxq666！
[考研] 286求调剂 +20 Faune 2026-04-06 20/1000 2026-04-07 11:33 by 诗与自由
[考研] 一志愿北交大材料工程总分358求调剂 +10 cs0106 2026-04-05 12/600 2026-04-06 19:41 by 无际的草原
[考研] 生物与医药求调剂 +7 heguanhua 2026-04-05 8/400 2026-04-06 18:41 by macy2011
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 296求调剂 +3 汪！？！ 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 329求调剂 +17 miaodesi 2026-04-02 20/1000 2026-04-05 18:33 by 蓝云思雨
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 309求调剂 +6 刘刘刘1231 2026-04-02 7/350 2026-04-04 13:41 by liucky
[考研] 一志愿沪985，326分求调剂 +3 刘墨墨 2026-04-03 3/150 2026-04-04 11:16 by 悲伤的芋头
[考研] 400分求调剂 +3 尴尬且挠头 2026-04-04 3/150 2026-04-04 08:41 by jp9609
[考研] 311求调剂 +11 勇敢的小吴 2026-04-02 11/550 2026-04-03 21:46 by qlm5820
[考研] 281求调剂 +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[考研] 289-求调剂 +4 这里是_ 2026-04-03 4/200 2026-04-03 14:23 by 1753564080
信息提示
关闭
请填处理意见
关闭