矩阵快速幂

最新推荐文章于 2024-06-02 15:09:13 发布

小刀刺大熊

最新推荐文章于 2024-06-02 15:09:13 发布

阅读量133

点赞数

分类专栏：数学文章标签： c++

本文链接：https://blog.csdn.net/seanbill/article/details/131437499

版权

数学专栏收录该内容

16 篇文章 0 订阅

订阅专栏

该代码示例展示了如何使用C++编程实现矩阵快速幂算法，用于高效地进行矩阵乘法。程序首先定义了一个Matrix结构体，包含了矩阵乘法、初始化、输入和输出等方法，然后在主函数中读取矩阵并进行快速幂运算，最后输出结果。

摘要由CSDN通过智能技术生成

#include <iostream>
#include <istream>
#include <sstream>
#include <vector>
#include <stack>
#include <list>
#include <map>
#include <set>
#include <deque>
#include <queue>
#include <cstring>
#include <unordered_map>
#include <unordered_set>
#include <algorithm>
#include <numeric>
#include <chrono>
#include <ctime>
#include <cmath>
#include <cctype>
#include <string>
#include <cstdio>
#include <iomanip>


#include <thread>
#include <mutex>
#include <condition_variable>
#include <functional>
#include <iterator>
using namespace std;

inline int read(int& x) {
	char ch = getchar();
	int f = 1; x = 0;
	while (ch > '9' || ch < '0') { if (ch == '-')f = -1; ch = getchar(); }
	while (ch >= '0' && ch <= '9') { x = (x << 1) + (x << 3) + ch - '0'; ch = getchar(); }
	return x * f;
}

const int MOD = 1e9 + 7;
long long int n,m;

struct Matrix {
	int a[202][202];
	Matrix() { memset(a, 0, sizeof(a)); }

	void Init() {
		for (int i = 0; i < n; i++)a[i][i] = 1;
	}

	Matrix operator*(const Matrix& input)const {
		Matrix ret;
		for (int i = 0; i < n; i++) {
			for (int j = 0; j < n; j++) {
				for (int k = 0; k < n; k++) {
					ret.a[i][j] = (ret.a[i][j] + (long long int)a[i][k] * input.a[k][j]) % MOD;
				}
			}
		}
		return ret;
	}
	void Input() {
		for (int i = 0; i < n; i++) {
			for (int j = 0; j < n; j++) {
				cin >> a[i][j];
			}
		}
	}
	void OutPut() {
		for (int i = 0; i < n; i++) {
			for (int j = 0; j < n; j++) {
				cout << a[i][j] << " ";
			}
			cout << endl;
		}
	}
}A, res;

void quickPow(long long int k) {
	res.Init();
	while (k) {
		if (k & 1) res = res * A;
		A = A * A;
		k >>= 1;
	}
	res.OutPut();
}

int main()
{
	cin >> n >> m;
	A.Input();
	quickPow(m);
	return 0;
}