#include <stdio.h>
#include <stdlib.h>
#define MAXSIZE 1250
typedef int ElemType ;
typedef struct
{
int i,j;
ElemType e;
}Triple;
typedef struct
{
Triple data[MAXSIZE+1];
int mu,nu,tu;
}TSMatrix;
/*三元组专置算法*/
void TransposeSMatrix(TSMatrix T,TSMatrix *M)
{
int i,col,p;
M->nu=T.mu;
M->mu=T.nu;
M->tu=T.tu;
for(p=col=1;col<=T.nu;col++)
{
for(i=1;i<=T.tu;i++)
{
printf("ok");
if(T.data .j==col)
{
M->data[p].i=T.data.j;
M->data[p].j=T.data.i;
M->data[p].e=T.data.e;
p++;
}
}
}
}
/*方法二*/
int FastTransposeSMatrix(TSMatrix T,TSMatrix *M)
{
int t,col,p,q;
int *num,*cpot;
M->nu=T.mu;
M->mu=T.nu;
M->tu=T.tu;
num=(int *)malloc(sizeof(int)*(T.nu+1));
cpot=(int *)malloc(sizeof(int)*(T.nu+1));
for(t=1;t<=T.nu;t++) num[t]=0;
for(p=1;p<=T.tu;p++)
{
++num[T.data[p].j];
}
cpot[1]=1;
for(col=2;col<=T.nu;col++)
cpot[col]=cpot[col-1]+num[col-1];
for(p=1;p<=T.tu;p++)
{
q=cpot[T.data[p].j];
M->data[q].i=T.data[p].j;
M->data[q].j=T.data[p].i;
M->data[q].e=T.data[p].e;
++cpot[T.data[p].j];
}
free(num);
free(cpot);
return 1;
}
int main()
{
int i,col,p;
Triple z[]={{1,2,12},{1,3,9},{3,1,-3},{3,6,12},{4,3,24},
{5,2,18},{6,1,15},{6,4,-7}};
TSMatrix T ,M;
T.mu=6;
T.nu=6;
T.tu=8;
for(i=0;i<8;i++)
{
T.data=z;
}
/*TransposeSMatrix(T,&M);*/
/*输出*/
FastTransposeSMatrix(T,&M);
printf("%d /n",M.tu);
for(i=1;i<M.tu;i++)
{
printf("%d %d %d/n",M.data.i,M.data.j,M.data.e);
}
printf("/n");
return 0;
}
#include <stdlib.h>
#define MAXSIZE 1250
typedef int ElemType ;
typedef struct
{
int i,j;
ElemType e;
}Triple;
typedef struct
{
Triple data[MAXSIZE+1];
int mu,nu,tu;
}TSMatrix;
/*三元组专置算法*/
void TransposeSMatrix(TSMatrix T,TSMatrix *M)
{
int i,col,p;
M->nu=T.mu;
M->mu=T.nu;
M->tu=T.tu;
for(p=col=1;col<=T.nu;col++)
{
for(i=1;i<=T.tu;i++)
{
printf("ok");
if(T.data .j==col)
{
M->data[p].i=T.data.j;
M->data[p].j=T.data.i;
M->data[p].e=T.data.e;
p++;
}
}
}
}
/*方法二*/
int FastTransposeSMatrix(TSMatrix T,TSMatrix *M)
{
int t,col,p,q;
int *num,*cpot;
M->nu=T.mu;
M->mu=T.nu;
M->tu=T.tu;
num=(int *)malloc(sizeof(int)*(T.nu+1));
cpot=(int *)malloc(sizeof(int)*(T.nu+1));
for(t=1;t<=T.nu;t++) num[t]=0;
for(p=1;p<=T.tu;p++)
{
++num[T.data[p].j];
}
cpot[1]=1;
for(col=2;col<=T.nu;col++)
cpot[col]=cpot[col-1]+num[col-1];
for(p=1;p<=T.tu;p++)
{
q=cpot[T.data[p].j];
M->data[q].i=T.data[p].j;
M->data[q].j=T.data[p].i;
M->data[q].e=T.data[p].e;
++cpot[T.data[p].j];
}
free(num);
free(cpot);
return 1;
}
int main()
{
int i,col,p;
Triple z[]={{1,2,12},{1,3,9},{3,1,-3},{3,6,12},{4,3,24},
{5,2,18},{6,1,15},{6,4,-7}};
TSMatrix T ,M;
T.mu=6;
T.nu=6;
T.tu=8;
for(i=0;i<8;i++)
{
T.data=z;
}
/*TransposeSMatrix(T,&M);*/
/*输出*/
FastTransposeSMatrix(T,&M);
printf("%d /n",M.tu);
for(i=1;i<M.tu;i++)
{
printf("%d %d %d/n",M.data.i,M.data.j,M.data.e);
}
printf("/n");
return 0;
}