一定有解的高斯消元

#include<bits/stdc++.h>
using namespace std;
const int N = 1005;
int n;
double a[N][N];
signed main()
{
    scanf("%d", &n);
    for(int i=1; i<=n; i++)
        for(int j=1; j<=n+1; j++)
            scanf("%lf", a[i]+j);
    for(int i=1; i<=n; i++)
    {
        int p = i;
        for(int j=i+1; j<=n; j++)
            if(fabs(a[j][i])>fabs(a[p][i])) p=j;
        swap(a[i], a[p]);
        for(int j=1; j<=n; j++)
        {
            if(j==i) continue;
            double r = a[j][i] / a[i][i];
            for(int k=i; k<=n+1; k++)
                a[j][k] -= a[i][k] * r;
        }
    }
    for(int i=1; i<=n; i++)
        printf("%d ", (int) round(a[i][n+1] / a[i][i]));
}
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容