当前位置: 首页 > news >正文

矩阵快速幂常用矩阵构造

斐波那契数列: \(F_i = F_{i-1} + F_{i-2}\)

\[\begin{bmatrix} F_i \\ F_{i-1} \end{bmatrix}= \begin{bmatrix} 1 & 1 \\ 1 & 0 \end{bmatrix} \times \begin{bmatrix} F_{i-1} \\ F_{i-2} \end{bmatrix} = \begin{bmatrix} 1 & 1 \\ 1 & 0 \end{bmatrix}^{i-1} \times \begin{bmatrix} F_1 \\ F_0 \end{bmatrix} \]

递推式 - 1: \(G_i = a \times G_{i-1} + b \times G_{i-2}\)

\[\begin{bmatrix} G_i \\ G_{i-1} \end{bmatrix} = \begin{bmatrix} a & b \\ 1 & 0 \end{bmatrix} \times \begin{bmatrix} G_{i-1} \\ G_{i-2} \end{bmatrix} = \begin{bmatrix} a & b \\ 1 & 0 \end{bmatrix}^{i-1} \times \begin{bmatrix} G_1 \\ G_0 \end{bmatrix} \]

递推式 - 2: \(G_i = a \times G_{i-1} + i^2\)

\[\begin{bmatrix} G_i \\ (i+1)^2 \\ i+1 \\ 1 \end{bmatrix} = \begin{bmatrix} a & 1 & 0 & 0 \\ 0 & 1 & 2 & 1 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & 0 & 1 \end{bmatrix} \times \begin{bmatrix} G_{i-1} \\ i^2 \\ i \\ 1 \end{bmatrix} = \begin{bmatrix} a & 1 & 0 & 0 \\ 0 & 1 & 2 & 1 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & 0 & 1 \end{bmatrix}^i \times \begin{bmatrix} G_0 \\ 1 \\ 1 \\ 1 \end{bmatrix} \]

递推式 - 3: \(G_i = a \times G_{i-1} + i^3\)

\[\begin{bmatrix} G_i \\ (i+1)^3 \\ (i+1)^2 \\ i+1 \\ 1 \end{bmatrix} = \begin{bmatrix} a & 1 & 0 & 0 & 0 \\ 0 & 1 & 3 & 3 & 1 \\ 0 & 0 & 1 & 2 & 1 \\ 0 & 0 & 0 & 1 & 1 \\ 0 & 0 & 0 & 0 & 1 \end{bmatrix} \times \begin{bmatrix} G_{i-1} \\ i^3 \\ i^2 \\ i \\ 1 \end{bmatrix} = \begin{bmatrix} a & 1 & 0 & 0 & 0 \\ 0 & 1 & 3 & 3 & 1 \\ 0 & 0 & 1 & 2 & 1 \\ 0 & 0 & 0 & 1 & 1 \\ 0 & 0 & 0 & 0 & 1 \end{bmatrix}^i \times \begin{bmatrix} G_0 \\ 1 \\ 1 \\ 1 \\ 1 \end{bmatrix} \]

递推式 - 4: \(G_i = a \times G_{i-1} + b^i\)

\[\begin{bmatrix} G_i \\ b^{i+1} \end{bmatrix} = \begin{bmatrix} a & 1 \\ 0 & b \end{bmatrix} \times \begin{bmatrix} G_{i-1} \\ b^i \end{bmatrix} = \begin{bmatrix} a & 1 \\ 0 & b \end{bmatrix}^i \times \begin{bmatrix} G_0 \\ b \end{bmatrix} \]

递推式 - 带常数:\(f(n) = a \times f(n-1) + b \times f(n-3) + c\)

\[\begin{bmatrix} f[n] \\ f[n-1] \\ f[n-2] \\ c \end{bmatrix} = \begin{bmatrix} a & 0 & b & 1 \\ 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 1 \end{bmatrix} \times \begin{bmatrix} f[n-1] \\ f[n-2] \\ f[n-3] \\ c \end{bmatrix} = \begin{bmatrix} a & 0 & b & 1 \\ 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 1 \end{bmatrix}^{n-2} \times \begin{bmatrix} f[2] \\ f[1] \\ f[0] \\ c \end{bmatrix} \]

