矩阵儿快速幂 - POJ 3233 矩阵力量系列
不要管上面的标题的bug
那是幂的意思,不是力量。。。
POJ 3233 Matrix Power Series
描述
Given a n × n matrix A and a positive integer k, find the sum $ S = A + A^2 + A^3 + … + A^k $.
给你个n×n大小的矩阵A和一个正整数k,求矩阵S = A + A^2 + A^3 + … + A^k。
输入
The input contains exactly one test case. The first line of input contains three positive integers n (n ≤ 30), k (k ≤ 109) and m (m < 104).
输入只包含一个测试点。 第一行输入包含三个正整数n(n≤30),k(k≤10^9)和m(m <10^4)。
Then follow n lines each containing n nonnegative integers below 32,768, giving A’s elements in row-major order.
接下来n行,每行有n个低于32,768的非负整数,描述A的组成。
输出
Output the elements of S modulo m in the same way as A is given.
把矩阵S的值按输入A的格式模上m输出。
样例
.in
2 2 4
0 1
1 1
.out
1 2
2 3
暴力求解的话要交给天河2号评测
吾等平民还是写点高效算法吧
先把那个式子拆成下面这样子:
$ \sum^k_{i=1} {A^i} = (\sum^{k'}_{i=1} {A^i}) \times {(1 + A^{k'+1})} + 2|k?A^k:0 $
然后DFS分一分合一下就好了。
代码蒯上
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
inline int gotcha()
{register int a=0;bool b=1;register char c=getchar();while(c>'9' || c<'0'){if(c=='-')b=0;c=getchar();}while(c>='0' && c<='9')a=a*10+c-48,c=getchar();return b?a:-a;
}
int n,tim,mo;
struct martix
{int a[32][32];martix(){memset(a,0,sizeof(a));}const martix operator + (const martix &b)const{martix c;register int i,j;for(i=1;i<=n;i++)for(j=1;j<=n;j++)c.a[i][j]=a[i][j]+b.a[i][j],c.a[i][j]%=mo;return c;}const martix operator * (const martix &b)const{martix c;register int i,j,k;for(k=1;k<=n;k++)for(i=1;i<=n;i++)for(j=1;j<=n;j++)c.a[i][j]+=a[i][k]*b.a[k][j],c.a[i][j]%=mo;return c;}
}ori;
martix powa(martix in,int tim)
{if(tim==1)return in;martix tmp;for(int i=0;i<=n;i++)tmp.a[i][i]=1;if(tim==0)return tmp;while(tim){if(tim&1)tmp=in*tmp;in=in*in;tim>>=1;}return tmp;
}
martix finder(martix now,int k)
{if(k==1)return now;if(k&1)return finder(now,k-1)+powa(now,k);return (powa(now,0)+powa(now,k>>1))*finder(now,k>>1);
}
int main()
{register int i,j;n=gotcha(),tim=gotcha(),mo=gotcha();for(i=1;i<=n;i++)for(j=1;j<=n;j++)ori.a[j][i]=gotcha()%mo;ori=finder(ori,tim);for(i=1;i<=n;i++){for(j=1;j<=n;j++)printf("%d ",ori.a[j][i]);printf("\n");}return 0;
}
转载于:https://www.cnblogs.com/finder-iot/p/7647906.html
矩阵儿快速幂 - POJ 3233 矩阵力量系列相关推荐
- I-Matrix Power Series POJ - 3233 矩阵快速幂+分治
I-Matrix Power Series POJ - 3233 矩阵快速幂+分治 Problem Description Given a n × n matrix A and a positive ...
- HDU 4565So Easy!2012长沙邀请赛A题(共轭构造+矩阵的快速幂)
So Easy! Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total S ...
- [矩阵乘法/快速幂专题]Arc of Dream,Recursive sequence,233 Matrix,Training little cats
矩阵快速幂习题 复习矩阵乘法及快速幂模板 乘法模板 快速幂模板 T1:Arc of Dream 题目 题解 code T2:Recursive sequence 题目 题解 code T3:233 M ...
- 求解斐波那契第n项的几种解法(含矩阵乘法+快速幂) Python实现
斐波那契数列 首先我们来定义一下斐波那契数列: f(n)={0n = 01n = 1f(n−1)+f(n−2)n > 1f(n)= \begin{cases} 0 & \text {n ...
- POJ 3233 矩阵快速幂 等比数列
//http://www.cnblogs.com/jiangjing/archive/2013/05/28/3103336.html/* 分析:求a^1+..a^n这是矩阵乘法中关于等比矩阵的求法:| ...
- 矩阵快速幂 POJ 3070 Fibonacci
题目传送门 1 /* 2 矩阵快速幂:求第n项的Fibonacci数,转置矩阵都给出,套个模板就可以了.效率很高啊 3 */ 4 #include <cstdio> 5 #include ...
- poj 3233 矩阵乘法(分块矩阵)
POJ 3233 题解:Sn为所求矩阵, 则 这样, 此题就变成了求矩阵幂和矩阵乘法, 分块矩阵乘法和普通矩阵一样的. code: /* adrui's submission Language : C ...
- 【BZOJ 2323】 2323: [ZJOI2011]细胞 (DP+矩阵乘法+快速幂*)
2323: [ZJOI2011]细胞 Description 2222年,人类在银河系外的某颗星球上发现了生命,并且携带了一个细胞回到了地球.经过反复研究,人类已经完全掌握了这类细胞的发展规律: 这种 ...
- bzoj 4417: [Shoi2013]超级跳马(矩阵合并+快速幂)
4417: [Shoi2013]超级跳马 Time Limit: 10 Sec Memory Limit: 256 MB Submit: 410 Solved: 252 [Submit][Stat ...
最新文章
- windows redis安装与配置
- 人活着系列之芳姐和芳姐的猪(Floyd)
- c语言循环8,C语言8 循环语句 | CN-SEC 中文网
- float相乘后的类型_Psychopy | 第1期:数据类型及运算符
- jQuery基础(未完待续)
- 通过Spring Data Neo4J操作您的图形数据库
- asp留言板源码-XYCMS留言板 v8.0
- JavaScript输入输出
- 为啥计算机课要带u盘,电脑课上,student 关掉以后会不会被老师发现,同时电脑机上,老师不知干了什么,U盘没法使用,求解...
- 在Oracle中写出性能优良的SQL语句
- 小米便签类间关系图绘制方法
- HTML5网页多媒体( 音频,视频,旧版本浏览器提示,格式兼容,多媒体标签属性)
- line划线计算机图像学,《计算机图形学基础》OpenGL中点画线法画直线(同时支持k1和k1的情况)...
- win7计算机管理无用户账户,win7系统控制面板中用户账户显示空白的解决方法
- 1、Centos之常见目录作用介绍
- Abp.vnext 使用VS2022+Docker 部署至阿里云
- 埋骨何须桑梓地,人生无处不青山!
- 分析波形信号解决开发板串口连PC后无输出或输出乱码的问题
- P4278 带插入区间K小值
- 2021版网络安全等级测评师培训教材(初级)
热门文章
- REST Framework 的用户认证组件
- 打通两台机器的ssh功能
- jdk1.8.0_45源码解读——Set接口和AbstractSet抽象类的实现
- 五大开源游戏引擎介绍
- mysql 判断是否已存在及插入表数据 的 简单存储过程
- linux的lsm--TCB组成部分的新视角
- Microsoft uaa bus driver for high definition audio
- 'SVN更新' has encountered a problem :An internal error occurred during: svn错误
- OpenCV删除面积小的区域 实现图像二值化分割 标记连通区域
- 视频会议H.323协议SIP协议