Умножение треугольных матриц«Методы обработки разреженных матриц»

#include "stdafx.h"
#include "targetver.h"
#include <stdlib.h>
#include <iostream>
#include <iomanip>
#include <time.h>
#define n 4
using namespace std;
int main()
{ 
        setlocale(LC_ALL,"Russian"); //русская кодировка
        int a[n][n]; //начальная матрица (1)
        int b[n][n]; //начальная матрица (2)
        const int m=n*(n+1)/2;//размерность свёрнутой матрицы
        //int a1[m];   //начальнай свёрнутая матрица (1)
        //int b1[m];   //начальнай свёрнутая матрица (2)
        //int c[n][n]; //произведение (1) и (2)
        //int c1[m];   //произведение свёрнутых (1) и (2)
        clock_t t1;  //переменная таймера
        clock_t t2;
 
        //Матрица 1 
        for(int i=0;i<n;i++)
        {
             for(int j=0;j<n;j++)
             {
                if(j<i)a[i][j]=0;
                else
                {
                    a[i][j]=rand()%10;
                }
             }
        }
        //-------------Свёрнутая матрица 1---------
        int a1[m];   //начальнай свёрнутая матрица (1)
        int k1=0;
        for(int j=0;j<n;j++)
        {
            for(int i=0;i<n;i++)
            {
                if(i<=j)
                {
                  a1[k1]=a[i][j];k1++;
                }
            
            }
        }
 
            //Вывод матрицы 1
        //--------------------------------
        cout<<"Матрица 1\n";
        cout<<"\n---------------------------------\n";
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                cout<<a[i][j]<<"  ";
            }
            cout<<endl;
        }
        cout<<"\n---------------------------------\n";
        //------------------------------------------
        //Вывод свёрнутой матрицы 1
        cout<<"\nСвёрнутая матрица 1\n";
        cout<<"\n---------------------------------\n";
        for(int i=0;i<m;i++)
        {
            cout<<a1[i]<<" ";
        }
            cout<<"\n---------------------------------\n";
            cout<<"\n\n\n";
        //------------------------------------------
            //Матрица 2
        for(int i=0;i<n;i++)
        {
             for(int j=0;j<n;j++)
             {
                if(j<i) b[i][j]=0;
                else
                {
                    b[i][j]=rand()%10;
                }
             }
        }
        //-------------Свёрнутая матрица 2---------
        int b1[m];   //начальнай свёрнутая матрица (2)
        int k2=0;
        for(int j=0;j<n;j++)
        {
            for(int i=0;i<n;i++)
            {
                if(i<=j)
                {
                  b1[k2]=b[i][j];k2++;
                }
            
            }
        }
            //Вывод матрицы 2   
        //-------------------------------------------
        cout<<"\nМатрица 2\n";
        cout<<"\n---------------------------------\n";
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                cout<<b[i][j]<<"  ";
            }
            cout<<endl;
        }
        cout<<"\n---------------------------------\n";
        //-------------------------------------------
        //Вывод свёрнутой матрицы 2
        cout<<"\nСвёрнутая матрица 2\n";
        cout<<"\n---------------------------------\n";
        for(int i=0;i<m;i++)
        {
            cout<<b1[i]<<" ";
        }
            cout<<"\n---------------------------------\n";
            cout<<"\n\n\n";
        //-----------------------------------------
 
 
///---------------------умножение простых матриц----------------------
 
            int c[n][n]; //произведение (1) и (2)
            for(int q = 0;q<n;q++)
            {
                    for(int w = 0; w<n;w++)
                    {
                        c[q][w]=0;
                    }
            }
    
 
   for (int i = 0; i < n; i++) 
    {
        for (int j = 0; j < n; j++)
        {
            // Multiply the row of A by the column of B to get the row, column of product.
            for (int inner = 0; inner < n; inner++)
            {
                c[i][j] += a[i][inner] * b[inner][j];
            }
        }
     }
 
 
 
  //--------------------Вывод результата умножения простых матриц------
 
 //-------------------------------------------
        cout<<"\nРезультат умножения простых матриц\n";
        cout<<"\n---------------------------------\n";
        for(int i = 0; i < n; i++)
                            {
                                for(int j = 0; j < n; j++)
                                    {
                                        cout << setw(3) << c[i][j];
                                    }
 
                                cout << "\n";
                            }
        cout<<"\n---------------------------------\n";
        //-------------------------------------------
 
        ///---------------------умножение свёрнутых матриц----------------------
 
        //-------------------------------------------k=j*(j-1)/2 +i
 
//     inn*(inn-1)/2 +i
 
//     j*(j-1)/2 +inn
int c1[m];   //произведение свёрнутых (1) и (2)
for(int q = 0;q<m;q++){c1[q]=0;}
   
 
  for (int i = 0; i < n; i++) 
   {
        for (int j = 0; j < n; j++) 
        {
            // Multiply the row of A by the column of B to get the row, column of product.
            for (int inn = 0; inn < n; inn++)
            {
                //if(i<=j)
                    c1[j*(j-1)/2 +i] += a1[inn*(inn-1)/2 +i] * b1[j*(j-1)/2 +inn];
            }
 
        }
        
    }
 
    //-------------------------------------------
//--------------------Вывод результата умножения свёрнутых матриц------
 
        cout<<"\nРезультат умножения свёрнутых матриц\n";
        cout<<"\n---------------------------------\n";
        for(int i=0;i<m;i++)
        {
            cout<<c1[i]<<" ";
        }
            cout<<"\n---------------------------------\n";
        //-------------------------------------------
 
        
 
        system("pause");
}