Перемножение матриц.

@Stas0n · Регистрация: 13.07.2011

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

Нужен код для перемножения матрицы и столбца за минимально возможное время.
Порядок матрицы ( и столбца ) огромен - около 100000. Помогите найти оптимальный код!
Ищу уже 2 день... что-то безрезультатно как то(((

@Mиxaил · 21.07.2011, 11:16

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
int multMatrixBlock(double *a, const int n, double *b, const int m, double *res, double *c, double *d, double*e){
    int n_block=n/m;
    int m_small;
    int i,j,k,o;
    m_small=n-n_block*m;
    clearMatrix(res, n*n);
    for ( i = 0; i <= n_block; i++)
    {
        for ( j = 0; j <= n_block; j++)
        {
            for ( k = 0; k <=n_block; k++)
            {
                if (i<n_block && j<n_block && k<n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixSq(c, m, d, e);
                    setSubMatrixAdded(res, n, e, m, i, j);
                }
                else if (i<n_block && j<n_block && k==n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m, d,m_small, m, e);
                    setSubMatrixAdded(res, n, e, m, i, j);  
                }
                else if (i==n_block && j<n_block && k<n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m_small, d,m, m, e);
                    setSubMatrixAdded(res, n, e, m, i, j);
                }
                else if (i==n_block && j<n_block && k==n_block) {   
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m_small, d,m_small, m, e);
                    setSubMatrixAdded(res, n, e, m, i, j);  
                }
                else if (i<n_block && j==n_block && k<n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m, d,m, m_small, e);
                    setSubMatrixAdded(res, n, e, m, i, j);
                }
                else if (i<n_block && j==n_block && k==n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m, d,m_small, m_small, e);
                    setSubMatrixAdded(res, n, e, m, i, j);
                }
                else if (i==n_block && j==n_block && k<n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixRect(c, m_small, d,m, m_small, e);
                    setSubMatrixAdded(res, n, e, m, i, j);  
                }
                else if (i==n_block && j==n_block && k==n_block) {
                    getSubMatrix(a, n, c, m, i, k);
                    getSubMatrix(b, n, d, m, k, j);
                    multMatrixSq(c, m_small, d, e);
                    setSubMatrixAdded(res, n, e, m, i, j);  
                }
            }
        }
    }
    return 0;
}

нашел в интернете

Добавлено через 1 минуту
матрица хранится в одномерном массиве.

@Stas0n · 21.07.2011, 11:21 **[ТС]**

Вот только там нет определния кучи ф-ций...

@Mиxaил · 21.07.2011, 11:38

Сообщение от hello19

Вот только там нет определния кучи ф-ций...

исходые функции не удалось найти, сорри... Зато есть еще вот это:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
int multMatrixSq(double *a, const int n, double *b, double *c){
    double s1=0, s2=0, s3=0, s4=0;
    double *f, *g, *h;
    int i, j, k;
    if (n%2==0){
        for (i = 0; i < n-1; i+=2)
        {
            f=a+i*n;
            for (j = 0; j < n-1; j+=2)
            {
                g=b+j;
                for (k=0;k<n;k++){
                    s1+=f[k]*g[n*k];    
                    s2+=f[n+k]*g[n*k];  
                    s3+=f[k]*g[n*k+1];  
                    s4+=f[n+k]*g[n*k+1];    
                }
                h=c+i*n+j;
                h[0]=s1; s1=0;
                h[n]=s2; s2=0;
                h[1]=s3; s3=0;
                h[n+1]=s4; s4=0;
            }
        }
    }
    else {
        for (i = 0; i < n-1; i+=2)
            {
                f=a+i*n;
                for (j = 0; j < n-1; j+=2)
                {
                    g=b+j;
                    for (k=0;k<n;k++){
                        s1+=f[k]*g[n*k];
                        s2+=f[n+k]*g[n*k];
                        s3+=f[k]*g[n*k+1];
                        s4+=f[n+k]*g[n*k+1];
                    }
                    h=c+i*n+j;
                    h[0]=s1; s1=0;
                    h[n]=s2; s2=0;
                    h[1]=s3; s3=0;
                    h[n+1]=s4; s4=0;
                }
                if  (j==n-1){
                    g=b+j;
                    for (k=0;k<n;k++){
                        s1+=f[k]*g[n*k];    
                        s2+=f[n+k]*g[n*k];  
                    }
                    h=c+i*n+j;
                    h[0]=s1; s1=0;
                    h[n]=s2; s2=0;
                }
            }
            if (i==n-1){
                f=a+i*n;
                for (j = 0; j < n-1; j+=2)
                {
                    g=b+j;
                    for (k=0;k<n;k++){
                        s1+=f[k]*g[n*k];    
                        s3+=f[k]*g[n*k+1];  
                    }
                    h=c+i*n+j;
                    h[0]=s1; s1=0;
                    h[1]=s3; s3=0;
                }
                if  (j==n-1){
                    g=b+j;
                    for (k=0;k<n;k++){
                        s1+=f[k]*g[n*k];    
                    }
                    h=c+i*n+j;
                    h[0]=s1; s1=0;
                }
            }
        }
    return 0;
}

тут, вроде, без сторонних функций

Jupiter · 21.07.2011, 11:42

Не по теме:

Сообщение от hello19

около 100000

лол) я представляю сколько веремени будут перемножаться такие матрицы кодом приведеным выше, а ТС рекомендую смирится с этим заданием, к сожалению это не студенческая задача

а по теме: вам нужно ленточное перемножение матриц

@asics · 21.07.2011, 11:49

Вот подобная тема -- Распараллеливание

Новые блоги и статьи Все статьи Все блоги /
PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Access VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .
Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .

@Stas0n 3 / 4 / 1 Регистрация: 13.07.2011 Сообщений: 313

	Перемножение матриц. 21.07.2011, 11:08. Показов 3585. Ответов 5 Метки нет (Все метки) Нужен код для перемножения матрицы и столбца за минимально возможное время. Порядок матрицы ( и столбца ) огромен - около 100000. Помогите найти оптимальный код! Ищу уже 2 день... что-то безрезультатно как то((( 0

@Stas0n 3 / 4 / 1 Регистрация: 13.07.2011 Сообщений: 313
	21.07.2011, 11:21 [ТС]
	Вот только там нет определния кучи ф-ций... 0

@asics Freelance 2891 / 1826 / 356 Регистрация: 09.09.2010 Сообщений: 3,841
	21.07.2011, 11:49
	Вот подобная тема -- Распараллеливание 0