#include <gmpxx.h>
#include <iostream.h>
#include <stdlib.h>

// **********************************************************************
//                  GLOBAL VARIABLES
// **********************************************************************
int i,j,k;  // global only for debbuging purpose

// **********************************************************************
//                  FUNCTION PROTOTYPE
// **********************************************************************
mpz_class deter(mpz_class *t,int n);

// **********************************************************************
//                 FUNCTION DEFINITION
// **********************************************************************
mpz_class deter(mpz_class *t,int n)
{
  // int i;
  static mpz_class D,a,b;
  mpz_class *LAUX,*L;

  D=1;
  for(i=1;i<n;i++)
    {
      for (j=0;j<i;j++)
	{
	  //  int k;
	  cout <<"// BEGIN i="<<i<<" j="<<j<<"  ";

	  a = t[i*n+j];
	  if (a==0)
	    continue;
	  b = t[j*n+j];
	  if ( b == 0)
	    {
	      cout << "# Pivot nul: can't continue\n";
	      exit(0);
	    }
	  // elinine a[i,j]: L(j) -> t(i,i)*L(j)-t(i,j)*L(i)
	  // elinine a[i,j]: L(j) -> b*L(j)-a*L(i)
	  
	  if (i<n-1)
	    D *= b;
	  L = t+i*n;
	  LAUX = t+j*n;
	  
	  for (k=j;k<n;k++)
	    L[k] = b*L[k]-a*LAUX[k];
	  cout <<" END i="<<i<<" j="<<j<<endl;
	}
    }

  if ( D == 0)
    cout << "# D nul!!\n";
  D =  t[n*n-1]/D;
  return D;
}                  // FIN deter()
// ***********************************************************************

int main(int na,char **para)
{
#define _Q 12
  mpz_class De;
  mpz_class *M;

  cerr << "# version gmp"<<gmp_version <<endl;

  M = new (mpz_class)[_Q*_Q];
  if (M==0)
    exit(0);

#include "inc.cpp" // lines like M[0] = "5767898980909...7667"; but HUGES !!
  De = deter(M,_Q);

  cout << De << endl;
  return 0;
}   //FIN main()
// **********************************************************************