Всем добрый день! Стоит задача решить слау, при этом число строк и столбцов не обязательно одинаковое (решаю задачу с тимуса 1580: http://acm.timus.ru/problem.aspx?space=1&num=1580). Матрица коэффициентов имеет значения только 1 или 0, а правые части могут быть любыми. Это завело меня в тупик связанный с тем, что на некоторые входные данные алгоритм дает NaN. Я разобрался что это из-за того что при приведении к треугольному виду используется деление на 0. Подскажите пожалуйста, как обойти это? Я нашел что есть т.н. метод Барейса, но литературы про него крайне мало. Помогите исправить код, буду очень благодарен:
using System;
class Program
{
static void Main(string[] args)
{
//INPUT
string[] inLine = Console.ReadLine().Split();
int N = int.Parse(inLine[0]); //N - число неизвестных (число студентов)
int M = int.Parse(inLine[1]); //M - число уравнений (число пар студентов пришедших к декану)
double[,] A = new double[M, N+1]; //A - матрица коэффициентов
for (int i = 0; i < M; i++)
{
inLine = Console.ReadLine().Split();
A[i, int.Parse(inLine[0]) - 1] = 1;
A[i, int.Parse(inLine[1]) - 1] = 1;
A[i, N] = int.Parse(inLine[2]);
}
if (M < N)
Console.WriteLine("IMPOSSIBLE");
else
{
double[] result = GaussMethod(A);
for (int i = 0; i < M; i++)
if (double.IsNaN(result[i]))
{
Console.WriteLine("IMPOSSIBLE");
return;
}
for (int j = 0; j < result.Length; j++)
Console.WriteLine("{0:0.00}", result[j]);
}
}
static double[] GaussMethod(double[,] a)
{
//найти размер матрицы
int width = a.GetLength(1);
int height = a.GetLength(0);
double[] x = new double[height]; //массив решений уравнения
//инициализация массива корней уравнения
for (int i = 0; i < x.Length; i++)
{
x[i] = a[i, width - 1];
}
//приведение матрицы к треугольному виду
double m;
for (int i = 0; i < width - 1; i++)
for (int j = i + 1; j < height; j++)
{
m = a[j, i] / a[i, i];
for (int k = 0; k < height; k++)
a[j, k] -= m * a[i, k];
x[j] -= m * x[i];
}
for (int i = 0; i < height; i++)
a[i, width - 1] = x[i];
for (int i = height - 1; i >= 0; i--)
{
for (int j = i + 1; j < width - 1; j++)
{
x[i] -= i != j ? a[i, j] * x[j] : 0;
}
x[i] = x[i] / a[i, i];
}
return x;
}
}