稀疏矩阵

#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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值