标签:space int 题解 三角 db cin 线性代数 高斯消 P3389
简要题意:
给定一个由 \(n\) 个方程组成的 \(n\) 元方程组。若有唯一解则输出,否则输出 No Solution
.
前置知识:线性代数相关知识。
很明显,这是线性代数中求解 \(Ax = B\) 的模板题。
考虑实现标准做法,即把 \(A \space | \space B\) 化为上三角的形式。
因为唯一解必须是 \(r(A) = r(A \space | \space B) = n\),因此直接利用上三角形式 \(A_{n,n}\) 是否为 \(0\) 就可以判断。为 \(0\),则有无穷解(因为至少最后一行全为零,即有一个基础解系);不为 \(0\),则必有唯一解。
于是对于化为上三角形式的矩阵,一波代入即可。
对于每一行,我们都要先除掉一个系数,然后把按照下面每一行对应列的系数消掉该元。
时间复杂度:\(\mathcal{O}(n^3)\).
#include <bits/stdc++.h>
using namespace std;
const int N = 1e2 + 1;
#define db double
db A[N][N],B[N],x[N];
int n;
int main() {
cin >> n;
for(int i = 1; i <= n; i++) {
for(int j = 1; j <= n; j++)
cin >> A[i][j];
cin >> B[i];
}
for(int i = 1; i < n; i++) {
int f = 0;
for(int j = i; j <= n; j++)
if(A[i][j]) {f = j; break;}
if(!f) continue;
for(int j = i + 1; j <= n; j++) A[f][j] /= A[f][i];
B[f] /= A[f][i]; A[f][i] = 1;
for(int j = f + 1; j <= n; j++) {
for(int k = i + 1; k <= n; k++) A[j][k] -= A[f][k] * A[j][i];
B[j] -= B[f] * A[j][i]; A[j][i] = 0;
}
}
if(!A[n][n]) return puts("No Solution"),0;
x[n] = B[n] / A[n][n];
for(int i = n - 1; i >= 1; i--) {
db p = 0;
for(int j = i + 1; j <= n; j++) p += A[i][j] * x[j];
x[i] = B[i] - p;
}
for(int i = 1; i <= n; i++)
cout << fixed << setprecision(2) << x[i] << endl;
return 0;
}
标签:space,int,题解,三角,db,cin,线性代数,高斯消,P3389 来源: https://www.cnblogs.com/bifanwen/p/15831571.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。