//算法5.3快速转置法
#include<iostream>
using namespace std;
#define MAXSIAE 1000
struct {
int row,rol;
int e;
};
struct {
int date[MAXSIAZE+1];
int m,n,len;
};
/*
部分代码
基于矩阵的三元组表示,采用一次定位快速转置法
三元组 只记录矩阵非零元素
*/
FastTransposeTSMatrix(TSMatrix A,TSMatrix *B)
{
int col,t,p,q;
/* num[MAXSIZE]表A第col列非零元素总个数
表B第col行非零元素总个数
position[MAXSIZE]表A第col列第一个非零元素的位置
表A第col行第一个非零元素的位置
第2行的即为第1行第一个非零元素位置+第1行第一行非零元素的个数
position[col]=position[col-1]+num[col-1];
*/
int num[MAXSIZE],position[MAXSIZE];
B->len=A.len;//len总非零元素个数
B->n=A.m;//A:m行n列 B:n行m列 B指针
B->m=A.n;
if(B->len)
{
//A:m行n列 B:n行m列
for(col=1;col<=A.n;col++)
num[col]=0;//初始化
//len总非零元素个数
for(t=0;t<A.len;t++)
num[A.date[t].col]++;
/*采用数组下标计数法,计算每一行的非零元素的个数
如果自己写date[t]数组不一定能写出来从0到len 非零元素总个数遍历
date[t].col 统计表A第col列非零元素总个数
例如 1只要遇到表A非零元素的列为 1就num[1]+1,为 2就num[2]+1,
遍历结束后num[1]到 num[MAXSIZE]就全部统计完
A.date[t].col 结构体存放
*/
position[1]=1;//初始化第一个位置
//依次计算第二个位置...
for(col=2;col<=A.n;col++)
position[col]=position[col-1]+num[col-1];
//将A从头扫描,实现矩阵转化
for(p=1;p<=A.len;p++){
//下标的处理 下标好多
col=A.date[p].col;
q=position[col];//正确位置
//行列互换,元素赋值
B->date[q].row=A.date[p].col;
B->date[q].col=A.date[p].row;
B->date[q].e=A.date[p].e;
// position[col]+1指向下一个列号为col的非零元素在三元组表B中的存放位置
position[col]++;
}
}
}
//A B 怎么定义??