2022春每日一题:Day 12
题目:[SDOI2006]线性方程组
显然,高斯消元模板题
代码:
#include
#include
#include
#include
#include
#define N 105
#define eps 1e-6
using namespace std;
int n;
double a[N][N];
int guass()
{
int r=0;
for(int c=0;cfabs(a[t][c]))
t=i;
if(fabs(a[t][c])=c;i--)
a[r][i]/=a[r][c];
for(int i=r+1;ieps)
for(int j=n;j>=c;j--)
a[i][j]-=a[r][j]*a[i][c];
++r;
}
if(reps)
return 2;
return 1;
}
for(int i=n-1;i>=0;i--)
for(int j=i+1;j