稀疏矩陣三元組（順序結構）

此處有前提

稀疏矩陣有以下特點：

    1.大量數值為0的元素

    2.非零元素分佈沒有規律

這兩點的描述都比較含糊，不夠理性，一般來說是自己把握的。也是因為這兩點特性，所以才能比較自然地理解這種資料結構。

先看執行結果如下：

原始碼如下：

#define ElemType int
#include<stdlib.h>
#include<stdio.h>
#include<iostream>
using namespace std;
typedef struct{
	int i,j;
	ElemType e;
}Triple;//第i行第j列的元素值為e 
typedef struct{
	Triple *data;
	int row,column,tu;//矩陣的行數、列數、和非零元個數 
}SMatrix;
void CreateSMatrix(SMatrix &M){	
	cout<<"輸入矩陣的行數、列數、和非零元個數："; 
	cin>>M.row>>M.column>>M.tu;
	M.data=(Triple *)malloc((M.tu+1)*sizeof(Triple));// 
	cout<<"輸入非零元的行下標、列下標、值："<<endl;
	for(int k=1;k<=M.tu;k++)
		cin>>M.data[k].i>>M.data[k].j>>M.data[k].e;
} 
bool TransposeSMatrix(SMatrix M,SMatrix &T){
	if(!M.tu) return false;
	T.row=M.row;T.column=M.column;T.tu=M.tu;
	int num[M.tu+1];	//儲存每一列非零元的個數 
	int cpot[M.tu+1];	//儲存每一列非零元在T.data中的下標 
	for(int k=1;k<=M.tu;k++)	//矩陣下標從1開始 
		num[k]=0;
	for(int k=1;k<=M.tu;k++)	//求M中每一列非零元個數 
		num[M.data[k].j]++;
	cpot[1]=1;
	for(int k=2;k<=M.tu;k++)	//第k列中第一個非零元在T.data中的序號 
		cpot[k]=cpot[k-1]+num[k-1];
	for(int k=1;k<=M.tu;k++){	//M矩陣中按行有序，自上而下轉置即可 
		int temp=M.data[k].j;
		int q=cpot[temp];
		T.data[q].i=M.data[k].j;
		T.data[q].j=M.data[k].i;
		T.data[q].e=M.data[k].e;
		cpot[temp]++;
	}
	return true;
}
void PrintSMatrix(SMatrix M){
	cout<<"\n輸出用正常矩陣表示：\n";
	ElemType A[M.row+1][M.column+1];
	for(int k=1;k<=M.row;k++){
		for(int l=1;l<=M.column;l++){
			A[k][l]=0;
		}
	}
	if(M.tu){
		for(int k=1;k<=M.tu;k++){
			A[M.data[k].i][M.data[k].j]=M.data[k].e;		
		}		
	}
	for(int k=1;k<=M.row;k++){
		for(int l=1;l<=M.column;l++){
			cout<<A[k][l]<<" ";
		}
		cout<<endl;
	}
}
int main(){
	SMatrix M,T;
	CreateSMatrix(M);
	cout<<"\n以下這個矩陣是用來存放倒置後的矩陣，其數值不影響倒置\n"; 
	CreateSMatrix(T);
	TransposeSMatrix(M,T);
	PrintSMatrix(T);
	return 0;
}