POJ 3233——Matrix Power Series（矩阵快速幂）

最新推荐文章于 2021-08-15 09:28:36 发布
最新推荐文章于 2021-08-15 09:28:36 发布
阅读量273 收藏 1
点赞数
分类专栏： 矩阵快速幂
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_39206089/article/details/81275233
版权
这篇博客首次分享POJ 3233题目的矩阵快速幂解法。题目要求计算S=(A+A^2+A^3+...+A^k)%m，其中涉及矩阵乘法和大数据范围。解决方案是利用2n*2n矩阵进行扩展，并应用快速幂优化，以处理矩阵套矩阵的问题。
摘要由CSDN通过智能技术生成

第一次在博客里发题解（被水淹没不知所措）。。。

题目：POJ 3233 Matrix Power Series

题目大意：这道题不用翻译也能看懂，有一个n*n的矩阵，给出一个k,求S=（A+A^2+A^3+…+A^k）%m。其中n ≤ 30，k ≤ 10^9，m<104。

题目分析：这道题首先要知道矩阵乘法（不知道的请移步百度），然后一看数据范围就知道必须要用快速幂。但是这道题又特坑，矩阵套矩阵！那么就开一个2n*2n的矩阵即可。
矩阵大概就是长这样的：| A | 0 |
……………………………| E | E |
emmmmm。。。其实也不太难啦

