POJ3233Matrix Power Series 矩阵快速幂(分块矩阵构造)

最新推荐文章于 2019-05-06 22:27:30 发布
最新推荐文章于 2019-05-06 22:27:30 发布
阅读量868 收藏
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.youkuaiyun.com/cccrown/article/details/20849013
版权
该博客主要解决POJ3233题目，通过构建分块矩阵来实现矩阵的快速幂运算，求解矩阵的前k项和。内容涉及矩阵乘法和快速幂算法在n*n矩阵上的应用，适用于解决大整数下的高效计算问题。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

/****************88

这题一看就找到感觉，构造分块矩阵 

  1   A    *   Sk-1     =  Sk                        其中A和Sk都是n*n的矩阵, 1是n*n的单位矩阵，0是n*n的零矩阵

  0   A        Ak-1         Ak

ad****************/

Matrix Power Series
Time Limit: 3000MS   Memory Limit: 131072K
Total Submissions: 13592   Accepted: 5875

Description

Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak.

Input

最低0.47元/天 解锁文章
yi淼菌君
关注
POJ3233---Matrix Power Series （矩阵快速幂（升级）：矩阵套矩阵）
Crazy
07-09 780
【题目来源】：https://vjudge.net/problem/POJ-3233 【题意】 题意让求Sn，其中呢，Sn= A + A2 + A3 + … + An. 【思路】 一开始就想到利用前n项和，找到对应的关系，最简单的关系是这样的: Sn=Sn-1+A(n-1)*A 写出来的矩阵关系式就是这样的（本人写法比较喜欢吧系数矩阵放在右边） 1 A 0 A 乘 Sn-1 0
POJ 3233 Matrix Power Series（矩阵快速幂）
qq_36552817的博客
01-21 341
题目：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]...
POJ 3233 Matrix Power Series (矩阵快速幂)
ACM之路，在希望的田野上。
04-04 1022
大致题意：简单题意就不解释了。    不过可以建议大家可以先做POJ 3070，先学会快速幂的基本思想。 没有做过的可以查看我的博客：点击打开链接 然后我们已经会使用矩阵快速幂求解A^k，则如何求解  A+A^2+A^3…… 这里还是运用了二分的思想 ，例如： 令S(N)=A+A^2+……+A^N;         则S(6)=A+A^2+A^3+A^4+A^5+A^6=(1+A^3)
Matrix Power Series POJ - 3233 (矩阵快速幂)
肘子的博客
10-11 297
传送门 题意： 题解： 附上代码： #include&lt;iostream&gt; #include&lt;cstdio&gt; using namespace std; typedef long long ll; const int MAXN=70; struct node{ ll a[MAXN][MAXN]; }; node shu,ans,mp; ll N; ...
POJ 3233 Matrix Power Series 矩阵快速幂求A+A2+A3+…+Ak
嗯。
04-03 1217
题意 ：给出n k m 和一个n*n的矩阵A 求A + A2 + A3 + … + Ak 参考http://blog.youkuaiyun.com/wangjian8006/article/details/7868864 构造矩阵很重要啊！！！ 弱菜不会啊 #include #include const int mod = 10000; const int maxn = 66; struct Ma
poj 3233 Matrix Power Series （构造分块矩阵）
Jiandong
11-25 299
题目链接：哆啦A梦传送门 题意：自己看。 参考博客：神犇 题解：分块矩阵：分块矩阵可以构造求和。 例如：我们可以这样构造， 还需注意一点的是：算完S(k+1)，取出右上角矩阵分块后，还需减掉单位矩阵E。   代码不是我写的，我就按自己习惯改了下变量，主要是博主写的太溜了，不搬都对不起他： #include&lt;cstdio&gt; #include&lt;algorithm...
POJ3233 Matrix Power Series构造矩阵解决问题
松子茶的专栏
05-07 5176
时间限制：1000 ms  |  内存限制：65535 KB  难度：4描述 Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak.输入The input contains exactly one test case. The first line of input contai
【poj3233】Matrix Power Series——矩阵快速幂
CFhM_R's Blog
02-26 515
题目： 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（矩阵快速幂）
Sporadic memory 的博客
07-29 288
第一次在博客里发题解（被水淹没不知所措）。。。 题目：POJ 3233 Matrix Power Series 题目大意：这道题不用翻译也能看懂，有一个n*n的矩阵，给出一个k,求S=（A+A^2+A^3+…+A^k）%m。其中n ≤ 30，k ≤ 10^9，m&lt;104。 题目分析：这道题首先要知道矩阵乘法（不知道的请移步百度），然后一看数据范围就知道必须要用快速幂。但是这道题又特坑，...
hdu 6395 Sequence 分块矩阵快速幂
大乐的博客
08-28 195
 题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6395 我真的服了，一直不知道T在哪里 最后改着改着变成wa了，发现是mod数写错了，改了就过了 #include&lt;bits/stdc++.h&gt; #define inf 0x3f3f3f3f #define mod 1000000007 #define For(i,m,n) fo...
