高斯消元 hihocoder 1195

参考博客：https://blog.csdn.net/pengwill97/article/details/77200372

 

做法：裸的高斯消元，自己敲模板是没过的，（答案都为整数，看了博客，发现居然要用double来写，而且我试了几组数据，验证是不对的，博客上的代码也不对。。。。。。）

 

///#include<bits/stdc++.h>
///#include<unordered_map>
///#include<unordered_set>
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<queue>
#include<set>
#include<stack>
#include<map>
#include<new>
#include<vector>
#define MT(a,b) memset(a,b,sizeof(a));
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const double pai=acos(-1.0);
const double E=2.718281828459;
const int mod=1e9+7;
const int INF=0x3f3f3f3f;
const double eps=1e-7;

int n,m;
double matrix[1005][1005];    ///增广矩阵
double x[1005];              ///解集
bool free_x[1005];       ///是否为不确定变元

int gauss(int equ,int var)/// equ 方程数量 ,var 变元数量
{
    ll max_row;         ///当前这列绝对值最大的行
    ll now_col=1;         ///当前列
    ll k=1;               ///当前行
    for(int i=1; i<=equ; i++)
    {
        x[i]=0;
        free_x[i]=1;
    }
    for(; k<=equ&&now_col<=var; k++,now_col++)
    {
        max_row=k;
        for(int i=k+1; i<=equ; i++)
        {
            if(fabs(matrix[i][now_col])>fabs(matrix[max_row][now_col]))
                max_row=i;
        }
        if(fabs(matrix[max_row][now_col])<eps)
        {
            k--;
            continue;
        }
        if(max_row!=k)
        {
            for(int i=1; i<=var+1; i++)
                swap(matrix[k][i],matrix[max_row][i]);
        }
        for(int i=k+1; i<=equ; i++)
        {
            if(fabs(matrix[i][now_col])>eps)
            {
                double s=matrix[i][now_col]/matrix[k][now_col];
                for(int j=now_col; j<=var+1; j++)
                    matrix[i][j]-=matrix[k][j]*s;
                matrix[i][now_col]=0;
            }
        }
    }
    ///1.无解的情况: 化简的增广阵中存在(0, 0, ..., a)这样的行(a != 0).
    for(int i=k; i<=equ; i++)
    {
        if(fabs(matrix[i][now_col])>eps)
            return -1;
    }
    ///2.无穷解的情况: 在var * (var + 1)的增广阵中出现(0, 0, ..., 0)这样的行，即说明没有形成严格的上三角阵.
    if(k<var+1)
        return var+1-k;
    /// 3. 唯一解的情况: 在var * (var + 1)的增广阵中形成严格的上三角阵.
    /// 计算出Xn, Xn-1 ... X1.
    for(int i=var; i>=1; i--)
    {
        double temp=matrix[i][var+1];
        for(int j=i+1; j<=var; j++)
            temp=temp-matrix[i][j]*x[j];
        ///if (temp % a[i][i] != 0) return -2; /// 说明有浮点数解，但无整数解.
        x[i]=temp/matrix[i][i];
    }
    return 0;
}
int main()
{
    while(~scanf("%d %d",&n,&m))
    {
        for(int i=1; i<=m; i++)
            for(int j=1; j<=n+1; j++)
                scanf("%lf",&matrix[i][j]);
        int flag=gauss(m,n);
        if(flag==-1)
            printf("No solutions\n");
        else if(flag>0)
            printf("Many solutions\n");
        else if(flag==-2)
            printf("No integer solution.\n");
        else
        {
            for(int i=1; i<=n; i++)
                printf("%d\n",(int)(x[i]+0.5));
        }
    }
    return 0;
}