代码：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
using<
最低0.47元/天 解锁文章
POJ 3233 Matrix Power Series（矩阵快速幂）
qq_36552817的博客
01-21 322
题目：3233 题意：给出一个方阵，求幂级数和，并对M取余 题解：采用矩阵快速幂，利用等比矩阵的性质 AC代码： #include<iostream> using namespace std; #define MAXN 61 int n, k, m; int A[MAXN][MAXN], B[MAXN][MAXN]; void mult(int a[MAXN][MAXN]...
【poj3233】Matrix Power Series——矩阵快速幂
CFhM_R's Blog
02-26 494
题目： Matrix Power Series Time Limit: 3000MS   Memory Limit: 131072K Total Submissions: 19029   Accepted: 8045 Description Given a n × n matrix A and a positive inte
【矩阵加速】[POJ3233]Matrix Power Series
weixin_30394981的博客
03-04 92
题目描述 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 posi...
【POJ 3233】Matrix Power Series（矩阵快速幂）
reverie_mjp的博客
10-26 379
既然难忘记，不如就让时光与回忆彼此相安
Matrix Power Series
xinxiaxindong的博客
08-16 269
直接暴力不行 因为每项底数都是A 可以用结合律 可以找到递推式 如果k是奇数 那么dfs(k) = dfs（k-1）+A^k 否则 dfs（k） = dfs(k/2)+A^k*(dfs(k/2))//先把dfs（k/2）算出来再带进去主要是 初始化一开始写成void函数的形式导致初始化失败 然后怎么改怎么wa 还有全局变量和局部变量冲突！！！！ 花式过样例 可你就是wa 为啥？为啥？哈
POJ3233Matrix Power Series 矩阵快速幂(分块矩阵构造)
cccrown的专栏
03-09 852
# include # include # include using namespace std; int mod; class Matrix{ public: int row; int col; int ** element; Matrix(); ~Matrix(); Matrix(const int , const
POJ 3233 Matrix Power Series(矩阵快速幂+二分)
scf0920
09-17 2649
题目地址：POJ 3233 题目大意：给定矩阵A，求A + A^2 + A^3 + … + A^k的结果（两个矩阵相加就是对应位置分别相加）。输出的数据mod m。k     这道题两次二分，相当经典。首先我们知道，A^i可以二分求出。然后我们需要对整个题目的数据规模k进行二分。比如，当k=6时，有：     A + A^2 + A^3 + A^4 + A^5 + A^6 =(A + A^2
变形斐波那契数列的和 I
欢迎您的光顾awa
12-12 160
变形斐波那契数列的和 I，矩阵乘法。 求 F[i]=F[i-2]+F[i-1]+n+1 这个函数前 n 个数的和。
Matrix Power Series（？）
yinghui_yht的博客
12-03 247
真的是，很多时候一个小错误看不出来，一觉睡醒就好了
Matrix Power Series(矩阵快速幂)
xdlove
02-09 1197
矩阵快速幂:http://www.cnblogs.com/kuangbin/archive/2012/08/17/2643347.html Matrix Power Series Time Limit: 3000MS   Memory Limit: 131072K Total Submissions: 16341   Accepted: 6966
[POJ3233] Matrix Power Series（矩阵快速幂）
weixin_30872157的博客
05-12 95
传送门 k <= 109 暴力肯定超时 根据矩阵性质，可以发现 S(4) = A1 + A2 + A2* (A1 + A2) S(5) = A1 + A2 +A2* (A1+ A2) + A5 所以可以递归二分求解，分别判断 Ak，k 为奇数和偶数的情况 ——代码 1 #include <cstdio> 2 #include &lt...
[矩阵加速] Matrix Power Series
烟尘的博客
04-16 377
文章目录题目描述解题思路参考代码 题目描述 Matrix Power Series 解题思路 先看题目，本题需要我们去求SSS的值，为了方便，不妨把A1+A2+…+AiA^1+A^2+…+A^iA1+A2+…+Ai记作S(i)S(i)S(i)（SSS和AAA都为矩阵） 再看数据范围，一个极大的变量，一个小变量，明显就是一道矩阵加速的题 很容易得到一个S(i)S(i)S(i)的递推式：S(i)=S...
Matrix Power Series(矩阵快速幂，二分法/构造矩阵套矩阵）
木静音yy的博客
08-15 254
传送门 题意： 思路: 方法一：（二分） 分析： 奇偶分类讨论： 1.代码： #include <iostream> #include <cstdio> #include <string> #include <cstring> #include <stdlib.h> #include <math.h> #include <ctype.h> #include <queue> #include <map&gt
用java解决POJ3233—矩阵幂序列问题
最新发布
11-19
以下是Java解决POJ3233—矩阵幂序列问题的代码和解释： ```java import java.util.Scanner; public class Main { static int n, k, m; static int[][] A, E; public static void main(String[] args) { Scanner sc = new Scanner(System.in); n = sc.nextInt(); k = sc.nextInt(); m = sc.nextInt(); A = new int[n][n]; E = new int[n][n]; for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { A[i][j] = sc.nextInt() % m; E[i][j] = (i == j) ? 1 : 0; } } int[][] res = matrixPow(A, k); int[][] ans = matrixAdd(res, E); printMatrix(ans); } // 矩阵乘法 public static int[][] matrixMul(int[][] a, int[][] b) { int[][] c = new int[n][n]; for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { for (int k = 0; k < n; k++) { c[i][j] = (c[i][j] + a[i][k] * b[k][j]) % m; } } } return c; } // 矩阵快速幂 public static int[][] matrixPow(int[][] a, int b) { int[][] res = E; while (b > 0) { if ((b & 1) == 1) { res = matrixMul(res, a); } a = matrixMul(a, a); b >>= 1; } return res; } // 矩阵加法 public static int[][] matrixAdd(int[][] a, int[][] b) { int[][] c = new int[n][n]; for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { c[i][j] = (a[i][j] + b[i][j]) % m; } } return c; } // 输出矩阵 public static void printMatrix(int[][] a) { for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { System.out.print(a[i][j] + " "); } System.out.println(); } } } ``` 解释： 1. 首先读入输入的n、k、m和矩阵A，同时初始化单位矩阵E。 2. 然后调用matrixPow函数求出A的k次幂矩阵res。 3. 最后将res和E相加得到结果ans，并输出。 4. matrixMul函数实现矩阵乘法，matrixPow函数实现矩阵快速幂，matrixAdd函数实现矩阵加法，printMatrix函数实现输出矩阵。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值