找回密码
 欢迎注册
楼主: gxqcn

[分享] 找一个含有你生日字串的回文素数

[复制链接]
发表于 2008-11-26 20:04:44 | 显示全部楼层


mathematica代码?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2008-11-26 20:22:46 | 显示全部楼层
原帖由 mathematica 于 2008-11-26 19:52 发表
郭先强先生的那个代码真长,以至于我都不想看一下,虽然有注释。
不知道我的输出结果为什么比郭先强多出两个,疑惑中...........


我的是每个不同的居中数字仅搜索一个回文素数,搜到后即再找下一个居中数字的;
而你的是找指定位数以内的,所以以6、7居中的各重复了一次。

代码长短并非关键,关键是算法和效率。
并从中体验到创造的乐趣。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2008-11-26 20:27:47 | 显示全部楼层
很高兴论坛来了位 mathematica,
想必是对 mathematica 非常熟悉,
以后可以给大家多介绍介绍,
做到知己知彼。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2008-11-26 20:28:36 | 显示全部楼层
呵呵

多零的数字可发掘出好多有意思的东西
不过,只求概率素数难度太低了
条件成熟时
大家求确定素数去
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2008-11-26 21:15:46 | 显示全部楼层
mathematica,你用你装的mathematica计算一下PrimePi(80000000000000),看用时多少,这个对我计算PFP数有比较大的帮助,还请你帮帮忙算一下.
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2008-11-27 19:35:16 | 显示全部楼层

帮你计算一下

In[1]:=Timing[PrimePi[80000000000000]]
Out[1]:={80.453 Second, 2582444113487}
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2008-11-27 20:37:01 | 显示全部楼层
下表中列出了103个素数的序号,能否算出对应的素数?
即:如果给素数的序号是5,那么第5个素数是11。
1003290624000
1015831756800
1024192512000
1031956070400
1044855521280
1045094400000
1053455155200
1070176665600
1074954240000
1083553873920
1100753141760
1114512556032
1114767360000
1128701952000
1146617856000
1160950579200
1170505728000
1185137049600
1189085184000
1194393600000
1203948748800
1218998108160
1223059046400
1238347284480
1248539443200
1253826625536
1254113280000
1274019840000
1284211998720
1289945088000
1316818944000
1320903770112
1337720832000
1354442342400
1375941427200
1393140695040
1404606873600
1426902220800
1433272320000
1444738498560
1463132160000
1467670855680
1486016741376
1504935936000
1528823808000
1536288768000
1547934105600
1560674304000
1567283281920
1580182732800
1592524800000
1605264998400
1612431360000
1625330810880
1651129712640
1664719257600
1671768834048
1672151040000
1712282664960
1719926784000
1741425868800
1755758592000
1761205026816
1783627776000
1791590400000
1805923123200
1828497162240
1834588569600
1857520926720
1872809164800
1881169920000
1911029760000
1926317998080
1934917632000
1981355655168
2006581248000
2031663513600
2048385024000
2063912140800
2089711042560
2090188800000
2106910310400
2140353331200
2149908480000
2167107747840
2194698240000
2201506283520
2229025112064
2229534720000
2257403904000
2293235712000
2321901158400
2341011456000
2370274099200
2388787200000
2407897497600
2437996216320
2476694568960
2497078886400
2508226560000
2548039680000
2568423997440
2579890176000
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2008-11-28 07:52:19 | 显示全部楼层
昨晚加班(平时基本不用加班的,因有个急于上市的产品需要一特定算法),
试着去计算了一下楼上需要的数据,
结果电脑“僵”了好长时间也没有出结果,
最后只好关闭了。

刚才仅算一个:
In[1]:=Timing[Prime[2579890176000]]
Out[1]:={77.859 Second, 79918242752791}

也就是说,130个全部算完估计得 2.8 小时。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-23 17:31 , Processed in 0.023035 second(s), 14 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表