#include <stdio.h>  // printf/scanf
#include <stdlib.h>  // malloc


/*! \brief Получение матрицы
 *
 *  \details Получение от пользователя матрицы размерности length
 *
 *  \param [out] matrix    Указатель для передачи матрицы
 *  \param length    	   Ожидаемая размерность матрицы
 */
void getMatrix(float** matrix, int length);


/*! \brief Получить алгебраическое дополнение
 *
 *  \details Возвращает алгебраическое дополнение элемента matrix[x][y]
 *
 *  \param  matrix    	Матрица
 *  \param length    	Размерность матрицы
 */
float ** getAlgemAddition(float** matrix, int length, int x, int y);


void printMatrix(float** matrix, int length){
	int i = 0;
	int j = 0;
	for(i = 0; i < length; i++){
		for(j = 0; j < length; j ++){
			printf("%f ", matrix[i][j]);
		}
		printf("\n");
	}
}


/*! \brief Функция Main
 *  \return Код завершения программы (0-OK)
 */
int main(){
	int length = 0;

	while(1){
		printf("\nВведите размерность матрицы: ");

		if(scanf("%d",&length) != 1){
			printf("\nНекорректный ввод\n");

			scanf("%*[^\n]");  // clearing STD from invalid input
			length = 0;  // just to be sure
		}
		else if(length > 0){
			break;
		}
		else{
			printf("\nНеверная размерность\n");
		}
	}

	printf("\n%d\n", length );

	float** matrix;

	getMatrix(matrix, length);

	printMatrix(matrix, length);

	return 0;
}


void getMatrix(float** p_matrix, int length){
	int i = 0;
	int j = 0;

	p_matrix = malloc( sizeof(float*) * length );

	for(i = 0; i < length; i++){
		p_matrix[i] = malloc( sizeof(float) * length );
		for(j = 0; j < length; j++){

			float elem;

			while(1){
				printf("\nMatrix[%d][%d]: ",i ,j);

				if(scanf("%f",&elem) != 1){
					printf("\nНекорректный ввод\n");
					scanf("%*[^\n]");  // clearing STD from invalid input
				}
				else{
					break;
				}
			}

			p_matrix[i][j] = elem;

		}
	}
}

float** getAlgemAddition(float** p_matrix, int length, int x, int y){
	// invalid arguments
	if(x >= length || x < 0 || y >= length || y < 0){
		return 0;
	}


	int i, j, ci, cj, di, dj;
					// 'c' stands for 'copy' (index of copy)
					// 'd' stands for 'delta'
					// (delta between copy index and original index)

	float** p_matrix_copy;
	p_matrix_copy = malloc( sizeof(float*) * length-1 ); // length-1
											// (w/o one line and row)
	di = 0;
	for(i = 0; i < length-1; i++){
		if(ci == i) di = 1;  // offset (skiping this line)

		p_matrix[i] = malloc( sizeof(float) * length-1 );
		dj = 0;
		for(j = 0; j < length-1; j++){
			if(cj == j) dj = 1;  // offset (skiping this row)
			p_matrix_copy[i][j] = p_matrix[ci + di][cj + dj];
 		}
	}
	return p_matrix_copy;
}