高斯消元模板（浮点数）

最新推荐文章于 2021-02-13 05:27:41 发布

maze_illusion

最新推荐文章于 2021-02-13 05:27:41 发布

阅读量283

点赞数

CC 4.0 BY-SA版权

分类专栏：数论模板

本文链接：https://blog.youkuaiyun.com/Game_Acm/article/details/82079292

数论模板专栏收录该内容

12 篇文章

订阅专栏

博客提供了浮点数高斯消元模板，高斯消元是重要的数值计算方法，在解线性方程组等方面有应用，该模板可助力相关计算。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

const double eps = 1e-9;
const int maxn = 210;
double a[maxn][maxn],x[maxn];
int equ,var;
int Gauss()
{
    int k,col,max_r;
    for ( k=0,col=0 ; k<equ&&col<var ; k++,col++ )
    {
        max_r = k;
        for ( int i=k+1 ; i<equ ; i++ )
            if ( fabs( a[i][col] )>fabs( a[max_r][col] ) )
                max_r = i;
        if ( fabs( a[max_r][col] )<eps ) return 0;
        if ( k!=max_r )
        {
            for ( int j=col ; j<var ; j++ )
                swap ( a[k][j] , a[max_r][j] );
            swap ( x[k] , x[max_r] );
        }
        x[k] /= a[k][col];
        for ( int j=col+1 ; j<var ; j++ )
            a[k][j] /= a[k][col];
        a[k][col] = 1;
        for ( int i=0 ; i<equ ; i++ )
            if ( i!=k )
            {
                x[i] -= x[k]*a[i][col];
                for ( int j=col+1 ; j<var ; j++ )
                    a[i][j] -= a[k][j]*a[i][col];
                a[i][col] = 0;
            }
    }
    return 1;
}