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

@8VAid8 · Регистрация: 21.04.2012

Author24 — интернет-сервис помощи студентам

Нужно перемножить треугольные матрицы в обычном виде и в свёрнутом. С обычным проблем нет.
Доступ к элементам свёрнутой матрицы осуществляется по формуле : k= j*(j-1)/2+i. Подстановкой этой формулы в алгоритм умножения обычных матриц пытаюсь добится результата, но он не совпадает. Возможно потому что обычные матрицы также ниже главной диагонали заполнены нулями...
Каким еще способом можно перемножить свёрнутые разряженные матрицы? Или где у меня ошибки?

C++

#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");
}

Вернуться к обсуждению:
Умножение треугольных матриц«Методы обработки разреженных матриц» C++

Следующий ответ

Другие темы раздела
C++ Вычислить количество элементов в массиве, равных нулю https://www.cyberforum.ru/ cpp-beginners/ thread709379.html Помогите пожалусто через 2 часа в институт пилить а мне зачет сделать нужно не в службу а в дружбу выручите	Ошибка при вызове конструктора с параметрами C++ подскажите, почему при вызове конструктора с параметрами выдает ошибку сегментации #include<math.h> #include<stdio.h> #include<stdlib.h> #include<conio.h> #include<string.h> #include <iostream> using namespace std;
C++ Новый масив из уникальных элементов другого массива Здравствуйте, помогите сделать функцию которая генерирует новый массив из элементов, которые не повторяются в другом массиву. Заранее спасибо. https://www.cyberforum.ru/ cpp-beginners/ thread709358.html	C++ Все элементы вектора, которые встречаются более одного раза, переписать в другой вектор https://www.cyberforum.ru/ cpp-beginners/ thread709355.html Все элементы вектора V,которые встречаються более одного раза ,переписать в другой вектор. Задача на С++ помогите решить.
C++ Работа с шрифтами Заносим в массив кнопки с диалогового окна: for(int index = 0; index < CONTROL_BUTTON_COUNT; index++) DDX_Control(pDX, index + 1000, ControlButtons); .... Устанавливаем текст кнопки и ее шрифт: CHAR Lang;	C++ найдите пожалуйста ошибку https://www.cyberforum.ru/ cpp-beginners/ thread709349.html почему то слепляет возраст и имя 2го и 3го рабочих #include <stdio.h> #include <stdlib.h> #include <conio.h> main() { int vA,vB,vC; char a; char b; char c;
C++ Как связать два cpp файла в одном проекте (visual studio 2010) https://www.cyberforum.ru/ cpp-beginners/ thread709325.html Есть два программы: одна создает файл.txt, вторая читает из него данные и выводит на экран. создает файл #include <iostream> #include <fstream> #include <string> using namespace std; void main () {	Найти ошибку в коде C++ пожалуйста найдите ошибку((( #include <stdio.h> #include <conio.h> #include <stdlib.h> #include <math.h> #define FALSE 0 #define TRUE 1 void menu(void);
C++ Получить матрицу из чисел, заполняющих матрицу по линиям, параллельным главной диагонали https://www.cyberforum.ru/ cpp-beginners/ thread709319.html Получить матрицу A=(aij), i,j=1..n, n<=100, элементами которой является ряд натуральных чисел, заполняющих матрицу по линиям параллельным главной диагонали, напр 1 2 6 7 3 5 8 13 4 9 12 14 10 11 15 16	C++ Битовые операции https://www.cyberforum.ru/ cpp-beginners/ thread709317.html Ребят, помогите кто чем может. Дана строка unsigned char. Необходимо через каждых 4 бита вставить 1, при этом ничего не затерев. Получается мы идем по битам элемента строки, заносим бит, сдвигаем. В результате у нас остается 7-й бит, который не уместился в байте + еще один бит, который необходимо вставить. Они переносятся в следующий элемент строки и т.д. С каждым шагом битов, которые не...
C++ Рекурсия: удаление из строки всех лишние пробелов Разработать рекурсивную функцию, удаляющую из строки все лишние пробелы. Пробелы считаются лишними, если их подряд идет более двух, если они стоят в конце строки после последней точки, если стоят после открывающегося парного знака препинания.(если можно . то на си , пожалуйста)	C++ Прошу! Уважаемые люди помогите с задачами по с++ вопрос жизни и смерти точнее отчисление из института вообщем сами задачи Семинар 2 и 3 вариант 14 и еще на фото будет прошу как можно скорее от этого зависит моя жизнь и жизнь моего рода За ранее гиганское спасибо вот https://www.cyberforum.ru/ cpp-beginners/ thread709306.html

	26.11.2012, 05:57