fork download
copy 
  1. #include <iostream>
  2. using namespace std;
  3.  
  4. template <typename T>
  5. class Matrix {
  6. protected:
  7.     T **a;
  8.     int size1, size2;
  9.  
  10. public:
  11.  
  12.     Matrix(int n, int m) {
  13.         if (n <= 0 || m <= 0) {
  14.             cout << "Помилка: негативні або нульові розміри матриці!" << endl;
  15.             size1 = size2 = 1; 
  16.         } else {
  17.             size1 = n;
  18.             size2 = m;
  19.         }
  20.  
  21.         a = new T*[size1];
  22.         for (int i = 0; i < size1; i++) {
  23.             a[i] = new T[size2];
  24.         }
  25.     }
  26.  
  27.     Matrix() : size1(1), size2(1) {
  28.         a = new T*[1];
  29.         a[0] = new T[1];
  30.     }
  31.  
  32.     ~Matrix() {
  33.         for (int i = 0; i < size1; i++) {
  34.             delete[] a[i];
  35.         }
  36.         delete[] a;
  37.     }
  38.  
  39.     T* operator [](const int i) {
  40.         return a[i];
  41.     }
  42.  
  43.     Matrix operator+(const Matrix& x) const {
  44.         if (size1 != x.size1 || size2 != x.size2) {
  45.             cout << "Розміри матриць не збігаються для додавання" << endl;
  46.         }
  47.         Matrix m3(size1, size2);
  48.         for (int i = 0; i < size1; i++) {
  49.             for (int j = 0; j < size2; j++) {
  50.                 m3.a[i][j] = a[i][j] + x.a[i][j];
  51.             }
  52.         }
  53.         return m3;
  54.     }
  55.  
  56.     Matrix operator-(const Matrix& x) const {
  57.         if (size1 != x.size1 || size2 != x.size2) {
  58.             cout  << "Розміри матриць не збігаються для віднімання"<< endl;
  59.         }
  60.         Matrix m3(size1, size2);
  61.         for (int i = 0; i < size1; i++) {
  62.             for (int j = 0; j < size2; j++) {
  63.                 m3.a[i][j] = a[i][j] - x.a[i][j];
  64.             }
  65.         }
  66.         return m3;
  67.     }
  68.  
  69.     Matrix operator*(const Matrix& x) const {
  70.         if (size2 != x.size1) {
  71.             cout << "Неможливо множити матриці, розміри не збігаються"<< endl;
  72.         }
  73.         Matrix c(size1, x.size2);
  74.         for (int i = 0; i < size1; i++) {
  75.             for (int j = 0; j < x.size2; j++) {
  76.                 c.a[i][j] = 0;
  77.                 for (int k = 0; k < size2; k++) {
  78.                     c.a[i][j] += a[i][k] * x.a[k][j];
  79.                 }
  80.             }
  81.         }
  82.         return c;
  83.     }
  84.  
  85.     friend ostream& operator<<(ostream& os, const Matrix& m) {
  86.         for (int i = 0; i < m.size1; i++) {
  87.             for (int j = 0; j < m.size2; j++) {
  88.                 os << m.a[i][j] << " ";
  89.             }
  90.             os << endl;
  91.         }
  92.         return os;
  93.     }
  94.  friend istream& operator>>(istream& is, Matrix& m) {
  95.         for (int i = 0; i < m.size1; i++) {
  96.             for (int j = 0; j < m.size2; j++) {
  97.                 is >> m.a[i][j];
  98.             }
  99.         }
  100.         return is;
  101.     }
  102.  
  103.     int getSize1() const { return size1; }
  104.     int getSize2() const { return size2; }
  105. };
  106.  
  107. template <typename T>
  108. class SquareMatrix : public Matrix<T> {
  109. public:
  110.     SquareMatrix(int n) : Matrix<T>(n, n) {}
  111.  
  112.     T calculateTrace() const {
  113.         if (this->getSize1() != this->getSize2()) {
  114.             cout << "Матриця не є квадратною" << endl;
  115.         }
  116.         T trace = 0;
  117.         for (int i = 0; i < this->getSize1(); i++) {
  118.             trace += this->a[i][i];
  119.         }
  120.         return trace;
  121.     }
  122. };
  123.  
  124. int main() {
  125.         Matrix<double> a(2, 2);
  126.         cin >> a;
  127.         Matrix<double> b(2, 2);
  128.         cin >> b;
  129.         Matrix<double> c = a + b;
  130.         cout << c;
  131.         Matrix<double> d = a - b;
  132.         cout << d;
  133.         Matrix<double> e = a * b;
  134.         cout << e;
  135.         SquareMatrix<double> sqMatrix(3);
  136.         cin >> sqMatrix;
  137.         cout << "Слід матриці: " << sqMatrix.calculateTrace() << endl;
  138.     return 0;
  139. }
  140.  
Success #stdin #stdout 0.01s 5288KB
comments (?)
stdin
copy 
1 2
1 2
4 5
5 6
4 5 6
6 7 8
2 3 4
stdout
copy 
5 7 
6 8 
-3 -3 
-4 -4 
14 17 
14 17 
Слід матриці: 15
https://ideone.com/Hqq0pa
language:
C++ (gcc 8.3)
created:
17 hours ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!