三元组表示的稀疏矩阵的相加

#include<iostream>
using namespace std;
const int Max=100;
const int Max_rows=4;
const int Max_cols=7;
   typedef class
        {
            public :
            int r;
            int c;
            int e;
        }TupNode;

    typedef class
        {
            public :
                int rows;
                int cols;
                int nums;
                TupNode date[Max];
        }TSMatrix;
    void display(TSMatrix &m); //输出压缩后的稀疏矩阵
    void show(TSMatrix &m);//输出完整的稀疏矩阵
    void init_matrix(TSMatrix &m,int r,int c,int a[][Max_cols]);//矩阵的初始化
    void add(TSMatrix m1,TSMatrix m2,TSMatrix &m3);
    int main()
        {
            int a1[Max_rows][Max_cols]=
            {
                {0,8,0,6,0,2,0},
                {0,0,34,54,0,0,2},
                {-7,0,0,0,0,0,0},
                {0,0,0,0,-12,0,0}

            };
            TSMatrix t1,t2,t3;
            init_matrix(t1,Max_rows,Max_cols,a1);
            show(t1);
            display(t1);
            cout<<endl;

            int a2[Max_rows][Max_cols]=
                {
                    {0,0,24,2,0,2,0 },
                    {0,8,0,7,14,2,0 },
                    {0,0,0,0,0,2,0  },
                    {0,-13,0,6,0,2,0},
                };
             init_matrix(t2,Max_rows,Max_cols,a2);
             show(t2);
             display(t2);
             add(t1,t2,t3);
             show(t3);
             display(t3);
            return 0;

        }


    void show(TSMatrix &m)//输出完整的稀疏矩阵
        {
            int i,j,k=0;
           for(i=0;i<m.rows;i++)
           {
               for(j=0;j<m.cols;j++)
               {
                   if(m.date[k].r==i&&m.date[k].c==j)
                   {
                       cout<<m.date[k].e<<"   ";
                       k++;
                   }

                   else
                   {
                      cout<<0<<"   ";
                   }
               }

               cout<<endl;
           }

           cout<<endl;

        }

        void display(TSMatrix &m) //输出压缩后的稀疏矩阵
            {
                int i;
                for(i=0;i<m.nums;i++)
                {
                    cout<<m.date[i].e<<" ";
                }

                cout<<endl;
            }

        void init_matrix(TSMatrix &m,int r,int c,int a[][Max_cols])//矩阵的初始化
            {
                    int i,j,counts=0;
                    for(i=0;i<r;i++)
                    {
                        for(j=0;j<c;j++)
                        {
                            if(a[i][j]!=0)
                            {
                                m.date[counts].r=i;
                                m.date[counts].c=j;
                                m.date[counts].e=a[i][j];
                                counts++;
                            }
                        }
                    }
                    m.nums=counts;
                    m.rows=r;
                    m.cols=c;
            }

        void add(TSMatrix m1,TSMatrix m2,TSMatrix &m3)
            {
                int i=0,j=0,nums=0;
                while(i<m1.nums&&j<m2.nums)
                {
                    if(m1.date[i].r==m2.date[j].r)
                    {
                        if(m1.date[i].c<m2.date[j].c)
                        {
                            m3.date[nums].c=m1.date[i].c;
                            m3.date[nums].r=m1.date[i].r;
                            m3.date[nums].e=m1.date[i].e;
                            nums++;
                            i++;
                        }

                        else if(m1.date[i].c>m2.date[j].c)
                        {
                            m3.date[nums].c=m2.date[j].c;
                            m3.date[nums].r=m2.date[j].r;
                            m3.date[nums].e=m2.date[j].e;
                            j++;
                            nums++;

                        }

                        else
                        {
                            if(m1.date[i].e+m2.date[j].e!=0)
                            {
                                m3.date[nums].c=m1.date[i].c;
                                m3.date[nums].r=m1.date[i].r;
                                m3.date[nums].e=m1.date[i].e+m2.date[j].e;
                                nums++;
                                i++;
                                j++;
                            }
                        }
                    }


                    else if(m1.date[i].r<m2.date[j].r)
                    {
                        m3.date[nums].r=m1.date[i].r;
                        m3.date[nums].c=m1.date[i].c;
                        m3.date[nums].e=m1.date[i].e;
                        nums++;
                        i++;
                    }

                    else
                    {
                        m3.date[nums].r=m2.date[j].r;
                        m3.date[nums].c=m2.date[j].c;
                        m3.date[nums].e=m2.date[j].e;
                        j++;
                        nums++;
                    }

                    if(i==m1.nums)
                        {
                            if(j<m2.nums)
                            {
                                while(j<m2.nums)
                                {
                                    m3.date[nums].r=m2.date[j].r;
                                    m3.date[nums].c=m2.date[j].c;
                                    m3.date[nums].e=m2.date[j].e;
                                    nums++;
                                    j++;
                                }
                                break;
                            }
                        }

                else if(j==m2.nums)
                    {
                        if(i<m1.nums)
                        {
                            while(i<m1.nums)
                            {

                                m3.date[i].r=m1.date[i].r;
                                m3.date[i].c=m1.date[i].c;
                                m3.date[i].e=m1.date[i].e;
                                i++;
                                nums++;
                            }
                            break;
                        }
                    }

                 else if(i==m1.nums&&j==m2.nums)
                    {
                        break;
                    }


                }

                m3.rows=m1.rows;
                m3.cols=m1.cols;
                m3.nums=nums;

            }
