矩阵的转置
创新互联服务项目包括清原网站建设、
清原网站制作、清原网页制作以及清原网络营销策划等。多年来，我们专注于互联网行业，利用自身积累的技术优势、行业经验、深度合作伙伴关系等，向广大中小型企业、政府机构等提供互联网行业的解决方案，
清原网站推广取得了明显的社会效益与经济效益。目前，我们服务的客户以成都为中心已经辐射到清原省份的部分城市，未来相信会继续扩大服务区域并继续获得客户的支持与信任！
将原矩阵的行、列对换，也就是将[i][j]和[j][i]位置上的数据对换。
程序代码：
#include //稀疏矩阵push pop operator[] 和顺序表一致
template
struct Triple //定义一个三元组 可以直接访问的定义成struct
{
size_t _row;
size_t _col;
T _value;
Triple(size_t row, size_t col, const T& value)
:_row(row)
, _col(col)
, _value(value)
{}
};
template
class SparseMatrix
{
public:
SparseMatrix(size_t M, size_t N, const T&invalid)
:_M(M)
, _N(N)
, invalid(invalid)
{
}
SparseMatrix(const T* a, size_t M, size_t N,const T& invalid)//const T& invalid表示哪个是无效数据
:_M(M)
, _N(N)
, invalid(invalid)
{
for (size_t i = 0; i < M; ++i)
{
for (size_t j = 0; j < N; ++j)
{
if (a[i*N + j] != invalid) //不等于无效值
{
Triple t(i, j, a[i*N + j]);
_a.push_back(t);
}
}
}
}
void Display()
{
size_t index = 0;
for (size_t i = 0; i < _M; ++i)
{
for (size_t j = 0; j < _N; ++j)
{
if (index<_a.size()&&
i == _a[index]._row && j == _a[index]._col)
{
cout << _a[index].value << " ";
++index;
}
else
{
cout << _invalid << " ";
}
}
cout << endl;
}
cout << endl;
}
SparseMatrix Transport() //转置
{
//时间复杂度 O(有效数据的个数*N（列））
SparseMatrix sm(_N,_M,_invalid);
for (size_t i = 0; i < N; ++i)
{
size_t index = 0;
while (index < _a.size())
{
if (_a[index].col == i)
{
Triple t(_a[index]._col, _a[index]._row, _a[index]._value);
sm._a.push_back(t);
}
++index;
}
}
return sm;
}
protected: //存三元组数组
//Triple* _a; 直接用动态顺序表
vector> _a;
size_t _M;
size_t _N;
T _invalid;
};
void Test2()
{
int a[6][5] = { { 1, 0, 3, 0, 5 },
{ 0, 0, 0, 0, 0 },
{ 0, 0, 0, 0, 0 },
{ 2, 0, 4, 0, 6 },
{ 0, 0, 0, 0, 0 },
{ 0, 0, 0, 0, 0 } };
SparseMatrix sm((int*)a,6,5,0); //强制转换成一维数组 数组 6行 5列 非法值0
sm.Display();
SparseMatrix sm2 = sm.Transport();
sm2.Display();
}
#include
using namespace std;
#include
#include"Matrix.h"
int main()
{
//Test1();
Test2();
system("pause");
return 0;
}
运行结果：
1 0 0 2 0 0
0 0 0 0 0 0
3 0 0 4 0 0
0 0 0 0 0 0
5 0 0 6 0 0
快速转置
程序代码：
#include //稀疏矩阵push pop operator[] 和顺序表一致
template
struct Triple //定义一个三元组 可以直接访问的定义成struct
{
size_t _row;
size_t _col;
T _value;
Triple(size_t row=0, size_t col=0, const T& value=T())//const临时对象具有常性
:_row(row)
, _col(col)
, _value(value)
{}
};
template
class SparseMatrix
{
public:
SparseMatrix(size_t M, size_t N, const T&invalid)
:_M(M)
, _N(N)
, invalid(invalid)
{
}
SparseMatrix(const T* a, size_t M, size_t N,const T& invalid)//const T& invalid表示哪个是无效数据
:_M(M)
, _N(N)
, invalid(invalid)
{
for (size_t i = 0; i < M; ++i)
{
for (size_t j = 0; j < N; ++j)
{
if (a[i*N + j] != invalid) //不等于无效值
{
Triple t(i, j, a[i*N + j]);
_a.push_back(t);
}
}
}
}
void Display()
{
size_t index = 0;
for (size_t i = 0; i < _M; ++i)
{
for (size_t j = 0; j < _N; ++j)
{
if (index<_a.size()&&
i == _a[index]._row && j == _a[index]._col)
{
cout << _a[index].value << " ";
++index;
}
else
{
cout << _invalid << " ";
}
}
cout << endl;
}
cout << endl;
}
SparseMatrix FastTransport() //快速转置
{
//时间复杂度 O(有效数据个数+N)
SparseMatrix sm(_N, _M, _invalid);
int* rowCounts = new int[_N];//统计转置后数据个数
memset(rowCounts, 0, sizeof(int)*_N);
size_t index = 0;
while (index < _a.size())
{
rowCounts[_a[index].col]++;
++index;
}
int rowStart = new int[_N];
rowStart[0] = 0;
for (size_t i = 1; i < _N; ++i)
{
rowStart[i] = rowStart[i - 1] + rowCounts[i - 1];
}
index = 0;
sm._a.resize(_a.size())';' //不能用pushback
while (index < _a.size()) //_a.size有效数据的个数
{
size_t row = _a[index]._col;
int& start = rowStart[row];
Triple t(_a[index]._col, _a[index]._row, _a[index]._value);
sm._a[start] = t;
++start; // rowStart[row]里的值++
++index；
}
delete[] rowCounts;
delete[] rowStart;
return sm;
}
protected: //存三元组数组
//Triple* _a; 直接用动态顺序表
vector> _a;
size_t _M;
size_t _N;
T _invalid;
};
void Test2()
{
int a[6][5] = { { 1, 0, 3, 0, 5 },
{ 0, 0, 0, 0, 0 },
{ 0, 0, 0, 0, 0 },
{ 2, 0, 4, 0, 6 },
{ 0, 0, 0, 0, 0 },
{ 0, 0, 0, 0, 0 } };
SparseMatrix sm((int*)a,6,5,0); //强制转换成一维数组 数组 6行 5列 非法值0
sm.Display();
SparseMatrix sm2 = sm.Transport();
sm2.Display();
SparseMatrix sm3 = sm.FastTransport();
sm3.Display();
}
#include
using namespace std;
#include
#include"Matrix.h"
int main()
{
//Test1();
Test2();
system("pause");
return 0;
}
运行结果：
1 0 0 2 0 0
0 0 0 0 0 0
3 0 0 4 0 0
0 0 0 0 0 0
5 0 0 6 0 0
创新互联www.cdcxhl.cn，专业提供香港、美国云服务器，动态BGP最优骨干路由自动选择，持续稳定高效的网络助力业务部署。公司持有工信部办法的idc、isp许可证， 机房独有T级流量清洗系统配攻击溯源，准确进行流量调度，确保服务器高可用性。佳节活动现已开启，新人活动云服务器买多久送多久。
本文标题：c++矩阵的转置和快速转置-创新互联
网址分享：
http://bjjierui.cn/article/dcicpi.html