递推式 - 带变量:\(f(1)=1\)\(f(2)=2\)\(f(n) = f(n-1) + 2f(n-2) + n^3\)

利用二项式定理:

\[n^3 = (n-1+1)^3 = (n-1)^3 + 3(n-1)^2 + 3(n-1) + 1 \]

\[\begin{bmatrix} f(n) \\ f(n-1) \\ n^3 \\ n^2 \\ n \\ 1 \end{bmatrix} = \begin{bmatrix} 1 & 2 & 1 & 3 & 3 & 1 \\ 1 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 3 & 3 & 1 \\ 0 & 0 & 0 & 1 & 2 & 1 \\ 0 & 0 & 0 & 0 & 1 & 1 \\ 0 & 0 & 0 & 0 & 0 & 1 \end{bmatrix} \times \begin{bmatrix} f(n-1) \\ f(n-2) \\ (n-1)^3 \\ (n-1)^2 \\ n-1 \\ 1 \end{bmatrix} = \begin{bmatrix} 1 & 2 & 1 & 3 & 3 & 1 \\ 1 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 3 & 3 & 1 \\ 0 & 0 & 0 & 1 & 2 & 1 \\ 0 & 0 & 0 & 0 & 1 & 1 \\ 0 & 0 & 0 & 0 & 0 & 1 \end{bmatrix}^{n-2} \times \begin{bmatrix} f(2) \\ f(1) \\ 8 \\ 4 \\ 2 \\ 1 \end{bmatrix} \]

递推式 - 前缀和:

递推式定义:

  • \( T[0]=T[1]=T[2]=1 \)
  • \( T[n]=T[n-1]+T[n-2]+T[n-3] \ (n \geq 3) \)
  • \( S[n] = T[0] + T[1] + \dots + T[n] = S[n-1] + T[n-1]+T[n-2]+T[n-3] \)

\[\begin{bmatrix} S[n] \\ T[n] \\ T[n-1] \\ T[n-2] \end{bmatrix} = \begin{bmatrix} 1 & 1 & 1 & 1 \\ 0 & 1 & 1 & 1 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \end{bmatrix} \cdot \begin{bmatrix} S[n-1] \\ T[n-1] \\ T[n-2] \\ T[n-3] \end{bmatrix} = \begin{bmatrix} 1 & 1 & 1 & 1 \\ 0 & 1 & 1 & 1 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \end{bmatrix}^{n-2} \cdot \begin{bmatrix} S[2] \\ T[2] \\ T[1] \\ T[0] \end{bmatrix} \]

http://www.jsqmd.com/news/26906/

相关文章:

  • 新学期每日总结(第17天)
  • 顶级CTF工具与资源大全
  • 第二章 数列极限
  • 小白也能看懂的RL-PPO
  • 第二十三天
  • ICPC2022南京 游记(VP)
  • [KaibaMath]1015 关于收敛数列迫敛性的证明
  • Manancher
  • 搜维尔科技:【技术分享】解析Xsens动捕与人形机器人的训练术语
  • Python while循环 _ 捕捉日落
  • 搜维尔科技:IROS 2025圆满落幕|MANUS手套展示世界级手部追踪技术,从遥操作到具身智能!
  • 2024 暑期模拟赛 #9
  • 三值纠缠模型:智能价值权衡的元能力与实现路径探索
  • 三值纠缠模型:智能价值权衡的元能力与实现路径探索
  • OceanBase系列---【如何拆分PMAX分区?】
  • AutoDL+Deepseek 7B
  • VLP平台与重组蛋白:新一代生物技术工具
  • 2025.10.30
  • 10/30
  • 实验任务3
  • 会计的职能 - 智慧园区
  • [CEOI 2020] 星际迷航
  • Chome插件Mathpix Snip对SDU信息服务平台的会话阻塞问题
  • 2025.10.30总结
  • AT_arc068_d [ARC068F] Solitaire 分析
  • 10/30观后感
  • 20251030周四日记
  • 手写汉字识别
  • Keil仿真条件断点10.30
  • 10.30 程序员的修炼之道:从小工到专家第三章 基本工具 - GENGAR