#include <iostream>
#include <time.h>


using namespace std;
const  GYOU =10;
const  RETU =GYOU+1;
//行列表示関数
void show(double *pA){
   for(int i=0;i<GYOU;i++){
      for(int j=0;j<RETU;j++){
         cout<<*(pA+(i*RETU)+j)<<" ";
      }
      cout<<"\n";
   }
cout<<"\n";
}
//行列表示関数

int main() {
double *pA;
double X[GYOU];//ｘの解を入れる

int flag=0; //A[p][p]が0の時のフラグ

double Atemp[GYOU][RETU];//一時行列置き場
double Temp=0;//一時数値置き場
double Wrareta=0;//-(Atemp[m+1][p]/Atemp[m][p])

srand((unsigned)time(NULL));//ランダム元


//################################
//対象行列
   //対象行列に数値を入れる

 static double A[GYOU][RETU];
   for(int i=0;i<GYOU;i++){
X[i]=0;
      for(int j=0;j<RETU;j++){
         A[i][j]=rand();
      Atemp[i][j]=0; 
      }

   }


/**
//2x*6y*4z=8
//3x*7y*2z=6
//1x*2y*3z=7
	static double A[GYOU][RETU]={ //constとあわせること( ..)φ注意
   {2,6,4,8 },{3,7,2,6},{1,2,3,7}
};

**/


pA=&A[0][0];
   show(pA);//行列表示

for( int p=0;p<GYOU-1;p++){//ppは 行列上の対角線を動く

   /*
      A[p][p]から[p+1][p]〜[p+n][p]までバブルソート
      して 行を入れ替える nは 端っこまでってこどだじぇ
   */
      for(int s=p;s<GYOU-1;s++) {
         for(int t=s+1;t<GYOU;t++){
            if(A[t][p]>A[s][p]) {
               for(int r=0;r<RETU;r++){
               
               Temp=A[t][r];
               A[t][r]=A[s][r];
               A[s][r]=Temp;

               }
            } 
         }
   
      } 
      
      if(A[p][p]==0){flag=1;break;}//[p][p]のチェック
      show(pA);
         
      //Atempに一時保存
         for( int i=0;i<GYOU;i++){
            for(int j=0;j<RETU;j++){
            Atemp[i][j]=A[i][j];
            }
         } 


      //A[p+1][p]〜A[p+n][p]までを 0にするための処理
      for(int m=p;m<GYOU-1;m++){
         Wrareta=- (Atemp[m+1][p]/Atemp[m][p]);
            
         for(int z=p;z<RETU;z++){ 
               
            A[m+1][z]=Atemp[m+1][z]+Atemp[m][z]*Wrareta;
         }
      }
      show(pA); 
}
if(flag==1) {cout<<"解は不定もしくは不能\n";}//不能、不定の表示

for(int m=GYOU-1; m>=0; --m) 
{
for(int n=m; n<GYOU-1; ++n) 
{
         A[m][GYOU]-=A[m][n+1]*X[n+1];
}
X[m]=A[m][GYOU]/(A[m][m]); 
}
   for (int t=0;t<GYOU;t++){ 
cout<<"X"<<t+1<<"="<<X[t]<<"\n";
   }
double check=0;
    for(int v=0;v<GYOU;v++){
      check+=X[v]*(A[0][v]);
   }
   cout<<Atemp[0][RETU-1]<<"=="<<check<<"\n";
return 0;
}