稀疏矩阵C 语言实现-豆柴文库

您所在位置：网站首页 / 稀疏矩阵C 语言实现.doc / 文档详情

免费试读已结束，剩余 1 页请下载文档后查看

16 金币

下载文档

/ 6

下载提示

如果您无法下载资料，请参考说明：

1、部分资料下载需要金币，请确保您的账户上有足够的金币

2、已购买过的文档，再次下载不重复扣费

3、资料包下载后请先用软件解压，在使用对应软件打开

文本预览

#include<iostream>
#include<iomanip>
usingnamespacestd;

template<classT>
//三元组
structTrituple
{
	introw;
	intcol;
	Tval;
};
//稀疏矩阵声明
template<classT>
classSparseMatrix
{
public:
	SparseMatrix(intmaxt=100);
	~SparseMatrix();
	boolTransposeTo(SparseMatrix&);
	boolAddTo(constSparseMatrix&);
	boolTransposeTo_Faster(SparseMatrix&);
	voidInput();
	voidOutput();
private:
	Trituple<T>*data;
	introws,cols,terms;
	intmaxterms;
};
template<classT>
SparseMatrix<T>::SparseMatrix(intmaxt)
{
	maxterms=maxt;
	data=newTrituple<T>[maxterms];
	terms=rows=cols=0;
}
template<classT>
SparseMatrix<T>::~SparseMatrix()
{
	if(data!=NULL)
	{
		delete[]data;
	}
}
//普通转置
template<classT>
boolSparseMatrix<T>::TransposeTo(SparseMatrix&B)
{
	if(terms>B.maxterms)
	{
		returnfalse;
	}
	B.rows=cols;
	B.cols=rows;
	B.terms=terms;
	if(terms>0)
	{
		intp=0;
		for(intj=1;j<=cols;j++)
		{
			for(intk=0;k<terms;k++)
			{
				if(data[k].col==j)
				{
					B.data[p].row=j;
					B.data[p].col=data[k].row;
					B.data[p].val=data[k].val;
					p++;
				}
			}
		}
	}
	returntrue;
}
//快速转置
template<classT>
boolSparseMatrix<T>::TransposeTo_Faster(SparseMatrix&B)
{
	if(terms>B.maxterms)
	{
		returnfalse;
	}
	B.rows=cols;
	B.cols=rows;
	B.terms=terms;
	if(terms>0)
	{
		int*num,*cpot;
		num=newint[cols];
		cpot=newint[cols];
		for(intj=0;j<cols;j++)
		{
			num[j]=0;
		}
		for(intk=0;k<terms;k++)
		{
			num[data[k].col-1]++;
		}
		//求出B中每一行的起始下标cpot[]
		cpot[0]=0;
		for(intj=1;j<cols;j++)
		{
			cpot[j]=cpot[j-1]+num[j-1];
		}
		//执行转置操作
		for(intp,k=0;k<terms;k++)
		{
			p=cpot[data[k].col-1]++;
			B.data[p].row=data[k].col;
			B.data[p].col=data[k].row;
			B.data[p].val=data[k].val;
		}
		delete[]num;
		delete[]cpot;
	}
	returntrue;
}
template<classT>
voidSparseMatrix<T>::Input()
{
	cout<<"intputthematrix'row:";
	introw1;
	cin>>row1;
	cout<<"intputthematrix'col:";
	intcol1;
	cin>>col1;
	cout<<"input"<<row1<<"*"<<col1<<"matrix"<<endl;
	intnumber;
	rows=row1;
	cols=col1;
	for(inti=0;i<rows;i++)
	{
		for(intj=0;j<cols