Lab_5/6 var (IIl semester)

//Lab_5.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 int operator ~ (Matrix&);
    friend Matrix operator -(Matrix&, Matrix&);
    friend Matrix operator >> (Matrix&, int k);
    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) = 1 + rand() % 100;
}

void Print(Matrix &Z)
{
    // Вывод матрицы на экран
    for (int i = 0; i<Z.n; i++)
    {
        for (int j = 0; j<Z.m; j++)
            cout << setw(4) << *(Z.M + i*Z.m + j);
        cout << endl;
    }
    cout << endl;
}

Matrix operator - (Matrix &M1, Matrix &M2)
{
    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 M3(M1.n, M1.m);
    for (int i = 0; i < M1.n; i++)
        for (int j = 0; j < M1.m; j++)
            *(M3.M + i*M3.m + j) = *(M1.M + i*M1.m + j) - *(M2.M + i *M2.m + j);
    return M3;
}

int operator ~ (Matrix &M3)
{
    int nstrmin = 0;
    int strsum = 0;
    int min = 0;
    for (int i = 0; i < M3.n; i++)
    {
        for (int j = 0; j < M3.m; j++)
            strsum = strsum + *(M3.M + i*M3.m + i);
        if (min > strsum)
        {

            min = strsum;
            nstrmin = i;
        }
    }

    Matrix Tmp(1, M3.m);
    for (int j = 0; j < Tmp.m; j++)
        *(Tmp.M + nstrmin *Tmp.m + j) = *(M3.M + nstrmin *M3.m + j);
    return nstrmin;
}



Matrix operator >> (Matrix &M1, int k)
{ //М2 = М1 >> К – циклический сдвиг строк матриц М1 и М2 вниз на К позиций;
    Matrix Tmp(M1.n + k, M1.m);
    Tmp.n = M1.n; //i
    Tmp.m = M1.m; //j
    int count = 0;
    cout << "Сдвиг на " << k << endl;
    for (int i = 0; i < Tmp.n; i++)
        for (int j = 0; j < Tmp.m; j++)
            *(Tmp.M + i*Tmp.m + j) = count++;
    for (int i = 0; i < M1.m; i++)
    {
        for (int j = 0; j < M1.n; j++)
        {
            int plusI = i + k;
            *(Tmp.M+ plusI* Tmp.m + j ) = *(M1.M + i* M1.m + j);
        }
        cout << 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, k = 0;
    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 I;
    I = B >> (~(A-C));
    Print(I);
    system("pause");
    return 0;
}