高斯消元(求解double,int型方程组)

#include<iostream>
#include<algorithm>
#include<vector>
#include<map>
#include<cstdio>
#include<cmath>
using namespace std;
const int maxn=507;
const double eps=1e-7;
double A[2*maxn][maxn];  //增广矩阵
double x[maxn];  //解集
//equ: 方程个数  var:变元个数
int Gauss(int equ,int var) {  //-1:无解   0:有唯一解  大于0:有多解,返回自由变元个数
    int row,col;
    for(row=0,col=0; col<var&&row<equ; col++,row++) {
        int max_r=row;
        for(int i=row+1; i<equ; i++) {
            if(fabs(A[i][col])>fabs(A[max_r][col])) {
                max_r=i;
            }
        }
        if(max_r!=row) {
            for(int j=col; j<var+1; j++)
                swap(A[row][j],A[max_r][j]);
        }
        if(fabs(A[row][col])<eps) {
            row--;
            continue;
        }
        for(int i=row+1; i<equ; i++) {
            if(fabs(A[i][col])<eps)continue;
            for(int j=col+1; j<var+1; j++) {
                A[i][j]=(A[i][j]*A[row][col]-A[i][col]*A[row][j])/A[row][col];
            }
            A[i][col]=0;
        }
    }
    for(int i=row; i<equ; i++) {
        if(fabs(A[i][col])>eps)return -1;
    }
    if(row<var)return var-row;
    for(int i=var-1; i>=0; i--) { //计算唯一解。
        double tmp=0;
        for(int j=i+1; j<var; j++) {
            tmp+=A[i][j]*x[j];
        }
        x[i]=(A[i][var]-tmp)/A[i][i];
    }
    return 0;
}
int main() {
//   freopen("in.txt","r",stdin);
    int equ,var;
    while(cin>>var>>equ) {
        for(int i=0; i<equ; i++) {
            for(int j=0; j<=var; j++) {
                cin>>A[i][j];
            }
        }
        int ans=Gauss(equ,var);
        if(ans==-1)
            cout<<"No solutions"<<endl;
        else if(ans>0)
            cout<<"Many solutions"<<endl;
        else {
            for(int i=0; i<var; i++)
                cout<<(int)(x[i]+0.5)<<endl;  //输出整数解记得四舍五入
        }
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值