poj3233 Matrix Power Series

最新推荐文章于 2023-12-30 17:20:46 发布
最新推荐文章于 2023-12-30 17:20:46 发布
阅读量218 收藏 1
点赞数 2
分类专栏： 矩阵快速幂 -------------------数论------------------
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_33362864/article/details/76945479
版权

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 ≤

最低0.47元/天 解锁文章
UMR小豪
关注
  • 2
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
POJ 3233 Matrix Power Series（矩阵快速幂）
qq_36552817的博客
01-21 288
题目：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
lifehappy
07-24 260
Matrix Power Series 思路 题意比较简单，就是要求S(n)=∑i=1nAiS(n) = \sum _{i = 1} ^{n} A^ {i}S(n)=∑i=1n​Ai，显然有S(n)=S(n−1)∗A+AS(n) = S(n - 1) * A + AS(n)=S(n−1)∗A+A，看到这里，那就简单了，递推式，加矩阵，矩阵快速幂无疑了嘛，所以我们开始构造矩阵。 显然有如下矩阵，EEE是单位矩阵，AAA是输入矩阵，OOO是零矩阵。 [EEOA]∗[OOAO]\begin{bmatrix} E
POJ 3233 Matrix Power Series
wlh1998的博客
07-20 145
//数据小于32768 #include #include #include using namespace std; const int maxn = 40; typedef long long ll; struct Mat { int mat[maxn][maxn]; int row,col; }; Mat mod_add(Mat a, Mat b, int p){ Mat ans
POJ3233 Matrix Power Series
dezhen7015的博客
05-12 70
Given an×nmatrixAand a positive integerk, find the sumS=A+A2+A3+ … +Ak. Input The input contains exactly one test case. The first line of input contains three positive int...
算法实验T8——POJ3233 Matrix Power Series
最新发布
qq_43010742的博客
12-30 351
但直接就和的话还是会有很多重复计算，可以分奇偶讨论k，对S二分：即k为偶数时，，由此可实现二分求解，重复利用子问题的解。已知一个 n * n 的矩阵 A，计算。对于单个矩阵的幂想到采用。，S 的值对 m 取模。
POJ3233 Matrix Power Series（矩阵快速幂）（矩阵快速幂+等比数列二分求和取模）
qq_43460290的博客
04-25 314
首先给出题目： Matrix Power Series 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 ...
POJ1459-Power Network
07-30
【标题】"POJ1459-Power Network"是一个编程竞赛题目，源自北京大学的在线判题系统POJ（Problem Set of Peking University）。这个题目主要涉及到图论和算法的应用，特别是广度优先搜索（BFS）以及一种优化技巧——...
POJ2109-Power of Cryptography
07-29
《POJ2109-Power of Cryptography：解题报告与AC代码解析》 在计算机科学领域，算法和编程竞赛是检验和提升编程技能的重要途径。POJ（Problem Set of Peking University）是由北京大学主办的一个在线编程竞赛平台，...
POJ1001-Precision power
07-29
【标题】"POJ1001-Precision power"是一个在线编程题目，源自北京大学的POJ（Problem Set of Peking University）平台。该题目主要考察的是算法设计与精度控制方面的知识，尤其是涉及到浮点数计算时的精度问题。 ...
poj 1459 Power Network.md
11-17
poj 1459 Power Network.md
poj.rar_poj
09-23
标题中的"poj.rar_poj"暗示了这是一个与POJ（Programming Online Judge）相关的压缩文件，POJ是一个在线编程竞赛平台，主要供程序员们练习和提交算法解决方案。在这个压缩包中，很可能包含了用户在POJ上参与编程挑战...
费马检测（判断一个数是不是素数）
樱花满地集于我心，楪舞纷飞祈愿相随
02-22 3565
#include using namespace std; typedef long long ll using namespace std; ll Quick_Mod(ll a, ll b, ll mod)//快速幂 { ll res = 1,term = a % mod; while(b) { if(b
HDU2643 Rank 第二类斯特林数
樱花满地集于我心，楪舞纷飞祈愿相随
08-09 1277
第二类斯特林数把一个包含n个元素的集合分成k个非空子集的方法数： 初始值： S2(n,0)=0,S2(n,k)=0(n<k)S_{2(n,0)}=0,S_{2(n,k)}=0(n<k) S2(n,1)=S2(n,n)=1S_{2(n,1)}=S_{2(n,n)}=1 递推式： S2(n,k)=S2(n−1,k−1)+S2(n−1,k)∗kS_{2(n,k)}=S_{2(n-1,k-1)}+
Cf 406 A. The Monster
樱花满地集于我心，楪舞纷飞祈愿相随
03-24 745
A. The Monster time limit per test 1 second memory limit per test 256 megabytes input standard input output standard output A monster is chasing after Rick and Morty on
poj2115 C Looooops(扩展欧几里得)
樱花满地集于我心，楪舞纷飞祈愿相随
08-23 427
C Looooops Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 23143   Accepted: 6357 Description A Compiler Mystery: We are given a C-language style for loop(
poj2635 The Embarrassed Cryptographer（高精度 同余取模）
樱花满地集于我心，楪舞纷飞祈愿相随
08-20 412
The Embarrassed Cryptographer Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 13831   Accepted: 3762 Description The young and very promising cryptographer
poj3292 Semi-prime H-numbers
樱花满地集于我心，楪舞纷飞祈愿相随
08-22 406
Semi-prime H-numbers Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 8845   Accepted: 3896 Description This problem is based on an exercise of David Hilber
Count primes
樱花满地集于我心，楪舞纷飞祈愿相随
09-19 403
Count primes Time Limit: 12000/6000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 521    Accepted Submission(s): 221 Problem Description Easy question!
用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、付费专栏及课程。

余额充值