现在的位置: 首页 > 综合 > 正文

求用一维数组表示的矩阵的乘积

2017年11月29日 ⁄ 综合 ⁄ 共 1554字 ⁄ 字号 评论关闭


版权所有,转载请注明出处!


写这篇博客并不是为了探讨算法效率的提升,而是尝试以一种更能够让人们读懂的方式进行书写程序,

这里我写了三个程序,分别使用三种编程方式进行书写,希望读者能够找到更加适合自己的方式。

1.直接按照一维数组的习惯操作:

#include <iostream>
#define N 3
using namespace std;
int main (int argc, const char * argv[])
{
    int A[] = {1,2,3,4,5,6,7,8,9};
    int B[] = {7,6,5,1,2,3,2,8,7};
    int C[9]= {0};
    memset(C, 0, sizeof(C));
    for(int i =0; i < N; i ++)
        for(int j = 0; j < N; j ++)
        {
            int ele = 0;
            for(int k =0; k < N; k ++)
            {
                ele += A[i*N + k]*B[k*N + j];
            }
            C[i*N +j] = ele;
        }
    
    for(int p =0; p < N; p++){
        for(int q = 0; q < N; q ++)
        {
            cout << C[p*N +q]<<" ";
        }
        cout<<endl;
    }
    return 0;
}


2.使用宏定义简化思维方式:

#include <iostream>
#define N 3
#define ELEMENT(A,i,j) A[i*N + j] 
using namespace std;
int main (int argc, const char * argv[])
{
    int A[] = {1,2,3,4,5,6,7,8,9};
    int B[] = {7,6,5,1,2,3,2,8,7};
    int C[9]= {0};
    memset(C, 0, sizeof(C));
    for(int i =0; i < N; i ++)
        for(int j = 0; j < N; j ++)
        {
            int ele = 0;
            for(int k =0; k < N; k ++)
            {
              
                ele += ELEMENT(A, i , k)*ELEMENT(B , k, j );
            }
            
            ELEMENT(C , i, j) = ele;
        }
    
    for(int p =0; p < N; p++){
        for(int q = 0; q < N; q ++)
        {
            cout << ELEMENT(C, p , q )<<" ";
        }
        cout<<endl;
    }

    return 0;
    
}


3.使用内联函数简化思维方式:

#include <iostream>
#define N 3
using namespace std;

inline int element(int *A, size_t i, size_t j)//实现下标操作
{
    return A[i*N + j];
}
inline void setElenent(int *A, size_t i, size_t j, int element)//设置指定的元素
{
    A[i*N + j] = element;
}

int main (int argc, const char * argv[])
{
    int A[] = {1,2,3,4,5,6,7,8,9};
    int B[] = {7,6,5,1,2,3,2,8,7};
    int C[9]= {0};
    for(int i =0; i < N; i ++)
        for(int j = 0; j < N; j ++)
        {
            int ele = 0;
            for(int k =0; k < N; k ++)
            {
                ele += element(A, i, k)*element(B , k, j);
            }
            setElenent(C,i,j,ele);
        }
    
    for(int p =0; p < N; p++){
        for(int q = 0; q < N; q ++)
        {
            cout << element(C, p , q )<<" ";
        }
        cout<<endl;
    }
    
    return 0;
    
}


从上面三个程序可以看出,他们使用的算法是完全一样的,仅仅是在书写的时候采用了不同的技巧,

第一个程序没有采用任何的技巧,直接进行操作,缺点就是容易思路混乱,而且很容易出错。

相比之下,第二和第三段程序采用了宏定义和内联函数来简化思维方式,使得我们能够像操作

二维数组一样对一维数组进行操作。


抱歉!评论已关闭.