高斯消元
https://www.luogu.org/problemnew/show/P3389
知识点:1.一个一个未知数的消
2.将未消的未知数的主式(用来消别的式子中的未知数的式子)的未消的未知数的系数先化为1,再与后面的式子加减消元
3.将未消的未知数的系数的绝对值最大的式子当主式
code:
#include <bits/stdc++.h> #define M 110 #define eps 1e-10 using namespace std; int n; double f[M][M]; double ans[M]; int main() { scanf("%d",&n); for(int i = 1;i <= n;i++) for(int j = 1;j <= n + 1;j++) scanf("%lf",&f[i][j]); for(int i = 1;i <= n;i++) { int r = i; for(int j = i + 1;j <= n;j++) if(fabs(f[r][i]) < fabs(f[j][i])) r = j; if(fabs(f[r][i]) < eps) { printf("No Solution"); return 0; } if(i != r)swap(f[i],f[r]); double div = f[i][i]; for(int j = i;j <= n + 1;j++) f[i][j] /= div; for(int j = i + 1;j <= n;j++) { div = f[j][i]; for(int k = i;k <= n + 1;k++) f[j][k] -= f[i][k] * div; } } ans[n] = f[n][n + 1]; for(int i = n - 1;i >= 1;i--) { ans[i] = f[i][n + 1]; for(int j = i + 1;j <= n;j++) ans[i] -= (f[i][j] * ans[j]); } for(int i = 1;i <= n;i++) printf("%.2lf\n",ans[i]); return 0; }