高斯消元小结

这里只是丢了一个板子,毕竟高斯消元这个东西原理说起来很简单,就是模拟了普通人手工解方程的过程,还是直接上代码来的方便

一道模板题:luogu2455

(在这里不推荐luogu的模板题,数据过水,此题数据强度还可以需要适当的和精度搏斗

主要提一下判断无解和无穷解的情况

一般的高斯消元是给出\(n\)个方程求\(n\)个未知数的解,但是有的时候会出现给出方程重复(\(x+y=2,2x+2y=4\))或者方程矛盾(\(x+y=2,2x+2y=2\)),这个时候我们进行高斯消元的话,最后的一个(或几个)方程会出现系数全部为\(0\)的情况,需要特殊处理

具体的,对于方程组中的某一方程\(k_1x_1+k_2x_2+\cdots+k_nx_n=y\)来说,若满足\(\forall i, k_i=0\),那么就需要特殊处理

1)若\(y\neq 0\),由于方程右边一定是\(0\),而它等于一个非零数,于是这个方程无解

2)若\(y=0\),方程也就是\(0=0\),这说明\(x\)取任意值均成立,而高斯消元最后其实每一个方程对应着一个未知数,所以有几个这样的方程就说明未知数中有多少个自有元,它们无论取何值均满足方程,也就是这个方程有无穷解

几点注意

1)最好设一个\(eps\),据说不设的话可以被卡

2)建议的处理方式是对于方程\(i\)仅保留含有\(x_i\)的那一项,方便最后判断无解和无穷解

#include<iostream>
#include<string.h>
#include<string>
#include<stdio.h>
#include<algorithm>
#include<math.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define lowbit(x) (x)&(-x)
#define fir first
#define sec second
#define rep(i,a,b) for (register int i=a;i<=b;i++)
#define per(i,a,b) for (register int i=a;i>=b;i--)
#define maxd 1000000007
#define eps 1e-9
typedef long long ll;
const int N=100000;
const double pi=acos(-1.0);
int n;
double ans[110],a[110][110];

int read()
{
    int x=0,f=1;char ch=getchar();
    while ((ch<'0') || (ch>'9')) {if (ch=='-') f=-1;ch=getchar();}
    while ((ch>='0') && (ch<='9')) {x=x*10+(ch-'0');ch=getchar();}
    return x*f;
}

double Fabs(double x)
{
    if (x>eps) return x;else return -x;
}

int main()
{
    n=read();
    rep(i,1,n)
    {
        rep(j,1,n+1) scanf("%lf",&a[i][j]);
    }
    rep(i,1,n)
    {
        int now=i;
        rep(j,i+1,n)
            if (Fabs(a[j][i])>Fabs(a[now][i])) now=j;
        if (now!=i) swap(a[now],a[i]);
        if (Fabs(a[i][i])<eps) continue;
        double div=a[i][i];
        rep(j,i,n+1) a[i][j]/=div;
        rep(j,1,n)
        {
            if (i==j) continue;
            div=a[j][i];
            rep(k,i,n+1) a[j][k]-=a[i][k]*div;
        }
    }
    bool nos=0,mors=0;
    rep(i,1,n)
    {
        int cnt=0;
        rep(j,1,n+1) cnt+=(Fabs(a[i][j])<eps);
        if ((cnt==n) && (Fabs(a[i][n+1])>eps)) nos=1;
        else if (cnt==n+1) mors=1;
    }
    if (nos) puts("-1");
    else if (mors) puts("0");
    else
    {
        ans[n]=a[n][n+1];
        per(i,n-1,1)
        {
            ans[i]=a[i][n+1];
            rep(j,i+1,n) ans[i]-=a[i][j]*ans[j];
        }
        rep(i,1,n)
            if (Fabs(ans[i])<=eps) printf("x%d=0\n",i);
            else printf("x%d=%.2lf\n",i,ans[i]);
    }
    return 0;
}

原文地址:https://www.cnblogs.com/encodetalker/p/10807152.html

时间: 2024-11-10 16:24:36

高斯消元小结的相关文章

小结:高斯消元

