Метод Гаусса - История изменений

Ctrlalt в 19:22, 4 июня 2023

2023-06-04T19:22:50Z

← Предыдущая версия		Версия от 19:22, 4 июня 2023
Строка 17:		Строка 17:
	int bestRow = row;		int bestRow = row;
	for (int i = row + 1; i < n; i++)		for (int i = row + 1; i < n; i++)
	if (~~fabs~~(a[i][col]) > ~~fabs~~(a[bestRow][col]))		if (abs(a[i][col]) > abs(a[bestRow][col]))
	bestRow = i;		bestRow = i;
	if (~~fabs~~(a[bestRow][col]) < 1e-9)		if (abs(a[bestRow][col]) < 1e-9)
	continue;		continue;
	for (int j = 0; j < m + 1; j++)		for (int j = 0; j < m + 1; j++)
Строка 39:		Строка 39:
	for (int col = 0; col < m; col++)		for (int col = 0; col < m; col++)
	s += a[row][col] * ans[col];		s += a[row][col] * ans[col];
	if (~~fabs~~(s - a[row][m]) > 1e-9) {		if (abs(s - a[row][m]) > 1e-9) {
	printf("NO");		printf("NO");
	return 0;		return 0;
Строка 60:		Строка 60:
	int bestRow = row;		int bestRow = row;
	for (int i = row + 1; i < eqCount; i++)		for (int i = row + 1; i < eqCount; i++)
	if (~~fabs~~(a[i][col]) > ~~fabs~~(a[bestRow][col]))		if (abs(a[i][col]) > abs(a[bestRow][col]))
	bestRow = i;		bestRow = i;

	if (~~fabs~~(a[bestRow][col]) < 1e-9)		if (abs(a[bestRow][col]) < 1e-9)
	continue;		continue;

Строка 92:		Строка 92:
	sum += a[row][col] * ans[col];		sum += a[row][col] * ans[col];

	if (~~fabs~~(sum - a[row].back()) > 1e-9)		if (abs(sum - a[row].back()) > 1e-9)
	return {};		return {};
	}		}

Ctrlalt в 20:14, 24 августа 2021

2021-08-24T20:14:28Z

@@ Строка 97: / Строка 97: @@
       return ans;
