#include<stdio.h>
#include<math.h>
#include<iostream>
#include<math.h>
#include<iomanip>
#include<iostream>
#include<math.h> using namespace std;
int main()
{
int n,i,j,l,k,t,s,k1,k2;
oat p;
cout<<"请?输 º?入¨?方¤?程¨¬个?数 ºyn"<<endl;
cin>>n;
oat A[n][n+1],X[n],Y[n];
cout<<"将?方¤?程¨¬组 Á¨¦用®?增?广?矩?阵¨®A 表 À¨ª 示 º?,ê?请?按 ã¡ä 行 D 依°¨¤次
ä?输 º?入¨?矩?阵¨®A"<<endl;
for(i=0;i<n;i++)
for(j=0;j<=n;j++)
cin>>A[i][j];
for(i=0;i<n;i++)Y[i]=i;
for(i=0;i<n-1;i++)
{
k1=i;
k2=i;
for(l=i;l<n;l++)
{
for(j=i;j<n;j++)
if(fabs(A[l][j])>fabs(A[k1][k2]))k1=l,k2=j;
}
//
if(A[k][i]==0)
break;
if(k1!=i)
for(j=0;j<=n;j++)
{
p=A[i][j];
A[i][j]=A[k1][j];
A[k1][j]=p;
}
if(k2!=i)
{
p=Y[i];
Y[i]=Y[k2];
Y[k2]=p;