C++稀疏矩阵的各种基本运算并实现加法乘法

 更新时间:2019年02月20日 15:14:17   作者:fanxingyue  
今天小编就为大家分享一篇关于C++稀疏矩阵的各种基本运算并实现加法乘法,小编觉得内容挺不错的,现在分享给大家,具有很好的参考价值,需要的朋友一起跟随小编来看看吧

代码:

#include <iostream>
#include<malloc.h>
#include<cstdio>
using namespace std;
#define M 4
#define N 4
#define MaxSize 100
typedef int ElemType;
typedef struct
{
  int r;
  int c;
  ElemType d;///元素值
} TupNode; ///三元组定义
typedef struct
{
  int rows;
  int cols;
  int nums;
  TupNode data[MaxSize];
} TSMatrix; ///三元组顺序表定义
void CreatMat(TSMatrix &t,ElemType A[M][N])
{
  t.rows=M;
  t.cols=N;
  t.nums=0;
  for(int i=0; i<M; i++)
    for(int j=0; j<N; j++)
      if(A[i][j]!=0)
      {
        t.data[t.nums].r=i;
        t.data[t.nums].c=j;
        t.data[t.nums].d=A[i][j];
        t.nums++;
      }
}
bool Value(TSMatrix &t,ElemType x,int i,int j)
{
  int k=0,k1;
  if(i>=t.rows||j>=t.cols)
    return false;
  while(k<t.nums&&i>t.data[k].r)k++;
  while(k<t.nums&&i==t.data[k].r&&j>t.data[k].c)k++;
  if(t.data[k].r==i&&t.data[k].c==j)
    t.data[k].d=x;
  else
  {
    for(k1=t.nums-1; k1>=k; k1--)
    {
      t.data[k1+1].r=t.data[k].r;
      t.data[k1+1].c=t.data[k].c;
      t.data[k1+1].d=t.data[k].d;
    }
    t.data[k].r=i;
    t.data[k].c=j;
    t.data[k].d=x;
    t.nums++;
  }
  return true;
}
bool Assign(TSMatrix t,ElemType &x,int i,int j)
{
  int k=0;
  if(i>=t.rows||j>=t.cols)
    return false;
  while(k<t.nums&&i>t.data[k].r)k++;
  while(k<t.nums&&i==t.data[k].r&&j>t.data[k].c)k++;
  if(t.data[k].r==i&&t.data[k].c==j)
    x=t.data[k].d;
  else
    x=0;
  return true;
}
void DispMat(TSMatrix t)
{
  if(t.nums<=0)
    return ;
  printf("\t%d\t%d\t%d\n",t.rows,t.cols,t.nums);
  printf("\t-----------------\n");
  for(int i=0; i<t.nums; i++)
    printf("\t%d\t%d\t%d\n",t.data[i].r,t.data[i].c,t.data[i].d);
}
void TranMat(TSMatrix t,TSMatrix &tb)
{
  int i,j,k=0;
  tb.rows=t.cols;
  tb.cols=t.rows;
  tb.nums=t.nums;
  if(t.nums!=0)
  {
    for(i=0; i<t.cols; i++)
      for(j=0; j<t.nums; j++)
        if(t.data[j].c==i)
        {
          tb.data[k].r=t.data[j].c;
          tb.data[k].c=t.data[j].r;
          tb.data[k].d=t.data[j].d;
          k++;
        }
  }
}
bool MatAdd(TSMatrix a,TSMatrix b,TSMatrix &c)
{
  int i=0,j=0,k=0;
  ElemType v;
  if(a.rows!=b.rows||a.cols!=b.cols)
    return false;
  c.rows=a.rows;
  c.cols=a.cols;
  while(i<a.nums&&j<b.nums)
  {
    if(a.data[i].r==b.data[j].r)///先控制行相等
    {
      if(a.data[i].c<b.data[j].c)
      {
        c.data[k].r=a.data[i].r;
        c.data[k].c=a.data[i].c;
        c.data[k].d=a.data[i].d;
        k++;
        i++;
      }
      else if(a.data[i].c>b.data[j].c)
      {
        c.data[k].r=b.data[j].r;
        c.data[k].c=b.data[j].c;
        c.data[k].d=b.data[j].d;
        k++;
        j++;
      }
      else
      {
        v=a.data[i].d+b.data[j].d;
        if(v!=0)
        {
          c.data[k].r=a.data[i].r;
          c.data[k].c=a.data[i].c;
          c.data[k].d=v;
          k++;
        }
        i++;
        j++;
      }
    }
    else if(a.data[i].r<b.data[j].r)
    {
      c.data[k].r=a.data[i].r;
      c.data[k].c=a.data[i].c;
      c.data[k].d=a.data[i].d;
      k++;
      i++;
    }
    else
    {
      c.data[k].r=b.data[j].r;
      c.data[k].c=b.data[j].c;
      c.data[k].d=b.data[j].d;
      k++;
      j++;
    }
    c.nums=k;
  }
  return true;
}
int getvalue(TSMatrix c,int i,int j)
{
  int k=0;
  while(k<c.nums&&(c.data[k].r!=i||c.data[k].c!=j))
    k++;
  if(k<c.nums)
    return (c.data[k].d);
  else
    return (0);
}
bool MatMul(TSMatrix a,TSMatrix b,TSMatrix &c)
{
  int i,j,k,p=0;
  ElemType s;
  if(a.cols!=b.rows)
    return false;
  for(i=0; i<a.rows; i++)
    for(j=0; j<b.cols; j++)
    {
      s=0;
      for(k=0; k<a.cols; k++)
        s+=getvalue(a,i,k)*getvalue(b,k,j);
      if(s!=0)
      {
        c.data[p].r=i;
        c.data[p].c=j;
        c.data[p].d=s;
        p++;
      }
    }
  c.rows=a.rows;
  c.cols=b.cols;
  c.nums=p;
  return true;
}
int main()
{
  ElemType a1[N][N]={{1,0,3,0},{0,1,0,0},{0,0,1,0},{0,0,1,1}};
  ElemType b1[M][M]={{3,0,0,0},{0,4,0,0},{0,0,1,0},{0,0,0,2}};
  TSMatrix a,b,c;
  CreatMat(a,a1);
  CreatMat(b,b1);
  printf("a的三元组:\n");
  DispMat(a);
  printf("b的三元组:\n");
  DispMat(b);
  printf("a转置为c\n");
  TranMat(a,c);
  printf("c的三元组\n");
  DispMat(c);
  printf("c=a+b\n");
  MatAdd(a,b,c);
  printf("c的三元组:\n");
  DispMat(c);
  printf("c=a*b\n");
  MatMul(a,b,c);
  printf("c的三元组:\n");
  DispMat(c);
  return 0;
}

