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

最新推荐文章于 2023-12-05 22:51:26 发布
原创 最新推荐文章于 2023-12-05 22:51:26 发布 · 343 阅读 · 1
这篇博客首次分享POJ 3233题目的矩阵快速幂解法。题目要求计算S=(A+A^2+A^3+...+A^k)%m，其中涉及矩阵乘法和大数据范围。解决方案是利用2n*2n矩阵进行扩展，并应用快速幂优化，以处理矩阵套矩阵的问题。

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

题目：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 namespace std;
const
最低0.47元/天 解锁文章
Windows计划任务自动化实战：从schtasks命令到运维脚本
最新发布
weixin_30709061的博客
05-13 239
本文详细解析了Windows计划任务自动化实战，重点介绍了schtasks命令的使用技巧和运维脚本编写方法。通过实际案例展示如何从基础单次任务进阶到企业级自动化方案，涵盖批处理脚本封装、PowerShell动态生成任务等高级技巧，并提供了权限管理、时间参数避坑等实用经验，帮助系统管理员提升工作效率。
POJ3233Matrix Power Series 矩阵快速幂(分块矩阵构造)
cccrown的专栏
03-09 992
# include # include # include using namespace std; int mod; class Matrix{ public: int row; int col; int ** element; Matrix(); ~Matrix(); Matrix(const int , const
Matrix Power Series(矩阵加速，矩阵套矩阵)
zyz_bz的博客
04-15 851
题目 描述 给定矩阵AAA，求矩阵S=A1+A2+…+AkS=A^1+A^2+…+A^kS=A1+A2+…+Ak，输出矩阵，SSS矩阵中每个元都要模mmm。 数据范围：n(n≤30),k(k≤109)，m(m&lt;104)n (n ≤ 30) , k (k ≤ 10^9) ，m (m &lt; 10^4)n(n≤30),k(k≤109)，m(m<104) 输入 输入三个正...
Matrix Power Series 题解
crisp_chifan的博客
08-16 210
类似于一般的快速幂求法，在求矩阵快速幂里的ans数组要先赋值成单位矩阵，其他的和快速幂一样，因为不能直接返回数组，所以这个用一个结构体表示矩阵（一开始形参不小心写了二维数组就错了，救命）。一般来说用等比数列求和来计算，但是这里是矩阵没法求1-A和1-A^k，所以用分治的方式去求，其中需要求矩阵的乘方，用矩阵快速幂求解。然后用分治的方式求等比数列的和。n为奇数：f[n] = f[n >> 1] + a ^ n;n为偶数：f[n] = f[n >> 1]题意是对一个n*n的矩阵，求上述式子的结果。
Matrix Power Series
lewissun1的博客
08-15 198
思路：矩阵快速幂，由于k的大小，如果直接快速幂会爆，所以需要通过二分缩小计算规模。题意：给出n*n的矩阵A，求和A + A2 + A3 + …
TZOJ7929: Matrix Power Series（矩阵k次幂前缀和模板、二分优化）
Pyciii的博客
08-15 146
很明显这是一道矩阵k次幂前缀和模板题目，然而k达到1e9，如果直接使用矩阵快速幂求解肯定会超时，所以考虑进行二分优化。
POJ 3233 Matrix Power Series(矩阵快速幂+二分)
scf0920
09-17 2796
题目地址：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
【poj3233】Matrix Power Series——矩阵快速幂
CFhM_R's Blog
02-26 562
题目： 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
POJ 3233 Matrix Power Series（矩阵快速幂）
qq_36552817的博客
01-21 418
题目：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]...
poj3233Matrix Power Series矩阵快速幂
When you go away
07-28 952
/* * 思路：这题数据k<10^9，很显然需要O(log n)的算法，自然就是矩阵 * 快速幂．矩阵快速幂和快速幂差不多就是求矩阵A的k次方时，每次折半 * 类似于二分的思想，比如要求A的１００次方，我们先求出A的５０次方， * 然后用得到的结果平方就好了．试想如果要算16次方，每次折半，我们 * 只需要做4次矩阵乘法，但是如果采用朴素的连乘则需要15次，效率对比 * 可想而知．
LintCode 947: Matrix Power Series (矩阵快速幂经典题)
roufoo的博客
09-17 238
Matrix Power Series 中文English Given a n × n matrix A and a positive integer k and a positive integer m, find the sum S = A + A^2 + A^3 + … + A^k,return the elements of S modulo m Example Example 1: ...
Matrix Power Series(矩阵快速幂，二分法/构造矩阵套矩阵）
木静音yy的博客
08-15 340
传送门 题意： 思路: 方法一：（二分） 分析： 奇偶分类讨论： 1.代码： #include <iostream> #include <cstdio> #include <string> #include <cstring> #include <stdlib.h> #include <math.h> #include <ctype.h> #include <queue> #include <map&gt
POJ 3233 Matrix Power Series 动态规划（矩阵的幂）
qq_53038151的博客
12-05 228
给出一个矩阵A，输出矩阵B的每一项对M取余数的值。
Matrix Power Series（矩阵快速幂-矩阵乘法）
画船听雨眠
09-09 986
Matrix Power Series（矩阵快速幂-矩阵乘法） judge：POJ 3233 source：POJ Monthly–2007.06.03, Huang, Jinsong Time Limit: 3000MS Memory Limit: 131072K Description Given a n×nn × nn×n matrix AAA and a positive intege...
第三章，矩阵，05-逆矩阵
qq_17065591的博客
10-11 969
第三章，矩阵，05-逆矩阵定义逆矩阵的唯一性非奇异矩阵矩阵可逆的判定伴随矩阵矩阵可逆的判定定理矩阵可逆的推论逆矩阵性质 定义 对于n阶矩阵A，如果有一个n阶矩阵B，使得 AB=BA=EAB=BA=EAB=BA=E 则称矩阵A为可逆矩阵，而矩阵B为A的逆矩阵，记B=A−1B=A^{-1}B=A−1 逆矩阵的唯一性 证明： 若B和C都是A的逆矩阵，则： AB=BA=E,AC=CA=EAB=BA=E,AC=CA=EAB=BA=E,AC=CA=E 那么： B=BE=B(AC)=(BA)C=EC=CB=BE=B(AC
Matrix Power Series（矩阵求求次方和）
xiaoxiede_wo的博客
04-26 310
Matrix Power SeriesTime Limit: 3000MS Memory Limit: 131072KTotal Submissions: 26209 Accepted: 10782DescriptionGiven a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak.Inp...
C++信息竞赛NOIP各种必备实用模板（逆元，费马小定理，扩展欧几里得，快速幂，卡常，欧拉函数，组合数。。。。。。）
C20200521的博客
03-28 941
为什么会有此博客？ 此博客是为了博主复习时方便地背模板。 博主在代码部分都是默写的，所以可能会有点小错误，见谅。 前半部分重点关注于数论。如果觉得很好的朋友可以收藏了慢慢看 卡常必备！快速读入 int read() { int s=0,f=1;char a=getchar(); while(a<'0' || a>'9') { if(a=='-') f...
Matrix Power Series----矩阵乘法(二分)
wk476855的专栏
07-15 1488
题目：http://poj.org/problem?id=3233 A+A^2+A^3……+A^k 如果k是偶数的话，原式=(I+A^k/2)(A+A^2……A^k/2)。 如果k是奇数的话，原式=(1+A^(k/2+1))(A+A^2……+A^k/2)+A^(k/2+1)。 我是用数组直接实现的，比较麻烦，速度比较慢，1000+ms。 下面是我的代码： #include int
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值