数据结构试验二

#include"iostream"

#include"stdio.h"

#include"algorithm"

#include"string.h"

#include"cmath"

#define mx 1005

using namespace std;

int matrix[mx][mx];

//三元组的结构定义

struct thNode //非零元素

{

int row,col;

int item;

};

struct thMatrix

{

thNode *data;

int m,n,t;//分别表示行数、列数和非零元素个数

};

//初始化三元组表

void initMatrix_Th(thMatrix &M)

{

M.data=new thNode[mx+1];

M.m=0;

M.n=0;

M.t=0;

}

//快速转置

void transMatrix_2(thMatrix M,thMatrix &T)

{

initMatrix_Th(T);

T.m=M.n;

T.n=M.m;

T.t=M.t;

int i,j,p,q;

int num[100];

int cpot[100];

if(T.t)

{

for(i=1;i<=M.n;i++)

num[i]=0;//M中每一列的非零元素个数初始化为0

for(i=1;i<=M.t;i++)

++num[M.data[i].col];//M中每一列的非零元素个数存储到一个数组中

cpot[1]=1;

for(i=2;i<=M.n;i++)

cpot[i]=cpot[i-1]+num[i-1];//每一列第一个非零元素的位置

for(p=1;p<=M.t;p++)

{

j=M.data[p].col;

q=cpot[j];//确定第j列第一个非零元素的位置

T.data[q].row=M.data[p].col;

T.data[q].col=M.data[p].row;

T.data[q].item=M.data[p].item;

++cpot[j];//指向下一个位置

}

}

}

//创建一个三元组

void creatMatrix(thMatrix &M)

{

initMatrix_Th(M);

int i,j,k=1;

cout<<"输入矩阵的行数和列数：";

cin>>M.m>>M.n;

for(i=1;i<=M.m;i++)

{

for(j=1;j<=M.n;j++)

{

cin>>matrix[i][j];

if(abs(matrix[i][j])>0)

{

M.data[k].row=i;

M.data[k].col=j;

M.data[k++].item=matrix[i][j];

}

}

}

M.t=k-1;

}

//输出一个三元组

void outputMatrix(thMatrix M)

{

int i;

cout<<"矩阵的行和列："<<M.m<<" "<<M.n<<endl;

cout<<"非零元素个数:"<<M.t<<endl;

for(i=1;i<=M.t;i++)

{

cout<<M.data[i].row<<" "<<M.data[i].col<<" "<<M.data[i].item<<endl;

}

}

int main()

{

thMatrix M,T;

initMatrix_Th(M);

cout<<"输入一个矩阵："<<endl;

creatMatrix(M);

cout<<"输出该矩阵的三元组表："<<endl;

outputMatrix(M);

transMatrix_2(M,T);

cout<<"输出该矩阵的转置矩阵的三元组表："<<endl;

outputMatrix(T);

return 0;

}