POJ 3233 Matrix Power Series 动态规划(矩阵的幂)

一、题目大意

给出一个矩阵A,

输出矩阵B的每一项对M取余数的值。

二、解题思路

以二维矩阵为例,首先计算K=2的情况,我们设结果矩阵为B

有如下表达式

那么不难看出,需要的矩阵其实就是以下的两个矩阵相乘后的左上角的N*N个

然后我们再来考虑K=3的情况,我们设结果矩阵为C

我们来考虑如何把C表示成矩阵B和A相乘的状态。

不难看出C矩阵就是以下两个矩阵相乘后的N*N的左上角

也是如下矩阵的左上角

综上,我们发现计算B和C时,乘号右边的矩阵是相同的,只是我们需要保证前N排的后N列必须始终为A[0][0]..A[0][N-1],A[1][0]..A[1][N-1],...

保留那些值我们可以通过给右边的矩阵右下角放置零矩阵即可。

那么最终题目的答案可以表示成下的表达式。

三、代码

#include <iostream>
#include <vector>
using namespace std;
typedef vector<int> vec;
typedef vector<vec> mat;
int M;
mat mul(mat &A, mat &B)
{mat C = mat(A.size(), vec(B[0].size()));for (int i = 0; i < A.size(); i++){for (int j = 0; j < B[0].size(); j++){for (int k = 0; k < B.size(); k++){C[i][j] = (C[i][j] + A[i][k] * B[k][j]) % M;}}}return C;
}
mat pow(mat &A, int N)
{mat B = mat(A.size(), vec(A[0].size()));for (int i = 0; i < B.size(); i++){B[i][i] = 1;}while (N > 0){if (N & 1){B = mul(B, A);}A = mul(A, A);N >>= 1;}return B;
}
void solve()
{int N, K;scanf("%d%d%d", &N, &K, &M);mat A = mat(2 * N, vec(2 * N));for (int i = 0; i < N; i++){for (int j = 0; j < N; j++){scanf("%d", &A[i][j]);}}for (int i = 0; i < N; i++){A[i + N][i] = 1;A[i + N][i + N] = 1;}mat B = mat(2 * N, vec(2 * N));for (int i = 0; i < N; i++){for (int j = 0; j < N; j++){B[i][j] = A[i][j];B[i][j + N] = A[i][j];}}A = pow(A, K - 1);A = mul(B, A);for (int i = 0; i < N; i++){for (int j = 0; j < N; j++){printf("%d%c", A[i][j] % M, j + 1 == N ? '\n' : ' ');}}
}
int main()
{solve();return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/200201.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

初识Linux——基本指令(详解)1

呀哈喽&#xff0c;我是结衣。 在学习数据结构的同时&#xff0c;也不要忘了Linux的学习啊。今天我们开始Linux的教学&#xff0c;在学习之前我们肯定要会搭建Linux的学习环境&#xff0c;在我们的以前的博客里是有讲解的&#xff0c;所以所以这里我们就不在多说&#xff0c;我…

UDP数据报套接字

文章目录 DatagramSocket APIDatagramPacket API示例一: 请求响应UDP服务端UDP客户端 DatagramSocket API Socket是操作系统中的一个概念&#xff0c;本质上是一种特殊的文件&#xff0c;Socket就属于把“网卡”这个设备给抽象成了文件。往 Socket 文件中写数据&#xff0c;就…

Usergolang 一些优质关于sip协议包

在Go语言中&#xff0c;有一些优质的SIP协议包&#xff0c;适用于构建SIP客户端和服务器。以下是其中一些常用的SIP库&#xff1a; 1. **github.com/cloudwebrtc/sip:** - GitHub 地址&#xff1a;[cloudwebrtc/sip](https://github.com/cloudwebrtc/sip) - 该库提供了用…

深入探讨Guava的缓存机制

第1章&#xff1a;引言 大家好&#xff0c;我是小黑&#xff0c;今天咱们聊聊Google Guava的缓存机制。缓存在现代编程中的作用非常大&#xff0c;它能提高应用性能&#xff0c;减少数据库压力&#xff0c;简直就是性能优化的利器。而Guava提供的缓存功能&#xff0c;不仅强大…

iptables入门

今天我的工作遇到了巡检网络配置的任务&#xff0c;这次巡检的主机都是运行十多年的机器&#xff0c;并不是新的firewalld&#xff0c;基本都是iptables&#xff0c;上学的时候以为这些都没人用&#xff0c;所以没有认真学习&#xff0c;现在需要用到了&#xff0c;所以写一篇文…

PostGIS学习教程九:空间连接

PostGIS学习教程九&#xff1a;空间连接 空间连接&#xff08;spatial joins&#xff09;是空间数据库的主要组成部分&#xff0c;它们允许你使用空间关系作为连接键&#xff08;join key&#xff09;来连接来自不同数据表的信息。我们认为“标准GIS分析”的大部分内容可以表示…

17.认识下Docker之docker的核心原理(2)

1.容器-我的小世界 不知道大家看没看过小说《完美时间》&#xff0c;里面石昊经常进入一个小世界在里面与世隔绝的修炼或者战斗&#xff0c;总之就是在一个完全封闭的空间里做他想做的事情而与外界隔离&#xff0c;不受侵扰。通过前面的分析我们知道&#xff0c;Namepace让应用…

SringBoot的启动原理,保姆级带你认识,让面试管对你刮目相看

SringBoot的启动原理&#xff0c;保姆级带你认识&#xff0c;让面试管对你刮目相看 一&#xff0c;介绍 graph TD; A[启动类] --> B[SpringApplication.run()] B --> C[创建SpringApplication实例] C --> D[初始化应用上下文] D --> E[加载应用配置] E --> F[…

树_二叉搜索树累加求和

//给出二叉 搜索 树的根节点&#xff0c;该树的节点值各不相同&#xff0c;请你将其转换为累加树&#xff08;Greater Sum Tree&#xff09;&#xff0c;使每个节点 node 的新值等于原树中大于或等于 // node.val 的值之和。 // // 提醒一下&#xff0c;二叉搜索树满足下列约束…

【Vulnhub 靶场】【hacksudo: FOG】【简单 - 中等】【20210514】

1、环境介绍 靶场介绍&#xff1a;https://www.vulnhub.com/entry/hacksudo-fog,697/ 靶场下载&#xff1a;https://download.vulnhub.com/hacksudo/hacksudo-FOG.zip 靶场难度&#xff1a;简单 - 中等 发布日期&#xff1a;2021年05月14日 文件大小&#xff1a;1.3 GB 靶场作…

解决 vite 中 import.meta.globEager is not function 的问题

本人正在重构两年前搭建到一半的博客网站&#xff0c;相关依赖都很陈旧&#xff0c;用到了 npm-check-updates 检测项目可升级依赖&#xff1a; 升级完成后解决完依赖发现控制台报错 import.meta.globEager is not function解决方案&#xff1a; vite版本降至4.3.0 亲测有效&…

忽略python运行出现的大量警告

添加以下代码即可 import warnings warnings.filterwarnings(ignore)

1107. 魔板(BFS,最小步数模型,unordered_map哈希)

1107. 魔板 - AcWing题库 Rubik 先生在发明了风靡全球的魔方之后&#xff0c;又发明了它的二维版本——魔板。 这是一张有 8 个大小相同的格子的魔板&#xff1a; 1 2 3 4 8 7 6 5我们知道魔板的每一个方格都有一种颜色。 这 8 种颜色用前 8 个正整数来表示。 可以用颜色的…

C语言面试之旅:掌握基础,探索深度(面试实战之ARM架构二)

读别人的代码&#xff0c;就像撕洋葱&#xff0c;每读一层&#xff0c;哭一遍。 引言 ARM 处理器是冯诺依曼存储结构&#xff0c;程序空间、RAM 空间及IO 映射空间统一编址&#xff0c;除对对RAM 操作以外&#xff0c;对外围IO、程序数据的访问均要通过加载/存储指令进…

介绍 Agency: 使AI与Go语言无缝对接

关注公众号【爱发白日梦的后端】分享技术干货、读书笔记、开源项目、实战经验、高效开发工具等&#xff0c;您的关注将是我的更新动力&#xff01; 在当今应用开发领域&#xff0c;类似OpenAI API等生成式AI技术的蓬勃发展正在彻底改变着应用开发的格局。Python和JavaScript等语…

App内存优化

一、内存优化介绍 1.背景介绍 内存是大问题但缺乏关注压实骆驼的最后一个稻草&#xff08;堆栈溢出&#xff09; 2.内存问题 内存抖动&#xff1a;锯齿状、GC导致卡顿内存泄露&#xff1a;可用内存减少、频繁GC内存溢出&#xff1a;OOM&#xff0c;程序异常 二、优化工具选…

小程序的分包加载具体流程

在小程序项目的 app.json 文件中配置分包信息。通过设置 subpackages 字段&#xff0c;指定每个分包的根目录路径和分包的名称。 "subpackages": [{"root": "subpackage1","name": "subpackage1"},{"root": &quo…

FreeRTOS简介

FreeRTOS简介 文章目录 FreeRTOS简介前言一、什么是FreeRTOS&#xff1f;二、FreeRTOS的特点三、FreeRTOS的版本和参考资料1、FreeRTOS版本2、FreeRTOS源码获取3、FreeRTOS参考资料 四、FreeRTOS源码简介 前言 FreeRTOS是一个可以基于ROM运行的、可裁剪的、抢占式、实时多任务…

Photoshop最新版PS2024安装使用 Ver25.0.0

Photoshop&#xff0c;这个是长红了几十年的软件&#xff0c;我大概从它的3.0版本开始用&#xff0c;目前已迭代到25.0&#xff0c;但一直还在用CS4/11.0版本&#xff0c;一直秉持着够用即可的原则&#xff0c;因为不是专业的平面设计人员&#xff0c;能够简单PP图片就行。&…

C++智能指针及简单实现

C智能指针 堆内存、栈内存与静态内存静态内存栈内存堆内存 动态内存管理new、delete运算符智能指针实现智能指针 shared_ptr智能指针的线程安全问题解决 unique_ptrweak_ptr循环引用 思维导图本模块思路 动态内存管理 - cppreference.com 堆内存、栈内存与静态内存 静态内存 …