邵哥日推題目使我產生學習高斯消元的想法qaq
學習的資料qwq
2. 小東西quq
**寶寶qaq
1 #include2 #include3 #include4using
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,帶入後面所有方程消元。考慮倒著往回進行代入消元,...