682|3

74

帖子

5

资源

一粒金砂(中级)

从《圣经密码》想到的数列问题及python解法

这几天,孩子的数学课外学习群里,老师不再发几何、组合、数论、应用题了,连续天天发数列题。

老师把好好的学习氛围搞得跟脑筋急转弯似的,咱也不敢说什么,只好到这里吐槽。当然,是有技术含量的吐。

今天的数列问题是 求 2、2、0、7、9、9的下一个数字是什么。

我们把标准答案和思路隐藏起来,您看完了本文再看:

 

游客,如果您要查看本帖隐藏内容请回复

 

很多年前,有一本书热过那么一阵,书名叫《圣经密码》,内容大概是说,经过逐一观察,发现圣经中预言了很多事件,形式是在圣经的某一页上,通过竖着、斜着、跳着,以及不知道怎么编出的规则数出的单个字母按顺序合在一起,就是某些历史事件的人名、地名等等。有兴趣的可以看下OSGeo中国中心的“圣经密码基本概况”(https://www.osgeo.cn/post/11250)。,这里不做详述。但是有一点是明确的——长度足够的随机数列,非常有可能获得你想要的某个较短数列(这个展开了得看哥德尔怎么搞垮数学,也不啰嗦了)。

于是,我准备找一个能够生成随机数列的“简单规则”来“猜测”任意脑筋急转弯数列的下一项。手头最顺手的语言是脚本语言python。
首先,我求出了1~1000的三角函数tan值的前12位,并对每一位上数字0~9进行了统计,发现出现的概率基本相同,这样就基本满足了分布要求。
然后我需要一个能够检测序列是否符合要求的工具,于是用状态机实现,_state表面上是统计现有的序列符合整体序列的长度,实际上当值等于标准序列长度时,表示发现了一个序列。
如果是测基因等序列,应当考虑同时检测等问题,但是吧我这只是要找出一个结果就可以,所以简化一下也能忍。
最终python类写成了这样:

class sq_auto_m():
    _seq=[]
    _full_flag= 0
    _stat =0
    _count=0
    _success=[]
    def __init__(self, seq):
        self._seq=seq[::]
        self.clean()
    
    def init(self,seq):
        self.__init__()
    
    def show_seq(self):
        print self._seq
    
    def show_cnt(self):
        print self._count
    
    def show_suc(self):
        return self._success
    
    def clean(self):
        self._full_flag= len(self._seq)
        self._stat =0
        self._count=0
        self._success=[]
    
    def push(self,val):
        self._count=self._count+1
        if self._seq[self._stat]==val:
            self._stat=self._stat+1
            if self._stat==self._full_flag:
                self._success.append(self._count)
                self._stat =0
                return True
        else:
            self._stat =0
            if self._seq[self._stat]==val:
                self._stat=self._stat+1
                #print self._stat
                if self._stat==self._full_flag:
                    self._success.append(self._count)
                    self._stat =0
                    return True
        
        return False

执行的时候,这样调用:

r=range(12)
qc=[ sq_auto_m([2,2,0,7,9,9]) for x in r]
i=1
while i< (1<<20):
    print i
    s=repr(abs(math.tan(i))).replace('.','')
    _=[ qc[x].push(int(s[x:x+1])) for x in r]
    i=i+1

for x in r:
    print x,qc[x].show_suc()

2的20次方超过一百万,我看着输出,跑到121516就不耐烦了,Ctrl+C强行终止,发现已经抛出两个结果

0 []
1 []
2 []
3 []
4 [115141]
5 []
6 [46758]
7 []
8 []
9 []
10 []
11 []

第6位发现的比较早,4万多就出现了。打印出来看看:

>>> for x in range(46753,46759+1):
...     repr(abs(math.tan(x)))[:15]
...
'0.1839028578233'
'1.0677023796878'
'3.9603504373298'
'0.3352375585022'
'0.8029491564140'
'9.4218595015405'
'0.5017617686735'

第6位(不算小数点)是2、2、0、7、9、9、1,所以答案是1?
很明显和正确答案不一样,但是也充分说明,没有真实环境可以检验的数列问题,只是脑筋急转弯。


此内容由EEWORLD论坛网友nemon原创,如需转载或用于商业用途需征得作者同意并注明出处

此帖出自创意市集论坛

回复

1万

帖子

133

资源

管理员

这个蛮有意思的~谢谢楼主分享~~

个人签名

玩板看这里:

http://bbs.eeworld.com.cn/elecplay.html

EEWorld测评频道众多好板等你来玩,还可以来频道许愿树许愿说说你想要玩的板子,我们都在努力为大家实现!


回复

185

帖子

0

资源

一粒金砂(中级)

脑壳疼

个人签名日拱一卒无有尽,功不唐捐终入海

回复

110

帖子

2

资源

一粒金砂(中级)

看看答案是啥 

 


回复
您需要登录后才可以回帖 登录 | 注册

关闭
站长推荐上一条 1/6 下一条

About Us 关于我们 客户服务 联系方式 器件索引 网站地图 最新更新 手机版

站点相关: 安防电子 汽车电子 手机便携 工业控制 家用电子 医疗电子 测试测量 网络通信 物联网

北京市海淀区知春路23号集成电路设计园量子银座1305 电话:(010)82350740 邮编:100191

电子工程世界版权所有 京ICP证060456号 京ICP备10001474号 电信业务审批[2006]字第258号函 京公海网安备110108001534 Copyright © 2005-2020 EEWORLD.com.cn, Inc. All rights reserved
快速回复 返回顶部 返回列表