高斯消元模板啦啦啦..

最新推荐文章于 2023-07-24 18:50:40 发布

Ra1nbow_Chan

最新推荐文章于 2023-07-24 18:50:40 发布

阅读量388

点赞数

分类专栏：高斯消元模板

本文链接：https://blog.csdn.net/Ra1nBow_Chan/article/details/51500840

版权

模板同时被 2 个专栏收录

3 篇文章 0 订阅

订阅专栏

高斯消元

1 篇文章 0 订阅

订阅专栏

放一下高斯消元模板啦啦啦..

    for ( i = 1; i <= n; i ++ ){
        for ( j = i; j <= n; j ++ ){
            if ( _abs ( fc[j][i] ) > eps ){
                for ( k = i; k <= n+1; k ++ ){
                    swap ( fc[j][k], fc[i][k] );
                }
                break;
            }
        }
        for ( j = 1; j <= n; j ++ ){
            if ( i != j && _abs (fc[j][i]) > eps ){
                double tmp = fc[j][i]/fc[i][i];
                for ( k = i; k <= n+1; k ++ ) fc[j][k] -= tmp*fc[i][k];
            }
        }
    }
    for ( i = 1; i <= n; i ++ ){
        for ( j = 1; j <= n; j ++ ){
            if ( _abs (fc[j][i]) > eps ) ans[i] = fc[j][n+1]/fc[j][i];
        }
    }