#include <stdio.h>
#include <cstdlib>
#include <ctime>
#include <math.h>
#include <iostream>
#include <omp.h>

using namespace std;

#define MATRIX_SIZE 200
#define MAX_RANDOM 100

int main()
{
	srand((unsigned)time(nullptr));

	float matrixA[MATRIX_SIZE][MATRIX_SIZE];
	float matrixB[MATRIX_SIZE][MATRIX_SIZE];
	float matrixC[MATRIX_SIZE][MATRIX_SIZE];
	float matrixD[MATRIX_SIZE][MATRIX_SIZE];
	float matrixAB[MATRIX_SIZE][MATRIX_SIZE];
	float matrixCD[MATRIX_SIZE][MATRIX_SIZE];

	for (int i = 0; i < MATRIX_SIZE; i++) {
		for (int j = 0; j < MATRIX_SIZE; j++) {
			matrixA[i][j] = (rand() % MAX_RANDOM) + 1;
		}
	}

	for (int i = 0; i < MATRIX_SIZE; i++) {
		for (int j = 0; j < MATRIX_SIZE; j++) {
			matrixB[i][j] = (rand() % MAX_RANDOM) + 1;
		}
	}


	int minB = matrixB[0][0];
	int maxA = matrixA[0][0];

	for (int i = 0; i < MATRIX_SIZE; i++) {
		for (int j = 0; j < MATRIX_SIZE; j++) {
			matrixC[i][j] = 1 / matrixA[i][j];
		}
	}
	for (int i = 0; i < MATRIX_SIZE; i++) {
		for (int j = 0; j < MATRIX_SIZE; j++) {
			matrixD[i][j] = sqrt(matrixB[i][j]);
		}
	}
	for (int i = 0; i < MATRIX_SIZE; i++) {
		for (int j = 0; j < MATRIX_SIZE; j++) {
			matrixAB[i][j] = matrixA[i][j] + matrixB[i][j];
		}
	}
	for (int i = 0; i < MATRIX_SIZE; i++) {
		for (int j = 0; j < MATRIX_SIZE; j++) {
			matrixCD[i][j] = matrixC[i][j] * matrixD[i][j];
		}
	}
	for (int i = 0; i < MATRIX_SIZE; i++) {
		for (int j = 0; j < MATRIX_SIZE; j++) {
			if (matrixA[i][j] > maxA) maxA = matrixA[i][j];
		}
	}
	for (int i = 0; i < MATRIX_SIZE; i++) {
		for (int j = 0; j < MATRIX_SIZE; j++) {
			if (matrixB[i][j] < minB) minB = matrixB[i][j];
		}
	}
	

	
	
	return 0;
}