poj 3233 Matrix Power Series

Matrix Power Series

思路

题意比较简单,就是要求S(n)=∑i=1nAiS(n) = \sum _{i = 1} ^{n} A^ {i}S(n)=i=1nAi,显然有S(n)=S(n−1)∗A+AS(n) = S(n - 1) * A + AS(n)=S(n1)A+A,看到这里,那就简单了,递推式,加矩阵,矩阵快速幂无疑了嘛,所以我们开始构造矩阵。

显然有如下矩阵,EEE是单位矩阵,AAA是输入矩阵,OOO是零矩阵。

[EEOA]∗[OOAO]\begin{bmatrix} E & E \\ O & A\end{bmatrix} * \begin{bmatrix} O & O\\ A & O \end{bmatrix}[EOEA][OAOO]

通过这个矩阵的递推,我们就可以通过快速幂,达到快速求解的目的。

我严重怀疑这道题目数据有问题,longlonglong\ longlong longwawawa,然后intintint就过了???

AC代码

/*Author : lifehappy
*/
// #pragma GCC optimize(2)
// #pragma GCC optimize(3)
// #include <bits/stdc++.h>#include <cstdio>
#include <iostream>
#include <stdlib.h>
#include <algorithm>
#include <cmath>#define mp make_pair
#define pb push_back
#define endl '\n'using namespace std;typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;const double pi = acos(-1.0);
const double eps = 1e-7;
const int inf = 0x3f3f3f3f;inline ll read() {ll f = 1, x = 0;char c = getchar();while(c < '0' || c > '9') {if(c == '-')    f = -1;c = getchar();}while(c >= '0' && c <= '9') {x = (x << 1) + (x << 3) + (c ^ 48);c = getchar();}return f * x;
}void print(ll x) {if(x < 10) {putchar(x + 48);return ;}print(x / 10);putchar(x % 10 + 48);
}const int N = 70;int n, k, mod;struct matrix {int a[N][N];matrix operator * (const matrix & t) const {matrix temp;for(int i = 1; i <= 2 * n; i++) {for(int j = 1; j <= 2 * n; j++) {temp.a[i][j] = 0;for(int k = 1; k <= 2 * n; k++) {temp.a[i][j] = (temp.a[i][j] + a[i][k] * t.a[k][j]) % mod;}}}return temp;}
}E, A, O;int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);n = read(), k = read(), mod = read();matrix fat, ans;for(int i = 1; i <= 2 * n; i++) {//先置零,for(int j = 1; j <= 2 * n; j++) {fat.a[i][j] = ans.a[i][j] = 0;}}for(int i = 1; i <= n; i++) {//读入的时候置A加上置E矩阵。for(int j = 1; j <= n; j++) {fat.a[i + n][j + n] = ans.a[i + n][j] = read();}fat.a[i][i] = fat.a[i][i + n] = 1;}while(k) {if(k & 1) ans = fat * ans;fat = fat * fat;k >>= 1;}for(int i = 1; i <= n; i++) {for(int j = 1; j <= n; j++) {printf("%d%c", ans.a[i][j], j == n ? '\n' : ' ');}}return 0;
}

调不出来的代码

写了一手逼格高一点的举证套矩阵的重载操作符的写法,可是太菜了,调不出来

/*Author : lifehappy
*/
// #pragma GCC optimize(2)
// #pragma GCC optimize(3)
// #include <bits/stdc++.h>#include <cstdio>
#include <iostream>
#include <stdlib.h>
#include <algorithm>
#include <cmath>#define mp make_pair
#define pb push_back
#define endl '\n'using namespace std;typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;const double pi = acos(-1.0);
const double eps = 1e-7;
const int inf = 0x3f3f3f3f;inline ll read() {ll f = 1, x = 0;char c = getchar();while(c < '0' || c > '9') {if(c == '-')    f = -1;c = getchar();}while(c >= '0' && c <= '9') {x = (x << 1) + (x << 3) + (c ^ 48);c = getchar();}return f * x;
}void print(ll x) {if(x < 10) {putchar(x + 48);return ;}print(x / 10);putchar(x % 10 + 48);
}const int N = 70;int n, k, mod;struct matrix {int a[N][N];matrix operator * (const matrix & t) const {matrix temp;for(int i = 1; i <= 2 * n; i++) {for(int j = 1; j <= 2 * n; j++) {temp.a[i][j] = 0;for(int k = 1; k <= 2 * n; k++) {temp.a[i][j] = (temp.a[i][j] + a[i][k] * t.a[k][j]) % mod;}}}return temp;}matrix operator + (const matrix & t) const {matrix temp;for(int i = 1; i <= n; i++) {for(int j = 1; j <= n; j++) {temp.a[i][j] = (a[i][j] + t.a[i][j]) % mod;}}return temp;}
}E, A, O;struct Matrix {matrix a[3][3];Matrix operator * (const Matrix & t) const {Matrix temp;for(int i = 1; i <= 2; i++) {for(int j = 1; j <= 2; j++) {temp.a[i][j] = O;for(int k = 1; k <= 2; k++) {temp.a[i][j] = (a[i][k] * t.a[k][j]) + temp.a[i][j];}}}}
};int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);n = read(), k = read(), mod = read();for(int i = 1; i <= n; i++) {for(int j = 1; j <= n; j++) {A.a[i][j] = read();E.a[i][j] = O.a[i][j] = 0;}E.a[i][i] = 1;}Matrix fat, ans;fat.a[1][1] = E, fat.a[1][2] = E, fat.a[2][1] = O, fat.a[2][2] = A;ans.a[1][1] = O, ans.a[1][2] = O, ans.a[2][1] = A, ans.a[2][2] = O;while(k) {if(k & 1) ans = ans * fat;fat = fat * fat;k >>= 1;}for(int i = 1; i <= n; i++) {for(int j = 1; j <= n; j++) {printf("%d%c", ans.a[1][1].a[i][j], j == n ? '\n' : ' ');}}return 0;
}

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

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

