題目鏈接：點擊查看

題目大意：給出一個線性方程組，求解答案

題目分析：模板題，掛個模板，時間復雜度是 O( n^3 ) 的，其中 a[ 1?][ 1 ] ~ a[ n ][ n ] 矩陣代表 x[ 1 ] ~ x[ n ] 在每個方程中的系數，a[ 1 ][ n + 1 ] ~ a[ n ][ n + 1 ] 代表的是等號右邊的 b

代碼：
?

#include<iostream> #include<cstdio> #include<string> #include<ctime> #include<cmath> #include<cstring> #include<algorithm> #include<stack> #include<climits> #include<queue> #include<map> #include<set> #include<sstream> #include<cassert> #include<bitset> using namespace std;typedef long long LL;typedef unsigned long long ull;const int inf=0x3f3f3f3f;const int N=110;double a[N][N];int n;bool Gauss() {for(int i=1;i<=n;++i)//枚舉列（項） {int max=i;for(int j=i+1;j<=n;++j)//選出該列最大系數 {if(fabs(a[j][i])>fabs(a[max][i]))//fabs是取浮點數的絕對值的函數{max=j;}}for(int j=1;j<=n+1;++j)//交換{swap(a[i][j],a[max][j]);}if(!a[i][i])//最大值等于0則說明該列都為0，肯定無解 return false;for(int j=1;j<=n;++j)//每一項都減去一個數（就是小學加減消元）{if(j!=i){double temp=a[j][i]/a[i][i];for(int k=i+1;k<=n+1;++k){a[j][k]-=a[i][k]*temp;//a[j][k]-=a[j][i]*a[i][k]/a[i][i];}}}}//上述操作結束后，矩陣會變成這樣/*k1*a=e1k2*b=e2k3*c=e3k4*d=e4*///所以輸出的結果要記得除以該項系數，消去常數for(int i=1;i<=n;++i)a[i][n+1]/=a[i][i];return true; }int main() { #ifndef ONLINE_JUDGE // freopen("data.in.txt","r",stdin); // freopen("data.out.txt","w",stdout); #endif // ios::sync_with_stdio(false);scanf("%d",&n);for(int i=1;i<=n;i++)for(int j=1;j<=n+1;j++)scanf("%lf",&a[i][j]);if(!Gauss())return 0*puts("No Solution");for(int i=1;i<=n;i++)printf("%.2f\n",a[i][n+1]);return 0; }

?

總結

以上是生活随笔為你收集整理的洛谷 - P3389 【模板】高斯消元法(高斯消元解方程)的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。