高斯消元

ACM模版

高斯消元

/* * 高斯消元(浮点数) */
#define eps 1e-9
const int MAXN = 220;

double a[MAXN][MAXN];   // 方程的左边的矩阵
double x[MAXN];         // 等式右边的值,求解后x存的就是结果
int equ, var;           // 方程数和未知数个数

// 返回0无解,返回1有解
int Gauss()
{
    int i, j, k, col, MAX;
    for (k = 0, col = 0; k < equ && col < var; k++, col++)
    {
        MAX = k;
        for (i = k + 1; i < equ; i++)
        {
            if (fabs(a[i][col]) > fabs(a[MAX][col]))
            {
                MAX = i;
            }
        }
        if (fabs(a[MAX][col]) < eps)
        {
            return 0;
        }
        if (k != MAX)
        {
            for (j = col; j < var; j++)
            {
                swap(a[k][j], a[MAX][j]);
            }
            swap(x[k], x[MAX]);
        }
        x[k] /= a[k][col];
        for (j = col + 1; j < var; j++)
        {
            a[k][j] /= a[k][col];
        }
        a[k][col] = 1;
        for (i = 0; i < equ; i++)
        {
            if (i != k)
            {
                x[i] -= x[k] * a[i][k];
                for (j = col + 1; j < var; j++)
                {
                    a[i][j] -= a[k][j] * a[i][col];
                }
                a[i][col] = 0;
            }
        }
    }
    return 1;
}
全部评论

相关推荐

不愿透露姓名的神秘牛友
07-07 18:05
哈哈哈哈哈感觉朋友找工作的已经疯掉了,直接上图
码农索隆:真老板娘:“我嘞个去,这不我当年的套路吗
点赞 评论 收藏
分享
牛客刘北:如果暑期实习是27届的话,你要晚一年才会毕业,企业为什么会等你呢?要搞清时间逻辑呀!27届现在实习只能是在暑假实习,这是日常实习,不是暑期实习。所以多去投日常实习吧,暑期实习肯定不会要你的
点赞 评论 收藏
分享
评论
点赞
收藏
分享

创作者周榜

更多
牛客网
牛客网在线编程
牛客网题解
牛客企业服务