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

简介: Best Solver Problem's Link   Mean:  给出x和M,求:(5+2√6)^(1+2x)的值。x

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

根据λ1=5+2√6,λ2=5-2√6是特征方程x^2-10x+1=0的解,可以求出:b=-10,c=1

再使用该特征方程反向推导出递推公式为: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,对于每一个模数M,都存在循环节,这样的话我们就不需要完整的计算n次幂运算了。

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

Time complexity: O(N+logx)

 

view code

 

目录
相关文章
|
8月前
迭代法求一元三次方程
迭代法求一元三次方程
97 0
|
8月前
PTA-求平方与倒数序列的部分和
求平方与倒数序列的部分和
70 1
【动态规划篇】斐波那契数列&&拆分词句&&三角矩阵
【动态规划篇】斐波那契数列&&拆分词句&&三角矩阵
【动态规划篇】斐波那契数列&&拆分词句&&三角矩阵
|
算法
【递归与递推 3】AcWing 717. 简单斐波那契(求解斐波那契数列的若干方法)
【递归与递推 3】AcWing 717. 简单斐波那契(求解斐波那契数列的若干方法)
107 0
【递归与递推 3】AcWing 717. 简单斐波那契(求解斐波那契数列的若干方法)
|
机器学习/深度学习 人工智能 算法
CF1446D Frequency Problem(思维 前缀和 根号分治 双指针)
CF1446D Frequency Problem(思维 前缀和 根号分治 双指针)
98 0
074.K阶斐波那契序列
074.K阶斐波那契序列
69 0
7-1 一元多项式求导 (10 分)
7-1 一元多项式求导 (10 分)
117 0
|
算法 C语言 C++
【算法练习】迭代法求平方根
【算法练习】迭代法求平方根
【算法练习】迭代法求平方根
|
机器学习/深度学习 算法
卡特兰数(Catalan Number) 算法、数论 组合~
卡特兰数(Catalan Number) 算法、数论 组合~
235 0
卡特兰数(Catalan Number) 算法、数论 组合~
|
知识图谱
【组合数学】生成函数 ( 使用生成函数求解不定方程解个数 )
【组合数学】生成函数 ( 使用生成函数求解不定方程解个数 )
313 0