码迷,mamicode.com
首页 > 其他好文 > 详细

特征根法求通项+广义Fibonacci数列找循环节 - HDU 5451 Best Solver

时间:2015-09-21 17:47:49      阅读:178      评论:0      收藏:0      [点我收藏+]

标签:

Best Solver

Problem‘s Link


 

Mean: 

给出x和M,求:(5+2√6)^(1+2x)的值。x<2^32,M<=46337.

analyse:

这题需要用到高中的数学知识点:特征根法求递推数列通项公式。

方法是这样的:

技术分享

技术分享

对于这题的解法:

记λ1=5+2√6,λ2=5-2√6,则λ1λ2=1,λ1+λ2=10

根据韦达定理可以推导出:λ1,λ2的特征方程为 x^2-10x+1=0

再使用该特征方程反向推导出递推公式为:a[n]=10*a[n-1]-a[n-2]

再由特征根法确定通项为:a[n]=(5+2√6)^n+(5-2√6)^n

观察通项,发现(5-2√6)^n<1,(5+2√6)^n>1。并且可以确定(5+2√6)^n的整数部分的值为a[n]-1

到这里,可以利用线性递推公式a[n]=10*a[n-1]-a[n-2],构造矩阵来找循环节。

为什么要找循环节呢?

因为n=2^x相当大,而模数M<=46337,对于每一个模数,都存在循环节。

而且我们发现这些循环节都比较小,所以我们可以直接暴力求矩阵循环节。

 

 

 

 

Time complexity: O(N)

 

Source code: 

 

特征根法在求递推数列通项中的运用

特征根法求通项+广义Fibonacci数列找循环节 - HDU 5451 Best Solver

标签:

原文地址:http://www.cnblogs.com/crazyacking/p/4826558.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!