#include <stdio.h>
#define N 100
int sgn (int x);
int read_dim (int *dim);
int matrix_adj (int matrix[N][N], int matrix_aux[N][N], int dim, int m, int n);
int det_matrix_N (int matrix[N][N], int dim);
int det_matrix_adj (int matrix[N][N], int dim, int m, int n);

main() {
	int i, j, dim, det, matrix[N][N];
	while (!read_dim(&dim))
		printf("Dimension no valida (<100).");
	printf("\nIntroduzca elementos matriz...\n");
	for (i=1; i<=dim; i++) {
		printf("Elementos fila %d: ", i);
		for (j=0; j<dim; j++)
			scanf("%d", &matrix[i-1][j]);
	}
	if (dim!=1)
		det=det_matrix_N(matrix,dim);
	else
		det=matrix[0][0];
	printf("\nDeterminante = %d.\n\n", det);
}

int sgn (int x) {
	if (x%2==0) return(1);
	return(-1);
}
int read_dim (int *dim) {
	int nl;
	char aux;
	printf("\nIntroduzca dimension matriz cuadrada (<100): ");
	nl=scanf("%d", &(*dim));
	if (nl!=1) {
		scanf("%c", &aux);
		while (aux!='\n')
			scanf("%c", &aux);
	}
	if ( (*dim>=100) || (*dim<=0) )
		return(0);
	return(1);
}
int matrix_adj (int matrix[N][N], int matrix_aux[N][N], int dim, int m, int n) {
	int i, j, p, q;
	for (j=0, q=0; j<dim; j++) 
		if (j!=n) {
			for (i=0, p=0; i<dim; i++)
				if (i!=m) {
					matrix_aux[p][q]=matrix[i][j];
					p++;
				}
			q++;
		}
	return(dim-1);
}
int det_matrix_N (int matrix[N][N], int dim) {
	int i, j, dim_aux, det;
	for (i=0, j=0, det=0, dim_aux=dim; i<dim; i++)
		det+=sgn(i+j)*matrix[i][j]*det_matrix_adj(matrix,dim_aux,i,j);
	return(det);
}
int det_matrix_adj (int matrix[N][N], int dim, int m, int n) {
	int matrix_aux[N][N];
	if ( (matrix_adj(matrix,matrix_aux,dim,m,n)) == 1 )
		return(matrix_aux[0][0]);
	return(det_matrix_N(matrix_aux,dim-1));
}