————————————————
版权声明:本文为CSDN博主「van之风」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/weixin_44077170/article/details/101480025

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
三元组稀疏矩阵是指矩阵中大部分元素为0,只有少部分元素非0。为了节约内存空间和降低计算成本,通常只将非零元素存储在矩阵中,采用三元组的方式。 三元组采用三个一维数组来表示,分别存储非零元素的行、列、和值,具体定义如下: ```c++ struct Triple { int row, col; // 行和列 ElemType value; // 元素值 }; ``` 矩阵的相加可以按照以下方法: ```c++ STATUS MatrixAdd(TSMatrix M, TSMatrix N, TSMatrix &Q) { // 判断两个矩阵是否可以相加 if (M.rows != N.rows || M.cols != N.cols) return ERROR; int m = 1, n = 1, k = 0; // 分别指向M、N、Q的当前元素 Q.rows = M.rows; Q.cols = M.cols; Q.nums = 0; // Q中实际存储的非零元素个数 while (m <= M.nums && n <= N.nums) { // 当前元素在M和N中的行和列不同时,按序号小的为当前元素 if (M.data[m].row < N.data[n].row || (M.data[m].row == N.data[n].row && M.data[m].col < N.data[n].col)) { Q.data[++k] = M.data[m]; m++; } else if (M.data[m].row > N.data[n].row || (M.data[m].row == N.data[n].row && M.data[m].col > N.data[n].col)) { Q.data[++k] = N.data[n]; n++; } else { // 当前元素在M和N中的行和列相同 Q.data[++k] = M.data[m]; Q.data[k].value += N.data[n].value; // 直接相加 m++; n++; } } // 将M或N中剩余的元素插入Q中 while (m <= M.nums) Q.data[++k] = M.data[m++]; while (n <= N.nums) Q.data[++k] = N.data[n++]; Q.nums = k; // 设置Q中实际存储的非零元素个数 return OK; } ``` 矩阵的转置可以按照以下方法: ```c++ STATUS MatrixTranspose(TSMatrix M, TSMatrix &T) { T.rows = M.cols; T.cols = M.rows; T.nums = M.nums; if (T.nums == 0) // M中没有非零元素 return OK; int num[M.cols + 1] = {0}; // 存储每一列的非零元素个数 for (int k = 1; k <= M.nums; k++) num[M.data[k].col]++; int cpot[M.cols + 1] = {0}; // 存储每一列在T中的起始位置 for (int i = 2; i <= M.cols; i++) cpot[i] = cpot[i - 1] + num[i - 1]; for (int k = 1; k <= M.nums; k++) { int j = M.data[k].col; // 获取M中的列 int q = cpot[j]; // 获取T中的位置 T.data[q].row = M.data[k].col; T.data[q].col = M.data[k].row; T.data[q].value = M.data[k].value; cpot[j]++; // 将T中当前列的位置向后移动一位 } return OK; } ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值