No Rules Rules

행렬 제곱 (feat. 백준, 10830번) 본문

생활/코테

행렬 제곱 (feat. 백준, 10830번)

개발하는 완두콩 2022. 8. 15. 23:04
728x90
반응형

행렬 제곱
https://www.acmicpc.net/problem/10830

 

10830번: 행렬 제곱

크기가 N*N인 행렬 A가 주어진다. 이때, A의 B제곱을 구하는 프로그램을 작성하시오. 수가 매우 커질 수 있으니, A^B의 각 원소를 1,000으로 나눈 나머지를 출력한다.

www.acmicpc.net

 

반응형

 

// woohyeon.kim
// kim519620.tistory.com
#include <iostream>
using namespace std;
long long tmp[5][5];
int N;
void solution(long long(*arr1)[5], long long(*arr2)[5]) {
	for (register int i = 0, j, k; i < N; ++i)
		for (j = 0; j < N; ++j) {
			tmp[i][j] = 0;
			for (k = 0; k < N; ++k)
				tmp[i][j] += arr1[i][k] * arr2[k][j];
			tmp[i][j] %= 1000;
		}
	for (register int i = 0, j; i < N; ++i)
		for (j = 0; j < N; ++j)
			arr1[i][j] = tmp[i][j];
}
int main() {
	ios::sync_with_stdio(false), cin.tie(NULL);
	register long long B, arr[5][5]{ 0 }, ans[5][5]{ 0 };
	cin >> N >> B;
	for (register int i = 0, j; i < N; ++i) {
		for (j = 0; j < N; ++j)
			cin >> arr[i][j];
		ans[i][i] = 1;
	}
	while (B) {
		if (B % 2)
			solution(ans, arr);
		solution(arr, arr);
		B /= 2;
	}
	for (register int i = 0, j; i < N; ++i) {
		for (j = 0; j < N; ++j)
			cout << ans[i][j] << " ";
		cout << "\n";
	}
	return 0;
}
// *&)*@*
728x90
반응형

'생활 > 코테' 카테고리의 다른 글

최소 힙 (feat. 백준, 1927번)  (0) 2022.08.16
최대 힙 (feat. 백준, 11279번)  (0) 2022.08.16
행렬 곱셈 (feat. 백준, 2740번)  (0) 2022.08.15
이항 계수 3 (feat. 백준, 11401번)  (0) 2022.08.15
곱셈 (feat. 백준, 1629번)  (0) 2022.08.15
Comments