#include <iostream>
#include <string>
#include <fstream>

using namespace std;

int inputValue(int min, int max);
int userInputFromConsole();
double** userInputFromFile(string path);
bool checkPath(string path);
string userInputPath();
int inputMethod();
void printInConsole(double opr);
string userOutputPath();
void printInFile(double opr, string path);
int outputMethod();
void start();
void printTask();

int inputValue(int min, int max) {
    int currentValue;
    bool isNotValid = true;
    do {
        cin >> currentValue;
        if (currentValue >= min && currentValue <= max)
            isNotValid = false;
        else
            cout << "Введите число в заданном диапазоне\n";
    } while (isNotValid);
    return currentValue;
}

int userInputFromConsole() {
    const int MIN_SIZE = 2;
    const int MAX_SIZE = 100;
    int num;
    cout << "Введите порядок матрицы в диапазоне " << MIN_SIZE << ".." << MAX_SIZE << ": ";
    num = inputValue(MIN_SIZE, MAX_SIZE);
    return num;
}

double** userInputMatrixFromConsole(int num) {
    const int MIN_VALUE = -10000;
    const int MAX_VALUE = 10000;
    cout << "Введите элементы матрицы в диапазоне " << MIN_VALUE << ".." << MAX_VALUE << ": " << endl;
    double** matrix = new double*[num];
    for (int i = 0; i < num; i++) {
        matrix[i] = new double[num];
        for (int j = 0; j < num; j++)
            matrix[i][j] = inputValue(MIN_VALUE, MAX_VALUE);
    }
    return matrix;
}

double** userInputFromFile(string path) {
    int num;
    ifstream file(path);
    file.open(path);
    file.clear();
    file >> num;
    double** matrix = new double*[num];
    for (int i = 0; i < num; i++) {
        matrix[i] = new double[num];
        for (int j = 0; j < num; j++)
            file >> matrix[i][j];
    }
    file.close();
    return matrix;
}

bool checkPath(string path) {
    ifstream file(path);
    if (file.is_open()) {
        cout << path << " найден" << endl;
        return true;
    }
    else {
        cout << path << " не найден" << endl;
        return false;
    }
}

string userInputPath() {
    string path;
    bool isNotValid = false;
    do {
        cout << "Введите абсолютный путь к файлу с входными данными" << endl;
        cin >> path;
    } while (!checkPath(path));
    return path;
}

int inputMethod() {
    int method;
    cout << "Каким способом хотите ввести данные?" << endl;
    cout << "1 - с помощью консоли" << endl;
    cout << "2 - с помощью файла" << endl;
    do {
        cin >> method;
        if (method != 1 && method != 2)
            cout << "Введите 1 или 2" << endl;
    } while (method != 2 && method != 1);
    return method;
}

void printInConsole(double opr) {
        cout << opr;
}

string userOutputPath() {
    string path;
    cout << "Введите абсолютный путь к файлу для вывода результата" << endl;
    cin >> path;
    cout << "Результат работа помещён в файл";
    return path;
}

void printInFile(double opr, string path) {
    ofstream file(path);
    file << opr;
}

int outputMethod() {
    int method;
    cout << "Куда хотите вывести результат?" << endl;
    cout << "1 - в консоль" << endl;
    cout << "2 - в файл" << endl;
    do {
        cin >> method;
        if (method != 1 && method != 2)
            cout << "Введите 1 или 2" << endl;
    } while (method != 2 && method != 1);
    return method;
}

double** userInput() {
    int num;
    double** matrix;
    short method = inputMethod();
    if (method == 1) {
        num = userInputFromConsole();
        matrix = userInputMatrixFromConsole(num);
    }
    else {
        string path = userInputPath();
        matrix = userInputFromFile(path);
    }
    return matrix;
}

void printResult(double opr) {
    short method = outputMethod();
    if (method == 1)
        printInConsole(opr);
    else {
        string path = userOutputPath();
        printInFile(opr, path);
    }
}

int calcZnak(int colPer) {
    if (colPer % 2 == 0)
        return 1;
    else
        return -1;
}

void per(int k, int n, double** matrix, int colPer) {
    double z = abs(matrix[k][k]);
    int i = k;
    colPer = 0;
    for (int j = k + 1; j < n; j++)
        if (abs(matrix[j][k]) > z) {
            z = abs(matrix[j][k]);
            i = j;
        }
    if (i > k) {
        colPer++;
        for (int j = k; j < n; j++) {
            z = matrix[i][j];
            matrix[i][j] = matrix[k][j];
            matrix[k][j] = z;
        }
    }
}

double calcOpr(int size, double** matrix, double det) {
    det = 1;
    int colPer = 0;
    for (int i = 0; i < size; i++) {
        if (matrix[i][i] == 0)
            per(i, size, matrix, colPer);
        det = calcZnak(colPer) * det * matrix[i][i];
        for (int j = i + 1; j < size; j++) {
            double r = matrix[j][i] / matrix[i][i];
            for (int k = i; k < size; k++)
                matrix[j][k] = matrix[j][k] - r * matrix[i][k];
        }
    }
    return det;
}

void printTask() {
    cout << "Данная программа вычисляет определитель матрицы" << endl;
}

void start() {
    printTask();
    double** matrix = userInput();
    double opr = calcOpr(_msize(matrix) / sizeof(matrix[0]), matrix, 0);
    printResult(opr);
}

int main()
{
    setlocale(LC_ALL, "Russian");
    start();
}