程序代码：

 Status CreateSMatrix(RLSMatrix *M)

 { /* 创建稀疏矩阵M */
   int i;
   Triple T;
   Status k;
   printf("请输入矩阵的行数,列数,非零元素数：");
   scanf("%d,%d,%d",&(*M).mu,&(*M).nu,&(*M).tu);
   (*M).data[0].i=0; /* 为以下比较做准备 */
   for(i=1;i<=(*M).tu;i++)
   {
     do
     {
       printf("请按行序顺序输入第%d个非零元素所在的行(1～%d),列(1～%d),元素值：",i,(*M).mu,(*M).nu);
       scanf("%d,%d,%d",&T.i,&T.j,&T.e);
       k=0;
       if(T.i<1||T.i>(*M).mu||T.j<1||T.j>(*M).nu) /* 行、列超出范围 */
         k=1;
       if(T.i<(*M).data[i-1].i||T.i==(*M).data[i-1].i&&T.j<=(*M).data[i-1].j) /* 没有按顺序输入非零元素 */
         k=1;
     }while(k); /* 当输入有误,重新输入 */
     (*M).data[i]=T;
   }
   for(i=1;i<=(*M).tu;i++) /* 计算rpos[] */
     if((*M).data[i].i>(*M).data[i-1].i)
       for(T.i=0;T.i<(*M).data[i].i-(*M).data[i-1].i;T.i++)
         (*M).rpos[(*M).data[i].i-T.i]=i;
   for(i=(*M).data[(*M).tu].i+1;i<=(*M).mu;i++) /* 给最后没有非零元素的几行赋值 */
     (*M).rpos[i]=(*M).tu+1;
   return OK;

 }

程序代码：

/*
看不懂说明，可以试运行改程序！
*/

#include <stdio.h>

#define MAXSIZE 100
#define MAXRC 20
#define OK 1

typedef int ElemType;
typedef int Status;

typedef struct
{
    int i, j;
    ElemType e;
}Triple;
typedef struct
{
    Triple data[MAXSIZE+1];
    int rpos[MAXRC+1];
    int mu, nu, tu;
}RLSMatrix;

Status CreateSMatrix(RLSMatrix *M)
{ /* 创建稀疏矩阵M */
   int i;
   Triple T;
   Status k;
   printf("请输入矩阵的行数,列数,非零元素数：");
   scanf("%d,%d,%d",&(*M).mu,&(*M).nu,&(*M).tu);
   (*M).data[0].i=0; /* 为以下比较做准备 */
   for(i=1;i<=(*M).tu;i++)
   {
     do
     {
       printf("请按行序顺序输入第%d个非零元素所在的行(1～%d),列(1～%d),元素值：",i,(*M).mu,(*M).nu);
       scanf("%d,%d,%d",&T.i,&T.j,&T.e);
       k=0;
       if(T.i<1||T.i>(*M).mu||T.j<1||T.j>(*M).nu) /* 行、列超出范围 */
         k=1;
       if(T.i<(*M).data[i-1].i||T.i==(*M).data[i-1].i&&T.j<=(*M).data[i-1].j) /* 没有按顺序输入非零元素 */
         k=1;
     }while(k); /* 当输入有误,重新输入 */
     (*M).data[i]=T;
   }
   for(i=1;i<=(*M).tu;i++) /* 计算rpos[] */
     if((*M).data[i].i>(*M).data[i-1].i) //比较非零数组中该非零元所在行与前一非零元所在行的大小
       for(T.i=0;T.i<(*M).data[i].i-(*M).data[i-1].i;T.i++)  //如果大于，则检查非零数组中相邻非零元所差的行数，并求出对应行首个非零元地址存入rpos[]中
       {
           (*M).rpos[(*M).data[i].i-T.i]=i;                    //比如：相邻非零数组中第三个元素data[3]与第四个元素data[4]的行(假设data[4].i=5，data[3].i=3)相差2
            printf("rpos[%d]=%d  ", (*M).data[i].i-T.i, i);    //则，rpos[4]=rpos[5]=4(即四、五行首个非零元在data[]数组中的位置为下标4)
       }
       for(i=(*M).data[(*M).tu].i+1;i<=(*M).mu;i++)      //最后一个非零元所在行下一行首非零元为data中最后一个非零元下一个地址，这对后面界定所在行非零元范围有用                
     (*M).rpos[i]=(*M).tu+1;                          //比如：存在非零元的最后一行界定范围从rpos[M.data[M.tu]]到M.tu+1即(rpos[M.data[M.tu]+1])
   return OK;                                         //设置成零的话，因为rpos[0]中是一个垃圾值可能会影响计算。
}

int main(void)
{
    RLSMatrix M;
    printf("请输入提供例子：本矩阵5行3列4个非零元(5,3,4)，\n第一行一个非零元(1,1,1)，\n第二行一个(2,2,2)，\n第三行一个(3,3,3)，\n第四行没有，第五行一个(5,2,4)\n");
    CreateSMatrix(&M);
    return 0;
}