高斯消元

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;
}
全部评论

相关推荐

不愿透露姓名的神秘牛友
06-29 17:30
点赞 评论 收藏
分享
05-09 12:23
已编辑
华南理工大学 Java
野猪不是猪🐗:给他装的,双九+有实习的能看的上这种厂我直接吃⑨✌们拿它练练面试愣是给他整出幻觉了
点赞 评论 收藏
分享
湫湫湫不会java:1.在校经历全删了2.。这些荣誉其实也没啥用只能说,要的是好的开发者不是好好学生3.项目五六点就行了,一个亮点一俩行,xxx技术解决,xxx问题带来xxx提升。第一页学历不行,然后啥有价值的信息也没有,到第二页看到项目了,第一个项目九点,第二个项目像凑数的俩点。总体给人又臭又长,一起加油吧兄弟
点赞 评论 收藏
分享
评论
点赞
收藏
分享

创作者周榜

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