PE440
一些证明,推荐复制入atom观看
首先我们考虑这个T(n)是什么,我们可以列出递归式:
(definition:T)
T(0) = 1
T(1) = 10
T(n) = 10*T(n-1) + T(n-2)
这个数列与Fibonacci数列有着类似的性质,我们考虑如何计算它某两项的gcd.对于Fibonacci数列我们有:
[此下为Fibonacci数列的性质,关系不大]{
(definition:F)
F(0) = 0
F(1) = 1
F(n) = F(n-1) + F(n-2)
(property:F)
1. F(x) Coprime F(x+1)
2. F(n) Divides F(k*n)
3. gcd(F(x),F(y)) = F(gcd(x,y))
我们来证明一下这些性质,首先是 [1.] 证明:辗转相除显然
这时我们可以发现,只要能证明 [3.] ,就可以推出 [2.] ,具体方法:
(proof:[2.])
假设命题 [3.] 为真:
gcd(F(x),F(k*x))=F(gcd(x,k*x))=F(x)->[2.]
我们考虑如何证明[3.], 我们发现有命题[4.]
(property:F)
4. F(a+b)=F(a)*F(b+1)+F(a-1)*F(b)
(proof:[4.])
考虑对b归纳,
basecase:
F(a+0)=F(a)
F(a+1)=F(a)+F(a-1)
assume forall 0<=b<=t is true, extend to b+1:
F(a+b )=F(a)*F(b+1)+F(a-1)*F(b)
F(a+b-1)=F(a)*F(b )+F(a-1)*F(b-1)
F(a+b+1)=F(a)*F(b+2)+F(a-1)*F(b+1) -> 符合命题
(proof:[3.])
gcd(F(a),F(a+t))=
gcd(F(a),F(a)*F(t+1)+F(a+1)*F(t))=
gcd(F(a),F(a+1)*F(t))=
gcd(F(a),F(t))
应用辗转相减法可得原命题
}
我们现在可以寻找这个T(n)的性质了,仿照寻找F(n)的性质(具体在上面),首先我们要找到T(n)关于很大的T(n+m)的关系.设
(equation:T1)
T(n+m)=T(n)*A(m)+T(n-1)*B(m)
考虑待定系数法
(derivation:T1)
T(n )=T(n)
A(0)=1
B(0)=0
T(n+1)=T(n)*A(1)+T(n-1)*B(1)
A(1)=10
B(1)=1
T(n+m+1)=
T(n+m)*10+T(n+m-1)=
T(n)*A(m)*10+T(n-1)*B(m)*10+T(n)*A(m-1)+T(n-1)*B(m-1)=
T(n)*(A(m)*10+A(m-1))+T(n-1)*(B(m)*10+B(m-1))
A(m+1)=A(m)*10+A(m-1)
B(m+1)=B(m)*10+B(m-1)
<implies that>
A(m)=T(m)
B(m)=T(m-1)
<reduce to>
T(n+m)=T(n)*T(m)+T(n-1)*T(m-1)
(equation:B1)
B(n+m)=T(n+m-1)=T(n)*T(m-1)+T(n-1)*T(m-2)=B(n+1)*B(m)+B(n)*B(m-1)
于是发现了与F(n)完全相同的形式,直接沿用F(n)的证明即可推出gcd(B(x),B(y))=B(gcd(x,y))也即gcd(T(n),T(m))=T(gcd(n+1,m+1)-1)
那么剩下的就是求gcd(c^a+1,c^b+1)(b>a), 我们可以考虑gcd(c^b-c^a,c^b+1)., c^b-c^a = c^(b-a)*(c^(b-a)-1) 显然 c^(b-a) Coprime c^b+1,
gcd(c^a+1,c^b+1)=gcd(c^b-c^a,c^a+1)=gcd(c^|b-2a|-1,c^a+1) 然后就变成了一个类似辗转相减的东西(我还是不很会整个分析),但是打个表就可以知道这东西值是
c^gcd(a,b)+1 a/gcd(a,b),b/gcd(a,b) odd otherwise
2 c even otherwise
1 c odd
然后这道题就做完了吧
PE440的更多相关文章
- ProjectEuler && Rosecode && Mathmash做题记录
退役选手打发时间的PE计划 挂在这里主要是dalao们看到有什么想交流的东西可以私聊哦(站内信或邮箱吧) 2017/8/11 PE595 :第一题QAQ 2017/8/12 PE598 2017/ ...
随机推荐
- Editplus常用快捷键
EditPlus 快捷键的使用 如果一个来你们公司面试程序员,连Ctrl + C 和Ctrl + V 都不用,而是使用“选中文本”→ 鼠标右键 → [复制],然后再鼠标右键→ [粘贴].你会不会录用他 ...
- C 语言学习 第五次作业总结
第五次作业,主要学习和复习的是几种循环结构的使用. 在前一次的课堂上,同学们已经学习了分支语句的使用.分支语句和循环语句配合使用,就可以写出更多的,逻辑功能丰富的代码了. 逻辑功能的丰富,也意味着学习 ...
- go 聊天室简单版总结
/* * 思路:在登录成功时将用户的id存进在线用户列表中的key value中链接的ws为空,并保存用户的信息. * 当跳转到聊天室时,将用户和聊天室链接的ws存进在线用户列表中的 * 问题:如何在 ...
- Android 线程池(转)
new Thread的弊端及Java四种线程池的使用: 1.new Thread的弊端 ); } catch (InterruptedException e) { // TODO Auto-gener ...
- jQuery链式操作[转]
用过jQuery的朋友都知道他强大的链式操作,方便,简洁,易于理解,如下 $("has_children").click(function(){ $(this).addClass( ...
- 关于解决asp.net mvc网站页面Banner图片即时更换css里背景图片url相对路径问题的新方案
最近在网站首页上想将Banner壁纸给做到后台上传随时更改的效果.遇到问题便是:将上传的图片路径动态添加到首页css代码中,结果尝试了网上提供的思路,更改相对路径,换为url中“../../Conte ...
- 解决:Win 10安装软件时提示:文件系统错误 (-1073740940)
1.win+R输入 gpedit.msc 2.左边计算机配置 windows设置——安全设置——本地策略——安全选项 3.在安全选项右边选择 用户账户控制:管理员批准模式中管理员的提升权限提示的行为, ...
- Leetcode 94. Binary Tree Inorder Traversal
Given a binary tree, return the inorder traversal of its nodes' values. For example:Given binary tre ...
- BZOJ4540 [Hnoi2016]序列
本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000作者博客:http://www.cnblogs.com/ljh2000-jump/转 ...
- please wait while windows configures microsoft visual studio professional 2013 [转载]
2016年5月30日 VS3013安装update 2以后,打开offie出现"please wait while windows configures microsoft visual s ...