相关文章

.NET Core 编写 Azure Function 并连接 GitHub 持续部署

点击上方蓝字关注“汪宇杰博客”导语Azure Function 是一个事件驱动型无服务器计算平台&#xff0c;可以解决复杂的业务流程问题&#xff0c;更加高效地进行开发。在本地构建和调试&#xff0c;而无需额外的设置&#xff0c;在云中大规模部署和操作&#xff0c;并使用触发器和绑…

Codeforces Round #766 (Div. 2) D. Not Adding 数学gcd

传送门 文章目录目录题意&#xff1a;思路&#xff1a;目录 题意&#xff1a; 给你一个长度为nnn的数组&#xff0c;你每次可以选择其中的两个数&#xff0c;如果他们的gcdgcdgcd在数组中没有出现那么就可以加在数组后面构成一个新的数组&#xff0c;问数组最长是多少。 2≤n…

「数据ETL」从数据民工到数据白领蜕变之旅(五)-使用dotNET脚本实现SSIS无限扩展...

在前面一文中&#xff0c;正式引出了SSIS专业数据ETL工具&#xff0c;笔者仅能作引路作用&#xff0c;未能使用文章的方式给大家写出更多的入门级的文章&#xff0c;希望读者们可以自行根据分享的学习资源自行完成入门及进阶的学习。同时也想给大家分享到SSIS的能力边界性&…

P3067 [USACO12OPEN]Balanced Cow Subsets G 折半搜索

传送门 文章目录目录题意&#xff1a;思路&#xff1a;目录 题意&#xff1a; 给你nnn个数&#xff0c;从中任意选出一组数&#xff0c;使这些数能分成和相等的两组&#xff0c;问有多少种选数方案。 2≤n≤20,1≤ai≤1e92\le n\le 20,1\le a_i\le 1e92≤n≤20,1≤ai​≤1e9 …

poj 2191 Mersenne Composite Numbers

思路 用pollard_rho miller_rabin来拆分数字&#xff0c;如果得到的质因子大于等于2的话就按照格式输出&#xff0c;否则就不是我们想要的梅森素数。 代码 /*Author : lifehappy */ // #pragma GCC optimize(2) // #pragma GCC optimize(3) // #include <bits/stdc.h>…

数据结构为什么那么难?

来源 | 异步 | 文末赠书2017年8月&#xff0c;本着让更多的人轻松学习算法的初心&#xff0c;我写作了第一本书《趣学算法》&#xff0c;该书在出版后受到广大读者一致好评&#xff0c;在一年内重印了10次&#xff0c;并输出了繁体版的版权。一位读者对我说&#xff0c;读这本书…

书籍推荐:《C#7.0本质论》

在dotNet平台中有多种开发语言可以使用&#xff0c;C#无疑是其中应用得最为广泛的。学习一门编程语言最好的方式就是找一本好书系统地学习&#xff0c;我读过的关于C#的书籍中&#xff0c;我认为下面三本最为经典&#xff1a;《C#本质论》&#xff1a;入门类&#xff0c;目前最…

乌龟棋(dp)

乌龟棋 思路 最优值问题&#xff0c;显然可以通过dpdpdp解决&#xff0c;我们定义dp[i][j][k][l]dp[i][j][k][l]dp[i][j][k][l]表示到达1i∗2∗j3∗k4∗l1 i * 2 * j 3 * k 4 * l1i∗2∗j3∗k4∗l这个点之前已经走过的价值最大的值(i,j,k,li, j, k, li,j,k,l分别是走一步&…

CF 1631 D. Range and Partition 尺取 + 前缀和

传送门 文章目录目录&#xff1a;题意&#xff1a;思路&#xff1a;目录&#xff1a; 题意&#xff1a; 给你一个长度为nnn的数组aaa以及kkk&#xff0c;让你选择一个值域[x,y][x,y][x,y]&#xff0c;满足能将该数组分成连续的kkk段并且每段中值域在[x,y][x,y][x,y]内的个数严…

