Линейное программирование. Транспортная задача. Метод минимального элемента

struct Pair {
    int i;
    int j;
};
int main()
{
    setlocale(LC_ALL,"RUS");
 
    int* a = nullptr; //Производство
    int m = 0;  //Кол-во машин
 
    int* b = nullptr; //потребность
    int n = 0;  //Кол-во работ
 
    int** c = nullptr; //стоимость
    int** x = nullptr; //решение
 
    Pair* found_elements = new Pair[m * n];
    int num_elements = 0; // Количество найденных элементов
 
    data_from_file(a, b, c, m, n, file_name);
    min_element_method(a, b, c, x, m, n, found_elements, num_elements);
    system("pause");
    return 0;
}

void data_from_file(int*& a, int*& b,int**& c, int& m, int& n, string name)
{
    ifstream file(name);
 
    // Считываем количество машин
    file.ignore(256, ':');
    file >> m;
 
    // Считываем количество работ
    file.ignore(256, ':');
    file >> n;
 
    // Выделяем память для массивов a и b
    a = new int[m];
    b = new int[n];
 
    //считывание параметров из файла
    file.ignore(256, ':');
    for (int i = 0; i < m; i++) {
        file >> a[i];
    }
 
    file.ignore(256, ':');
    for (int i = 0; i < n; i++) {
        file >> b[i];
    }
 
    c = new int*[m];
    for (int i = 0; i < m; i++) {
        c[i] = new int[n];
    }
 
    // Считываем стоимость перевозок
    file.ignore(256, ':');
    for (int i = 0; i < m; i++) {
        for (int j = 0; j < n; j++) {
            file >> c[i][j];
        }
    }
}

void print_matrix_with_marks(int** x, int m, int n, Pair* found_elements, int num_elements) {
    for (int i = 0; i < m; i++) {
        for (int j = 0; j < n; j++) {
            bool is_found = false;
            for (int k = 0; k < num_elements; k++) {
                if (found_elements[k].i == i && found_elements[k].j == j) {
                    is_found = true;
                    break;
                }
            }
            
            if (is_found) {
                cout << x[i][j]<<"* ";
            } else {
                cout << x[i][j]<< " ";
            }
        }
        cout << endl;
    }
}

void min_element_method(int*& a, int*& b, int**& c, int**& x, int& m, int& n, Pair*& found_elements, int& num_elements)
{
    int count_freeze_rows = 0;
    int count_freeze_coloms = 0;
 
    int z = 0;
    x = new int*[m];
    for (int i = 0; i < m; i++) {
        x[i] = new int[n];
    }
 
    num_elements = 0; // Инициализируем количество найденных элементов
 
    for (int i = 0; i < m; i++) {
        for (int j = 0; j < n; j++) {
            x[i][j] = 0;
        }
    }
 
    // Создаем массивы для "замороженных" строк и столбцов
    bool* frozen_rows = new bool[m];
    bool* frozen_cols = new bool[n];
    for (int i = 0; i < m; i++) {
        frozen_rows[i] = false;
    }
    for (int j = 0; j < n; j++) {
        frozen_cols[j] = false;
    }
 
    while (num_elements < m + n - 1) {
        int min_cost = INT_MAX;
        int min_i = -1;
        int min_j = -1;
        // Поиск минимального элемента
        for (int i = 0; i < m; i++) {
            if (!frozen_rows[i]) 
            { // Проверяем, не "заморожена" ли строка
                for (int j = 0; j < n; j++) 
                {
                    if (!frozen_cols[j]) { // Проверяем, не "заморожен" ли столбец
                        if (a[i] > 0 && b[j] > 0 && x[i][j] == 0 && c[i][j] < min_cost) 
                        {
                            min_cost = c[i][j];
                            min_i = i;
                            min_j = j;
                        }
                    }
                }
            }
        }
 
        if (min_i == -1 || min_j == -1) {
            break; // Нет больше минимальных элементов
        }
 
        // Выбираем минимум
        int work = min(a[min_i], b[min_j]);
        x[min_i][min_j] = work;
        a[min_i] -= work;
        b[min_j] -= work;
 
        z += x[min_i][min_j] * c[min_i][min_j];
 
        // Сохраняем найденный элемент в массиве found_elements
        found_elements[num_elements].i = min_i;
        found_elements[num_elements].j = min_j;
        num_elements++;
 
         
        // Замораживаем строку и столбец
            frozen_rows[min_i] = true;
            frozen_cols[min_j] = true;
        
    }
    //cout<< count_freeze_rows;
 
    cout << "Количество элементов в базисе: " << num_elements << endl;
 
    cout << "Целевая функция равна: " << z << endl;
 
    cout << "Полученный допустимый базис: {";
    for (int k = 0; k < num_elements; k++) {
        cout << "(" << found_elements[k].i + 1 << "," << found_elements[k].j + 1 << ")";
        if (k < num_elements - 1) {
            cout << ",";
        }
    }
    cout << "}" << endl;
 
    cout << "Базисное допустимое решение T-задачи: " << endl;
    print_matrix_with_marks(x, m, n, found_elements, num_elements);
}