글 읽기 - 채점 90%에서 틀리다고 나오는데 혹시 맞으신 분 가르침 부탁드립니다.

(追記) (追記ここまで)

채점 90%에서 틀리다고 나오는데 혹시 맞으신 분 가르침 부탁드립니다.

9254번 - 역행렬

jcy9503 2년 전 0

테스트 케이스가 상당히 한정적이라 구글링하여 역행렬 관련 예제를 모조리 넣어보아도 정상적으로 출력됨을 확인할 수 있었습니다.

오차 범위가 문제인가? 해서 long double 자료형을 바꾸어보았더니 시간 초과가 나서 패스했구요.

아무래도 n의 값이 엄청 큰 케이스에서 문제가 생기는 거 같은데 혹시 코드 상에 어떤 점이 문제일 수 있는지 힌트라도 주시면 감사하겠습니다.

#include <iostream>
using namespace std;
long num[500][500];
long result[500][500];
void Swap(const int n, const int row1, const int row2)
{
	long tp1[500];
	long tp2[500];
	for (int i = 0; i < n; ++i)
	{
		tp1[i] = num[row1][i];
		tp2[i] = result[row1][i];
		num[row1][i] = num[row2][i];
		result[row1][i] = result[row2][i];
		num[row2][i] = tp1[i];
		result[row2][i] = tp2[i];
	}
}
void Multiply(const int n, const int row, const long k)
{
	for (int i = 0; i < n; ++i)
	{
		num[row][i] *= k;
		result[row][i] *= k;
	}
}
void Divide(const int n, const int row, const long k)
{
	for (int i = 0; i < n; ++i)
	{
		num[row][i] /= k;
		result[row][i] /= k;
	}
}
void AddMult(const int n, const int from, const int to, const long k)
{
	for (int i = 0; i < n; ++i)
	{
		num[to][i] += num[from][i] * k;
		result[to][i] += result[from][i] * k;
	}
}
void DiagonalOne(const int n)
{
	for (int i = 0; i < n; ++i)
	{
		if (num[i][i] != 1.0)
			Divide(n, i, num[i][i]);
	}
}
bool CheckInverse(const int n)
{
	for (int i = 0; i < n; ++i)
	{
		if (num[i][i] != 1.0) return false;
	}
	return true;
}
int main(int argc, char* argv[])
{
	int n;
	cin >> n;
	for (int i = 0; i < n; ++i)
	{
		for (int j = 0; j < n; ++j)
			cin >> num[i][j];
	}
	for (int i = 0; i < n; ++i)
		result[i][i] = 1;
	bool flag = true;
	// Sorting Rows to Prevent Containing Diagonal Elements having 0 value.
	for (int i = 0; i < n; ++i)
	{
		if (num[i][i] == 0.0)
		{
			for (int j = 0; j < n; ++j)
			{
				if (j == i) continue;
				if (num[i][j] != 0.0 && num[j][i] != 0.0)
				{
					Swap(n, i, j);
					break;
				}
			}
		}
	}
	DiagonalOne(n);
	// Transformation to Reduced Row Echelon Form
	for (int i = 0; i < n - 1; ++i)
	{
		for (int j = i + 1; j < n; ++j)
		{
			AddMult(n, i, j, -num[j][i]);
			Divide(n, j, num[j][j]);
		}
	}
	for (int i = n - 1; i >= 0; --i)
	{
		for (int j = i - 1; j >= 0; --j)
		{
			AddMult(n, i, j, -num[j][i]);
		}
	}
	if (!CheckInverse(n))
	{
		flag = false;
	}
	if (flag)
	{
		for (int i = 0; i < n; ++i)
		{
			for (int j = 0; j < n; ++j)
				cout << result[i][j] << ' ';
			cout << '\n';
		}
	}
	else
		cout << "no inverse\n";
	return 0;
}

댓글을 작성하려면 로그인해야 합니다.

(追記) (追記ここまで)

Baekjoon Online Judge

채점 현황

채점 현황

문제

유저 대회

고등학교 대회

출처

대학교 대회

도움말

사업자 등록 번호: 541-88-00682
대표자명: 최백준
주소: 서울시 서초구 서초대로74길 29 서초파라곤 412호
전화번호: 02-521-0487 (이메일로 연락 주세요)
이메일: contacts@startlink.io
통신판매신고번호: 제 2017-서울서초-2193 호

한국어 | English (Beta)

문제

문제

출처

ICPC