对于一组多项式方程(增广矩阵中,x[i, n+1]表示式子的值:x[i,j]表示第i个方程第j项的系数,在这里,增广矩阵可能不一定是n个,可能多可能少:opt表示运算规则): (x[1,1]*a[1]) opt (x[1,2]*a[2]) opt ... opt (x[1,n]*a[n])=x[1, n+1] (x[2,1]*a[1]) opt (x[2,2]*a[2]) opt ... opt (x[2,n]*a[n])=x[2, n+1] ... (x[n,1]*a[1]) opt (x[n

高斯消元小小结

所谓高斯消元,就是一种解线性方程组的算法. 学过线性代数的同学都知道,线性方程组本质就是一个向量X1左乘一个系数矩阵A得到另一个向量X2,我们要求解的就是所有未知数构成的向量X1. 设一个n元一次方程组,我们把所有未知数的系数以及等号右边的常数在保持相对位置不变的情况下组成一个n行n+1列的矩阵, a11 a12 a13 ...... a1n b1a21 a22 a23 ...... a2n b2......an1 an2 an3 ...... ann bn 然后我们通过矩阵的初等变换,每个方程

poj_1222_高斯消元

第一次学习使用高斯消元,将灯板化为线性方程组,进行求解. /*######################################################################### # File Name: poj_1222.cpp # Author: CaoLei # Created Time: 2015/7/20 15:48:04 ###################################################################

HDU 4870 Rating(高斯消元)

HDU 4870 Rating 题目链接 题意:一个人注册两个账号,初始rating都是0,他每次拿低分的那个号去打比赛,赢了加50分,输了扣100分,胜率为p,他会打到直到一个号有1000分为止,问比赛场次的期望 思路:f(i, j)表示i >= j,第一个号i分,第二个号j分时候,达到目标的期望,那么可以列出转移为f(i, j) = p f(i', j') + (1 - p) f(i'' + j'') + 1 f(i', j')对应的是赢了加分的状态,f(i'', j'')对应输的扣分的状态

【BZOJ 4171】 4171: Rhl的游戏 (高斯消元)

4171: Rhl的游戏 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 74  Solved: 33[Submit][Status][Discuss] Description RHL最近迷上一个小游戏:Flip it.游戏的规则很简单,在一个N*M的格子上,有一些格子是黑色,有一些是白色 .每选择一个格子按一次,格子以及周围边相邻的格子都会翻转颜色(边相邻指至少与该格子有一条公共边的格子 ),黑变白,白变黑.RHL希望把所有格子都变成白色的.不幸

POJ 1830 开关问题 高斯消元,自由变量个数

http://poj.org/problem?id=1830 如果开关s1操作一次,则会有s1(记住自己也会变).和s1连接的开关都会做一次操作. 那么设矩阵a[i][j]表示按下了开关j,开关i会被操作一次,记得a[i][i] = 1是必须的,因为开关i操作一次,本身肯定会变化一次. 所以有n个开关,就有n条方程, 每个开关的操作次数总和是:a[i][1] + a[i][2] + ... + a[i][n] 那么sum % 2就代表它的状态,需要和(en[i] - be[i] + 2) % 2

BZOJ 3105: [cqoi2013]新Nim游戏 [高斯消元XOR 线性基]

以后我也要用传送门! 题意:一些数,选择一个权值最大的异或和不为0的集合 终于有点明白线性基是什么了...等会再整理 求一个权值最大的线性无关子集 线性无关子集满足拟阵的性质,贪心选择权值最大的,用高斯消元判断是否和已选择的线性相关 每一位记录pivot[i]为i用到的行 枚举要加入的数字的每一个二进制为1的位,如果有pivot[i]那么就异或一下(消元),否则pivot[i]=这个数并退出 如果最后异或成0了就说明线性相关... #include <iostream> #include &l

[bzoj1013][JSOI2008]球形空间产生器sphere-题解[高斯消元]

Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体.现在,你被困在了这个n维球体中,你只知道球面上n+1个点的坐标,你需要以最快的速度确定这个n维球体的球心坐标,以便于摧毁这个球形空间产生器. Input 第一行是一个整数n(1<=N=10).接下来的n+1行,每行有n个实数,表示球面上一点的n维坐标.每一个实数精确到小数点后6位,且其绝对值都不超过20000. Output 有且只有一行,依次给出球心的n维坐标(n个实数),两个实数之间用一个空格隔开.每个实数精确到

[spoj104][Highways] (生成树计数+矩阵树定理+高斯消元)

In some countries building highways takes a lot of time... Maybe that's because there are many possiblities to construct a network of highways and engineers can't make up their minds which one to choose. Suppose we have a list of cities that can be c