3行3列の行列 A に対して、Aの2乗, Aの3乗 を求めよ。

質問

3行3列の行列 A に対して、Aの2乗, Aの3乗 を求めよ。
という問題で

#include <stdio.h>

#define N 3

int A[N][N];
int A_2[N][N]; /* 行列Aを2乗したもの */
int A_3[N][N]; /* 行列Aを3乗したもの */

int main()
{
int i,j,k,l;

/* 3行3列の係数行列Aの成分を入力 */
printf("係数行列Aを%d行%d列で入力してください\n", N, N);

for( i=0; i<N; i++)
{
for( j=0; j<N; j++)
{
printf("A[%d][%d]=", i+1, j+1);
scanf("%d", &A[i][j]);
}
}

ここまでがあたえられえているのですが、単純にA_2＝A×Aとしてみたのですが、セグメントエラーが出てしまいうまく出力できませんでした。
A_2とA_3の計算方法とあたえられているプログラムの動作を教えてください。
よろしくお願いします

回答No.3

asuncion · ベストアンサー

| A B C | | a b c | | Aa+Bd+Gc Ab+Be+Ch Ac+Bf+Ci |
| D E F | | d e f | | Da+Ed+Fc Db+Ee+Fh Dc+Ef+Fi |
| G H I | | g h i | | Ga+Hd+Ic Gb+He+Ih Gc+Hf+Ii |

本当ですか？
積のところで、gがどこにも登場していません。

まあ、それはさておき、
X = A * B
のとき、
積Xのi行j列要素は、Aのi行目とBのj行目の積和です。

いきなり3行3列の積を求めようとするより、
まずは2行2列あたりでじゅうぶん考え方をマスターなさるのが
よいのではないでしょうか。

お礼

返信が遅れてしまい申し訳ありません。

＞本当ですか？
＞積のところで、gがどこにも登場していません。

間違った1行目をコピーしたためｇが出ていませんでした。
2行2列から勉強しようと思います。
こんな私に付き合っていただき、ありがとうございました。

回答No.4

asuncion

＞積Xのi行j列要素は、Aのi行目とBのj行目の積和です。

おっと失礼しました。

誤：Aのi行目とBのj行目の積和
正：Aのi行目とBのj列目の積和

回答No.2

asuncion

＞A_2[i][j] = A[i][j] * A[i][j];

間違っています。
行列の掛け算の方法を勉強し直してください。

補足

ご指摘ありがとうございます。

例として、3行3列の行列の積は

| A B C | | a b c | | Aa+Bd+Gc Ab+Be+Ch Ac+Bf+Ci |
| D E F | | d e f | | Da+Ed+Fc Db+Ee+Fh Dc+Ef+Fi |
| G H I | | g h i | | Ga+Hd+Ic Gb+He+Ih Gc+Hf+Ii |

となるので、このプログラムでその計算をさせるには
for(i=0;i<3;i++)
for(j=0;j<3;j++)
A_2[i][J] += A[I][j]*A[j][i];
と計算すればよいのでしょうか？

回答No.1

asuncion

＞単純にA_2＝A×Aとしてみた

考え方は正しいです。
うまくいかない、というコードを貼ってください。

補足

回答ありがとうございます。
A_2[i][j] = A[i][j] * A[i][j];
としました。
あたえられているプログラムだと
int i,j,k,l;
となっていますが、k,lはどのようにつかうのでしょうか？