总结

以上就是这篇文章的全部内容了,希望本文的内容对大家的学习或者工作具有一定的参考学习价值,谢谢大家对脚本之家的支持。如果你想了解更多相关内容请查看下面相关链接

相关文章

  • C++小知识:不要节约代码行数

    C++小知识:不要节约代码行数

    今天小编就为大家分享一篇关于C++小知识:不要节约代码行数,小编觉得内容挺不错的,现在分享给大家,具有很好的参考价值,需要的朋友一起跟随小编来看看吧
    2019-01-01
  • C++实现LeetCode(134.加油站问题)

    C++实现LeetCode(134.加油站问题)

    这篇文章主要介绍了C++实现LeetCode(134.加油站问题),本篇文章通过简要的案例,讲解了该项技术的了解与使用,以下就是详细内容,需要的朋友可以参考下
    2021-07-07
  • C++实现Go的defer功能(示例代码)

    C++实现Go的defer功能(示例代码)

    defer和go一样都是Go语言提供的关键字。defer用于资源的释放,会在函数返回之前进行调用。接下来通过本文给大家介绍C++实现Go的defer功能,感兴趣的朋友跟随小编一起看看吧
    2021-07-07
  • c++隐式类型转换示例分享

    c++隐式类型转换示例分享

    这篇文章主要介绍了c++隐式类型转换的二个示例,需要的朋友可以参考下
    2014-03-03
  • 解析C/C++ Capstone 引擎源码编译问题

    解析C/C++ Capstone 引擎源码编译问题

    Capstone的编译非常简单只需要一步即可轻松得到对应的Lib库文件,如下将介绍该引擎如何被编译,以及简单的测试编译,这篇文章主要介绍了C/C++ Capstone 引擎源码编译,需要的朋友可以参考下
    2022-09-09
  • C语言中数据结构之链式基数排序

    C语言中数据结构之链式基数排序

    这篇文章主要介绍了C语言中数据结构之链式基数排序的相关资料,希望通过本文能帮助到大家,需要的朋友可以参考下
    2017-09-09
  • 详解C语言#define预处理宏定义

    详解C语言#define预处理宏定义

    本文主要介绍了C语言#define预处理宏定义,文中通过示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2021-09-09
  • C程序中唯一序列号的生成实例详解

    C程序中唯一序列号的生成实例详解

    这篇文章主要介绍了C程序中唯一序列号的生成实例详解的相关资料,需要的朋友可以参考下
    2017-06-06
  • C++使用WideCharToMultiByte函数生成UTF-8编码文件的方法

    C++使用WideCharToMultiByte函数生成UTF-8编码文件的方法

    用来映射Unicode字符串的WideCharToMultiByte函数经常被用来进行UTF-8编码的转换,以下我们将看到C++使用WideCharToMultiByte函数生成UTF-8编码文件的方法,首先先来对WideCharToMultiByte作一个详细的了解:
    2016-06-06
  • C语言实现24点问题详解

    C语言实现24点问题详解

    24点问题就是在屏幕上输入1〜10范围内的4个整数(可以有重复),对它们进行加、减、乘、除四则运算后(可以任意的加括号限定计算的优先级),寻找计算结果等于24的表达式。本文将通过C语言实现24点问题的求解,需要的可以参考一下
    2021-12-12

最新评论