高斯消元模板

最新推荐文章于 2022-04-12 19:24:58 发布

HownoneHe

最新推荐文章于 2022-04-12 19:24:58 发布

阅读量527

点赞数 1

分类专栏： gauss 模板模板与算法文章标签： gauss matrix

本文链接：https://blog.csdn.net/HownoneHe/article/details/51355500

版权

模板与算法同时被 3 个专栏收录

25 篇文章 3 订阅

订阅专栏

模板

6 篇文章 0 订阅

订阅专栏

gauss

1 篇文章 0 订阅

订阅专栏

忘了（懒癌发作）可以提醒一下自己，嘿嘿嘿。。。

void gauss() {
    fo(i,1,n-1) {
        fo(j,i+1,n) if (abs(matrix[j].a[i])>abs(matrix[i].a[i])) 
        swap(matrix[i],matrix[j]);
        fo(j,i+1,n) {
            db f=matrix[j].a[i]/matrix[i].a[i];
            fo(k,i,n) matrix[j].a[k]=matrix[i].a[k]*f-matrix[j].a[k];
            matrix[j].b=matrix[i].b*f-matrix[j].b;
        }
    }
    fd(i,n,2) {
        fo(j,1,i-1) {
            db f=matrix[j].a[i]/matrix[i].a[i];
            fo(k,j,n) matrix[j].a[k]=matrix[i].a[k]*f-matrix[j].a[k];
            matrix[j].b=matrix[i].b*f-matrix[j].b;
        }
    }
    fo(i,1,n) an[i]=matrix[i].b/matrix[i].a[i];
}