POJ - 3233 —— 矩阵快速幂（分块矩阵与递归）
lucius-liu.cn
05-19 996
题干： 求A + A^2 + A^3 +...+ A^k 思路： 第一种做法：（递归） 从小栗子开始思考 —— S = A + A^2 + A^3 + ... + A^6 则 S = （A + A^2 + A^3）*（E + A^3）= 【 A * （E + A）+ A^3】* （E + A^3） 发现可以将长度不断减半，可以进行递归，然后写出递归公式 S(k)...
Sequence( 分块+矩阵快速幂 )
weixin_30693183的博客
11-01 111
题目链接 1 #include<bits/stdc++.h> 2 using namespace std; 3 #define e exp(1) 4 #define pi acos(-1) 5 #define mod 1000000007 6 #define inf 0x3f3f3f3f 7 #define ll long long 8 ...
HDU 6395 Sequence (分块矩阵快速幂)
NPU_SXY的博客
09-03 265
给出F的定义求F(n) % 1e9+7。 如果没有这个p/n，就是一个简单的矩阵快速幂。 向下取整的值在1-n中，显然每一个数出现的时候都是连续的一段，这样就可以分块地去求矩阵快速幂。 对于i，可以得出，p/i最后一次出现就是p/(p/i)。 以及为什么G++比C++提交快这么多啊嘤嘤嘤   #pragma once #include &lt;cstdio&gt; #include ...
POJ 3233 矩阵快速幂
java开发指南博客 【转载】
04-16 192
其实就是考擦构造矩阵。 我瞎构造也能构造出来 我构造的矩阵是 Let B = A A 0 I 然后看discuss里这么构造了一个矩阵 Let B= A I 0 I B^(k+1) = A^k I+A+...+A^k 0 I 最后输出的时候还要剪掉一个单位矩阵 然后我觉得我构造的挺好，比...
分块+矩阵快速幂
周周的博客
08-14 2330
http://acm.hdu.edu.cn/showproblem.php?pid=6395 P/i会有相同的一段，所以想到分块，这个也经常用 #include&lt;algorithm&gt; #include &lt;iostream&gt; #include &lt;stdlib.h&gt; #include &lt;string.h&gt; #include &lt;std...
分块加矩阵快速幂
weixin_34289454的博客
08-14 233
#include<stdio.h> #include<string.h> #include<algorithm> using namespace std; typedef long long ll; const int mod=1e9+7; int a,b,c,d,p,n; struct node { int m[3][3]; ...
poj 3233 矩阵乘法（分块矩阵）
退役ACMer
10-03 1251
POJ 3233 题解：Sn为所求矩阵， 则 这样， 此题就变成了求矩阵幂和矩阵乘法， 分块矩阵乘法和普通矩阵一样的。 code： /* adrui's submission Language : C++ Result : Accepted Love : ll Favorite : Dragon Balls Standing in the Ha
P5343 【XR-1】分块 矩阵快速幂
Wind_bow的博客
05-06 860
题目描述 有一个长度为 nnn 的序列，xht37 现在想分块维护它。 PinkRabbit 要求他只准将序列分成 PRPRPR 种长度的块。 NaCly_Fish 要求他只准将序列分成 NFNFNF 种长度的块。 同一个人可能会要求 xht37 多次相同的块长。 xht37 想同时满足 PinkRabbit 和 NaCly_Fish 要求，只好使用两个人都允许的长度分块。 xht37...
hdu6395 /// 分块矩阵快速幂
babimei370256276的博客
02-11 376
题目大意： F(1)=A， F(2)=B, F(i)=C*F(i-2)+D*F(i-1)+p/i(向下取整) 给定A B C D p n 求F(n) 构造 　　　矩阵A　　　　*　　 矩阵B = 矩阵C ┌ F(n-2) F(n-1) 1 ┐　　┌ 0 C 0 ┐ ┌ F(n-1) F(n) 1 ┐...
用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、付费专栏及课程。

余额充值