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

最新推荐文章于 2019-05-06 22:27:30 发布
原创 最新推荐文章于 2019-05-06 22:27:30 发布 · 902 阅读 · 0 ·
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

该博客主要解决POJ3233题目，通过构建分块矩阵来实现矩阵的快速幂运算，求解矩阵的前k项和。内容涉及矩阵乘法和快速幂算法在n*n矩阵上的应用，适用于解决大整数下的高效计算问题。

/****************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&nbs

最低0.47元/天 解锁文章
yi淼菌君
关注
POJ3233---Matrix Power Series （矩阵快速幂（升级）：矩阵套矩阵）
Crazy
07-09 803
【题目来源】：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 366
题目：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 1049
大致题意：简单题意就不解释了。    不过可以建议大家可以先做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 320
传送门 题意： 题解： 附上代码： #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 1242
题意 ：给出n k m 和一个n*n的矩阵A 求A + A2 + A3 + … + Ak 参考http://blog.csdn.net/wangjian8006/article/details/7868864 构造矩阵很重要啊！！！ 弱菜不会啊 #include #include const int mod = 10000; const int maxn = 66; struct Ma
poj 3233 Matrix Power Series （构造分块矩阵）
Jiandong
11-25 319
题目链接：哆啦A梦传送门 题意：自己看。 参考博客：神犇 题解：分块矩阵：分块矩阵可以构造求和。 例如：我们可以这样构造， 还需注意一点的是：算完S(k+1)，取出右上角矩阵分块后，还需减掉单位矩阵E。   代码不是我写的，我就按自己习惯改了下变量，主要是博主写的太溜了，不搬都对不起他： #include&lt;cstdio&gt; #include&lt;algorithm...
POJ3233 Matrix Power Series构造矩阵解决问题
松子茶的专栏
05-07 5289
时间限制：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
POJ 3233——Matrix Power Series（矩阵快速幂）
Sporadic memory 的博客
07-29 309
第一次在博客里发题解（被水淹没不知所措）。。。 题目：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。 题目分析：这道题首先要知道矩阵乘法（不知道的请移步百度），然后一看数据范围就知道必须要用快速幂。但是这道题又特坑，...
POJ 3233 Matrix Power Series(矩阵快速幂+二分)
scf0920
09-17 2744
题目地址：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
hdu 6395 Sequence 分块矩阵快速幂
大乐的博客
08-28 218
 题目链接：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 1091
题干： 求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 126
题目链接 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 290
给出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 211
其实就是考擦构造矩阵。 我瞎构造也能构造出来 我构造的矩阵是 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 2362
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 256
#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 1279
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 908
题目描述 有一个长度为 nnn 的序列，xht37 现在想分块维护它。 PinkRabbit 要求他只准将序列分成 PRPRPR 种长度的块。 NaCly_Fish 要求他只准将序列分成 NFNFNF 种长度的块。 同一个人可能会要求 xht37 多次相同的块长。 xht37 想同时满足 PinkRabbit 和 NaCly_Fish 要求，只好使用两个人都允许的长度分块。 xht37...
hdu6395 /// 分块矩阵快速幂
babimei370256276的博客
02-11 413
题目大意： 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) { ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值