Lab_5/8 var (IIl semester)

//Lab_5_var8.cpp
#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <ctime>

using namespace std;
class Matrix
{
private:
    int *M; // указатель на матрицу
    int n, m; // размер матрицы
public:
    Matrix() {}; // конструктор по умолчанию
    Matrix(int, int); //  конструктор с параметрами
    ~Matrix(); // деструктор
    Matrix(Matrix&); // Конструктор копирования
    int get_X(int i, int j) { return *(M + i*m + j); }
    int get_n() { return n; }
    int get_m() { return m; }
    // Дружественные  функции
    // инициализация матрицы случайными числами
    friend void Formation(Matrix&);
    // вывод матрицы на экран
    friend void Print(Matrix&);
    friend Matrix operator -(Matrix&, Matrix&);
    friend Matrix operator /(Matrix&, Matrix&);
    friend Matrix operator <<(Matrix&, Matrix&);
    friend void normN(Matrix &A, int N);
    friend void normM(Matrix &A, int M);
    //Переопределенные операции
    Matrix& operator =(Matrix const& M2)
    {
        int i, j;
        if (this == &M2)return *this;
        else
        {
            n = M2.n;
            m = M2.m;
            M = new int[n*m];
            for (i = 0; i < n; i++) {
                for (j = 0; j < m; j++)
                    *(M + i*m + j) = *(M2.M + i*M2.m + j);
            }
            return *this;
        }
    };
};
Matrix::Matrix(int _n, int _m)
{// конструктор с параметрами
    n = _n;
    m = _m;
    M = new int[n*m];
    cout << "Конструктор с параметрами класса Matrix отработал" << "\n\n";
}
Matrix::~Matrix()
{// деструктор
    delete[] M;
    cout << "Деструктор класса Matrix" << endl;
}
Matrix::Matrix(Matrix &Z)
{// Конструктор копирования
    n = Z.get_n();
    m = Z.get_m();
    M = new int[n*m];
    for (int i = 0; i<n; i++)
        for (int j = 0; j<m; j++)
            *(M + i*m + j) = Z.get_X(i, j);
    cout << "Конструктор копирования класса Matrix\n" << endl;
}
void Formation(Matrix &Z)
{// Инициализация матрицы случайными числами
    for (int i = 0; i<Z.n; i++)
        for (int j = 0; j<Z.m; j++)
            *(Z.M + i*Z.m + j) = 0 + rand() % 20;
}
void Print(Matrix &Z)
{// Вывод матрицы на экран
    for (int i = 0; i<Z.n; i++)
    {
        for (int j = 0; j<Z.m; j++)
            cout << setw(5) << *(Z.M + i*Z.m + j);
        cout << endl;
    }
    cout << endl;
}
Matrix operator << (Matrix &M1, Matrix &M2)
{   /*V = М1 << М2 – из М1 и М2 выбрать строку с наименьшей суммой элементов*/
    float buf;
    // приведение матриц к квадратному виду
    if (M1.m > M1.n) normN(M1, M1.m);
    if (M1.n > M1.m) normM(M1, M1.n);
    if (M2.m > M2.n) normN(M2, M2.m);
    if (M2.n > M2.m) normM(M2, M2.n);
    Matrix Tmp(M1.n, M1.m);
    Tmp.n = M1.n;
    Tmp.m = M1.m;
    for (int i = 0; i < M1.n; i++)
    {
        int sumM1 = 0;
        int sumM2 = 0;
        for (int j = 0; j < M1.m; j++)
        {
            sumM1 += *(M1.M + i*M1.m + j);
            sumM2 += *(M2.M + i*M2.m + j);
        }
        if (sumM1 > sumM2)
        {
            for (int j = 0; j < M1.m; j++)
            {
                *(Tmp.M + i*Tmp.m + j) = *(M2.M + i*M2.m + j);
            }
        }
        else
        {
            for (int j = 0; j < M1.m; j++)
            {
                *(Tmp.M + i*Tmp.m + j) = *(M1.M + i*M1.m + j);
            }
        }
    }
    return Tmp;
}
Matrix operator /(Matrix &M1, Matrix &M2)
{
    //V = М1 / М2 – из минимальных элементов столбцов М1 и М2 составить последовательность V;
    Matrix Tmp(1, M1.m);
    Tmp.n = 1;
    Tmp.m = M1.m;
    int *minM1 = new int[M1.n];
    for (int j = 0; j < M1.n; j++)
        minM1[j] = *(M1.M + 0 * M1.m + j);
    int *minM2 = new int[M2.n];
    for (int j = 0; j < M2.n; j++)
        minM2[j] = *(M2.M + 0 * M2.m + j);
    for (int j = 0; j < M1.m; j++)
        for (int i = 0; i < M1.n; i++)
            if (minM1[j] > *(M1.M + i * M1.m + j))
                minM1[j] = *(M1.M + i * M1.m + j);
    for (int j = 0; j < M2.m; j++)
        for (int i = 0; i < M2.n; i++)
            if (minM2[j] > *(M2.M + i * M2.m + j))
                minM2[j] = *(M2.M + i * M2.m + j);

    for (int j = 0; j < M2.n; j++)
        if (minM1[j] > minM2[j])
            *(Tmp.M + 0*Tmp.m + j) = minM2[j];
        else
            *(Tmp.M + 0 * Tmp.m + j) = minM1[j];
    return Tmp;
}
Matrix operator -(Matrix &M1, Matrix &V)
{
    //М2=М1-V – из М1 удалить строки, сумма элементов которых больше суммы элементов V.
    Matrix Tmp(M1.n, M1.m);
    Tmp.n = M1.n;
    Tmp.m = M1.m;
    for (int i = 0; i < M1.n; i++)
        for (int j = 0; j < M1.m; j++)
            *(Tmp.M + i*Tmp.m + j) = 0;
    int sum = 0;
    for (int i = 0; i < V.m; i++)
        sum += *(V.M + 0 * V.m + i);
    cout << "Sum V =" << sum << endl;
    for (int i = 0; i < M1.n; i++)
    {
        int sumM1Row = 0;
        for (int j = 0; j < M1.m; j++)
            sumM1Row += *(M1.M + i*M1.m + j);
    }
    int countRow = 0;
    for (int i = 0; i < M1.n; i++)
    {
        int sumM1Row = 0;
        for (int j = 0; j < M1.m; j++)
            sumM1Row += *(M1.M + i*M1.m + j);
        cout << "Sum M1["<<i<<"]" << sumM1Row << endl;
        if (sumM1Row < sum)
        {
            for (int j = 0; j < M1.m; j++)
                *(Tmp.M + countRow*Tmp.m + j) = *(M1.M + i*M1.m + j);
            countRow++;
        }
        else
            cout << "Row ["<<i<<"] delete" << endl;
    }
    return Tmp;
}
void normN(Matrix &A, int N)
{// доопределение матрицы по N до максимального размера нулями
    int i, j;
    Matrix Tmp(A);
    A.n = N;
    A.M = new int[A.n*A.m];
    for (i = 0; i < A.n; i++)
    {
        if (i >= Tmp.n)
        {
            for (j = 0; j < A.m; j++)
                *(A.M + i*A.m + j) = 0;
        }
        else {
            for (j = 0; j < A.m; j++)
                *(A.M + i*A.m + j) = *(Tmp.M + i*Tmp.m + j);
        }
    }
}
void normM(Matrix &A, int M)
{// доопределение матрицы по M до максимального размера нулями
    int i, j;
    Matrix Tmp(A);
    A.m = M;
    A.M = new int[A.n*A.m];
    for (j = 0; j < A.m; j++)
    {
        if (j >= Tmp.m)
        {
            for (i = 0; i < A.n; i++)
                *(A.M + i*A.m + j) = 0;
        }
        else {
            for (i = 0; i < A.n; i++)
                *(A.M + i*A.m + j) = *(Tmp.M + i*Tmp.m + j);
        }
    }
};
int main()
{
    setlocale(LC_ALL, "Rus");
    srand(time(NULL));
    int n, m;
    cout << "Введите размерность матрицы A:\n";
    cout << "Введите количество строк: ";
    cin >> n;
    cout << "Введите количество столбцов: ";
    cin >> m;
    cout << "\n";
    Matrix A(n, m);
    Formation(A);
    Print(A);
    cout << "Введите размерность матрицы B:\n";
    cout << "Введите количество строк: ";
    cin >> n;
    cout << "Введите количество столбцов: ";
    cin >> m;
    cout << "\n";
    Matrix B(n, m);
    Formation(B);
    Print(B);
    cout << "Введите размерность матрицы C:\n";
    cout << "Введите количество строк: ";
    cin >> n;
    cout << "Введите количество столбцов: ";
    cin >> m;
    cout << "\n";
    Matrix C(n, m);
    Formation(C);
    Print(C);

    Matrix D, F, G, I;
    D = A<<B;
    Print(D);
    F = A / B;
    Print(F);
    G = A - F;
    Print(G);
    cout << "I = (C - (A << B)) - (B / A)" << endl;
    I = (C - (A << B)) - (B / A);
    Print(I);
    system("pause");
}