+  }
@@ Строка 102: / Строка 131: @@
 Теория:
 * [http://e-maxx.ru/algo/linear_systems_gauss e-maxx.ru — Метод Гаусса решения системы линейных уравнений]
 Код:
 * [https://github.com/indy256/codelibrary/blob/master/cpp/linearalgebra/gauss.cpp indy256/codelibrary/cpp/linearalgebra/gauss.cpp]
 * [http://github.com/ADJA/algos/blob/master/NumberTheory/Gauss.cpp ADJA/algos/NumberTheory/Gauss.cpp]

Ctrlalt в 23:05, 6 августа 2021

2021-08-06T23:05:42Z

@@ Строка 51: / Строка 51: @@
       for (int i = 0; i < m; i++)
           printf("%.10lf ", ans[i]);
+  }
 == Ссылки ==
 Теория:
-* [http://e-maxx.ru/algo/linear_systems_gauss e-maxx.ru &mdash; Метод Гаусса решения системы линейных уравнений]
+* [http://e-maxx.ru/algo/linear_systems_gauss e-maxx.ru — Метод Гаусса решения системы линейных уравнений]
 Код:
-* [http://github.com/indy256/codelibrary/blob/master/java/src/Gauss.java CodeLibrary &mdash; Gaussian elimination algorithm in O(N^3)]
+* [https://github.com/indy256/codelibrary/blob/master/cpp/linearalgebra/gauss.cpp indy256/codelibrary/cpp/linearalgebra/gauss.cpp]
-* [http://github.com/ADJA/algos/blob/master/NumberTheory/Gauss.cpp Algos &mdash; Gauss method of solving systems of linear algebraic equation.]
+* [http://github.com/ADJA/algos/blob/master/NumberTheory/Gauss.cpp ADJA/algos/NumberTheory/Gauss.cpp]

Ctrlalt: /* Ссылки */

2015-07-16T03:49:01Z

Ссылки

← Предыдущая версия		Версия от 03:49, 16 июля 2015
Строка 56:		Строка 56:
	Теория:		Теория:
	* [http://e-maxx.ru/algo/linear_systems_gauss e-maxx.ru — Метод Гаусса решения системы линейных уравнений]		* [http://e-maxx.ru/algo/linear_systems_gauss e-maxx.ru — Метод Гаусса решения системы линейных уравнений]
			Код:
	* [http://github.com/indy256/codelibrary/blob/master/java/src/Gauss.java CodeLibrary — Gaussian elimination algorithm in O(N^3)]		* [http://github.com/indy256/codelibrary/blob/master/java/src/Gauss.java CodeLibrary — Gaussian elimination algorithm in O(N^3)]
	~~Код:~~
	* [http://github.com/ADJA/algos/blob/master/NumberTheory/Gauss.cpp Algos — Gauss method of solving systems of linear algebraic equation.]		* [http://github.com/ADJA/algos/blob/master/NumberTheory/Gauss.cpp Algos — Gauss method of solving systems of linear algebraic equation.]

Ctrlalt: /* Ссылки */

2015-07-16T03:48:40Z

Ссылки

← Предыдущая версия		Версия от 03:48, 16 июля 2015
Строка 54:		Строка 54:

	== Ссылки ==		== Ссылки ==
			Теория:
	* [http://e-maxx.ru/algo/linear_systems_gauss e-maxx.ru — Метод Гаусса решения системы линейных уравнений]		* [http://e-maxx.ru/algo/linear_systems_gauss e-maxx.ru — Метод Гаусса решения системы линейных уравнений]
	* [http://github.com/indy256/codelibrary/blob/master/java/src/Gauss.java CodeLibrary — Gaussian elimination algorithm in O(N^3)]		* [http://github.com/indy256/codelibrary/blob/master/java/src/Gauss.java CodeLibrary — Gaussian elimination algorithm in O(N^3)]
			Код:
	* [http://github.com/ADJA/algos/blob/master/NumberTheory/Gauss.cpp Algos — Gauss method of solving systems of linear algebraic equation.]		* [http://github.com/ADJA/algos/blob/master/NumberTheory/Gauss.cpp Algos — Gauss method of solving systems of linear algebraic equation.]

Ctrlalt: Новая страница: « #include #include #include using namespace std; int n, m, ansRow[110]; double a[110][110], ans[110]; int main() { s…»

2015-03-14T04:22:08Z

Новая страница: « #include <stdio.h> #include <math.h> #include <algorithm> using namespace std; int n, m, ansRow[110]; double a[110][110], ans[110]; int main() { s…»

Новая страница

#include <stdio.h>
#include <math.h>
#include <algorithm>
using namespace std;

int n, m, ansRow[110];
double a[110][110], ans[110];

int main() {
scanf("%d%d", &n, &m);
for (int i = 0; i < n; i++)
for (int j = 0; j < m + 1; j++)
scanf("%lf", &a[i][j]);

fill(ansRow, ansRow + m, -1);
for (int col = 0, row = 0; col < m && row < n; col++) {
int bestRow = row;
for (int i = row + 1; i < n; i++)
if (fabs(a[i][col]) > fabs(a[bestRow][col]))
bestRow = i;
if (fabs(a[bestRow][col]) < 1e-9)
continue;
for (int j = 0; j < m + 1; j++)
swap(a[row][j], a[bestRow][j]);
ansRow[col] = row;
for (int i = 0; i < n; i++) {
if (i == row)
continue;
double k = -a[i][col] / a[row][col];
for (int j = col; j < m + 1; j++)
a[i][j] += k * a[row][j];
}
row++;
}
for (int col = 0; col < m; col++)
ans[col] = (ansRow[col] != -1 ? a[ansRow[col]][m] / a[ansRow[col]][col] : 0);
for (int row = 0; row < n; row++) {
double s = 0;
for (int col = 0; col < m; col++)
s += a[row][col] * ans[col];
if (fabs(s - a[row][m]) > 1e-9) {
printf("NO");
return 0;
}
}
if (count(ansRow, ansRow + m, -1)) {
printf("INF");
return 0;
}
printf("YES\n");
for (int i = 0; i < m; i++)
printf("%.10lf ", ans[i]);
}

== Ссылки ==
* [http://e-maxx.ru/algo/linear_systems_gauss e-maxx.ru — Метод Гаусса решения системы линейных уравнений]
* [http://github.com/indy256/codelibrary/blob/master/java/src/Gauss.java CodeLibrary — Gaussian elimination algorithm in O(N^3)]
* [http://github.com/ADJA/algos/blob/master/NumberTheory/Gauss.cpp Algos — Gauss method of solving systems of linear algebraic equation.]