[C/C++][矩陣內積]

/*
* File Name: Innerproduct.c
* Author: MH
* Since 2011/03/21
* Toolkit: Dev C++
*/
 
# include <stdlib.h>
# include <stdio.h>
 
int main(){
 
    int i, j, k;
    int array1[3][3] = {{1, 2, 3}, {4, 5, 6}, {7, 8, 9}};    //array[row][column]
    int array2[3][3] = {{2, 3, 4}, {5, 6, 7}, {8, 9, 1}};
    int array3[3][3] = {0};
 
    printf("A = \n");
 
    for(i=0; i<3; i++){
 
        for(j=0; j<3; j++){
            printf("%-3d", array1[i][j]);
        }
 
        printf("\n");
 
    }
 
    printf("\nB = \n");
 
    for(i=0; i<3; i++){
 
        for(j=0; j<3; j++){
            printf("%-3d", array2[i][j]);
        }
 
        printf("\n");
 
    }
 
    printf("\n\nA・B =\n");
 
    for(i=0; i<3; i++){
 
        for(j=1; j<3; j++){
 
            for(k=1; k<=3; k++){
                array3[i][j] = array3[i][j] + array1[i][k] * array2[k][j];
                // array1: from left to right, from top to bottom
                // array2: from top to bottom, from left to right
            }
 
            printf("%-5d", array3[i][j]);
 
        }
 
        printf("\n\n");
 
    }
 
    printf("\n");
 
    system("Pause");
    return 0;
}