一、矩陣乘法

https://blog.csdn.net/weixin_43272781/article/details/82899737

二、矩陣構造方法?

Fibonacci數列：F(0)=1 , F(1)=1 , F(n)=F(n-1)+F(n-2)

我們以前快速求Fibonacci數列第n項的方法是 構造常系數矩陣

（一）?? Fibonacci數列f[n]=f[n-1]+f[n-2],f[1]=f[2]=1的第n項快速求法（不考慮高精度）

解法：

考慮1×2的矩陣【f[n-2],f[n-1]】。根據Fibonacci數列的遞推關系，我們可以通過乘以一個2×2的矩陣A，得到矩陣：【f[n-1],f[n]】。

即：【f[n-2],f[n-1]】*A = 【f[n-1],f[n]】=【f[n-1],f[n-1]+f[n-2]】

很容易構造出這個2×2矩陣A，即：
０ １?
１ １

所以，有【f[1],f[2]】×A＝【f[2],f[3]】
又因為矩陣乘法滿足結合律，故有：
【f[1],f[2]】×A ^(n-1) =【f[n],f[n+1]】
這個矩陣的第一個元素f[n]即為所求。

?

（二）?? 數列f[n]=f[n-1]+f[n-2]+1,f[1]=f[2]=1的第n項的快速求法（不考慮高精度）

解法：
仿照前例，考慮1×3的矩陣【f[n-2],f[n-1],1】，希望求得某3×3的矩陣A，使得此1×3的矩陣乘以A得到矩陣：【f[n-1],f[n],1】

即：【f[n-2],f[n-1],1】* A ＝【f[n-1],f[n],1】=【f[n-1],f[n-1]+f[n-2]+1,1】

容易構造出這個3×3的矩陣A，即：
０ １ ０?
１ １ ０?
０ １ １

故：【f[1],f[2],1】* A^(n-1) = 【f[n],f[n+1],1】

?

（三）數列f[n]=f[n-1]+f[n-2]+n+1,f[1]=f[2]=1的第n項的快速求法（不考慮高精度）.
解法：
仿照前例，考慮1×4的矩陣【f[n-2],f[n-1],n,1】，希望求得某4×4的矩陣A，使得此1×4的矩陣乘以A得到矩陣：【f[n-1],f[n],n+1,1】
即：【f[n-2],f[n-1],n,1】* A? = 【f[n-1],f[n],n+1,1】＝【f[n-1],f[n-1]+f[n-2]+n+1,n+1,1】
容易構造出這個4×4的矩陣A，即：
０ １ ０ ０?
１ １ ０ ０?
０ １ １ ０?
０ １ １ １

故：【f[1],f[2],3,1】* A^(n-1) = 【f[n],f[n+1],n+2,1】

?

?

（四）?? 數列f[n]=f[n-1]+f[n-2],f[1]=f[2]=1的前n項和s[n]=f[1]+f[2]+……+f[n]的快速求法（不考慮高精度）.

解法：

仿照之前的思路，考慮1×3的矩陣【f[n-2],f[n-1],s[n-2]】，我們希望通過乘以一個3×3的矩陣A，得到1×3的矩陣：【f[n-1],f[n],s[n-1]】
即：【f[n-2],f[n-1],s[n-2]】 * A? = 【f[n-1],f[n],s[n-1]】=【f[n-1],f[n-1]+f[n-2],s[n-2]+f[n-1]】
容易得到這個3×3的矩陣A是：
０ １ ０?
１ １ １?
０ ０ １

這種方法的矩陣規模是(r+1)*(r+1)

f(1)=f(2)=s(1)=1 ，所以，有

【f(1),f(2),s(1)】* A? = 【f(2),f(3),s(2)】

故：【f(1),f(2),s(1)】* A^(n-1)? = 【f(n),f(n+1),s(n)】

?

（五）?? 數列f[n]=f[n-1]+f[n-2]+n+1,f[1]=f[2]=1的前n項和s[n]=f[1]+f[2]+……+f[n]的快速求法（不考慮高精度）.

解法：

考慮1×5的矩陣【f[n-2],f[n-1],s[n-2],n,1】,
我們需要找到一個5×5的矩陣A，使得它乘以A得到如下1×5的矩陣【f[n-1],f[n],s[n-1],n+1,1】
即：【f[n-2],f[n-1],s[n-2],n,1】* A? =【f[n-1],f[n],s[n-1],n+1,1】

