學習 模板 高斯消元

2022-04-28 17:30:09 字數 1020 閱讀 8646

邵哥日推題目使我產生學習高斯消元的想法qaq

學習的資料qwq

2. 小東西quq

**寶寶qaq

1 #include2 #include3 #include4

using

namespace

std;

5const

int sz = 111;6

intn;

7double eps = 1e-7;8

double

plat[sz][sz], ans[sz];

9int

main()

21if(fabs(plat[flag][i]) < eps)

25if(i !=flag) swap(plat[i], plat[flag]);

26//

對換一行,找最大當前係數, 這樣可以只處理當前行的係數

27double div = plat[i][i];//

當前行的主元

28for(int j = i; j <= n+1; j++)

29 plat[i][j] /= div;//

30for(int j = i+1; j <= n; j++) //

手玩一下就好啦quq

35 }//

消元部分

36 ans[n] = plat[n][n+1

];37

//經過上述操作, 最後乙個(設為z)係數為一,答案為plat[n][n+1]

38for(int i = n-1; i >= 1; i--) //

回帶操作qaq依舊需要手玩

43for(int i = 1; i <= n; i++)

44 printf("

%.2lf\n

", ans[i]);//

終於完結, 沒勁撒花qaq

45return0;

46}47//

希望下次看自己還能看懂orz

高斯消元模板

include include include include include include include include typedef long long ll const int n 1008 高斯消元模板 const double eps 1e 12 double aug n n 增廣矩...

高斯消元模板

用迭代的辦法打會簡潔一些。有些精度上的細節需要注意。多次消元要清空use和cho陣列。實數高斯消元 int use maxn cho maxn void solve equation int n,int m break fd i,n,1 if cho i 自由元的個數就是cho為0的個數 無解的話去...

高斯消元模板

y,z.x,y,z.x,y,z.可以得到乙個上三角矩陣。這時已經可以直接讀出最後乙個方程的解了,再回代到上面的方程就行了。具體實現 對於x ixi xi,找到xixi xi係數最大的乙個方程,以減少精度誤差。然後,把該方程的xixi xi係數轉化為1,帶入後面所有方程消元。考慮倒著往回進行代入消元,...