gRPC的简单使用

前言八月初的时候&#xff0c;在公司内部做了一个主题为《gRPC的简单使用》的分享&#xff0c;其实就是和小伙伴们扯扯淡&#xff0c;现在抽空回忆一下&#xff0c;也算是一个小小的总结吧。现在市面上耳熟能详的RPC框架也很多&#xff0c;下面列举几个遇到比较多的。谷歌的gRP…

P6282 [USACO20OPEN] Cereal S 思维

传送门 文章目录目录&#xff1a;题意&#xff1a;思路&#xff1a;目录&#xff1a; 题意&#xff1a; 你有nnn头奶牛&#xff0c;mmm中不同种类的麦片&#xff0c;每个麦片只有一箱&#xff0c;给你每个奶牛第一和第二喜欢的麦片&#xff0c;奶牛会先看第一喜欢的是否有&am…

HDU 4417 Super Mario(莫队 + 树状数组 + 离散化)

Super Mario 思路 区间查找问题&#xff0c;容易想到离线莫队&#xff0c;确实这题就是莫队&#xff0c;接下来我们考虑如何维护区间高度值问题。 既然是离线嘛&#xff0c;我们容易想到离散化和他挂钩&#xff0c;想想这题是否需要离散化&#xff0c;高度的最大值是1000000…

生命周期结束,Spring Boot 1.x退役

一年前 Spring 官方宣布 Spring Boot 1.x 生命周期将于今年 8 月 1 日结束&#xff0c;如今时间已到&#xff0c;在发布 Spring Boot 1.5.22 的同时&#xff0c;Spring 确认将不再为 1.x 系列发布维护版本。官方希望用户尽快迁移到 Spring Boot 2.x 上&#xff0c;为此还制作了…

P3085 [USACO13OPEN]Yin and Yang G 点分治

文章目录题意&#xff1a;思路&#xff1a;传送门 题意&#xff1a; 给你一颗nnn个点的树&#xff0c;每条边为黑色或者白色&#xff0c;问满足以下条件的路径条数&#xff1a;路径上存在一个不是端点的点&#xff0c;使得两端点到该点的两条路径上两种颜色的边数相等。 1≤n…

G. Xor-MST(异或最小生成树)

G. Xor-MST 思路 异或最小生成树&#xff0c;这里采用了一种分治的方法来贪心求解最值&#xff1a; 首先我们对所有的点权值从小到大排个序&#xff0c;从高位开始在中间找到一个这个位置上的0&#xff0c;10&#xff0c;10&#xff0c;1分界点分成两个集合&#xff0c;然后…

CF 1638 E. Colorful Operations set 区间平推

文章目录题意&#xff1a;思路&#xff1a;传送门 题意&#xff1a; 给你一个数组aaa&#xff0c;初始价值全为000&#xff0c;颜色全为111&#xff0c;让后让你实现以下三个操作&#xff1a; 将[l,r][l,r][l,r]区间内的颜色都染成ccc。将所有颜色为ccc的位置价值都加上xxx。…

B Graph(异或最小生成树)

Graph 思路 图是联通的&#xff0c;并且加边的时候要保证环一定是异或值为0&#xff0c;所以我们可以保证从一个点到另一个点的路径异或值是不变的&#xff0c;这个时候就简单了&#xff0c;不就是一个异或最小生成树了嘛。 我们只要预处理一下&#xff0c;任选一个点作为根…

Apollo 配置中心:分布式部署

Apollo&#xff08;阿波罗&#xff09;是携程框架部门研发的分布式配置中心&#xff0c;能够集中化管理应用不同环境、不同集群的配置&#xff0c;配置修改后能够实时推送到应用端&#xff0c;并且具备规范的权限、流程治理等特性&#xff0c;适用于微服务配置管理场景。服务端…

CF 1637 D. Yet Another Minimization Problem dp + 思维

文章目录题意&#xff1a;思路&#xff1a;传送门 题意&#xff1a; 给你两个长度为nnn的数组a,ba,ba,b&#xff0c;每次交换可以选择一个iii&#xff0c;交换ai,bia_i,b_iai​,bi​&#xff0c;最小化∑i1n∑ji1n(aiaj)2∑i1n∑ji1n(bibj)2\sum_{i1}^n \sum_{ji1}^n(a_ia_j)…

各种逆元推导

逆元 求解一&#xff08;费马小定理&#xff09; ppp是一个质数&#xff0c;并且a%p̸0a \% p \not 0a%p​0&#xff0c;则有ap−1≡1(modp)a ^ {p - 1} \equiv 1 \pmod pap−1≡1(modp)&#xff0c;ap−2≡a−1a ^ {p - 2} \equiv a ^ {-1}ap−2≡a−1&#xff0c;即可得到…