行列の積の求めかた

質問

行列の積を求めるプログラムを組んだのですがうまく表示されないです。
x=１,２,３
　４,５,６
y=１,５
　５,３
　８,１
としの行列の積です。
#include <stdio.h>
int main(void)
{
int i,j;
int ma[2][3] = {{1,2,3},{4,5,6}};
int mb[3][2] = {{1,5},{5,3},{8,1}};
int mc[2][2] = {0};

for(i=0;i<2;i++)
for(j=0;j<3;j++)
mc[i][i] += ma[i][j]*mb[j][i];
mc[1-i][i] += ma[1-i][j]*mb[j][i];

for(i=0;i<2;i++){
for(j=0;j<2;j++)
printf("%3d",mc[i][j]);
putchar('\n');
}

return(0);

}
結果は↓のようになります。
38 0
0 41
まだプログラミングをはじめて日が浅いので基本てきなとこが抜けてる
かもしれません。よろしくお願いします。

回答No.5

N-Ishikawa

行列の積ですか。もう２０年前のことなのでどう計算すればいいのか
サッパリ忘れてしまいましたが．．．．．
大括弧で括っていないためにループの範囲が間違っているように見えます。
※for文やif文で大括弧を使用しないと直後の1文しか実行の対象にしてくれません。

for(i=0;i<2;i++)
for(j=0;j<3;j++)
mc[i][i] += ma[i][j]*mb[j][i]; ←ここしか実行されない
mc[1-i][i] += ma[1-i][j]*mb[j][i]; ←i=2,j=3で一回しか実行されない。ここもループさせるなら上の行と大括弧で括る必要あり

ループ範囲を可視化しやすいように以下のように書いた方がいいかと思います。

for(i=0;i<2;i++)
for(j=0;j<3;j++) {
mc[i][i] += ma[i][j]*mb[j][i];
mc[1-i][i] += ma[1-i][j]*mb[j][i];
}

ちなみに行列の積のサンプルは以下URLでも記載されていますね。
参考にされては如何でしょうか？

http://mitv2.net/matrix/multi_matrix.html

回答No.4

Tacosan

今は最終的に 2×2行列になるからいいけど, 10×15行列になるときにも 10行 (あるいは 15行) ずらっと並べるつもり?
行や列の数が不定だった場合にどうすればいいと思う?

回答No.3

Tacosan

読み返してみたら #1 の最初の文は誤解を生みそうな気がするので変更.
上の for の部分は忘れてしまって, #2 でいわれるような 3重ループに直しちゃってください (補足ありがとうございます＞#2).

お礼

回答ありがとうございます。
まだ一般的になぜダメかすぐにわからないです　汗
まずは三重ループをがんばってみます！。

回答No.2

asuncion

#1さんの回答に補足します。
#1さんが３重ループと書かれているのは、

ｉ行目の処理
　ｊ列目の処理
　　積のｉ行ｊ列の要素は、「掛けられる行列のｉ行ｋ列　＊　掛ける行列のｋ行ｊ列」の和

という意味です。

お礼

回答ありがとうございます。三重ループで作ってみます！

回答No.1

Tacosan

mc に積を求めている上の 2つの for文と, 求めた mc を表示している下の 2つの for 文の違いをじっくりと見てください.
でも, 「行列の積」を求めるときには普通 3重ループにするものだよ. 今の場合は「たまたま」 2行 2列だからいいかもしれないけど, 一般的にはダメなことは理解できるよね?