51 Nod1113 矩阵快速幂

题目：给出一个N * N的矩阵，其中的元素均为正整数。求这个矩阵的M次方。由于M次方的计算结果太大，只需要输出每个元素Mod (10^9 + 7）的结果。

思路：这道题是快速幂乘法的升级版，普通快速幂是整数，这个只是转换成整数乘法而已。

#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;
typedef long long ll;
#define mod (1000000000+7)
struct node {ll a[105][105];
};    //构造矩阵结构体
int n;
node mulpity(node x, node y)//矩阵乘法
{node ans;memset(ans.a, 0, sizeof(ans.a));for(int i=0;i<n;i++)for(int j=0;j<n;j++)for (int k = 0; k < n; k++)//k表示x的行和b的列{ans.a[i][j] =(ans.a[i][j]+ (x.a[i][k] * y.a[k][j])%mod)%mod;}return ans;
}
node quikmul(node x, int m)
{node ans;memset(ans.a, 0, sizeof(ans.a));for (int i = 0; i < n; i++)ans.a[i][i] = 1;//构造单位矩阵while (m ){if (m & 1)ans = mulpity(ans,x);x = mulpity(x,x);m >>= 1;}return ans;}
int main()
{int m;node ans;while (~scanf("%d%d", &n, &m)){for (int i = 0; i < n; i++)for (int j = 0; j < n; j++)scanf("%lld", &ans.a[i][j]);               ans = quikmul(ans, m); for (int i = 0; i < n; i++){for (int j = 0; j < n; j++){if (j == 0)cout << ans.a[i][j];elsecout << " " << ans.a[i][j];}printf("\n");}}return 0;
}

51 Nod1113 矩阵快速幂相关推荐

矩阵快速幂 HDU3483
1 #include <iostream> 2 #include <cstring> 3 4 using namespace std; 5 6 //矩阵大小上限 7 const ...
快速幂 + 矩阵快速幂
快速幂 1 #include<iostream> 2 #include<algorithm> 3 #include<cstring> 4 #define LL lo ...
HDU 2243考研路茫茫——单词情结（AC自动机+矩阵快速幂）
背单词,始终是复习英语的重要环节.在荒废了3年大学生涯后,Lele也终于要开始背单词了. 一天,Lele在某本单词书上看到了一个根据词根来背单词的方法.比如"ab",放在单词前一般 ...
矩阵快速幂---BestCoder Round#8 1002
当要求递推数列的第n项且n很大时,怎么快速求得第n项呢? 可以用矩阵快速幂来加速计算. 我们可以用矩阵来表示数列递推公式比如fibonacci数列可以表示为 [f(n) f(n-1)] = [ ...
HDU4549 M斐波那契数列 —— 斐波那契、费马小定理、矩阵快速幂
题目链接:https://vjudge.net/problem/HDU-4549 M斐波那契数列 Time Limit: 3000/1000 MS (Java/Others) Memory Li ...
又见斐波那契~矩阵快速幂入门题
链接:https://www.nowcoder.com/acm/contest/105/G 来源:牛客网题目描述这是一个加强版的斐波那契数列. 给定递推式求F(n)的值,由于这个值可能太大,请对 ...
BZOJ4818 [SDOI2017] 序列计数【矩阵快速幂】
题目分析: 一个很显然的同类项合并.注意到p的大小最大为100,考虑把模p意义下相同的求出来最后所有的减去没有质数的做矩阵快速幂即可. 代码: 1 #include<bits/stdc++.h& ...
矩阵快速幂的一份小结
矩阵真是个好东西!虽然矩乘的复杂度有点难看... ... 这几天也做了不少矩阵题目,还是有几道好题目的.不过我打算从入门开始. 矩阵乘法:A[i][k]*B[k][j]=C[i][j];(A的第i行的 ...
[BZOJ3583]杰杰的女性朋友(矩阵快速幂)
杰杰的女性朋友时间限制:10s 空间限制:256MB 题目描述杰杰是魔法界的一名传奇人物.他对魔法具有深刻的洞察力,惊人的领悟力,以及令人叹为观止的创造力.自从他从事魔法竞赛以来,短短几 ...

51 Nod1113 矩阵快速幂

51 Nod1113 矩阵快速幂相关推荐

最新文章

热门文章