POJ3233 Matrix Power Series 矩阵乘法】的更多相关文章

http://poj.org/problem?id=3233 挺有意思的..学习到结构体作为变量的转移, 题意 : 给定矩阵A,求A + A^2 + A^3 + ... + A^k的结果(两个矩阵相加就是对应位置分别相加).输出的数据mod m.k<=10^9. http://blog.csdn.net/rowanhaoa/article/details/21024093 代码: #include<cstdio> #include<cstring> #include<i…
Matrix Power Series Time Limit: 3000MS Memory Limit: 131072K Description Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + - + Ak. Input The input contains exactly one test case. The first line of input contains three po…
Matrix Power Series Time Limit: 3000MS   Memory Limit: 131072K Total Submissions: 27277   Accepted: 11143 Description Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak. Input The input contains exactly one test ca…
解题思路 题目里要求\(\sum_{i=1}^kA^i\),我们不妨再加上一个单位矩阵,求\(\sum_{i=0}^kA^i\).然后我们发现这个式子可以写成这样的形式:\(A(A(A...)+E)+E)+E\)于是,我们可以将\(*A+E\)看做一次变换,然后尝试构造一个矩阵.我们发现: \[ (\left[ \begin{matrix} A & E \\ 0 & E \end{matrix} \right])^n= \left[ \begin{matrix} A^{n+1} &…
http://poj.org/problem?id=3233 题目大意:给定矩阵A,求A + A^2 + A^3 + … + A^k的结果(两个矩阵相加就是对应位置分别相加).输出的数据mod m.k<=10^9.这道题两次二分,相当经典.首先我们知道,A^i可以二分求出.然后我们需要对整个题目的数据规模k进行二分.比如,当k=6时,有:A + A^2 + A^3 + A^4 + A^5 + A^6 =(A + A^2 + A^3) + A^3*(A + A^2 + A^3)应用这个式子后,规模…
Description Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak. Input The input contains exactly one test case. The first line of input contains three positive integers n (n ≤ 30), k (k ≤ 109) and m (m < 104). Then…
传送门 题意 给出n,m,k,求 \[\sum_{i=1}^kA^i\] A是矩阵 分析 我们首先会想到等比公式,然后得到这样一个式子: \[\frac{A^{k+1}-E}{A-E}\] 发现要用矩阵除法,可以用求矩阵逆来做,现在我们换一种做法,我们发现有这样一个性质: \[\left[ \begin{matrix} A & E \\ 0 & E \\ \end{matrix} \right]^n= \left[ \begin{matrix} A^{n} & \sum_{i=0}…
S(k)=A^1+A^2...+A^k. 保利求解就超时了,我们考虑一下当k为偶数的情况,A^1+A^2+A^3+A^4...+A^k,取其中前一半A^1+A^2...A^k/2,后一半提取公共矩阵A^k/2后可以发现也是前一半A^1+A^2...A^k/2.因此我们可以考虑只算其中一半,然后A^k/2用矩阵快速幂处理.对于k为奇数,只要转化为k-1+A^k即可.n为矩阵数量,m为矩阵大小,复杂度O[(logn*logn)*m^3] #include <iostream> #include &…
本文为博主原创文章,欢迎转载,请注明出处 www.cnblogs.com/yangyaojia [POJ3233]Matrix Power Series 分治+矩阵 题目大意 A为n×n(n<=30)的矩阵,让你求 \(\sum\limits_{i=1}^{k}A^i\) 并将答案对取模p 输入格式: 有多组测试数据,其中第一行有3个正整数,为n,k(k<=\(10^9\)),p(p<=\(10^4\)) 后面有n行,每行n个数. 输出格式: 输出最后答案的矩阵. 输入输出样例 inpu…
Matrix Power Series Time Limit: 3000MS   Memory Limit: 131072K Total Submissions: 23187   Accepted: 9662 Description Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak. Input The input contains exactly one test cas…