=【f[n-1], f[n-1]+f[n-2]+n+1,s[n-2]+f[n-1],n+1,1】
容易構造出A為：
０ １ ０ ０ ０?
１ １ １ ０ ０?
０ ０ １ ０ ０?
０ １ ０ １ ０?
０ １ ０ １ １

故：【f(1),f(2),s(1),3,1】* A^(n-1)? = 【f(n),f(n+1),s(n),n+2,1】

?

一般地，如果有f[n]=p*f[n-1]+q*f[n-2]+r*n+s
可以構造矩陣A為：
0? q ?0 ?0 ?0?
1 ?p ?1 ?0 ?0?
0 ?0 ?1 ?0 ?0?
0 ?r ?0 ?1 ?0?
0 ?s ?0 ?1 ?1

?

更一般的，對于f[n]=Sigma(a[n-i]*f[n-i])+Poly(n)，其中0<i<=某常數c, Poly (n)表示n的多項式，我們依然可以構造類似的矩陣A來解決問題。
設Degree(Poly(n))=d, 并規定Poly(n)=0時，d=-1，此時對應于常系數線性齊次遞推關系。則本方法求前n項和的復雜度為：
((c+1)+(d+1))3*logns

?

例如：A(0) = 1 , A(1) = 1 , A(N) = X * A(N - 1) + Y * A(N - 2) (N >= 2)；給定三個值N，X，Y求S(N):S(N) = A(0)2 +A(1)2+……+A(n)2。

解：

考慮1*4 的矩陣【s[n-2],a[n-1]^2,a[n-2]^2,a[n-1]*a[n-2]】

我們需要找到一個4×4的矩陣A，使得它乘以A得到1×4的矩陣

【s[n-1],a[n]^2,a[n-1]^2,a[n]*a[n-1]】

即：【s[n-2],a[n-1]^2,a[n-2]^2,a[n-1]*a[n-2]】* A = 【s[n-1],a[n]^2,a[n-1]^2,a[n]*a[n-1]】

= 【s[n-2]+a[n-1]^2 , x^2 * a[n-1]^2 + y^2 * a[n-2]^2 + 2*x*y*a[n-1]*a[n-2] ,

a[n-1]^2 , x*a[n-1]^2 + y*a[n-2]a[n-1]】

可以構造矩陣A為：

1???? 0??? 0??? 0

1??? x^2?? 1??? x

0??? y^2?? 0??? 0

0??? 2xy?? 0??? y

故：【S[0],a[1]^2,a[0]^2,a[1]*a[0]】 * A^(n-1) = 【s[n-1],a[n]^2,a[n-1]^2,a[n]*a[n-1]】

所以：【S[0],a[1]^2,a[0]^2,a[1]*a[0]】 * A^(n) = 【s[n],a[n+1]^2,a[n]^2,a[n+1]*a[n]】

若A = (B * C ) 則AT?= ( B * C )T?= CT?* BT

故

三、矩陣快速冪

https://blog.csdn.net/weixin_43272781/article/details/85939539

四、例題

http://poj.org/problem?id=3070

http://poj.org/problem?id=3233

http://acm.hdu.edu.cn/showproblem.php?pid=2276

http://acm.hdu.edu.cn/showproblem.php?pid=5015

https://ac.nowcoder.com/acm/contest/338/L（題解：https://blog.csdn.net/weixin_43272781/article/details/85937715）

http://acm.hdu.edu.cn/showproblem.php?pid=3306?

五、參考文章

http://www.matrix67.com/blog/archives/276

?http://www.cnblogs.com/frog112111/archive/2013/05/19/3087648.html?tdsourcetag=s_pcqq_aiomsg

https://www.cnblogs.com/frog112111/archive/2013/05/15/3080678.html?

https://www.cnblogs.com/frog112111/archive/2013/05/16/3081493.html?

https://www.cnblogs.com/frog112111/archive/2013/05/16/3082416.html

https://www.cnblogs.com/frog112111/archive/2013/05/20/3089358.html

https://www.cnblogs.com/frog112111/archive/2013/05/20/3089687.html

https://www.cnblogs.com/frog112111/archive/2013/05/21/3089953.html

https://www.cnblogs.com/frog112111/archive/2013/05/22/3092012.html

https://www.cnblogs.com/frog112111/archive/2013/05/20/3089802.html

https://www.cnblogs.com/frog112111/archive/2013/05/22/3093697.html

https://www.cnblogs.com/frog112111/archive/2013/05/22/3093757.html?

總結

以上是生活随笔為你收集整理的矩阵构造方法的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。