神奇代码——高斯消元求调
查看原帖
神奇代码——高斯消元求调
1265632
xhxxwcr楼主2025/1/21 16:04
#include<bits/stdc++.h>

using namespace std;
double a[20][20], b[20], c[20][20];

int main() {
    int n;
    cin >> n;
    for (int i = 1; i <= n + 1; i++) {
        for (int j = 1; j <= n; j++) {
            cin >> a[i][j];
        }
    }
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= n; j++) {
            c[i][j] = 2 * (a[i][j] - a[i + 1][j]);
            b[i] += a[i][j] * a[i][j] - a[i - 1][j] * a[i - 1][j];
        }
    }
    for (int i = 1; i <= n; i++) {
        for (int j = i; j <= n; j++) {
            if (fabs(c[j][i]) > 1e-8) {
                for (int k = 1; k <= n; k++) {
                    swap(c[i][k], c[j][k]);
                }
                swap(b[i], b[j]);
            }
        }
        for (int j = 1; j <= n; j++) {
            if (i == j) continue;
            double _ = c[j][i] / c[i][i];
            for (int k = 1; k <= n; k++) {
                c[j][k] -= c[i][k] * _;
            }
            b[j] -= b[i] * _;
        }
    }
    for (int i = 1; i <= n; i++) printf("%.3lf ", b[i] / c[i][i]);
    printf("\n");
    return 0;
}

蒙了

2025/1/21 16:04
加载中...