稀疏矩阵

 /*
思路：将稀疏矩阵a转置为b，先扫描一遍a三元表（行列值），
获得a各列元素个数->得到b各行元素个数->知道a各元素在b中的插入位置->再扫描一遍a即可
T（n）=O(Cols+Terms);
*/
1 const int MaxTerms = 100;
class SparseMatrix
{public:
    SparseMatrix  FastTranspose();
    SparseMatrix Add(SparseMatrix b);
    SparseMatrix Multiply(SparseMatrix b);
private:
    int Rows, Cols, Terms;//矩阵的行数 列数 非零项数
    MatrixTerm smArray[MaxTerms];
};
class MatrixTerm
{
    friend SparseMatrix;
private:
    int row, col, value;//非零项的位置和值。三元表行优先
};
//稀疏矩阵快速转置
SparseMatrix SparseMatrix::FastTranspose()//将矩阵a(自身)转置为b
{
    int* RowSize = new int[Cols];//下标为行
    int* RowStart = new int[Cols];
    SparseMatrix b;
    b.Rows = Cols;b.Cols = Rows;b.Terms = Terms;
    if (Terms > 0)
    {
        for (int i = 0;i < Cols;++i)
            RowSize[i] = 0;
        for (int i = 0;i < Cols;++i)
            RowSize[smArray[i].Cols]++;
        RowStart[0] = 0;
        for (int i = 1;i < Cols;++i)
            RowStart[i] = RowStart[i - 1] + RowSize[i - 1];//RowStart目前存放矩阵b各行元素起始位置
        for (int i = 1;i < Terms;++i)
        {
            int j = RowStart[smArray[i].col];//smArray[i]在b三元表中的位置j
            b.smArray[j].row = smArray[i].col;
            b.smArray[j].col= smArray[i].row;
            b.smArray[j].value = smArray[i].value;
            RowStart[smArray[i].col]++;//RowStart存放本行下一个元素进来的位置
        };
    };
    delete[]RowSize; delete[] RowStart;
    return b;
};