稀疏矩阵三元组乘法进行输出c语言

假设稀疏矩阵用三个数组来表示，分别为行下标数组row，列下标数组col和值数组val。代码如下：

void sparse_matrix_multiply(int rowA[], int colA[], int valA[], int nA, int mA, int kA,
                            int rowB[], int colB[], int valB[], int nB, int mB, int kB) {
    if(mA != nB) {
        printf("Error: The number of columns of A is not equal to the number of rows of B.\n");
        return;
    }
    // 初始化稀疏矩阵C
    int max_size = nA * mB;
    int rowC[max_size+1], colC[max_size], valC[max_size];
    rowC[0] = 0;
    int i, j, k, p, q;
    for(i = 0; i < nA; i++) {
        int row_start = rowC[i];
        for(j = 0; j < mB; j++) {
            int v = 0;
            for(p = rowA[i]; p < rowA[i+1]; p++) {
                // 在B矩阵中查找对应元素
                for(q = rowB[colA[p]]; q < rowB[colA[p]+1]; q++) {
                    if(colB[q] == j) {
                        v += valA[p] * valB[q];
                        break;
                    }
                }
            }
            if(v != 0) {
                colC[row_start] = j;
                valC[row_start] = v;
                row_start++;
            }
        }
        rowC[i+1] = row_start;
    }
    // 输出稀疏矩阵C
    for(i = 0; i < nA; i++) {
        for(j = 0; j < mB; j++) {
            int v = 0;
            for(k = rowC[i]; k < rowC[i+1]; k++) {
                if(colC[k] == j) {
                    v = valC[k];
                    break;
                }
            }
            printf("%d ", v);
        }
        printf("\n");
    }
}

其中，nA和mA分别表示矩阵A的行数和列数，kA表示A矩阵中的非零元素个数；nB和mB分别表示矩阵B的行数和列数，kB表示B矩阵中的非零元素个数。函数通过遍历稀疏矩阵A和B的非零元素，计算出稀疏矩阵C的每个元素，并将其存储在rowC、colC和valC数组中。最后再按照普通矩阵的方式输出稀疏矩阵C。注意，需要判断每个位置是否有非零元素，如果有则输出其值，否则输出0。