Найти произведение двух симметричных матриц А и В

[Alexandr-]

Найти произведение двух симметричных матриц А и В. Матрицы хранятся в одномерных массивах, где построчно записаны элементы, стоящие не ниже главной диагонали.

В общем не понимаю вообще, как реализовать 2 предложение в задачи.
К примеру из 2 массивов(допустим заполненных случайными числами) получается так(ниже код) матрица насколько я понимаю, значит нужна ещё 1 такая же и умножить их соответственно. Вот только, как реализовать это: "где построчно записаны элементы, стоящие не ниже главной диагонали".

Помогите написать задачу, кому не сложно.

#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <time.h>
int tr(int *C, int *D, int n){//функция создания 2 массивов(не знаю зачем создал)
    int i;
 for (i=1;i<=n;i++) {
C[i] = rand()%10+1; 
 printf ("%d ",C[i]); }
 printf("\n");
   for (i=1;i<=n;i++) {
D[i] = rand()%10+1; 
 printf ("%d ",D[i]); }
 printf("\n");
}
 
 
int main () {  
 int i,j,n; int *C,*D,**A;
srand(time(NULL));
printf ("vvedite razmer massiva\n");
scanf ("%d",&n);
C = (int *)malloc(n * sizeof(int));
D = (int *)malloc(n * sizeof(int)); 
A = (int **)malloc(n * sizeof(int *)); 
   
tr(C,D,n);
for (i=1;i<=n;i++) { A[i]=(int *)malloc(sizeof(int)*n); printf ("\n");
for (j=1;j<=n;j++) {
 
 if (i==1) A[i][j]=C[j]; //условие умножения массивов
 if (j==1) A[i][j]=D[i];
 if(i!=1 && j!=1) A[i][j]=C[j]*D[i];
A[1][1]=C[1]*D[1];
printf ("%3d",A[i][j]); }}
free(C);
free(D);
free(A);
getch (); return 0; }

[Alanel']

Ну если ты изначально вводишь матрицу, которая потом преобразуется в массив, то я бы так делала:

for(i=0; i<n; i++) //идем по строке матрицы
{
   for(j=0; j<i; j++)
      scanf("%d", &a[i*n+i]);//можно любую переменную, просто при считывании пропускаем числа, идем до диагонали
   for(j=i; j<n; j++)
      scanf("%d", &a[i*n+j]);//формируем одномерный массив
}

Если матрица записана в проге, то первый цикл по j убираем, а во втором будет просто присвоение a[i*n+j]=A[i][j];
Ну а перемножение в таких массивах думать надо. Надеюсь, что правильно поняла задание.

[BDA]

#include <iostream>
#include <cmath>
#include <cstdlib>
#include <ctime>

using namespace std;

int
pos(int s, int row, int col)
{
    if (row > col) {
        int tmp = row;
        row = col;
        col = tmp;
    }
    return (col - row + (s + s - row + 1) * row / 2);
}

void
print_matrix1(int a[], int n)
{
    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < n; ++j) {
            cout << a[pos(n, i, j)] << " ";
        }
        cout << endl;
    }
}

void
print_matrix2(int a[], int n)
{
    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < n; ++j) {
            cout << a[n * i + j] << " ";
        }
        cout << endl;
    }
}

int
main()
{
    srand(time(NULL));
    int count;
    cin >> count;
    int n = (sqrt(8 * count + 1) - 1) / 2;
    int a[count], b[count], c[n * n];
    for (int i = 0; i < count; ++i) {
        a[i] = rand() % 41 - 20;
        b[i] = rand() % 41 - 20;
    }
    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < n; ++j) {
            int tmp = 0;
            for (int k = 0; k < n; ++k) {
                tmp += a[pos(n, i, k)] * b[pos(n, k, j)];
            }
            c[n * i + j] = tmp;
        }
    }
    cout << "A:" << endl;
    print_matrix1(a, n);
    cout << "B:" << endl;
    print_matrix1(b, n);
    cout << "C:" << endl;
    print_matrix2(c, n);
}

Ваша задача, разобраться в программе и переписать ее с С++ на Си.
У меня матрица хранится в виде одномерного массива, причем только элементы, расположенные не ниже главной диагонали, т.е.:

1 2 3
2 4 5
3 5 6

хранится:

1 2 3 4 5 6

[Sciv]

Объясняю: для симметричных матриц нет необходимости хранить всю матрицу в двумерном массиве, так как:

1 2 3
2 5 6
3 6 4

Как видим, элемент [2,1]=[1,2], [3,1]=[1,3] и тд - в общем виде: [i,j]=[j,i]. Таким образом, все элементы ниже главной диагонали будут всего лишь дублировать элементы выше нее.

Для данной в примере матрицы одномерный массив запишется так: [1,2,3,5,6,4] при порядке матрицы n=3.

[Alexandr-]

Что нужно сделать теперь понимаю, но все равно не получается из одномерного массива получить такую матрицу(просто не понимаю, что с ним делать то).
Помогите с частью данного кода пожалуйста.

[BDA]

Чем Вас не устраивает код из моего поста выше?

Кстати, программа работала неправильно.
Исправлена.
http://ideone.com/adRU5x

[Alexandr-]

Там С++, а я его ещё даже не начинал проходить, поэтому на С переписать не могу. Не могли бы вы написать часть кода формирования матрицы из одномерного массива на С, если не сложно.

[BDA]

#include <math.h>
#include <stdlib.h>
#include <time.h>

int
pos(int s, int row, int col)
{
    if (row > col) {
        int tmp = row;
        row = col;
        col = tmp;
    }
    return (col - row + (s + s - row + 1) * row / 2);
}

void
print_matrix1(int a[], int n)
{
    int i, j;
    for (i = 0; i < n; ++i) {
        for (j = 0; j < n; ++j) {
            printf("%d ", a[pos(n, i, j)]);
        }
        printf("\n");
    }
}

void
print_matrix2(int a[], int n)
{
    int i, j;
    for (i = 0; i < n; ++i) {
        for (j = 0; j < n; ++j) {
            printf("%d ", a[n * i + j]);
        }
        printf("\n");
    }
}

int
main()
{
    srand(time(NULL));
    int count;
    scanf("%d", &count);
    int n = (sqrt(8 * count + 1) - 1) / 2;
    int a[count], b[count], c[n * n];
    int i, j, k;
    for (i = 0; i < count; ++i) {
        a[i] = rand() % 41 - 20;
        b[i] = rand() % 41 - 20;
    }
    for (i = 0; i < n; ++i) {
        for (j = 0; j < n; ++j) {
            int tmp = 0;
            for (k = 0; k < n; ++k) {
                tmp += a[pos(n, i, k)] * b[pos(n, k, j)];
            }
            c[n * i + j] = tmp;
        }
    }
    printf("A:\n");
    print_matrix1(a, n);
    printf("B:\n");
    print_matrix1(b, n);
    printf("C:\n");
    print_matrix2(c, n);
    return 0;
}