Pagini recente » Cod sursa (job #3305239) | Cod sursa (job #1403696) | Cod sursa (job #2832312) | Cod sursa (job #355567) | Cod sursa (job #3343563)
#include <bits/stdc++.h>
using namespace std;
ifstream fin("gauss.in");
ofstream fout("gauss.out");
const int MAX=300;
double EPS=1e-4;
int n,m,i,j,where[MAX+5];
double a[MAX+5][MAX+5],sol[MAX+5],init_a[MAX+5][MAX+5];
bool gauss()
{
for (j=1; j<=m; j++)
{
sol[j]=0;
where[j]=0;
}
for (i=1,j=1; i<=n && j<=m; j++)
{
int row_nonzero=i;
double val=fabs(a[i][j]);
for (int row=i; row<=n; row++)
if (fabs(val)<fabs(a[row][j]))
{
row_nonzero=row;
val=fabs(a[row][j]);
}
if (val<EPS)
continue;
where[j]=i;
for (int col=j; col<=m+1; col++)
swap(a[i],a[row_nonzero]);
for (int row=1; row<=n; row++)
{
if (row!=i)
{
double c=a[row][j]/a[i][j];
for (int col=j; col<=m+1; col++)
a[row][col]-=a[i][col]*c;
}
}
i++;
}
for (j=1; j<=m; j++)
{
if (where[j]!=0)
sol[j]=a[where[j]][m+1]/a[where[j]][j];
}
for (i=1; i<=n; i++)
{
double sum=0;
for (j=1; j<=m; j++)
sum+=sol[j]*init_a[i][j];
sum-=init_a[i][m+1];
if (fabs(sum)>EPS)
return false;
}
return true;
}
int main()
{
ios_base::sync_with_stdio(false);
fin.tie(nullptr); fout.tie(nullptr);
fin>>n>>m;
for (i=1; i<=n; i++)
for (j=1; j<=m+1; j++)
fin>>a[i][j];
for (i=1; i<=n; i++)
for (j=1; j<=m+1; j++)
init_a[i][j]=a[i][j];
if (gauss())
{
for (i=1; i<=m; i++)
fout<<fixed<<setprecision(10)<<sol[i]<<" ";
fout<<'\n';
}
else
fout<<"Imposibil";
return 0;
}