首页 > 其他分享 >常系数齐次线性递推学习笔记

常系数齐次线性递推学习笔记

时间:2023-02-24 17:35:28浏览次数:40  
标签:系数 frac log 笔记 齐次 线性 递推

求一个满足 \(k\) 阶齐次线性递推数列 \(a_i\) 的第 \(n\) 项,即:

\[f_n=\sum_{i=1}^ka_if_{n-i} \]

\(n \leq 10^{18},k\leq 32000\) 。

使用矩阵乘法加速可以做到 \(O(k^3\log n)\) ,运用这一项科技,可以做到 \(O(k\log n)\) 。

数列转化为生成函数

将该线性递推序列看作一个生成函数:

\[F(x)=\sum_{i=0}^{n}f_ix^{i} \]

将递推式也视为一个生成函数:

\[A(x)=\sum_{i=0}^ka_ix^i \]

显然有:

\[F(x)=F(x)A(x)+R(x) \]

其中 \(R(x)\) 为加以修正的常数项。进一步化式子,将其表示为分式形式。

\[F(x)=\frac{R(x)}{1-{A(x)}} \]

然后就可以用常系数齐次线性递推来优化了。\(\text{ps.}\) 分母那一坨东西也被叫做数列的特征方程。

常系数齐次线性递推

现在已经将所求数列转化成了一个生成函数的分式形式,所求即为:

\[[x^n]F(x)=\frac{P(x)}{G(x)} \]

接下来的一步非常巧妙!将其系数按照奇偶分开:

\[[x^n]F(x)=\frac{P(x)G(-x)}{G(x)G(-x)} \]

先看这个分母,将其视为一个整体显然是一个偶函数,因此其奇数项系数为 \(0\) 。将分子系数按照奇偶分开:

\[[x^n]F(x)=\frac{xO(x^2)+E(x^2)}{G(x^2)} \]

对 \(n\) 按照奇偶性分类讨论,如果是奇数的话保留左边,否则保留右边不断递归下去,直到 \(n=0\) 也就是仅剩常数项为止。

多项式乘法使用 \(\text{FFT}\) 优化复杂度可以做到 \(O(k \log n \log k)\) 。

标签:系数,frac,log,笔记,齐次,线性,递推
From: https://www.cnblogs.com/oscaryangzj/p/17152289.html

相关文章

  • 【笔记】springboot使用Spring-data-jpa
    Spring-data-jpaSpring-data-jpa依赖于Hibernate。通过整合Hibernate之后,我们以操作Java实体的⽅式最终将数据改变映射到数据库表中。添加依赖:<dependency<groupId......
  • 笔记
     判断datatable指定列是否重复varquery=fromtindtTemp.AsEnumerable()grouptbynew{t1=t.Field<string>("St......
  • Python学习笔记--网络通信--socket
    1.socket里面的,AF_INET和AF_UNIX有什么区别?AF_INET用于真实的两台机器进行通信。AF_UNIX用于本地自己跟自己通信。参考资料:http://www.langdebuqing.com/  2.soc......
  • MySQL数据库学习笔记1
    MySQL数据库学习笔记1MySQL服务器启动与连接#启动mysql.serverstart#连接mysql-uroot-pMySQL数据库的数据模型客户端访问MySQL数据库,是与数据库管理系统交......
  • SA 学习笔记
    前言这是我发布的第一篇博客,从学习SA开始。基本概念一些规定字符串的下标从1开始,\(n=|S|\),记\(suf(x)\)表示以\(s[x]\)开头的后缀,\(pre(x)\)表示以\(s[x]\)为结尾的前......
  • 《分布式技术原理与算法解析》学习笔记Day21
    分布式数据存储三要素什么是分布式数据存储系统?分布式存储系统的核心逻辑,就是将用户需要存储的数据根据某种规则存储到不同的机器上,当用户想要获取指定数据时,再按照规则......
  • uni-app学习笔记之----getCurrentPages()的使用
    1、判断是否是首页如果得到数组元素只有一个,说明是首页2、得到页面中的信息得到数组中的第一个元素代表首页,最后一个元素代表当前页 ......
  • 0003001第三章-灰度变换与空间滤波笔记
    目录第三章-灰度变换与空间滤波(空间域处理)第三章-灰度变换与空间滤波(空间域处理)空间域处理指的是对像素进行操作,用一个映射函数\(T(原像素)\)得到一个新像素,即:\(s=T(r)......
  • 【笔记】IDEA中maven导入依赖提示证书错误解决方法
    先是提示:一定要备份配置文件!!! 一定要备份配置文件!!! 一定要备份配置文件!!!先说原因:idea内置了jre,与你开发用的jre不是同一个软件,你通过命令修改的是开发用的jre的证书库,导入m......
  • 纸质笔记和电子笔记哪个好纸质笔记和电子笔记哪个好
    大家在日常的学习、工作中少不了要做笔记,把一些注意事项、重要的内容、容易忘记的事情记录下来,防止自己忘记。不过有一些网友习惯于使用纸质笔记,有的网友则更倾向于使用电......